Add_cell with quadratic element

10 months ago by

I developed a script that allows to pass from a geometry made with Salome ( into a mesh for Fenics. However, this program only works for linear elements. I would like update my script in order to make a quadratic mesh. I use add_cell function from MeshEditor for linear triangle, linear edge, and linear tetrahedron. I tried to use this function for tetrahedrons with ten nodes using this line : add_cell(0,0,1,2,3,4,5,6,7,8,9).  I have an error when i compile my program because there  are too many arguments. Can you help me? is there another command that allows to create this type of cell?

Thank you in advance,

Best regards,
Community: FEniCS Project

1 Answer

10 months ago by
Take a look at the code for UnitDiscMesh.cpp in dolfin library - you need to create the vertex points first, then attach the midpoint data later.
Please login to add an answer/comment or follow this question.

Similar posts:
Search »