On Wed, Sep 23, 2009 at 7:47 AM, Dave May <span dir="ltr"><<a href="mailto:dave.mayhem23@gmail.com">dave.mayhem23@gmail.com</a>></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'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't think what I'm doing should cause a problem, but I'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"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></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"><<a href="mailto:dave.mayhem23@gmail.com" target="_blank">dave.mayhem23@gmail.com</a>></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's via DAGetInterpolation()?<br>
<br>Under what circumstance will DAGetInterpolation() fail when used between DA'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>