Made ReactionDiffusionFE solver behave similarly to DiffusionSolveFE (i.e. position dependent diff/decay constants secretion interleaved with diffusive steps). Old ReactionDiffusionSolverFE is available under name FlexibleReactionDiffusionSolverFE - because it behaved like FlexibleDiffusionSolver having same limitations .