liam/geometric_algebra
1
0
Fork 0
Code Issues Pull requests Projects Releases Wiki Activity

Uses SIMD for division.

Browse source
This commit is contained in:
Alexander Meißner 2022-10-16 12:47:05 +02:00
parent 6e56b66ec9
commit b721a48313
1 changed files with 7 additions and 12 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

19
src/lib.rs
View file

@ -86,7 +86,7 @@ impl Exp for ppga2d::IdealPoint {
ppga2d::Translator { ppga2d::Translator {
g0: simd::Simd32x3 { g0: simd::Simd32x3 {
f32x3: [1.0, self.g0[0], self.g0[1]], f32x3: [1.0, self.g0[0], self.g0[1]],
} },
} }
} }
} }
@ -95,11 +95,8 @@ impl Ln for ppga2d::Translator {
type Output = ppga2d::IdealPoint; type Output = ppga2d::IdealPoint;
fn ln(self) -> ppga2d::IdealPoint { fn ln(self) -> ppga2d::IdealPoint {
ppga2d::IdealPoint { let result: ppga2d::IdealPoint = self.into();
g0: simd::Simd32x2 { result / ppga2d::Scalar { g0: self.g0[0] }
f32x2: [self.g0[1] / self.g0[0], self.g0[2] / self.g0[0]],
}
}
} }
} }
@ -110,6 +107,7 @@ impl Powf for ppga2d::Translator {
(ppga2d::Scalar { g0: exponent } * self.ln()).exp() (ppga2d::Scalar { g0: exponent } * self.ln()).exp()
} }
} }
impl Exp for ppga2d::Point { impl Exp for ppga2d::Point {
type Output = ppga2d::Motor; type Output = ppga2d::Motor;
@ -244,7 +242,7 @@ impl Exp for ppga3d::IdealPoint {
ppga3d::Translator { ppga3d::Translator {
g0: simd::Simd32x4 { g0: simd::Simd32x4 {
f32x4: [1.0, self.g0[0], self.g0[1], self.g0[2]], f32x4: [1.0, self.g0[0], self.g0[1], self.g0[2]],
} },
} }
} }
} }
@ -253,11 +251,8 @@ impl Ln for ppga3d::Translator {
type Output = ppga3d::IdealPoint; type Output = ppga3d::IdealPoint;
fn ln(self) -> ppga3d::IdealPoint { fn ln(self) -> ppga3d::IdealPoint {
ppga3d::IdealPoint { let result: ppga3d::IdealPoint = self.into();
g0: simd::Simd32x3 { result / ppga3d::Scalar { g0: self.g0[0] }
f32x3: [self.g0[1] / self.g0[0], self.g0[2] / self.g0[0], self.g0[3] / self.g0[0]],
}
}
} }
} }