CFD Online Discussion Forums - print coefficent of the linear system

Dear foamer,

I am modifying the pimpleFoam solver to achieve a LES nesting (common operation in NWP models (numerical weather perdition)). For debug purpose, I would like to print the matrix of the linear system Ax=b.

According this page http://openfoamwiki.net/index.php/Op...es_in_OpenFOAM, those coefficients should be accessible with

Code:

scalarField *lowerPtr_; scalarField *diagPtr_; scalarField *upperPtr_;

But I'dont really know the object which holds those pointer. Any idea?

Once those scalarField available, how to know the location of their conponant in the matrix (from scalarField(1D datas) to a matrix (2D data))?

Cheers,
Marcel