[jlr]

Hi,

I am trying to use a pointScalarField in my solver for foam-extend 4.0, which compiles with no issue. However, while running the solver crashes with a segmentation fault at output timesteps.

Debugging I know the crashing line is GeometricBoundaryField.C:530

this->operator[](patchi).patch().name()

It crashes when calling the name() method. patch() returns a pointPatch object.

The name() method in the pointPatch class is a pure virtual function (see pointPatch.H:127), so that's why calling it crashes the solver. Debugging other foam-extend's default solvers that use a pointTypeField it seems most call the name() method from facePointPatch, which is a base class for most pointPatch types.

So, I guess the problem is that the patches of my pointScalarField are not of any type. Shouldn't they know the type from boundary definitions at polymesh/boundary?. All of my patches are of type patch. I just don't see what I'm doing differently from the default solvers...

Thank you

[hjasak]

Hi,

My guess is that the pointMesh is not constructed correctly.

Just callin a normal constructor won’t do: you need a reference to a meshObject

Const pointMesh& pm = pointMesh::New(mesh);

Hrv

[jlr]

Turns out my (silly) mistake was creating the pointMesh object locally on a function, thus erasing it from memory once the function ended.

Thank you for your help!

[hjasak]

What can I say debugging service without looking at the code!

I think I’ll call myself Professor for the rest of the day...