“…Alternatively, most previous works have considered semi-implicit (IMEX) time discretizations and/or operator splitting schemes, where the reaction and diffusion terms are treated separately, see e.g. [5,6,7,8,9,10,11,12,13,14,15]. The advantage of IMEX and operator splitting schemes is that they only require the solution of linear systems at each time step.…”