True, Paraview does draw them in as "flat" edges, but this is less of a concern for me; the important part is if I want to do a simulation on a curved geometric model, I need to be able to achieve higher order numerical accuracy because it reduces my simulation times by orders of magnitude.
On Wed, Aug 24, 2016 at 10:22 AM, Matthew Knepley <knep...@gmail.com> wrote: > On Wed, Aug 24, 2016 at 12:57 AM, Andrew Ho <andre...@uw.edu> wrote: > >> 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). >> > > Okay, I will look at getting this to read in alright. Notice that Paraview > does exactly the wrong thing here in that it has > straight lines connecting the midpoints and corners of the triangles. They > should be smooth quadratic curves, and this > misunderstanding comes from the horrible format choice which pretends that > there are "extra" vertices. > > Matt > > >> 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-ma...@mcs.anl.gov---------- >>> >> >> >> -- >> Andrew Ho >> > > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > -- Andrew Ho