Sorry Barry, I should have been more specific. I get a segmentation violation. I think writing a simple test to isolate the bug is the best option to debug this error. Do you agree?
Cheers, Dave On Wed, Sep 23, 2009 at 3:00 PM, Barry Smith <bsmith at mcs.anl.gov> wrote: > > On Sep 23, 2009, at 7:47 AM, Dave May wrote: > > Hey Matt, >> In one piece of code I have, yes the call to DAGetInterpolation() does >> seem to cause a nasty crash. >> > > 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. > > Barry > > > 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. >> >> I think I will have to write a stand alone test case to see if I can >> reproduce the error in a simpler code. >> 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. >> >> Any hints would be appreciated. :) >> >> Cheers, >> Dave >> >> >> >> On Wed, Sep 23, 2009 at 2:33 PM, Matthew Knepley <knepley at gmail.com> >> wrote: >> Since it just inserts a point on every edge and face (in 2D), I do not see >> why it would fail. Does it? >> >> Matt >> >> >> On Wed, Sep 23, 2009 at 12:04 AM, Dave May <dave.mayhem23 at gmail.com> >> wrote: >> Hello, >> 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()? >> >> Under what circumstance will DAGetInterpolation() fail when used between >> DA's generated in this manner? >> >> Cheers, >> Dave >> >> >> >> >> -- >> What most experimenters take for granted before they begin their >> experiments is infinitely more interesting than any results to which their >> experiments lead. >> -- Norbert Wiener >> >> > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090923/19c1ea6d/attachment.htm>
