-pass along errors returned by cycle reduction algorithms;
-logarithmic reduction fails on error
-replace expensive and random condest() by call to linsolve()
-uses linsolve() to compute hx