<div class="gmail_quote">On Tue, Apr 19, 2011 at 16:56, Ethan Coon <span dir="ltr"><<a href="mailto:ecoon@lanl.gov">ecoon@lanl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
This is in 3.1 I assume?</blockquote></div><br><div>Yeah, it's a problem in petsc-dev too (only tested with all periodic BCs). I get a seg-fault in</div><div><br></div><div><div>[0]PETSC ERROR: [0] VecScatterBegin_1 line 22 /home/jed/petsc/include/../src/vec/vec/utils/vpscat.h</div>
<div>[0]PETSC ERROR: [0] VecScatterBegin line 1546 /home/jed/petsc/src/vec/vec/utils/vscat.c</div><div>[0]PETSC ERROR: [0] DMGlobalToLocalBegin_DA line 15 /home/jed/petsc/src/dm/impls/da/dagtol.c</div><div>[0]PETSC ERROR: [0] DMGlobalToLocalBegin line 630 /home/jed/petsc/src/dm/interface/dm.c</div>
<div>[0]PETSC ERROR: [0] pci_function line 58 "unknowndirectory/"divclean.cpp</div><div>[0]PETSC ERROR: [0] SNES user function line 0 unknownunknown</div><div>[0]PETSC ERROR: [0] SNESComputeFunction line 1152 /home/jed/petsc/src/snes/interface/snes.c</div>
<div>[0]PETSC ERROR: [0] SNESSolve_LS line 142 /home/jed/petsc/src/snes/impls/ls/ls.c</div><div>[0]PETSC ERROR: [0] SNESSolve line 2385 /home/jed/petsc/src/snes/interface/snes.c</div></div><div><br></div><div><br></div><div>
earlier valgrind issues</div><div><br></div><div><div>==14160== Invalid read of size 2</div><div>==14160==    at 0x4C26C2C: memcpy (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)</div><div>==14160==    by 0x7589D01: opal_convertor_pack (in /usr/lib/openmpi/libmpi.so.1.0.1)</div>
<div>==14160==    by 0x114B90C8: ??? (in /usr/lib/openmpi/openmpi/mca_btl_sm.so)</div><div>==14160==    by 0x10C9B1B7: ??? (in /usr/lib/openmpi/openmpi/mca_pml_ob1.so)</div><div>==14160==    by 0x10C89C0C: ??? (in /usr/lib/openmpi/openmpi/mca_pml_ob1.so)</div>
<div>==14160==    by 0x74D8650: PMPI_Isend (in /usr/lib/openmpi/libmpi.so.1.0.1)</div><div>==14160==    by 0x700B103: PetscGatherMessageLengths (mpimesg.c:127)</div><div>==14160==    by 0x6D2B90A: VecScatterCreate_PtoS (vpscat.c:1623)</div>
<div>==14160==    by 0x6D074FD: VecScatterCreate (vscat.c:1334)</div><div>==14160==    by 0x5C5939A: DMSetUp_DA_3D (da3.c:953)</div><div>==14160==    by 0x5CAA4BF: DMSetUp_DA (dareg.c:22)</div><div>==14160==    by 0x5C1CA6C: DMSetUp (dm.c:189)</div>
<div>==14160==  Address 0x12addbb2 is 2 bytes inside a block of size 16 alloc'd</div><div>==14160==    at 0x4C23918: memalign (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)</div><div>==14160==    by 0x707ADEF: PetscMallocAlign (mal.c:30)</div>
<div>==14160==    by 0x6D2B5EA: VecScatterCreate_PtoS (vpscat.c:1603)</div><div>==14160==    by 0x6D074FD: VecScatterCreate (vscat.c:1334)</div><div>==14160==    by 0x5C5939A: DMSetUp_DA_3D (da3.c:953)</div><div>==14160==    by 0x5CAA4BF: DMSetUp_DA (dareg.c:22)</div>
<div>==14160==    by 0x5C1CA6C: DMSetUp (dm.c:189)</div><div>==14160==    by 0x5C5CD8D: DMDACreate3d (da3.c:1366)</div><div>==14160==    by 0x406608: multi_init(int, char**) (multi_petsc.cpp:107)</div><div>==14160==    by 0x403A67: main (mhd_main.cpp:351)</div>
</div><div><br></div><div><div>==14160== Invalid read of size 8</div><div>==14160==    at 0x6D15283: Pack_1 (vpscat.c:498)</div><div>==14160==    by 0x6D1BFC2: VecScatterBegin_1 (vpscat.h:100)</div><div>==14160==    by 0x6D0A483: VecScatterBegin (vscat.c:1575)</div>
<div>==14160==    by 0x5C62E24: DMGlobalToLocalBegin_DA (dagtol.c:19)</div><div>==14160==    by 0x5C1FFE3: DMGlobalToLocalBegin (dm.c:631)</div><div>==14160==    by 0x40B408: pci_function(_p_SNES*, _p_Vec*, _p_Vec*, void*) (divclean.cpp:61)</div>
<div>==14160==    by 0x4E8D9B5: SNESComputeFunction (snes.c:1162)</div><div>==14160==    by 0x4E74BB1: SNESSolve_LS (ls.c:158)</div><div>==14160==    by 0x4E9A27C: SNESSolve (snes.c:2407)</div><div>==14160==    by 0x40AB6F: divclean() (divclean.cpp:27)</div>
<div>==14160==    by 0x403AFD: main (mhd_main.cpp:365)</div><div>==14160==  Address 0x12af1570 is 0 bytes after a block of size 800 alloc'd</div><div>==14160==    at 0x4C23918: memalign (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)</div>
<div>==14160==    by 0x707ADEF: PetscMallocAlign (mal.c:30)</div><div>==14160==    by 0x6CD46A4: VecCreate_MPI_Private (pbvec.c:409)</div><div>==14160==    by 0x6CD4BF7: VecCreate_MPI (pbvec.c:453)</div><div>==14160==    by 0x6C8FB00: VecSetType (vecreg.c:53)</div>
<div>==14160==    by 0x6CD4F46: VecCreate_Standard (pbvec.c:482)</div><div>==14160==    by 0x6C8FB00: VecSetType (vecreg.c:53)</div><div>==14160==    by 0x5C738AB: DMCreateGlobalVector_DA (dadist.c:34)</div><div>==14160==    by 0x5C1D2D4: DMCreateGlobalVector (dm.c:277)</div>
<div>==14160==    by 0x5C732FE: VecDuplicate_MPI_DA (dadist.c:17)</div><div>==14160==    by 0x6C81665: VecDuplicate (vector.c:485)</div><div>==14160==    by 0x40C565: pci_init() (divclean.cpp:177)</div></div><div><br></div>