<div class="gmail_quote">On Thu, Sep 22, 2011 at 17:14, Milan Mitrovic <span dir="ltr">&lt;<a href="mailto:milan.v.mitrovic@gmail.com">milan.v.mitrovic@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":22g">I use snes for an implicit time stepping scheme... I call it in a loop<br>
for every time step...<br>
<br>
the jacobian is recomputed at every step, and the preallocation can<br>
change from time to time... what is the proper way to do it then?<br>
destroy the matrix, then recreate it?</div></blockquote><div><br></div><div>Don&#39;t destroy the matrix, just call the preallocation routines when the nonzero pattern changes.</div><div><br></div><div>You can actually destroy the matrices from C, but I don&#39;t think there is a pure Fortran way to do it (it needs another pointer indirection) and you should never need to anyway.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":22g"> what should I do with the<br>
preconditioner matrix? (petsc complains if I dont do P = J)</div></blockquote></div><br><div>What do you mean &quot;complains&quot;? You can certainly use a different preconditioning matrix.</div>