|
From: | Konstantinos Poulios |
Subject: | Re: point numbering when reading two meshes |
Date: | Wed, 24 Nov 2021 15:36:50 +0100 |
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 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.
[Prev in Thread] | Current Thread | [Next in Thread] |