TSTHETA

Notes

This method is cast as a 1-stage implicit Runge-Kutta method.

  Theta | Theta
  -------------
        |  1

For the default Theta=0.5, this is also known as the implicit midpoint rule.

When the endpoint variant is chosen, the method becomes a 2-stage method with first stage explicit

  0 | 0         0
  1 | 1-Theta   Theta
  -------------------
    | 1-Theta   Theta

For the default Theta=0.5, this is the trapezoid rule (also known as Crank-Nicolson, see TSCN).

To apply a diagonally implicit RK method to DAE, the stage formula

 Y_i = X + h sum_j a_ij Y'_j

is interpreted as a formula for Y'_i in terms of Y_i and known stuff (Y'_j, j<i)

See Also

Level:beginner
Location:src/ts/impls/implicit/theta/theta.c
Index of all TS routines
Table of Contents for all manual pages
Index of all manual pages

Examples