set Additional shared libraries = ./libexponential_decay.so
set Dimension                              = 2
set Start time                             = 0
set End time                               = 100

# We use a new solver scheme option that enables the operator split.
set Nonlinear solver scheme                = single Advection, single Stokes
set Use operator splitting                 = true
set Maximum time step                      = 10

subsection Solver parameters
  subsection Operator splitting parameters
    set Reaction time step                 = 0.0005
  end
end
