“…We point out that the latter property was already implemented in preconditioning techniques to reduce the number of steps of Conjugate Gradient (CG) methods for pd linear systems Bx = v [4], [12], [13], [22], [28], [29], [36], [40]. Moreover, the class of algebras of L-type, for which the above property holds, has been also used to define more efficient direct methods solving the system Bx = v. By utilizing simple displacement decompositions of B −1 in terms of matrices of L-type, it is, in fact, possible to compute B −1 v via a finite number of fast unitary transforms [5], [8], [17], [21], [30].…”