Dear myself,
On Thu, 13 Nov 2008, Tim Kroeger wrote:
> *system.solution += *system.old_local_solution;
> system.update();
>
> Unfortunately, this doesn't work in parallel since solution is a
> parallel vector and old_local_solution is a serial vector.
I have now found the following solution that
Dear all,
For a transient linear system in my application, the equation is quite
complicted, and it seemed easier to me to solve for (u_{n+1}-u_n) than
for u_{n+1} (that is, I solve for the update rather than for the
solution). To get the new solution, I then do:
*system.solution += *system.o
