In this section we will describe the implementation of the TD-ADGA algorithm as it has been implemented in the Molecular Interactions, Dynamics And Simulations Chemistry Program Package (MidasCpp). The implementation relies on the previous implementations in MidasCpp for the integration and propagation of ordinary differential equations for timedependent nuclear wave functions as well as the ADGA framework for handling SP calculations, grid placement, and evaluation of integrals. A high level description of how the TD-ADGA procedure is implemented as a part the ordinary differential equation integrator is shown in Algorithm 1. The ordinary differential equation integrator is on input given the initial and end times (t 0 and t end ), the initial wave funtion parameters (y 0 ), and the initial PES (V 0 ). Further TD-ADGA inputs are also given, these include a boolean indicating if a full restarts should be utilized, the time interval between running TD-ADGA convergence checks (t TD-ADGA ), the time interval at which the one-mode densities are calculated (t dens ), and the time interval between reset points (t reset ). These intervals will be explained further below. In the initial 5 lines of Algorithm 1 the integration is prepared and in the while loop from line 6 to line 38 the actual time-propagation takes place until the dynamics have been propagated to the end time.