Automatic continuation on time steps

Add new syntax for `solve` with list of grid sizes, such that ```solve(..., grid_size=[N1, N2, N3])```
will solve in sequence for N1, N2, N3 steps, reusing solutions as initial guesses.

Add method with list for grid_size, that will call the usual solve in sequence with warm start.
