[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.
>

Yes, I do not support that since I think its a crazy way to talk about
things. All the topological information is in the Tri3 mesh, and
Cubit has no business telling me about the function space.

   Matt


> 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
>



-- 
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/d9efc71f/attachment.html>