extern/kalker
1
0
Fork 0
mirror of https://github.com/PaddiM8/kalker.git synced 2024-12-12 17:40:52 +01:00
Code Issues Packages Projects Releases Wiki Activity

Division with complex numbers

Browse Source
This commit is contained in:
bakk 2021-05-20 19:50:59 +02:00
parent ee3082ad66
commit 328810d29d
1 changed files with 25 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
kalk/src/kalk_num/mod.rs
View File

@ -292,8 +292,23 @@ impl KalkNum {
}
pub(crate) fn div(self, context: &mut crate::interpreter::Context, rhs: KalkNum) -> KalkNum {
let right = calculate_unit(context, &self, rhs.clone()).unwrap_or(rhs);
KalkNum::new(self.value / right.value, &right.unit)
let right = calculate_unit(context, &self, rhs.clone()).unwrap_or(rhs.clone());
// Avoid unecessary calculations
if self.imaginary_value == 0f64 && right.imaginary_value == 0f64 {
KalkNum::new(self.value / right.value, &right.unit)
} else {
// Multiply both the numerator and denominator
// with the conjugate of the denominator, and divide.
let conjugate = rhs.get_conjugate();
let numerator = self.clone().mul(context, conjugate.clone());
let denominator = rhs.clone().mul(context, conjugate);
KalkNum::new_with_imaginary(
numerator.value / denominator.value.clone(),
&right.unit,
numerator.imaginary_value / denominator.value,
)
}
}
pub(crate) fn rem(self, context: &mut crate::interpreter::Context, rhs: KalkNum) -> KalkNum {
@ -301,6 +316,14 @@ impl KalkNum {
KalkNum::new(self.value % right.value, &right.unit)
}
pub fn get_conjugate(&self) -> KalkNum {
KalkNum::new_with_imaginary(
self.value.clone(),
&self.unit,
self.imaginary_value.clone() * (-1f64),
)
}
/// Get an estimate of what the number is, eg. 3.141592 => π
pub fn estimate(&self) -> Option<String> {
let fract = self.value.clone().fract().abs();