OnPositiveSideOf fails if only one element between positive and negative side
Hi,
This is related to Gmsh issue #1750 regarding the usage of the Crack plugin in spiral geometries.
If one instead of cracking the mesh wants to use basis functions with support limited to one side of the turn-to-turn surfaces by e.g. using OnPositiveSideOf
(or equivalently OnNegativeSideOf
), the behavior as described in the above issue is apparent again.
That is, the classification into positive and negative side is not done correctly if there is a mesh element touches both the positive side of one turn and the negative side of another turn with two different nodes (see Gmsh issue #1750 for visualization). I suppose that a fix similar to that of #1750 would also fix this issue. Furthermore, a manual approach using OnOneSideOf
seems feasible and will serve as a hot fix.
Once I have the time, I will locate the corresponding source code and try to fix it. I'd be happy about any pointers on where to start. Until then, this issue serves to make other users aware of this problem.
Many thanks,
Erik