I have a TransientNonlinearImplicitSystem that I need help solving. I
solution technique I am trying to use is the multistep predictor corrector
method, which is summarized as follows.
(1) The predictor step which initializes by unknown h and the time rate of
change of h, h_dot.
(2) Solve Mx = R for x. Note, R is a function of both h and h_dot from
(3) The corrector step updates h and h_dot based on x that was solved in
step (2). Then, repeat staring at (2) until convergence.
If anyone is familiar with this method and may know how to solve it with
libmesh I would greatly appreciate the help and would be happy to provide