DA with imposed parallel decomposition

[Matthew Knepley]

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>