[petsc-users] DMPlex higher order elements

[Matthew Knepley]

On Wed, Aug 24, 2016 at 12:12 PM, Jed Brown <jed at jedbrown.org> wrote:

> Matthew Knepley <knepley at gmail.com> writes:
> >> 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.
>
> The context is clearly that the mesh generator needs to express the
> curved elements.  DMPlex doesn't have a geometric model available, so it
> doesn't know how to make the Tri3 elements curve to conform more
> accurately to the boundary.  The mesh generator has no business telling
> you what function space to use for your solution, but it'd be a shame to
> prevent it from expressing element geometry.
>

Oh, so you want me to accept this Tri6 format for specifying quadratic
triangle
geometry. I guess I could stomach that, but I really hate that you would
mix the
topological information with the function space information, so that I have
to peel
them apart. They should just have a separate thing that indicates a
function space
and values for the geometry instead of screwing up a perfectly good
topological
definition.

   Matt

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