|From:||Lesage,Anne Cecile J|
|Subject:||point numbering when reading two meshes|
|Date:||Tue, 23 Nov 2021 19:27:03 +0000|
I am reading two meshes in a python script for getfem
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('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)
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?
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.
|[Prev in Thread]||Current Thread||[Next in Thread]|