yes, GetFEM mesh regions work with element faces, not points. BR Kostas
On Wed, Nov 24, 2021 at 12:11 AM Lesage,Anne Cecile J < [email protected]> wrote: > Sorry the problem was not with node numbering > > Correct me if I am wrong but my conclusion is that you shall respect the > following ruke > > “When you give getfem a list of pointid to impose a boundary conditions, > make sure that the set of points defines faces” > > There were only disconnected surface points in my list > > > > *From:* Lesage,Anne Cecile J > *Sent:* Tuesday, November 23, 2021 1:27 PM > *To:* [email protected] > *Subject:* point numbering when reading two meshes > > > > Dear all > > > > I am reading two meshes in a python script for getfem > > > > > > meshb=gf.Mesh('import','gid','brainr25tet4h5.GiD.msh') > > meshh=gf.Mesh('import','gid','headr25tet4h5.GiD.msh') > > > > The I want to select points for Dirichlet boundary conditions with a list > that I am reading from files > > > > # selection of faces for dirichlet boundary condition > > DUB = 1 # dirichlet on displacement brain > > DPBI = 2 # dirichlet on displacement non homogeneous > > DUH = 3 > > > > print('Read DIRICHLET boundary condition brainstem'); > > # read dirichlet conditions bstem > > filebstem = open('bstem.txt', 'r') # 'r' = read > > pidbstem = np.loadtxt(filebstem,dtype=int) > > #print(pidbstem) > > filebstem.close() > > fbstem=meshb.faces_from_pid(pidbstem) > > meshb.set_region(DUB,fbstem) > > > > > > print('Read DIRICHLET boundary condition head fix\n'); > > # read dirichlet conditions on dipslacement > > filehead = open('headdir.txt', 'r') # 'r' = read > > pidhead = np.loadtxt(filehead,dtype=int) > > #print(pidhead) > > filehead.close() > > fhead=meshh.faces_from_pid(pidhead) > > meshh.set_region(DUH,fhead) > > > > the file I read consider that for each mesh part the numbering starts from > 1 > > I can observe that the point selection is working well for the first mesh > > However it does not work for the second mesh > > I think this is because getfem renumbers the nodes of the 2nd mesh. > > How can I be sure how it is done? > > The first mesh here has 3647 nodes. Does getfem start numbering the points > of the second mesh at 3648? > > Do I just need to edit my second list and add 3647? > > > > Thank you > > Regards > > 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. >
