[frobaux]

Hello Foamers,
I need help on a PtrList object,
My goal is to create a solver that has the capabilities of interFoam with multiregion physics.



So I convert interFoam

   volScalarField& alpha1(mixture.alpha1());

To

    PtrList<volScalarField> l_alpha1(fluidRegions.size());
    l_alpha1.set
    (
        i,
        l_mixture[i].alpha1()
    );

The problem is, when I update l_mixture[i].alpha1()=something (modify the internal field),
The l_alpha1[i] does not update and both fields are thus different.

Reciprocally, when I modify l_alpha1[i], l_mixture[i].alpha1() is not updated.

Maybe I'm wrong, but I order to mimic the interFoam behavior, both fields are in fact one, and should behave as such?

Thank you a lot for your help

[adhiraj]

Please correct me if I'm wrong, but from what I understand,

volScalarField& alpha1(mixture.alpha1());

creates a reference to the volScalarField returned by mixture.alpha1(), and thus if you change either, the other is updated.

When you write

l_alpha1.set
(
    i,
    l_mixture[i].alpha1()
);

you are initializing l_alpha1[i] using the "value" of the volScalarField returned by mixture.alpha1(). Thus they are not the same object.

[frobaux]

Originally Posted by adhiraj

When you write

l_alpha1.set
(
    i,
    l_mixture[i].alpha1()
);

you are initializing l_alpha1[i] using the "value" of the volScalarField returned by mixture.alpha1(). Thus they are not the same object.

Yes, This is what happens, but as l_alpha1 is a list of pointers (PtrList), I thought that l_alpha1[i] would point to the l_mixture[i].alpha1(). Maybe it is the normal behavior, but then I don't understand the name of the structure : "PtrList"...


Otherwise, is there a way to obtain l_alpha1[i] => l_mixture[i].alpha1().



Thanks a lot for your help!