In line 22 the actual integration step is taken. The integration is done using the Dormand-Prince 8(5,3) explicit Runge-Kutta method as implemented in MidasCpp. After the integration step is completed, the current time has been updated and the wave function parameters and the time-derivative of these have also been updated. In line 24 the one-mode densities for the specific wave function method are calculated using the updated wave function parameters and the relevant space-and time-smoothening is carried out. In order to not evaluate the wave function densities on a grid in each time step, the one-mode densities are only calculated at intervals spaced by t dens . Note that t dens should be chosen smaller than or equal to t TD-ADGA since it makes no sense to check the TD-ADGA convergence if the one-mode densities have not been updated.