search

FreeFem / FreeFem-sources

FreeFEM source code
https://freefem.org/
Other
770 stars 191 forks source link

Trunc not working anymore #221

Closed ffeppon closed 2 years ago

ffeppon commented 2 years ago

Dear developpers,

I encounter again some bug on some meshes when I trunc them with FreeFEM. 

mesh Th=readmesh("Th_0000.mesh");
mesh Thf=trunc(Th,region==3);
mesh Ths=trunc(Th,region==2);

mesh.zip

This results in 

Assertion fail : (iseuil)
    line :1819, in file ../femlib/fem.cpp
Assertion fail : (iseuil)
    line :1819, in file ../femlib/fem.cpp
 err code 6 ,  mpirank 0

Many thanks if you can fix this bug.

Florian

frederichecht commented 2 years ago

I have correct in last push on develop https://github.com/FreeFem/FreeFem-sources/compare/develop?expand=1

frederichecht commented 2 years ago

Thank for find this mistake,