Re: [deal.II] Higher order Tetrahedron support for gmsh
On 3/21/23 11:41, Kumar Saurabh wrote: Is there any documentation that I can take a look at to understand the order of the vertices in second order tetrahedron. I mean how the mapping is done to the reference element for second order tetrahedron. You mean what the order of nodes is that gmsh uses? The place for that would be the gmsh documentation. Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/f746c012-bd9a-ef0f-0487-a481bf81eadd%40colostate.edu.
Re: [deal.II] Higher order Tetrahedron support for gmsh
Hi Dr. Bangerth, Thanks for your message. Is there any documentation that I can take a look at to understand the order of the vertices in second order tetrahedron. I mean how the mapping is done to the reference element for second order tetrahedron. I just want to use volumetric integration for assembly with Dirichlet boundary computation. Thanks, Kumar Saurabh On Monday, March 20, 2023 at 1:28:32 PM UTC-7 Wolfgang Bangerth wrote: > On 3/20/23 12:26, Kumar Saurabh wrote: > > > > I wanted to check if deal.ii has a support for second order tetrahedron > > meshes. I tried to load the file from GMSH and it gave the error for > > unsupported element type. > > > > Is it that deal.ii do not support second order tetrahedron from GMSH or > > do not support second order tetrahedron at all? > > Kumar, > for triangles and tetrahedra, mappings are generally implemented by > relying in MappingFE. As a consequence, if there is a Lagrange element > of polynomial degree p, there is a corresponding mapping of polynomial > degree p. Right now, FE_SimplexP is available up to polynomial degree > p=2, so quadratic elements should be ok. > > Whether or not GridIn can read such meshes is a different question. If > you get an error, then the answer is apparently no. But that shouldn't > be too hard to fix, and we would be quite happy to accept a patch (and > point you in the right direction regarding what needs to be > implemented). My suspicion is that it would not take more then 100 lines > of code to read these meshes. > > The bigger problem with quadratic meshes (which also exists for the case > of quadratic quadrilaterals/hexahedra) is that we do not currently have > a good way to translate the information encoded by mid-edge nodes into a > mapping. If we can read such meshes for quadrilaterals, I suspect that > we just ignore the mid-edge nodes and rely on the user's ability to > provide a mapping that restores this kind of information. > > Best > W. > > > -- > > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/7b08cce9-2693-4373-920d-771f2411b0a9n%40googlegroups.com.
Re: [deal.II] Higher order Tetrahedron support for gmsh
On 3/20/23 12:26, Kumar Saurabh wrote: I wanted to check if deal.ii has a support for second order tetrahedron meshes. I tried to load the file from GMSH and it gave the error for unsupported element type. Is it that deal.ii do not support second order tetrahedron from GMSH or do not support second order tetrahedron at all? Kumar, for triangles and tetrahedra, mappings are generally implemented by relying in MappingFE. As a consequence, if there is a Lagrange element of polynomial degree p, there is a corresponding mapping of polynomial degree p. Right now, FE_SimplexP is available up to polynomial degree p=2, so quadratic elements should be ok. Whether or not GridIn can read such meshes is a different question. If you get an error, then the answer is apparently no. But that shouldn't be too hard to fix, and we would be quite happy to accept a patch (and point you in the right direction regarding what needs to be implemented). My suspicion is that it would not take more then 100 lines of code to read these meshes. The bigger problem with quadratic meshes (which also exists for the case of quadratic quadrilaterals/hexahedra) is that we do not currently have a good way to translate the information encoded by mid-edge nodes into a mapping. If we can read such meshes for quadrilaterals, I suspect that we just ignore the mid-edge nodes and rely on the user's ability to provide a mapping that restores this kind of information. Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/b9fae7fd-ca33-80be-9a0b-b638475276e4%40colostate.edu.
[deal.II] Higher order Tetrahedron support for gmsh
Hi, I wanted to check if deal.ii has a support for second order tetrahedron meshes. I tried to load the file from GMSH and it gave the error for unsupported element type. Is it that deal.ii do not support second order tetrahedron from GMSH or do not support second order tetrahedron at all? Thanks, Kumar Saurabh -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/c2ff0147-dfa8-4b0a-96a8-e4594f5c0274n%40googlegroups.com.
