Dear Anne Cecile Yes, I can confirm this. Calling set_classical_fem(2) makes tertrahedral elements into 10-node tetras, and hexahedral elements into 27-node hexas, independent of whether they have the same number of nodes in the mesh that you provide.
BR Kostas On Thu, Mar 24, 2022 at 4:56 PM Lesage,Anne Cecile J < ajles...@mdanderson.org> wrote: > Dear all > > > > I am writing poroelastic python script for two equations coupling pressure > and displacement > > If I write the following, does it mean that I am using a second order > tetrahedra approximation for my displacement? > > > > > ########################################################################################################## > > # > > # Definition of finite elements methods and integration method > > # > > > ######################################################################################################### > > > > press_fem_order = 1 # pressure finite element order > > disp_fem_order = 2 # displacements finite element order > > > > > > mfub = gf.MeshFem(meshb, 3) > > mfub.set_classical_fem(disp_fem_order) > > > > > > > > mimu3b=gf.MeshIm(meshb, 5) # integration displacement on tetrahedra brain > > > > > > mfpb= gf.MeshFem(meshb, 1) > > mfpb.set_classical_fem(press_fem_order) > 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. >