[AliE]

Dear Foamers,

I am actually trying to print to screen the coefficients related to the pressure- correction equation in icoFoam (diagonal and off-diagonal). Until now, I have followed the suggestions in this post:

Accessing Coefficient Matrix

Here is the code:

// -------

fvScalarMatrix pEqn
(
fvm::laplacian(rAU, p) == fvc::div(phiHbyA)
);


label Nc = mesh.nCells(); //Total number of cells

//Diagonal coefficients
for(label i=0; i<Nc; i++)
{
Info<< pEqn.diag()[i] << endl;
}

//Off-diagonal coefficients
for(label faceI=0; faceI<UEqn.lduAddr().lowerAddr().size(); faceI++)
{
Info<< pEqn.upper()[faceI] << endl; //Lower coefficients
Info<< pEqn.lower()[faceI] << endl; //Upper coefficients
}

// ------

It seems to work, but for some reasons in the diagonal coefficents printed out the contribution of the boundary conditions is not included.

Do you have any suggestion to fix my problem?

Thank you very much for the help,

AliE