[petsc-users] DMPlex higher order elements

[Andrew Ho]

I created an unstructured tri6 mesh in Cubit and am trying to read it into
PETSc using DMPlexCreateFromFile. However, when I do PETSc gives me an
error that it doesn't support this type of element.

I know my PETSc install has Exodus II support because by giving a different
Exodus mesh file with Tri3 elements works just fine. I've attached the
Exodus mesh generated by Cubit. I can open up the mesh file in Paraview
just fine and it does appear to be valid (it has 2 elements approximating a
quarter circle wedge).

Here's the error message:

[0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
>
> [0]PETSC ERROR: Argument out of range
> [0]PETSC ERROR: Cone size 6 not supported for dimension 2
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [0]PETSC ERROR: Petsc Development GIT revision: v3.7.3-1257-ge04de0c  GIT
> Date: 2016-08-24 00:14:37 -0500
>
> [0]PETSC ERROR: bin/mesh_testing on a arch-linux2-c-opt named ahome by
> andrew Tue Aug 23 22:48:57 2016
>
> [0]PETSC ERROR: Configure options --with-debugging=0 --COPTFLAGS="-O3
> -march=native" --CXXOPTFLAGS="-O3 -march=native" --FOPTFLAGS="-O3
> -march=native" --download-exodusii=yes --download-hdf5=yes
> --download-netcdf=yes
> [0]PETSC ERROR: #8 DMPlexGetRawFaces_Internal() line 93 in
> petsc/src/dm/impls/plex/plexinterpolate.c
>
> [0]PETSC ERROR: #9 DMPlexGetFaces_Internal() line 20 in
> petsc/src/dm/impls/plex/plexinterpolate.c
>
> [0]PETSC ERROR: #10 DMPlexInterpolateFaces_Internal() line 172 in
> petsc/src/dm/impls/plex/plexinterpolate.c
>
> [0]PETSC ERROR: #11 DMPlexInterpolate() line 532 in
> petsc/src/dm/impls/plex/plexinterpolate.c
>
> [0]PETSC ERROR: #12 DMPlexCreateExodus() line 168 in
> petsc/src/dm/impls/plex/plexexodusii.c
>
> [0]PETSC ERROR: #13 DMPlexCreateExodusFromFile() line 46 in
> petsc/src/dm/impls/plex/plexexodusii.c
>
> [0]PETSC ERROR: #14 DMPlexCreateFromFile() line 1967 in
> petsc/src/dm/impls/plex/plexcreate.c
>
> [0]PETSC ERROR: #15 main() line 11 in mesh_testing.cpp
>
> [0]PETSC ERROR: No PETSc Option Table entries
> [0]PETSC ERROR: ----------------End of Error Message -------send entire
> error message to petsc-maint at mcs.anl.gov----------
>


-- 
Andrew Ho
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160823/13947d68/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: tri6_wedge.exo
Type: application/octet-stream
Size: 1360 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160823/13947d68/attachment.obj>