  The Wilson fermion determinant can be written as product of the determinants
of two hermitian positive definite matrices. This formulation allows to
simulate non-degenerate quark flavors by means of the hybrid Monte Carlo
algorithm. A major numerical difficulty is the occurrence of nested inversions.
We construct a Uzawa iteration scheme which treats the nested system within one
iterative process.
