<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 20 mar 2019 alle ore 23:40 Derek Gaston via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Trying to track down some memory corruption I'm seeing on larger scale runs (3.5B+ unknowns).</div></div></div></div></div></div></blockquote><div><br></div><div>Uhm.... are you using 32bit indices? is it possible there's integer overflow somewhere?</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Was able to run Valgrind on it... and I'm seeing quite a lot of uninitialized value errors coming from ghost updating.  Here are some of the traces:</div><div><br></div><div>==87695== Conditional jump or move depends on uninitialised value(s)<br>==87695==    at 0x73236D3: PetscMallocAlign (mal.c:28)<br>==87695==    by 0x7323C70: PetscMallocA (mal.c:390)<br>==87695==    by 0x739048E: VecScatterMemcpyPlanCreate_Index (vscat.c:284)<br>==87695==    by 0x73A5D97: VecScatterMemcpyPlanCreate_PtoP (vpscat_mpi1.c:312)<br>==64730==    by 0x7393E8A: VecScatterSetUp_vectype_private (vscat.c:857)<br>==64730==    by 0x7395E5D: VecScatterSetUp_MPI1 (vpscat_mpi1.c:2543)<br>==64730==    by 0x73DDD39: VecScatterSetUp (vscatfce.c:212)<br>==64730==    by 0x73DCD73: VecScatterCreateWithData (vscreate.c:333)<br>==64730==    by 0x7444232: VecCreateGhostWithArray (pbvec.c:685)<br>==64730==    by 0x744490D: VecCreateGhost (pbvec.c:741)<br></div><div><br></div><div>==133582== Conditional jump or move depends on uninitialised value(s)<br>==133582==    at 0x4030384: memcpy@@GLIBC_2.14 (vg_replace_strmem.c:1034)<br>==133582==    by 0x739E4F9: PetscMemcpy (petscsys.h:1649)<br>==133582==    by 0x739E4F9: VecScatterMemcpyPlanExecute_Pack (vecscatterimpl.h:150)<br>==133582==    by 0x739E4F9: VecScatterBeginMPI1_1 (vpscat_mpi1.h:69)<br>==133582==    by 0x73DD964: VecScatterBegin (vscatfce.c:110)<br>==133582==    by 0x744E195: VecGhostUpdateBegin (commonmpvec.c:225)<br></div><div><br></div><div>This is from a Git checkout of PETSc... the hash I branched from is: 0e667e8fea4aa from December 23rd (updating would be really hard at this point as I've completed 90% of my dissertation with this version... and changing PETSc now would be pretty painful!).</div><div><br></div><div>Any ideas?  Is it possible it's in my code?  Is it possible that there are later PETSc commits that already fix this?</div><div><br></div><div>Thanks for any help,</div><div>Derek<br></div><div><br></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Stefano</div></div>