[frobaux]

Hey,


I have results from another run that I want to import, (on a different dt as the current one)


I would like to keep a field with the previous values and one for the next value (to interpolate between both of them). For now I have something like

PtrList<pointScalarField> potPhiPrevious(externalPhiRegions.size());
PtrList<pointScalarField> potPhiNext(externalPhiRegions.size());

I set them reading files in 0/ and 0+outerDt/

forAll(externalPhiRegions, i)
{
 potPhiPrevious.set
    (
        i,
         new pointScalarField
       (
            IOobject
            (
                "Phi",
                pTime,
                externalPhiRegions[i],
                IOobject::MUST_READ,
                IOobject::NO_WRITE
            ),
            pMesh
        )
  );

 }

where pTime is of type word and stand for the previous time where this field is availae. And same for potPhiNext
The problem is at some times, I need to send

potPhiNext[i] into potPhiPrevious[i], and then read the new potPhiNext[i].



I've tried everything I could think of (

potPhiPrevious.set(i, potPhiNext[i]);
potPhiPrevious[i]=potPhiNext[i];

Always have an error at run that I dont understand...

#3  Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::Boundary::Boundary(Foam::DimensionedField<double, Foam::pointMesh> const&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::Boundary const&) at ??:?
#4  Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::GeometricField(Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh> const&) at ??:?
#5  Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::clone() const at ??:?
#6  Foam::tmp<Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh> >::ptr() const at ??:?

How could I do it properly, without re-reading the field or hard copying it.

Thanks a lot!

[frobaux]

After a entire day of following error, debugging and trying new stuffs, the only solution I came with was rereading it completely.



But the same problem appeared for the interpolated in time (I use a simple linear interpolation)


I set this one to be written at each time step:

    potPhi.set
    (
        i,
        new pointScalarField
        (
            IOobject
            (
                "Phi_interp",
                runTime.timeName(),
                externalPhiRegions[i],
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            potPhiPrevious[i]*(1-alph)+potPhiNext[i]*alph
        )
    );

But the neither runTime.write() nor potPhi[i].write() works, always the same memory problem.

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? in "/lib/x86_64-linux-gnu/libc.so.6"

Does anyone have an Idea?



Also, when I try to print potPhiPrevious[0].internalField() it always send an error like:

...
-2.24018e-09
-2.53537e-09
-2.85198e-09
)
;

    boundaryField
    {
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::Boundary::writeEntries(Foam::Ostream&) const at ~/These/OpenFOAM/v1712/OpenFOAM-v1712/src/OpenFOAM/lnInclude/GeometricBoundaryField.C:586
#4  Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::Boundary::writeEntry(Foam::word const&, Foam::Ostream&) const at ~/These/OpenFOAM/v1712/OpenFOAM-v1712/src/OpenFOAM/lnInclude/GeometricBoundaryField.C:574 (discriminator 2)
#5  Foam::Ostream& Foam::operator<< <double, Foam::pointPatchField, Foam::pointMesh>(Foam::Ostream&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh> const&) at ~/These/OpenFOAM/v1712/OpenFOAM-v1712/src/OpenFOAM/lnInclude/GeometricField.C:1337 (discriminator 2)
#6  Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::writeData(Foam::Ostream&) const at ~/These/OpenFOAM/v1712/OpenFOAM-v1712/src/OpenFOAM/lnInclude/GeometricField.C:1036
#7  ? at ~/These/OpenFOAM/v1712/OpenFOAM-v1712/src/OpenFOAM/lnInclude/DimensionedFieldIO.C:151
#8  ? at ~/OpenFOAM/robaux-v1712/applications/solver/myCoupledHPCSolver/./externalPhi/createExternalPhiFields.H:115 (discriminator 7)
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? at ??:?

If someone has any idea of where to look of understand this, I would be very grateful =)

Thanks a lot !

[frobaux]

well, after a while and reading all treads containing PtrList and pointScalarField, I ended up finding the error, and it was my mistake..
The Sigrev error meant that there was empty memory somewhere and I tought it was due to the PtrList stuffs, but I was actually due to the fact that the pointMesh(es) associated with the pointScalarField(s) were not stored in memory.

I endend up adding a

PtrList<pointMesh> pMeshes(type1Regions.size())