Transform a nonzero initial velocity into a radial momentum scaling with 1 / r^(dim -1).
Member ryujin::Interface::ParabolicModule< dim, Number >::backward_euler_step (const StateVector< Number > &old_state_vector, StageVectors stage_state_vectors, StageWeights stage_weights, StateVector< Number > &new_state_vector, Number t, Number tau, Number tau_max=std::numeric_limits< Number >::max()) const =0
This function should always operate on the selected memory space.
For HyperbolicSystem/ ParabolicSystem / InitialValues realizations that are already MPIEnsemble "aware" we should fall back to only creating one instance - possibly by checking for and calling a modified constructor that takes an mpi_ensemble as an argument.
Member ryujin::quadratic_newton_step (Number &p_1, Number &p_2, const Number phi_p_1, const Number phi_p_2, const Number dphi_p_1, const Number dphi_p_2, const Number sign=Number(1.0))
Write out the quadratic Newton step in more detail.