Recently, global and local relaxation Runge–Kutta methods have been developed for guaranteeing the conservation, dissipation, or other solution properties general convex functionals whose dynamics are crucial an ordinary differential equation solution. These novel time integration procedures application in a wide range of problems that require dynamics-consistent stable numerical methods. The s...