“…For this reason, time dependent implementations of RD always feature some form of implicit time integration [12,13,14,32,33], or a fully coupled space-time formulation [17,34,35]. Moreover, positivity preservation always requires the satisfaction of time step constraints [36,14], unless some form of nonlinear time (or space-time) discretization is used [37,34,35].…”