The starting step-size is computed following the algorithm in [1]}. In the case of simple ODEs (Van der Pol and Brusselator), we used the computed optimal starting step-size. For the PDEs, we take the minimum of the starting step-size returned from the algorithm and the step-size restriction based on the CFL constraints to avoid any stability issues from influencing the numerical study.
