[petsc-dev] 6-node wedge cells in DMPlex

[Matthew Knepley]

On Tue, Apr 15, 2014 at 5:50 PM, Adrian Croucher
<a.croucher at auckland.ac.nz>wrote:

>
> On 16/04/14 06:06, Matthew Knepley wrote:
>
>   There is certainly a workaround. You give me the interpolated mesh. The
> only thing this code does is
> allow you to input just cells and vertices. Since you are using these
> elements now, you know the faces,
> so just use something like DMPlexCreateFromDAG() to input the fully
> interpolated mesh.
>
>
> Yeah, I was afraid that might be the workaround ;-) It seems to me that
> for a large unstructured mesh, creating the whole DAG would basically
> amount to duplicating what the very handy DMPlexInterpolate() function does.
>
>
>  The problem here is that if I allow faces of different shapes, I have to
> keep track of the shape of each
> face when I make it. This is possible,  and I have added it to the TODO
> list, but I think it will be a while
> before I get to it. If you want to try doing it, I would help.
>
>
> I'd have to get more familiar with how it works first, and it probably
> isn't top of my priority list either at present- I can continue getting my
> code together without support for wedge cells at the moment. But at some
> point I will need it.
>

Okay, when you need it, we will make it work.

  Thanks,

     Matt


> Cheers, Adrian
>
> --
> Dr Adrian Croucher
> Department of Engineering Science
> University of Auckland
> New Zealand
> tel 64-9-373-7599 ext 84611
>
>


-- 
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-dev/attachments/20140416/7b47d6ec/attachment.html>