A factorization algorithm to compute Pfaffians Original Research Article

We describe an explicit algorithm to factorize an even antisymmetric image matrix into triangular and trivial factors. The construction resembles the Crout algorithm for LU factorization. It allows for a straight forward computation of Pfaffians (including their signs) at the cost of image flops.