seemps.evolution.radau#

seemps.evolution.radau(H: MPO, time: float | tuple[float | int | floating | integer, float | int | floating | integer] | Sequence[float | int | floating | integer] | ndarray[tuple[Any, ...], dtype[_ScalarT]], state: MPS, steps: int = 1000, stages: int = 3, inv_tol: float = 1e-7, strategy: Strategy = DEFAULT_STRATEGY, callback: Callable[[float, MPS], Any] | None = None, itime: bool = False) TypeAliasForwardRef('~seemps.state.MPS') | list[Any][source]#

Solve a Schrödinger equation 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.

Parameters:
stagesint, default = 3

Number of Radau IIA stages (3 or 5).

inv_tolfloat, default = 1e-7

Tolerance for the GMRES solver.