[petsc-users] DMPlex higher order elements

[Matthew Knepley]

On Wed, Aug 24, 2016 at 12:57 AM, Andrew Ho <andrewh0 at 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-maint at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160824/d48e2dc9/attachment-0001.html>