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

[Adrian Croucher]

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.

Cheers, Adrian

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140416/99b8f742/attachment.html>