Stiffness Matrix

Try:

mypde=LienarPDE(...)
mypde.setValue(..)
mypde.getOperator().saveMM("mymatrix.mm")

to write the stiffness matrix to the file mymatrix.mm in the Matrix Market File Format, see http://math.nist.gov/MatrixMarket/formats.html . This does not work for more than one MPI rank.

Thanks Lutz Gross, that solved my question.

Thanks for this Lutz, is there a similar way to extract the right hand side i.e f of KU = f.

Many Thanks T

Try:

mypde=LienarPDE(...)
mypde.setValue(..)
f=mypde.getRightHandSide()
for i in xrange(f.getNumDataPoints()):
          print f.getValueOfDataPointAsTuple(i)

This gives the following error

Attribute Error: 'Data' object has no attribute 'getDataPointRank'

Sorry meant to write 'getNumDataPoints' instead of 'getDataPointRank'

Thanks Lutz Gross, that solved my question.