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 etc. As Jed says it doesn't make sense to save "partial" solutions within the linear solver. Just save it every several time-steps. Barry Also the code should not be "crashing" at seemingly long times (after hours) with a Segmentation fault. Send us the full error message and we'll see if there is some way we can help you fix this problem. > On Feb 18, 2019, at 9:47 AM, Sal Am via petsc-users <petsc-users@mcs.anl.gov> > wrote: > > 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 the results until that point seems to have been lost. > > > >
