
<html>

  <head>


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

  </head>

  <body text="#000000" bgcolor="#FFFFFF">

    Dear all,<br>

    <div class="moz-forward-container"> <br>

      I am trying to simulate a time-dependent PDE on a fairly large 2D

      domain. <br>

      <br>

      I use the BCGS_AMG_SSOR linear solver for parallel computations,

      in combination with Newton method,      <br>

          Dune::PDELab::ISTLBackend_BCGS_AMG_SSOR<IGO> LS;<br>

          LS ls(gfs,5000,0, <u><i>false</i></u>, true);<br>

      <br>

      I am using dune 2.4 version.<br>

      I noticed that if I do not reuse the matrix in LS then I

      experience massive memory leaks. For reference, I have 64 GB RAM

      and a swap of 128 GB, and even then I run out of memory within

      just a quarter way through my simulation!<br>

      But if I reuse the matrix (i.e., LS ls(gfs,5000,0, <u><i>true</i></u>,

      true)), I no longer get the memory leak issue, but the simulation

      becomes painfully slow. It starts out at normal speed, but gets

      progressively slower, to a point where it is essentially useless.<br>

      <br>

      Is this normal? I am wondering if anyone can suggest a fix or an

      alternative for my problem. <br>

      <br>

      Thanks in advance,<br>

      warm wishes, Shubhangi<br>

      <pre class="moz-signature" cols="72">-- 

Dr. Shubhangi Gupta

Marine Geosystems

GEOMAR Helmholtz Center for Ocean Research

Wischhofstraße 1-3, 

D-24148 Kiel


Room: 12-206

Phone: +49 431 600-1402

Email: <a class="moz-txt-link-abbreviated" href="mailto:sgupta@geomar.de" moz-do-not-send="true">sgupta@geomar.de</a></pre>

    </div>

  </body>

</html>