[luca]

Hello,

supposing that the fvMatrix.residual() member give us the residual of the matrix in the form R = Ax - B, I've noticed a strange behavior.

In some way .residual() member applied to fvMatrix<vector> modifies the extra diagonal elements of the matrix:

the operator .H() returns different results for fvMatrix<vector> before and after calling the fvMatrix.residual() member.

Attached you find a simple code to reproduce the "phenomena" , you can run it under tutorial directory of the scalarTransportFoam solver

Tested on 1.7.x version yesterday update

scalarTransportFoam.C

Thanks in advance