In this paper, we discuss the numerical solution of systems of differential equations with both linear (stiff) and nonlinear components that arise from the semi-discretization of certain partial differential equations. When faced with the task of solving problems with discontinuous solutions numerically, the linear stability properties are insufficient for convergence. So, to overcome this drawback, strong stability preserving (SSP) methods are introduced. While using SSP implicit-explicit methods lead to severe constraints on the allowed time-stepping, an integrating factor approach reduces this limitation by solving the linear part exactly. In this work, we develop integrating factor general linear methods (IFGLMs) with strong stability properties. Sufficient conditions for IFGLMs to be SSP are discussed. The construction of IFGLMs of order p ≤ 6 and stage order 1 ≤ q ≤ p with r = p external stages and 2 ≤ s ≤ 10 internal stages is presented, which have larger effective SSP coefficients than the class of Runge–Kutta and two–step Runge–Kutta methods. Results are verified numerically on several representative test cases.
