On Wed, Sep 23, 2009 at 7:47 AM, Dave May <span dir="ltr">&lt;<a href="mailto:dave.mayhem23@gmail.com">dave.mayhem23@gmail.com</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hey Matt,<br>    In one piece of code I have, yes the call to DAGetInterpolation() does seem to cause a nasty crash.<br>It doesn&#39;t occur all the time, just with certain processor sizes (64) and certain mesh sizes (80x80x40). I was wondering if there was some pathological cases I did not know about. <br>

<br>I think I will have to write a stand alone test case to see if I can reproduce the error in a simpler code. <br>I don&#39;t think what I&#39;m doing should cause a problem, but I&#39;m not sure how best to debug the problem I have.<br>

<br>Any hints would be appreciated.  :)<br></blockquote><div><br>I guess it might be possible for you to specify a partition that breaks it, but I have a hard<br>time envisioning it. I would print out the DA sizes (happens with DAView()) for each level<br>
every time.<br><br>  Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Cheers,<br><font color="#888888">  Dave</font><div><div>
</div><div class="h5"><br><br><br><br><div class="gmail_quote">On Wed, Sep 23, 2009 at 2:33 PM, Matthew Knepley <span dir="ltr">&lt;<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Since it just inserts a point on every edge and face (in 2D), I do not see why it would fail. Does it?<br>

<br>  Matt<div><div></div><div><br><br><div class="gmail_quote">On Wed, Sep 23, 2009 at 12:04 AM, Dave May <span dir="ltr">&lt;<a href="mailto:dave.mayhem23@gmail.com" target="_blank">dave.mayhem23@gmail.com</a>&gt;</span> wrote:<br>


<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hello,<br>  Suppose I have a DA and I enforce the parallel decomposition during creation by specifying the arrays lx[], ly[], lz[] in DACreate3d(). If I now create a second DA using DARefine(), am I alays also able to obtain an interpolation operator between the two DA&#39;s via DAGetInterpolation()?<br>



<br>Under what circumstance will DAGetInterpolation() fail  when used between DA&#39;s generated in this manner?<br><br>Cheers,<br><font color="#888888">  Dave<br><br>
</font></blockquote></div><br><br clear="all"><br></div></div><font color="#888888">-- <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<br>

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