complex_algebra 0.1.8

A crate supporting complex number algebra
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83

use super::c;
use std::ops::{Add, Div, Mul, Sub};

macro_rules! complex_div {
  ($LHS:ty, $RHS:ty, $T:tt ) => {
    impl<$T: Add<Output = $T> + Mul<Output = $T> + Sub<Output = $T> + Div<Output=T> + Copy + PartialEq> Div<$RHS>
      for $LHS
    {
      type Output = c<$T>;
      fn div(self, rhs: $RHS) -> Self::Output {
        let r = rhs.0*rhs.0+rhs.1*rhs.1;
        c(
          (self.0 * rhs.0 + self.1 * rhs.1)/r,
          (self.1 * rhs.0 - self.0 * rhs.1)/r,
        )
      }
    }
  };
}
complex_div!(c<T>, c<T>, T);
complex_div!(&c<T>, &c<T>, T);
complex_div!(c<T>, &c<T>, T);
complex_div!(&c<T>, c<T>, T);

macro_rules! complex_real_div {
  ($LHS:ty, $T:tt ) => {
    impl<$T: Add<Output = $T> + Mul<Output = $T> + Sub<Output = $T> + Div<Output=T> + Copy + PartialEq> Div<$T>
      for $LHS
    {
      type Output = c<$T>;
      fn div(self, rhs: $T) -> Self::Output {
        c(self.0/rhs, self.1/rhs)
      }
    }
  };
}
complex_real_div!(c<T>, T);
complex_real_div!(&c<T>, T);

#[cfg(test)]
mod tests {
  use super::*;

  #[test]
  fn test_div_1() {
    let z1 = c(3f64, 4f64);
    let z2 = c(8f64, -2f64);
    assert_eq!(z1 / z2, c(4.0 / 17.0, 19.0 / 34.0));
  }

  #[test]
  fn test_div_2() {
    let z1 = c(3f64, 4f64);
    let z2 = c(8f64, -2f64);
    assert_eq!(z1 / &z2, c(4.0 / 17.0, 19.0 / 34.0));
  }

  #[test]
  fn test_div_3() {
    let z1 = c(3f64, 4f64);
    let z2 = c(8f64, -2f64);
    assert_eq!(&z1 / z2, c(4.0 / 17.0, 19.0 / 34.0));
  }

  #[test]
  fn test_div_4() {
    let z1 = c(3f64, 4f64);
    let z2 = c(8f64, -2f64);
    assert_eq!(&z1 / &z2, c(4.0 / 17.0, 19.0 / 34.0));
  }

  #[test]
  fn test_div_5() {
    let z = c(3f64, 4f64);
    assert_eq!(z / 3.0, c(1.0, 4.0 / 3.0));
  }

  #[test]
  fn test_div_6() {
    let z = c(3f64, 4f64);
    assert_eq!(&z / 3.0, c(1.0, 4.0 / 3.0));
  }
}