On Wed, Sep 23, 2009 at 7:47 AM, Dave May <dave.mayhem23 at gmail.com> wrote:
> Hey Matt, > In one piece of code I have, yes the call to DAGetInterpolation() does > seem to cause a nasty crash. > 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. :) > I guess it might be possible for you to specify a partition that breaks it, but I have a hard time envisioning it. I would print out the DA sizes (happens with DAView()) for each level every time. Matt > 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 >> > > -- 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/4eb0bc66/attachment.htm>
