seemps.evolution.euler2#

seemps.evolution.euler2(L: MPO, time: float | tuple[float | int | floating | integer, float | int | floating | integer] | Sequence[float | int | floating | integer] | NDArray, state: MPS, steps: int = 1000, strategy: Strategy = DEFAULT_STRATEGY, callback: Callable[[float, MPS], Any] | None = None)[source]#

Solve d|state>/dt = L|state> using the 2nd order Euler method.

Implements a two-step integration method,

\[\xi(t_{n}) = 2 \psi(t_{n}) + \delta t L \psi(t_{n}) \psi(t_{n+1}) = \psi(t_{n}) + (\delta t/2) L \xi(t_{n})\]

The Euler algorithm is a very bad integrator and is offered only for illustrative purposes. See euler() to understand the parameters and the function’s output.