2This is done by choosing parameters that minimizes the largest eigenvalue of the amplification matrix for a simple 1D advection equation with a penalty method implementation of periodic boundary conditions as suggested in [2].