This is part of my continuing my work on making sympy capable of generating
code for evaluating matrix functions. So far the codeprinters and codegen
is done(ish); all that's left is autowrap.
Questions:
1. Should functions created by autowrap create matrices that would be
in-out parameters
1. Maybe if mat is NULL then we malloc? But probably the answer is what
you suggest with an explicit kwarg
2. Theano has a solution for this. It may be worth looking at how they
handle broadcasting. While this is accessible through SymPy it looks like
the SymPy documentation on this feature