Hello,
I am working on POD decomposition and for obtaining the interpolate solution I need the stiffness matrix and load vector but doing assemble and applying the boundaries conditions I get a wrong stiffness matrix. I don’t if it because the Dolfin stiffness matrix is sorted out using a local DOF or why but I obtaining a wrong all the time. I want the stiffness matrix for operating using numpy matrices. How can I get the stiffness matrix in global coordinates?
Thank you in advance