Numerical method for unitary systems (Q1077130)

The Lanczos method is modified to construct probability amplitudes when the corresponding propagator can be expressed in terms of the Cayley representation. The method is the most appropriate when the Hermitian matrix in the Cayley matrix representation is sparse. This algorithm requires much less storage compared with the conventional method. The algorithm is monitored by means of the magnitude of the residual vector.