[martel]

Hi all,

I am trying to se the effect of the boundary conditions on the matrix diagonal coefficients for a very simple laplacian operator.

I use this simple code:

// Define a matrix that contains the coefficients of the Laplacian
    tmp<fvScalarMatrix> TEqnLaplaciano 
	  ( 
	   fvm::laplacian(DT,T) 
	    ); 

// Add BC contribution to matrix diagonal 
     scalarField& diag    = TEqnLaplaciano().diag();
     TEqnLaplaciano().addBoundaryDiag(diag,0);

But it does not compile. I get the compilation error:

In file included from /opt/openfoam240/src/finiteVolume/lnInclude/fvMatrix.H:1058:0,
                 from /opt/openfoam240/src/finiteVolume/lnInclude/ddtScheme.C:29,
                 from /opt/openfoam240/src/finiteVolume/lnInclude/ddtScheme.H:325,
                 from /opt/openfoam240/src/finiteVolume/lnInclude/fvcDdt.C:28,
                 from /opt/openfoam240/src/finiteVolume/lnInclude/fvcDdt.H:199,
                 from /opt/openfoam240/src/finiteVolume/lnInclude/fvc.H:44,
                 from /opt/openfoam240/src/finiteVolume/lnInclude/fvCFD.H:8,
                 from laplacianFoam.C:32:
/opt/openfoam240/src/finiteVolume/lnInclude/fvMatrix.C: In function ‘int main(int, char**)’:
/opt/openfoam240/src/finiteVolume/lnInclude/fvMatrix.C:116:6: error: ‘void Foam::fvMatrix<Type>::addBoundaryDiag(Foam::scalarField&, Foam::direction) const [with Type = double; Foam::scalarField = Foam::Field<double>; Foam::direction = unsigned char]’ is protected
 void Foam::fvMatrix<Type>::addBoundaryDiag
      ^
laplacianFoam.C:54:45: error: within this context
      TEqnLaplaciano().addBoundaryDiag(diag,0);
                                             ^

Related to some variables protected (??). And I don't know how to solve this protection problem.

Any help would be greatly welcome

Regards,

--CME