<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, 9 May 2020 at 23:44, Mohammed Ashour <<a href="mailto:ashour.msc@gmail.com">ashour.msc@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I'm using TSALPHA1, as in ierr = TSSetType(ts,TSALPHA1); CHKERRQ(ierr);<div><br></div></div></blockquote><div><br></div><div>Well, TSALPHA1 uses some vectors keeping intenal state, as you know, generalized-\alpha is not really a one-step method. That intenal state should be saved and loaded if you want to implement a truly correct checkpoint/restart. But PETSc does not currently provide access to this vector in its public API, so you will have to hack things around. IIRC, just saving and restoring the internal V1 vector is all what you need, maybe also set ts->steprestart to FALSE and properly set the ts->steps counter.</div><div><br></div><div>Or you can just ignore these detail, and restart the easy way from just the saved solution, at the price of a very small "glitch" at the restart time. But note that this glitch just means that you temporarily change the scheme to reinitialize, and it is exactly the same thing the implementation does at the initial time to initialize the method (as you surely know, the method is also not trully self-starting). The restarting procesure as implemented in TSALPHA1 is not in the literature, it is of my own cooking but based on rather trivial relations performing two steps with dt/2 time step size, the fist of those inner steps using backward-Euler.</div><div><br></div><div><br></div></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div>Lisandro Dalcin<br>============<br>Research Scientist<br>Extreme Computing Research Center (ECRC)<br>King Abdullah University of Science and Technology (KAUST)<br><a href="http://ecrc.kaust.edu.sa/" target="_blank">http://ecrc.kaust.edu.sa/</a><br></div></div></div></div>