seemps.evolution.radau#
- seemps.evolution.radau(L: MPO, time: float | tuple[float | int | floating | integer, float | int | floating | integer] | Sequence[float | int | floating | integer] | NDArray, state: MPS, steps: int = 1000, stages: int = 3, inv_tol: float = 1e-7, strategy: Strategy = DEFAULT_STRATEGY, callback: Callable[[float, MPS], Any] | None = None) TypeAliasForwardRef('~seemps.state.MPS') | list[Any][source]#
Solve
d|state>/dt = L|state>using an implicit Radau IIA method with either 3 or 5 stages (order 5 or 9, respectively).See
seemps.evolution.euler()for a description of the missing function arguments and the function’s output.