Hi guys, I am trying to programatically generate a Julia function that calculates a 6x6 matrix. I have written a short LaTeX description of the problem at:
http://nbviewer.ipython.org/github/berceanu/notebooks/blob/master/OPODrag/block_matrix.ipynb The basic idea is, I have an analytical expression for the matrix elements of the 3x3 matrices M and Q, M_{m,n} representing the element on row m, column n. The problem is that these depend on sums over Kronecker symbols, which need to be contracted. Once I build the M and Q matrices, I need to join them in the block-form given in the above notebook. A working strategy that I have so far is to generate the matrix L in Mathematica, and copy-paste the resulting (huge) expression to Julia. I would now like to bypass Mathematica completely, and thought of using Sympy. While it has the capability of contracting the sums over the delta symbols, it is not clear to me how to then export the resulting function to Julia. If this can be done without Sympy, i.e. in pure Julia, that would be awesome! On a last note, performance of the resulting function is important, because I will need to digonalize a lot of these 6x6 matrices, calculated for different parameter values (eg. varying k_p, k_s, omega_p, etc).
