Abstract
We consider the mixed-variable finite-element discretization of the biharmonic problem. The numerical solution of the associated linear equations by preconditioned conjugate-gradient methods is studied. The role of squaring matrices for preconditioning is investigated.