I have run into an issue, which I am not sure how to solve myself. I am solving multiphase (Navier-Stokes + phase-field) equations and marching in time. I use “adaptmesh” function to tailor my FE mesh to the phase function. The code I developed was working fine until recently I got the error:
“Exec error : periodic: Sorry one vertex of edge is losted”
I have found the error in the mailing list here, but adding “fixeborder=true” to “buildmesh” did not do anything. I also tried to make “adaptmesh” use old vertices with “keepbackvertices=true”, but that also did not help.
I created minimal example to reproduce the error (~ 8 MB ZIP file here). When I run it, I get:
minExample> FreeFem++ -ne -nw tryAdapt.edp -- FreeFem++ v4.6 (Tue 16 Jun 2020 06:53:31 AM CEST - git v4.6-113-g243cbb52) Load: lg_fem lg_mesh lg_mesh3 eigenvalue load: iovtk sizestack + 1024 =3856 ( 2832 ) I am here 1!!! I am here 2!!! Abscisse: s0 = 0.399285 <==> s1 0.315121 current line = 78 Exec error : periodic: Sorry one vertex of edge is losted -- number :1 Exec error : periodic: Sorry one vertex of edge is losted -- number :1 err code 8 , mpirank 0
Any idea regarding what I am doing wrong? I have also checked the “c” function I am using as an input for adaptmesh and did not see any weird behaviour.
My FreeFEM version is v4.6 (Tue 16 Jun 2020 06:53:31 AM CEST - git v4.6-113-g243cbb52).