<div dir="ltr">I am trying to replace the non-linear solver in a flow simulation problem where the matrix sparsity can change during the iterations. I tried successfully to create the matrix within the FormJacobian function but I have a memory leak which I don't know how to fix.<div><br></div><div>My FormJacobian function is something like this:</div><div><br></div><div>.. FormJacobian( ..){</div><div>     matrixconstruction(); // returns csr matrix in c++ regular vectors</div><div>     MatSeqAIJSetPreallocationCSR ...;</div><div>     //loop over the entries and </div><div>     MatSetValues;</div><div>     MatAssemblyBegin();</div><div>     MatAssemblyEnd();</div><div>     return 0;</div><div>}</div><div><br></div><div>This approach works for small problems but, when dealing with large ones I observe a memory leak. I am pretty sure that It is related to the fact that I am creating new matrices every time FormJacobian is called but I don't know how to destroy it between iterations.</div><div><br></div><div>Any suggestions?</div><div><br></div><div>Regards,</div><div><br></div><div>Adolfo</div></div>