[jlr]

Hi,

I'm programming a solver in foam-extend 4.0 that needs a bool value defined in every point of a finite area mesh (faMesh)

At first I used a simple boolIOList, which works but doesn't get decomposed / reconstructed when using decomposePar / reconstructPar. So, now I'm trying to use a GeometricField, namely a pointBoolField (it is not defined in the source code, I define it myself in the solver)

However, there doesn't seem to be a way to define a pointTypeField on the faMesh points only. pointTypeFields require a pointMesh in their constructor. For instance, in $FOAM_SOLVERS/solidMechanics/stressFemFoam/calculateStress.H :

// Create pointMesh for field post-processing
const pointMesh& pMesh = pointMesh::New(mesh);

// U
pointVectorField Up
(
    IOobject
    (
        "Up",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    pMesh,
    U.dimensions()
);

The pointMesh class only has one constructor, which calls on a polyMesh, which can only represent 3D volume meshes and not a faMesh

I've tried defining the pointBoolField on the fvMesh tied to the faMesh, and then only work with the points included in the faMesh. I haven't been successful though, even when using a type defined in the source code such pointScalarField instead, but for the sake of not making this post too long I won't go into details about that for now. I first want to confirm that there is no other way to work with a pointTypeField on a faMesh, since defining the field in the whole fvMesh seems like a waste of memory to me, and much less intuitive for the user.

Thank you