[klausb]

Hello,

is it possible to access interfaceBouCoeffs in a custom function for a custom PCG solver version at the location shown below?

I have the impression, it's protected.

Code section:

namespace Foam
{ 
    defineTypeNameAndDebug(myPCG, 0);

    lduMatrix::solver::addsymMatrixConstructorToTable<myPCG>
        addmyPCGSymMatrixConstructorToTable_;
}

// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

Foam::myPCG::myPCG
(
    const word& fieldName,
    const lduMatrix& matrix,
    const FieldField<Field, scalar>& interfaceBouCoeffs,
    const FieldField<Field, scalar>& interfaceIntCoeffs,
    const lduInterfaceFieldPtrsList& interfaces,
    const dictionary& solverControls
)
:
    lduMatrix::solver
    (
        fieldName,
        matrix,
        interfaceBouCoeffs,
        interfaceIntCoeffs,
        interfaces,
        solverControls
    )
{}
                     
Foam::solverPerformance Foam::myPCG::solve 
(
    scalarField& psi,
    const scalarField& source,
    const direction cmpt
) const
{
    word precond_name = lduMatrix::preconditioner::getName(controlDict_);

// --- Setup class containing solver performance data
    solverPerformance solverPerf
    (
        lduMatrix::preconditioner::getName(controlDict_) + typeName,
        fieldName_
    );
          
    // How to get access to interfaceBouCoeffs_ for this custom function?
    myWork::custom_function(... interfaceBouCoeffs_); 

    return solverPerf;    
}

Klaus

[mAlletto]

Hello all,

I was just asking my how the interfaceBouCoeffs are calculated and were. I understand it is an argument to the constructor but I didn't found how they are calculated

[klausb]

Have a look at section 4/15 of: https://www.linkedin.com/pulse/openf...menico-lahaye/


This should give you some hints but it's not the complete answer to your question.