On Wed, Feb 20, 2019 at 4:43 AM Sal Am wrote:
> Hi Matthew you were right,
>
> The matrix I have is very ill conditioned and my supervisor gave it for
> testing purposes. Having said that, I was able to solve it previously
> however, for some reason it said convergence reached at e-3 even though
I'm guessing you mange your own time stepping (and nonlinear solver if there
is one)?
You can save the solution with a call to VecView() and then reload the
solution with a VecLoad() but you need to manage any other restart data that
you may need, like the value of the current time
What kind of problems are you solving? Running for a Krylov method for
tens of thousands of iterations is very rarely recommended.
Regarding storage, it's significantly more expensive to store the Krylov
basis (even when it's a recurrence) than the current approximation.
Some methods require
What kind of solver are you using and how often do you want to write?
Sal Am via petsc-users writes:
> Is there a function/command line option to save the solution as it is
> solving (and read in the file from where it crashed and keep iterating from
> there perhaps)?
> Had a seg fault and all
