Hello,
I want to implement a nonlinear-variational time-dependent problem.
For time stepping I currently do
// Time-stepping
dolfin::Progress p("Time-stepping");
while (t < T) {
// Assemble vector and apply boundary conditions
assemble(b, L);
// bc.apply(b);
// Solve the linear system (re-use the already factorized matrix A)
lu.solve(u->vector()), b);
L.u0 = u;
// Move to next interval
p = t / T;
t += dt;
}
Is there any possibility to use adaptive time-stepping or calculate the maximum time-step to reach a certain tolerance for a norm ||u_t - u_t+1|| < tol?
