Before revision:  These approximations make the scheme to be unconditionally stable in time and preserve positivity of the solution which is computed  via a matrix exponential . The paper demonstrates that the proposed method is computationally efficient, accurate and simple to implement } .
After revision:  These approximations make the scheme to be unconditionally stable in time and preserve positivity of the solution which is computed via a matrix exponential , or via P } .