I want to implement my mapping between OpenFOAM and my structure solver. The mapping should be based on edges/faces.
As direct access to edges/faces is not supported by PreCICE, I want to implement it by accessing Fluid-Mesh-Nodes and Fluid-Mesh-Centers provided by OpenFOAM. I compute the edges/faces from both sets of vertices and recalculate data to my surface mesh.
I have an initial implementation. However, there is probably a bug in PreCICE, as it is impossible to call ‘setMeshAccessRegion’ for both Nodes and Centers. On the other hand, if I call it only once, I cannot access the second mesh.
If I remove the following line in ‘src/precice/impl/ParticipantImpl.cpp’
PRECICE_CHECK(!_accessRegionDefined, “setMeshAccessRegion may only be called once.”);
It works as expected (at least sequentially, but there should not be any difference in parallel runs).
Is it possible to fix it so that this check will be per mesh and not per application?
Or is there any problem with it?