I encounter an assertion fails with some 2D meshes that I didn’t have before. Again these problems with the trunc function (which should be combinatorial, not ?)!
Th.o.mesh (197.0 KB)
The following code fails:
mesh Th=readmesh("Th.o.mesh"); mesh Thf=trunc(Th,region==3);
Many thanks in advance for your help.