by combining a function transform approach with the scalar auxiliary variable (SAV)
approach. A suitable function transform ensures that the schemes are bound/positivity
preserving, while the SAV approach enables us to construct unconditionally stable schemes.
The first class of schemes requires solving a coupled second-order linear systems with
variable coefficients at each time step, while the second class of schemes only requires …