[dlahaye]

Greetings,

The PETSC4FOAM solve() function calls the scalarSolve() with tpsi.ref() as argument. See L396 of https://develop.openfoam.com/modules.../petscSolver.C .

Why is tpsi.ref() used instead of tpsi? What is .ref() of tpsi ?

Thx!

[dlahaye]

I imagine .ref() to return a non-constant reference to tpsi as explained at
Difference between boundaryField() and boundaryFieldRef()

Correct?

[olesen]

Originally Posted by dlahaye

Greetings,

The PETSC4FOAM solve() function calls the scalarSolve() with tpsi.ref() as argument. See L396 of https://develop.openfoam.com/modules.../petscSolver.C .

Why is tpsi.ref() used instead of tpsi? What is .ref() of tpsi ?

Thx!

Take a look at tmp.H, refPtr.H and PrecisionAdaptor.H for where ref() is being used. In short, you are correct - it is for a modifiable reference to the underlying object. Unlike boundaryFieldRef() however, there is no event counter to track changes.

[dlahaye]

Thx!

What tmp.H, refPtr.H and PrecisionAdaptor.H do you mean? I do not find it immediately in PETSC4Foam.

Best. D.

[olesen]

Originally Posted by dlahaye

Thx!

What tmp.H, refPtr.H and PrecisionAdaptor.H do you mean? I do not find it immediately in PETSC4Foam.

Best. D.

They are all located here:

https://develop.openfoam.com/Develop...penFOAM/memory

[dlahaye]

Thx! Much appreciated!

[dlahaye]

I am trying to overwrite computed values of TEqn.psi() within laplacianFoam (or any other solver). I am struggling with the syntax. To modify the first component of TEqn.psi(), I tried the syntax

TEqn.psi().ref()[0] = 1.0;

and variants thereof. These variants result in various compiler error messages.

Help is appreciated.

[dlahaye]

Below for v1906 and v1912

GeometricField<scalar, fvPatchField, volMesh>& psi =
    const_cast<GeometricField<scalar, fvPatchField, volMesh>&>(TEqn.psi());
psi.primitiveFieldRef()[0] = 1.0;