Automatic continuation on time steps

[PierreMartinon]

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.