Dear Lorenzo, It would be very useful if you could provide a minimal test case that demonstrates this issue. Strange that md.set_time() helps updating the model correctly.
Best regards Kostas On Mon, Jul 18, 2022 at 4:36 PM Lorenzo Ferro <lorenzo.i...@gmail.com> wrote: > Dear All, > > I'm experimenting with mesh changes (a sequence like del_convex -> > del_point -> add_point -> add_convex -> set_region) and relevant MeshFEM > automatic adaptation with the Python interface. > > I found out that as soon as I update mesh and regions, if after that I try > to run a test assembly method for example just to calculate the total area > of the mesh elements (via gf.asm_generic) I achieve a completely wrong > result, > but, > after just running a Model method like "md.set_time(0.)" and then running > immediately "gf.asm_generic" then I achieve the correct result! > > Is this a bug? or is there some logic behind this behaviour? > > Thank you, > Lorenzo >