Abstract
We obtain an algorithm for the skew-Takagi factorization based on the singular value decomposition. As an example application of such an algorithm, we consider the computation of Heisenberg evolution of the fermionic creation and annihilation operators in the case of the quadratic Hamiltonian of a special form
Disclosure statement
No potential conflict of interest was reported by the author(s).