printing global structural stiffness matrix

[cenk tort]

Hello ,

I added a new function as getA to the fullgenlinSOE class to print out the global structural stiffness matrix. I am calling this function in the newtonraphson.cpp right after the stiffness matrix is assembled. But I found that all entries of my global stiffness matrix is equal to my penalty constraint value which is 1e12. Which does not seem correct. Where should I print my global stiffness matrix to to see it before solving for the displacements ?

Thanks.
Cenk

[fmk]

you could put it at the top of the solve method.