At each stage, elements above the diagonal are eliminated as well as
those eliminated below as in standard Gaussian elimination. No
back-substitution is required, but only limited pivoting is possible, so
that accuracy may be lost. The operational count is
.
John Gilbert
1999-02-25