I am using the PDE Toolbox to perform a FEM simulation. I am using the following toolbox function to aquire the matrices I am after:
FEM = assembleFEMatrices(model,'nullspace')
However, the only matrices I require are the global stiffness matrix (FEM.Kc) and nodal force vector (FEM.Fc).
The issue is that the assembleFEMatrices function is quite costly. I was wondering if there is a way improve the efficiency of this? Could the stiffness matrix and force vector be calculated another way, i.e. could they be calculated from the results of the 'solve' function that is less time consuming:
Or if it is possible to assemble ONLY the matrices you desire to reduce the computational expense?
Thank you for any advices/solutions that can be offered