[petsc-users] DMPlex higher order elements

[Andrew Ho]

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 <knepley at gmail.com> wrote:

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



-- 
Andrew Ho
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160824/0900c6b5/attachment.html>