Dear Anne-Cecile To be more precise, an "anatomical area" should be defined as a set of element faces and not as a set of discrete points. This is why GetFEM honours the definition of regions in terms of element volumes (entire elements) and element surfaces (faces).
You can apply BCs on points in GetFEM but you should in general avoid it. So I recommend you to either identify the "anatomic area" with a set of faces in your external program or to find inside GetFEM which faces include the given set of points. Then then define a standard mesh region for that. The second option you can implement with fc = mesh.faces_from_pid(point_ids) [image: image.png] Best regards Kostas On Tue, Nov 16, 2021 at 11:20 PM Lesage,Anne Cecile J < ajles...@mdanderson.org> wrote: > Dear all, > > > > I have a 3d mesh and a list of point id which I would like to read in a > python script. This is to define a set of points with Dirichlet boundary > conditions on displacement > > I see in the python examples the use of geometrical criteria, my mesh is a > complex organ and I would like to be able to build set of points that > corresponds to a specific anatomical area. > > > > Thank you > > Anne-Cecile > The information contained in this e-mail message may be privileged, > confidential, and/or protected from disclosure. This e-mail message may > contain protected health information (PHI); dissemination of PHI should > comply with applicable federal and state laws. If you are not the intended > recipient, or an authorized representative of the intended recipient, any > further review, disclosure, use, dissemination, distribution, or copying of > this message or any attachment (or the information contained therein) is > strictly prohibited. If you think that you have received this e-mail > message in error, please notify the sender by return e-mail and delete all > references to it and its contents from your systems. >
