Time vector is defined by the nb_time_steps + 1.

Control are resampled (interpolated to the state variable grid). It notably simplifies the code. Avoid to specify for each discretization scheme.