Krylov implicit integration factor (IIF) methods were developed in Chen and Zhang (J Comput Phys 230:4336–4352, 2011) for solving stiff reaction–diffusion equations on high dimensional unstructured meshes. The methods were further extended to solve stiff advection–diffusion–reaction equations in Jiang and Zhang (J Comput Phys 253:368–388, 2013). Recently we studied the computational power of Kr...