SSP methods are used heavily for numerically integrating IVP resulting from discretized hyperbolic problems since they offer properties not present in non-SSP methods. Modern robust IVP solvers include many important features such as error estimation and automatic step-size control [1]}, [2]}, [3]}. Most of these important features have not yet been developed for existing higher-order optimal explicit SSP methods. The current work provides these important features for the optimal SSP methods. The numerical results provide evidence of the robustness of these pairs.
