gbrochar
/
computorv1
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

feat(maths): complex div

This commit is contained in:
gbrochar 2023-08-06 13:55:18 +02:00
parent 58a0ec983d
commit 3500de6cb0
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/maths.rs
View File

@ -31,6 +31,11 @@ impl Rational {
pub fn inverse(&self) -> Self { pub fn inverse(&self) -> Self {
Rational::new(self.denominator, self.numerator) Rational::new(self.denominator, self.numerator)
} }
pub fn simplify(lhs: Rational, rhs: Rational) -> Self {
let res = Rational::new(lhs.numerator * rhs.denominator, lhs.denominator * rhs.numerator);
res.reduce()
}
} }
impl ops::Add<Rational> for Rational { impl ops::Add<Rational> for Rational {
@ -172,4 +177,12 @@ impl ops::Mul<GaussianRational> for GaussianRational {
GaussianRational::new(self.real * rhs.real - self.imaginary * rhs.imaginary, self.real * rhs.imaginary + self.imaginary * rhs.real) GaussianRational::new(self.real * rhs.real - self.imaginary * rhs.imaginary, self.real * rhs.imaginary + self.imaginary * rhs.real)
} }
} }
impl ops::Div<GaussianRational> for GaussianRational {
type Output = GaussianRational;
fn div(self, rhs: GaussianRational) -> GaussianRational {
let cd_squared = rhs.real * rhs.real + rhs.imaginary * rhs.imaginary;
GaussianRational::new(Rational::simplify(self.real * rhs.real + self.imaginary * rhs.imaginary, cd_squared), Rational::simplify(self.imaginary * rhs.real - self.real * rhs.imaginary, cd_squared))
}
} }