Dear all,

I have a simple question regarding a matrix assembly in Firedrake which I believe is not computed correctly.

The relevant code can be found here and the issue is with the B_mu matrix in solvers.py, line 51. I compare the result with my own Matlab code and the two don't agree. Also, the rest of the matrices M, A, C_lam agree in both sets of code.

Can someone help with what the problem might be?

Thanks,

Anna.

-- 


 Dr Anna Kalogirou
 Research Fellow
 School of Mathematics
 University of Leeds

 http://www1.maths.leeds.ac.uk/~matak/