
<html>


<head>


<meta http-equiv="Content-Type" content="text/html; charset=utf-8">


</head>


<body>


<div dir="ltr">


<div dir="ltr">Do you have a test example?<br clear="all">


<div>


<div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">


<div dir="ltr">--Junchao Zhang</div>


</div>


</div>


</div>


<br>


<div class="gmail_quote">


<div dir="ltr" class="gmail_attr">On Tue, Jan 14, 2020 at 4:44 AM Y. Shidi <<a href="mailto:ys453@cam.ac.uk">ys453@cam.ac.uk</a>> wrote:<br>


</div>


<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">


Dear developers,<br>


<br>


I have a 2x2 nested matrix and the corresponding nested vector.<br>


When I running the code with field splitting, it gets the following<br>


errors:<br>


<br>


[0]PETSC ERROR: PetscTrFreeDefault() called from VecRestoreArray_Nest() <br>


line 678 in /home/ys453/Sources/petsc/src/vec/vec/impls/nest/vecnest.c<br>


[0]PETSC ERROR: Block at address 0x3f95f60 is corrupted; cannot free;<br>


may be block not allocated with PetscMalloc()<br>


[0]PETSC ERROR: --------------------- Error Message <br>


--------------------------------------------------------------<br>


[0]PETSC ERROR: Memory corruption: <br>


<a href="http://www.mcs.anl.gov/petsc/documentation/installation.html#valgrind" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/documentation/installation.html#valgrind</a><br>


[0]PETSC ERROR: Bad location or corrupted memory<br>


[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">


http://www.mcs.anl.gov/petsc/documentation/faq.html</a> <br>


for trouble shooting.<br>


[0]PETSC ERROR: Petsc Release Version 3.9.3, unknown<br>


[0]PETSC ERROR: 2DPetscSpuriousTest on a arch-linux2-c-debug named <br>


merlin by ys453 Tue Jan 14 10:36:53 2020<br>


[0]PETSC ERROR: Configure options --download-scalapack --download-mumps <br>


--download-parmetis --download-metis --download-ptscotch <br>


--download-superlu_dist --download-hypre<br>


[0]PETSC ERROR: #1 PetscTrFreeDefault() line 269 in <br>


/home/ys453/Sources/petsc/src/sys/memory/mtr.c<br>


[0]PETSC ERROR: #2 VecRestoreArray_Nest() line 678 in <br>


/home/ys453/Sources/petsc/src/vec/vec/impls/nest/vecnest.c<br>


[0]PETSC ERROR: #3 VecRestoreArrayRead() line 1835 in <br>


/home/ys453/Sources/petsc/src/vec/vec/interface/rvector.c<br>


[0]PETSC ERROR: #4 VecRestoreArrayPair() line 511 in <br>


/home/ys453/Sources/petsc/include/petscvec.h<br>


[0]PETSC ERROR: #5 VecScatterBegin_SSToSS() line 671 in <br>


/home/ys453/Sources/petsc/src/vec/vscat/impls/vscat.c<br>


[0]PETSC ERROR: #6 VecScatterBegin() line 1779 in <br>


/home/ys453/Sources/petsc/src/vec/vscat/impls/vscat.c<br>


[0]PETSC ERROR: #7 PCApply_FieldSplit() line 1010 in <br>


/home/ys453/Sources/petsc/src/ksp/pc/impls/fieldsplit/fieldsplit.c<br>


[0]PETSC ERROR: #8 PCApply() line 457 in <br>


/home/ys453/Sources/petsc/src/ksp/pc/interface/precon.c<br>


[0]PETSC ERROR: #9 KSP_PCApply() line 276 in <br>


/home/ys453/Sources/petsc/include/petsc/private/kspimpl.h<br>


[0]PETSC ERROR: #10 KSPFGMRESCycle() line 166 in <br>


/home/ys453/Sources/petsc/src/ksp/ksp/impls/gmres/fgmres/fgmres.c<br>


[0]PETSC ERROR: #11 KSPSolve_FGMRES() line 291 in <br>


/home/ys453/Sources/petsc/src/ksp/ksp/impls/gmres/fgmres/fgmres.c<br>


[0]PETSC ERROR: #12 KSPSolve() line 669 in <br>


/home/ys453/Sources/petsc/src/ksp/ksp/interface/itfunc.c<br>


<br>


I am not sure why it happens.<br>


<br>


Thank you for your time.<br>


<br>


Kind Regards,<br>


Shidi<br>


</blockquote>


</div>


</div>


</body>


</html>