iter wrong code but it has no red squiggle
This commit is contained in:
parent
ee7df7787a
commit
f7ee18076a
@ -2,31 +2,31 @@
|
|||||||
use crate::integer::Planar64;
|
use crate::integer::Planar64;
|
||||||
|
|
||||||
#[inline]
|
#[inline]
|
||||||
pub fn zeroes2<'a>(a0:Planar64,a1:Planar64,a2:Planar64)->&'a [Planar64]{
|
pub fn zeroes2(a0:Planar64,a1:Planar64,a2:Planar64)->impl std::iter::IntoIterator{
|
||||||
if a2==Planar64::ZERO{
|
// if a2==Planar64::ZERO{
|
||||||
return zeroes1(a0,a1);
|
// return zeroes1(a0,a1);
|
||||||
}
|
// }
|
||||||
let radicand=a1.get() as i128*a1.get() as i128-a2.get() as i128*a0.get() as i128*4;
|
let radicand=a1.get() as i128*a1.get() as i128-a2.get() as i128*a0.get() as i128*4;
|
||||||
if 0<radicand{
|
if 0<radicand{
|
||||||
//start with f64 sqrt
|
//start with f64 sqrt
|
||||||
let planar_radicand=Planar64::raw(unsafe{(radicand as f64).sqrt().to_int_unchecked()});
|
let planar_radicand=Planar64::raw(unsafe{(radicand as f64).sqrt().to_int_unchecked()});
|
||||||
//TODO: one or two newtons
|
//TODO: one or two newtons
|
||||||
if Planar64::ZERO<a2{
|
if Planar64::ZERO<a2{
|
||||||
[(-a1-planar_radicand)/(a2*2),(-a1+planar_radicand)/(a2*2)].as_slice()
|
[(-a1-planar_radicand)/(a2*2),(-a1+planar_radicand)/(a2*2)].into_iter()
|
||||||
}else{
|
}else{
|
||||||
[(-a1+planar_radicand)/(a2*2),(-a1-planar_radicand)/(a2*2)].as_slice()
|
[(-a1+planar_radicand)/(a2*2),(-a1-planar_radicand)/(a2*2)].into_iter()
|
||||||
}
|
}
|
||||||
}else if radicand==0{
|
}else if radicand==0{
|
||||||
[a1/(a2*-2)].as_slice()
|
[a1/(a2*-2),Planar64::ZERO].into_iter()
|
||||||
}else{
|
}else{
|
||||||
&[]
|
[Planar64::ZERO,Planar64::ZERO].into_iter()
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
#[inline]
|
#[inline]
|
||||||
pub fn zeroes1<'a>(a0:Planar64,a1:Planar64)->&'a [Planar64]{
|
pub fn zeroes1(a0:Planar64,a1:Planar64)->impl std::iter::IntoIterator{
|
||||||
if a1==Planar64::ZERO{
|
if a1==Planar64::ZERO{
|
||||||
&[]
|
return [Planar64::ZERO].into_iter();
|
||||||
}else{
|
}else{
|
||||||
[-a0/a1].as_slice()
|
return [-a0/a1].into_iter();
|
||||||
}
|
}
|
||||||
}
|
}
|
Loading…
Reference in New Issue
Block a user