In this paper, we suggest a technique to avoid order reduction in time when integrating reaction–diffusion boundary value problems under non-homogeneous boundary conditions with exponential splitting methods. More precisely, we consider Lie–Trotter and Strang splitting methods and Dirichlet, Neumann and Robin boundary conditions. Beginning from an abstract framework in Banach spaces, a thorough error analysis after full discretization is performed and some numerical results are shown which corroborate the theoretical results.
