<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Jan 25, 2017 at 11:12 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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
  

    
  
  <div bgcolor="#FFFFFF">
    <p>hi Matt,</p>
    <p>A couple of years ago we discussed adding support for 6-node
      wedge cells in DMPlex- the hurdle being that at present DMPlex
      doesn't support cells with faces of different shapes (see below)-
      and you said it could be done, if and when I needed it.<br>
    </p>
    <p>People are starting to want to use my code on real models now,
      many of which have 6-node wedge cells. So it would be very handy
      to have support for those in DMPlex fairly soon if possible.</p>
    <p>Any chance you would be able to take a look at it? In principle I
      would be happy to help if I can. I have started looking at
      plexinterpolate.c to try and see how it works, but at this stage
      I'm not at all confident that I understand the subtleties.</p></div></blockquote><div>Okay. First I need to understand exactly what shape we want to handle. Can you draw it or explain precisely?</div><div><br></div><div>Second, we need to prescribe a total order on the vertices for a wedge, and then a set of faces as ordered arrays of vertices.</div><div>That information goes here:</div><div><br></div><div>  <a href="https://bitbucket.org/petsc/petsc/src/e8e596daf24561903cea0839b7b68a432a6283e4/src/dm/impls/plex/plexinterpolate.c?at=master&fileviewer=file-view-default#plexinterpolate.c-37">https://bitbucket.org/petsc/petsc/src/e8e596daf24561903cea0839b7b68a432a6283e4/src/dm/impls/plex/plexinterpolate.c?at=master&fileviewer=file-view-default#plexinterpolate.c-37</a></div><div><br></div><div>I think that might be all we have to do, but until you try something, secret assumptions lurking in the code are not obvious.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div bgcolor="#FFFFFF">
    <p>Cheers, Adrian<br>
    </p>
    <p><br>
      On 16/04/14 08:42, Matthew Knepley wrote:</p>
    <p>
      </p><blockquote type="cite">
        <p> </p>
        <blockquote type="cite"> <br>
          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.<br>
          <blockquote type="cite">
            <div dir="ltr">
              <div class="gmail_extra">
                <div class="gmail_quote">
                  <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>
              </div>
            </div>
          </blockquote>
          <br>
          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.</blockquote>
        <br>
        Okay, when you need it, we will make it work.<br>
        <br>
      </blockquote>
      <br>
    <p></p>
    --
    <pre class="gmail-m_3547875481669262486moz-signature" cols="72">Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: <a class="gmail-m_3547875481669262486moz-txt-link-abbreviated" href="mailto:a.croucher@auckland.ac.nz" target="_blank">a.croucher@auckland.ac.nz</a>
tel: <a href="tel:+64%209-923%204611" value="+6499234611" target="_blank">+64 (0)9 923 4611</a>
</pre>
  </div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>