Joris Gillis jgillis

CasADi dev, postdoc mech engineering @ KULeuven, freelancer through yacoda

jgillis / mfunctionparser.py

Created March 14, 2013 12:58

Parse m-functions in casadi

	import casadi
	import re

	class MFunctionParser:
	def parse(self,filename):
	f = open(filename,"r")

	signature = f.readline()
	m = re.search("function \[?([\w,]+)\]? = \w+\((.*)\)",signature)
	if not m:

jgillis / multisens

Created May 24, 2013 07:49

casadi multiple sensitivities

jgillis / gist:5683620

Last active December 17, 2015 22:39

Method to print the class inheritance hierarchy

	def print_subclasses(myclass, depth=0):
	print (" " * depth) + " - " + myclass.__name__
	for s in myclass.__subclasses__():
	print_subclasses(s,depth=depth+1)

jgillis / qpsolve.py

Created August 4, 2014 11:29 — forked from jaeandersson/qpsolve.py

	import numpy as NP
	import casadi as C
	def qpsolve(H,g,lbx,ubx,A=NP.zeros((0,0)),lba=NP.zeros(0),uba=NP.zeros(0)):
	# Convert to CasADi types
	H = C.DMatrix(H)
	g = C.DMatrix(g)
	lbx = C.DMatrix(lbx)
	ubx = C.DMatrix(ubx)
	A = C.DMatrix(A)
	A = A.reshape((A.size1(),H.size1())) # Make sure matching dimensions

jgillis / HackyFunction.py

Last active September 6, 2021 13:37

using vertcat(MX,MX) as inputs to Function

	from casadi import *

	def HackyFunction(name,ins,outs):
	# Obtain all symbolic primitives present in the inputs
	Vs = symvar(veccat(*ins))

	# Construct a helper function with these symbols as inputs
	h = Function("helper",Vs, ins)

	# Assert dense inputs

jgillis / gist:c9c381e9c65b059cbe7b

Last active August 29, 2015 14:15

treemap

	<!DOCTYPE html PUBLIC '-//W3C//DTD XHTML 1.0 Strict//EN' 'http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd'>
	<html xmlns='http://www.w3.org/1999/xhtml'>
	<head><script src='http://d3js.org/d3.v3.min.js'></script><script src='http://cpettitt.github.io/project/dagre-d3/v0.1.5/dagre-d3.min.js'></script><style>svg { overflow: hidden;}.node rect { stroke: #333; stroke-width: 1.5px; fill: #fff;}.edgeLabel rect { fill: #fff;}.edgePath { stroke: #333; stroke-width: 1.5px; fill: none;}.outer { width: 1024px; height: 960px; overflow: auto;}.inner { width: 8000px; height: 6000px;}svg { display: block; width: 100%; height: 100%;}


	.link {
	stroke: #999;
	stroke-width: 2px;
	}

jgillis / ipynbdoctest.py

Last active August 29, 2015 14:17 — forked from akimd/ipynbdoctest,py

	#!/usr/bin/env python
	"""
	simple example script for running and testing notebooks.

	Usage: `ipnbdoctest.py foo.ipynb [bar.ipynb [...]]`

	Each cell is submitted to the kernel, and the outputs are compared with those stored in the notebook.
	"""

	import os,sys,time

jgillis / index.py

Created June 24, 2015 06:18

	import os
	import subprocess
	import io
	import shutil


	dir = '/home/frs/project/casadi/CasADi/'

	shutil.rmtree(os.path.join(dir,'branches'))
	os.mkdir(os.path.join(dir,'branches'))

jgillis / sysid_singleshooting.m

Created August 16, 2015 14:37

	import casadi.*
	load SilverboxSimulated
	u_data = u;
	y_data = y;

	% simulation horizon
	N = size(u_data, 1);

	% declare symbols for states and controls
	y = MX.sym('y');

jgillis / sysid_multipleshooting.m

Created August 18, 2015 06:07