RE: point numbering when reading two meshes

[Lesage,Anne Cecile J]

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: getfem-users@nongnu.org
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 2^nd 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.