Dear all
Now I understand why it prints three times, it is because mfub (displacements)
is three dimensions
Thank you
AC
From: Lesage,Anne Cecile J
Sent: Wednesday, January 19, 2022 10:18 AM
To: getfem-users@nongnu.org
Subject: FW: [EXT] Re: field output size
Dear Kostas
I tried this to test the nodes output and its transpose
The problem is that the nodes are printed sequentially three times in a row
np.savetxt("brain_displacements_nodeslines.txt",
mfub.basic_dof_nodes(),delimiter=' ', fmt='%f')
a = np.empty([meshb.nbpts(),3])
for i in range(meshb.nbpts()):
Pt=mfub.basic_dof_nodes(i)
a[i][0]=Pt[0]
a[i][1]=Pt[1]
a[i][2]=Pt[2]
np.savetxt("brain_displacements_nodes.txt", a, delimiter=' ', fmt='%f')
thank you
Anne-Cecile
From: Konstantinos Poulios
<logar...@googlemail.com<mailto:logar...@googlemail.com>>
Sent: Wednesday, January 19, 2022 4:49 AM
To: Lesage,Anne Cecile J
<ajles...@mdanderson.org<mailto:ajles...@mdanderson.org>>
Cc: getfem-users@nongnu.org<mailto:getfem-users@nongnu.org>
Subject: [EXT] Re: field output size
WARNING: This email originated from outside of MD Anderson. Please validate the
sender's email address before clicking on links or attachments as they may not
be safe.
Dear Anne-Cecilie
You can just save the nodal values to text files, e.g.
numpy.save_txt("brain_displacements.txt",md.variable("ub"))
numpy.save_txt("brain_pressure.txt",md.variable("pb"))
and so on.
If you want to know the coordinates of each degree of freedom, you can also
save these in files
numpy.save_txt("brain_displacements_nodes.txt", mfub.basic_dof_nodes())
numpy.save_txt("brain_pressure_nodes.txt", mfpb.basic_dof_nodes())
For displacements and node coordinates you can also reshape your arrays to
having 3 columns before exporting.
I hope it helps.
Best regards
Kostas
On Tue, Jan 18, 2022 at 10:56 PM Lesage,Anne Cecile J
<ajles...@mdanderson.org<mailto:ajles...@mdanderson.org>> wrote:
Dear all
I am using python scripting with getfem
I have a mesh with 54026
I output the results in a vtu format with the following script lines
mfoutb = gf.MeshFem(meshb)
mfoutb.set_classical_discontinuous_fem(2)
mfoutb.export_to_vtu("PoroEResect25brain_%i.vtu" % it,
mfub, md.variable("ub"), "Displacements",
mfpb, md.variable("pb"), "Pressure")
when I save the spreadsheet of the results with paraview, I can the results on
2646060
Ideally I would like to have the results only on the 54026 nodes
Is there a way to correct my definition of mfoutb that corrects that?
I am preforming an inverse method on the material parameter typically on 700+
simulations, I need my results files to be smaller
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.
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.
