[petsc-dev] DACreate3d should work for 1 node in Z

[Ethan Coon]

This is in 3.1 I assume?

Ethan

On Tue, 2011-04-19 at 16:13 +0200, Jed Brown wrote:
> This would make it easier to solve 2D problems with a 3D code. It
> currently crashes or produces NaN because the scatter is set up
> incorrectly:
> 
> 
> ==10402== Invalid read of size 2
> ==10402==    at 0x4C26C2C: memcpy
> (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
> ==10402==    by 0x6D87D01: opal_convertor_pack
> (in /usr/lib/openmpi/libmpi.so.1.0.1)
> ==10402==    by 0xF33D0C8: ???
> (in /usr/lib/openmpi/openmpi/mca_btl_sm.so)
> ==10402==    by 0xEB1F1B7: ???
> (in /usr/lib/openmpi/openmpi/mca_pml_ob1.so)
> ==10402==    by 0xEB0DC0C: ???
> (in /usr/lib/openmpi/openmpi/mca_pml_ob1.so)
> ==10402==    by 0x6CD6650: PMPI_Isend
> (in /usr/lib/openmpi/libmpi.so.1.0.1)
> ==10402==    by 0x551F85A: PetscGatherMessageLengths (mpimesg.c:128)
> ==10402==    by 0x580D5AC: VecScatterCreate_PtoS (vpscat.c:1621)
> ==10402==    by 0x5823C9F: VecScatterCreate (vscat.c:1316)
> ==10402==    by 0x518A547: DACreate_3D (da3.c:957)
> ==10402==    by 0x51D1ACA: DASetType (dareg.c:48)
> ==10402==    by 0x5194FA7: DASetTypeFromOptions_Private
> (dacreate.c:65)
> ==10402==  Address 0xbc33562 is 498 bytes inside a block of size 516
> alloc'd
> ==10402==    at 0x4C23918: memalign
> (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
> ==10402==    by 0x53E3D3B: PetscMallocAlign (mal.c:30)
> ==10402==    by 0x5566432: PetscTrMallocDefault (mtr.c:192)
> ==10402==    by 0x580D2B4: VecScatterCreate_PtoS (vpscat.c:1601)
> ==10402==    by 0x5823C9F: VecScatterCreate (vscat.c:1316)
> ==10402==    by 0x518A547: DACreate_3D (da3.c:957)
> ==10402==    by 0x51D1ACA: DASetType (dareg.c:48)
> ==10402==    by 0x5194FA7: DASetTypeFromOptions_Private
> (dacreate.c:65)
> ==10402==    by 0x5195A81: DASetFromOptions (dacreate.c:131)
> ==10402==    by 0x518FFF5: DACreate3d (da3.c:1742)
> ==10402==    by 0x40669A: multi_init(int, char**)
> (multi_petsc.cpp:115)
> ==10402==    by 0x403AF5: main (mhd_main.cpp:351)
> 
> 

-- 
------------------------------------
Ethan Coon
Post-Doctoral Researcher
Applied Mathematics - T-5
Los Alamos National Laboratory
505-665-8289

http://www.ldeo.columbia.edu/~ecoon/
------------------------------------