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.Comment: 11 pages, to appear in proceedings of the workshop "Numerical
Challenges in Lattice QCD", Springer Verla
