<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Apr 14, 2014 at 8:44 PM, Adrian Croucher <span dir="ltr"><<a href="mailto:a.croucher@auckland.ac.nz" target="_blank">a.croucher@auckland.ac.nz</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">hi<br>
<br>
I just discovered that DMPlexInterpolate() complains if I try to interpolate a mesh containing 6-node wedge cells in 3-D:<br>
<br>
[0]PETSC ERROR: Argument out of range<br>
[0]PETSC ERROR: Cone size 6 not supported for dimension 3<br>
<br>
Looking in the code, it appears that DMPlexGetRawFaces_Internal() only supports cells with 3, 4 or 8 nodes in 3-D.<br>
<br>
At first I thought it would be simple enough to add another case in there for cone size 6. But then I noticed that DMPlexInterpolateFaces_<u></u>Internal() will throw an error for any type of cell with faces that don't all have the same number of nodes - and there's a comment in there saying "I see no way to do this if we admit faces of different shapes".<br>

<br>
I don't know the routine well enough to understand what the difficulty is, but are there any plans to generalise DMPlexInterpolate() so that it can handle commonly-used cell types (like 6-node wedge) that have faces with different numbers of nodes? We use wedge cells quite a bit, so it would limit the usefulness of DMPlex if it can't support them. Unless there is some kind of workaround?<br>
</blockquote><div><br></div><div>There is certainly a workaround. You give me the interpolated mesh. The only thing this code does is</div><div>allow you to input just cells and vertices. Since you are using these elements now, you know the faces,</div>
<div>so just use something like DMPlexCreateFromDAG() to input the fully interpolated mesh.</div><div><br></div><div>The problem here is that if I allow faces of different shapes, I have to keep track of the shape of each</div>
<div>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</div><div>before I get to it. If you want to try doing it, I would help.</div><div><br></div><div>  Thanks,</div>
<div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Cheers, Adrian<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Dr Adrian Croucher<br>
Department of Engineering Science<br>
University of Auckland<br>
New Zealand<br>
tel 64-9-373-7599 ext 84611<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>