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.
