Sorry Barry, <br> I should have been more specific. I get a segmentation violation.<br>I think writing a simple test to isolate the bug is the best option to debug this error.<br>Do you agree?<br><br>Cheers,<br> Dave<br>
<br><br><div class="gmail_quote">On Wed, Sep 23, 2009 at 3:00 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</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;">
<div class="im"><br>
On Sep 23, 2009, at 7:47 AM, Dave May wrote:<br>
<br>
<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>
</blockquote>
<br></div>
A crash? Segmentation violation etc? Or an error message saying that the decomposition is not supported? They are very different things: a crash is a bug we need to deal with, unsupported is just because it is too hard to support refinement with all decompositions.<br>
<font color="#888888">
<br>
Barry</font><div><div></div><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
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>
<br>
Cheers,<br>
Dave<br>
<br>
<br>
<br>
On Wed, Sep 23, 2009 at 2:33 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br>
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<br>
<br>
<br>
On Wed, Sep 23, 2009 at 12:04 AM, Dave May <<a href="mailto:dave.mayhem23@gmail.com" target="_blank">dave.mayhem23@gmail.com</a>> wrote:<br>
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>
Dave<br>
<br>
<br>
<br>
<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>
<br>
</blockquote>
<br>
</div></div></blockquote></div><br>