<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div>Hi, I have a solver for time dependent problem and in each forward time step, I do </div><br><div>//<em> after update values in Amat and Pmat  /</em> /</div><div>KSPSetOperators(ksp, Amat, Pmat)</div><div>KSPSolve(ksp, rhs, solution)</div><br><div>The program runs and gives correct solution. But I got following memory leak errors from valgrind. </div><div>Do I need to create a new ksp in order to avoid the memory leak? Or is it due to the fieldsplit preconditioner I used?</div><br><div>==1340== 675,330 (1,232 direct, 674,098 indirect) bytes in 7 blocks are definitely lost in loss record 3,123 of 3,155</div><div>==1340==    at 0x4C320A6: memalign (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)</div><div>==1340==    by 0x4FDDA26: PetscMallocAlign (mal.c:49)</div><div>==1340==    by 0x4FDF3AD: PetscMallocA (mal.c:422)</div><div>==1340==    by 0x6564982: PCCreate_LU (lu.c:241)</div><div>==1340==    by 0x66BD26A: PCSetType (pcset.c:84)</div><div>==1340==    by 0x66BDDF8: PCSetFromOptions (pcset.c:154)</div><div>==1340==    by 0x680DAD0: KSPSetFromOptions (itcl.c:335)</div><div>==1340==    by 0x6606EA7: PCSetUp_FieldSplit (fieldsplit.c:944)</div><div>==1340==    by 0x66B29C4: PCSetUp (precon.c:894)</div><div>==1340==    by 0x68174DF: KSPSetUp (itfunc.c:376)</div><div>==1340==    by 0x681B07D: KSPSolve (itfunc.c:703)</div><div>==1340==    by 0x118F87: stokesSolver::forwardStep() (transient.cpp:1256)</div><div>==1340== </div><div>==1340== 851,586 (6,840 direct, 844,746 indirect) bytes in 9 blocks are definitely lost in loss record 3,126 of 3,155</div><div>==1340==    at 0x4C320A6: memalign (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)</div><div>==1340==    by 0x4FDDA26: PetscMallocAlign (mal.c:49)</div><div>==1340==    by 0x4FDF3AD: PetscMallocA (mal.c:422)</div><div>==1340==    by 0x55C6C46: ISCreate (isreg.c:37)</div><div>==1340==    by 0x55B5983: ISCreateGeneral (general.c:627)</div><div>==1340==    by 0x55B1C01: ISInvertPermutation_General (general.c:174)</div><div>==1340==    by 0x55BD8D5: ISInvertPermutation (index.c:481)</div><div>==1340==    by 0x5C275DC: MatLUFactorSymbolic_SeqAIJ (aijfact.c:296)</div><div>==1340==    by 0x57824CF: MatLUFactorSymbolic (matrix.c:3005)</div><div>==1340==    by 0x656311C: PCSetUp_LU (lu.c:109)</div><div>==1340==    by 0x66B29C4: PCSetUp (precon.c:894)</div><div>==1340==    by 0x68174DF: KSPSetUp (itfunc.c:376)</div><div>==1340== </div><div>==1340== 979,830 (25,488 direct, 954,342 indirect) bytes in 9 blocks are definitely lost in loss record 3,129 of 3,155</div><div>==1340==    at 0x4C320A6: memalign (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)</div><div>==1340==    by 0x4FDDA26: PetscMallocAlign (mal.c:49)</div><div>==1340==    by 0x4FDF3AD: PetscMallocA (mal.c:422)</div><div>==1340==    by 0x5DC7F7D: MatCreate (gcreate.c:83)</div><div>==1340==    by 0x5C24F29: MatGetFactor_seqaij_petsc (aijfact.c:101)</div><div>==1340==    by 0x5797454: MatGetFactor (matrix.c:4426)</div><div>==1340==    by 0x656308C: PCSetUp_LU (lu.c:108)</div><div>==1340==    by 0x66B29C4: PCSetUp (precon.c:894)</div><div>==1340==    by 0x68174DF: KSPSetUp (itfunc.c:376)</div><div>==1340==    by 0x681B07D: KSPSolve (itfunc.c:703)</div><div>==1340==    by 0x660BDA1: PCApply_FieldSplit_Schur (fieldsplit.c:1180)</div><div>==1340==    by 0x66AC096: PCApply (precon.c:444)</div><div>==1340== </div><div>==1340== 1,070,621 (5,320 direct, 1,065,301 indirect) bytes in 7 blocks are definitely lost in loss record 3,134 of 3,155</div><div>==1340==    at 0x4C320A6: memalign (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)</div><div>==1340==    by 0x4FDDA26: PetscMallocAlign (mal.c:49)</div><div>==1340==    by 0x4FDF3AD: PetscMallocA (mal.c:422)</div><div>==1340==    by 0x55C6C46: ISCreate (isreg.c:37)</div><div>==1340==    by 0x55B5983: ISCreateGeneral (general.c:627)</div><div>==1340==    by 0x5DA3869: MatGetOrdering_ND (spnd.c:27)</div><div>==1340==    by 0x5DA7712: MatGetOrdering (sorder.c:247)</div><div>==1340==    by 0x6562E19: PCSetUp_LU (lu.c:98)</div><div>==1340==    by 0x66B29C4: PCSetUp (precon.c:894)</div><div>==1340==    by 0x68174DF: KSPSetUp (itfunc.c:376)</div><div>==1340==    by 0x681B07D: KSPSolve (itfunc.c:703)</div><div>==1340==    by 0x660BDA1: PCApply_FieldSplit_Schur (fieldsplit.c:1180)</div><br><img class="mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border:0; width:0; height:0;" src="https://link.getmailspring.com/open/1CACEAE1-56BD-451A-8915-2D16DF1BB565@getmailspring.com?me=05842e25&recipient=cGV0c2MtdXNlcnNAbWNzLmFubC5nb3Y%3D">