diff --git a/fixed_wide/src/zeroes.rs b/fixed_wide/src/zeroes.rs index 97da6c9..865d7d3 100644 --- a/fixed_wide/src/zeroes.rs +++ b/fixed_wide/src/zeroes.rs @@ -17,10 +17,7 @@ macro_rules! impl_zeroes{ let radicand=a1.[](a1)-a2.[](a0)*4; match radicand.cmp(&Fixed::<{$n*2},{$n*2*32}>::ZERO){ Ordering::Greater=>{ - //start with f64 sqrt - //failure case: 2^63 < sqrt(2^127) let planar_radicand=radicand.sqrt().halve_precision(); - //TODO: one or two newtons //sort roots ascending and avoid taking the difference of large numbers match (a2pos,Self::ZERO[Ratio::new(-a1-planar_radicand,a2*2),Ratio::new(a0*2,-a1-planar_radicand)].into(),