Selon Andriy Andreykiv <[EMAIL PROTECTED]>: > Dear Yves, > > When I create a new level set on a getfem_mesh, I have to also set > values of > > the level set function by setting > > ls.values( )[i] = .. > > for all the nodes of my original mesh, that I want to cut by the level set. > I browsed through the examples, located in getfem instalation, and I see that > > all of them have quite simple procedure of doing this. Namelly, in every > example this procedure strongly depends on the geometry of the problem and I > > didn't see some generic way of doing this. This is particularly a problem if > > your level set (geometry of your interface) is quite arbitrary, as in my > case. > Is there a way to set those values in some generic way? > By, say, computing a signed distance between your interface (possibly defined > > as a getfem::mesh itself) and the nodes where you have to define the value of > > your level set? > > Thank you in advance, > Andriy
Dear Andriy, In order to work without problems, it is indeed important for the level set function to be closed to (or not too far from) a signed distance. So if you have a signed distance to a surface, you can interpolate it on the mesh_fem of the level set (of course, there is an approximation there). But may be I did not understand well your question. Best regards, Yves. _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
