[petsc-users] DMPlex higher order elements

[Matthew Knepley]

On Wed, Aug 24, 2016 at 10:45 AM, Boyce Griffith <griffith at cims.nyu.edu>
wrote:

>
>
> On Aug 24, 2016, at 11:02 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.
>>
>
> 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.
>
>
> Do you support / plan to support curved elements?
>

I had "support" in there, but there were bugs. Toby and Mark discovered
these, and Toby has fixed them. I think
all of the fixes are in master now.

We currently support isoparametric elements, but it will not be hard (I
think) to support super/sub-parametrics
by attaching an independent FEM space to the coordinate DM. Right now it
just defaults to the one in the
master DM. I don't think there should be a problem with that, but of course
we need to try.

  Matt

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


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