Multiple ways to solve diffusion reaction systems (fully-implicit reduced approach, operator- splitting approach, domain-decomposition approach, exploitation of block structure) are available. Criteria for the selection of methods have been determined.