Re: [Libmesh-users] AMR question

[Roy S. <roy...@ic...>]

On Thu, 15 Apr 2010, Boyce Griffith wrote:

>
>
> On 4/15/10 11:27 AM, Roy Stogner wrote:
>> What version and configuration of PETSc?
>
> PETSc 3.0.0-p9 running on OS X 10.5 using gcc 4.4.2 and OpenMPI 1.3.3.
>
> PETSc was compiled at optimization level -O3 using --with-debugging=0.

So the problem is OS, compiler, and MPI implementation independent.
Still a slim chance it's PETSc (I've only reproduced it on 3.0.0,
couldn't get it to crop up with a brief test using 3.1) but I'd bet
it's us.

>> I'd appreciate it. Also if you do manage to reproduce the problem in
>> a simple test case, let me know. ex21 is currently hard-coded to use
>> a 3D L-shaped domain, the error doesn't occur until after a few
>> refinement steps, and by that time the size of the send_list (where
>> the problem mostly likely is) is up in the thousands and hard to
>> examine manually.
>
> Where is the relevant send_list --- within EquationSystems::reinit()?

It's constructed by a DofMap function that gets called from there, and
it's used by the failing localize() that also gets called from in
there.

> This error is cropping up after several refinement steps for me too (usually 
> many), but the total number of nodes in the problem is not too huge, so it 
> might be a little easier to diagnose.

Thanks.  I'm going to add some more parameters to ex21.in, see if this
still crops up in 2D or 1D.  It was a low priority when I thought it
might just be a bug in my PETSc install, but if it's not?  I'd been
hoping to release a libMesh 0.6.5 soon and this needs to be fixed
first.
---
Roy