burbomath 0.0.1

Burbokop's rust math library
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
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121

use super::{Angle, Cos, Point, Sin};
use crate::math::{Sq, Vector};
use core::ops::{Add, Div, Mul, Neg, Not, Sub};

#[derive(Debug, Clone, Copy)]
pub struct Complex<T> {
    real: T,
    imag: T,
}

impl<T> From<(T, T)> for Complex<T> {
    fn from(value: (T, T)) -> Self {
        Self {
            real: value.0,
            imag: value.1,
        }
    }
}

impl<T> Complex<T> {
    pub fn from_cartesian(real: T, imag: T) -> Self {
        Self { real, imag }
    }

    pub fn from_polar(r: T, a: Angle<T>) -> Self
    where
        T: Clone + Cos<Output = T> + Sin<Output = T> + Mul<Output = T>,
    {
        Self {
            real: a.clone().cos() * r.clone(),
            imag: a.sin() * r,
        }
    }

    pub fn from_uneven_polar(r: Vector<T>, a: Angle<T>) -> Self
    where
        T: Clone + Cos<Output = T> + Sin<Output = T> + Mul<Output = T>,
    {
        let (x, y) = r.into();
        Self {
            real: a.clone().cos() * x,
            imag: a.sin() * y,
        }
    }

    pub fn real(&self) -> &T {
        &self.real
    }

    pub fn imag(&self) -> &T {
        &self.imag
    }

    pub fn into_cartesian(self) -> Point<T> {
        (self.real, self.imag).into()
    }

    pub fn div(v0: Vector<T>, v1: Vector<T>) -> Self
    where
        T: Sq<Output = T>,
        T: Add<Output = T>,
        T: Sub<Output = T>,
        T: Mul<Output = T>,
        T: Div<Output = T>,
        T: Neg<Output = T>,
        T: Clone,
    {
        let (a, b) = v0.into();
        let (c, d) = v1.into();
        let len_sq = c.clone().sq() + d.clone().sq();

        (
            (a.clone() * c.clone() + b.clone() * d.clone()) / len_sq.clone(),
            (b * c - a * d) / len_sq,
        )
            .into()
    }
}

impl<T> Add for Complex<T>
where
    T: Add,
{
    type Output = Complex<<T as Add>::Output>;

    fn add(self, rhs: Self) -> Self::Output {
        Self::Output {
            real: self.real + rhs.real,
            imag: self.imag + rhs.imag,
        }
    }
}

impl<T> Mul for Complex<T>
where
    T: Mul<Output = T> + Clone + Add<Output = T> + Sub<Output = T>,
{
    type Output = Complex<T>;

    fn mul(self, rhs: Self) -> Self::Output {
        let (a, b) = self.into_cartesian().into();
        let (c, d) = (rhs.real, rhs.imag);
        (a.clone() * c.clone() - b.clone() * d.clone(), a * d + b * c).into()
    }
}

impl<T> Not for Complex<T>
where
    T: Sq<Output = T>,
    T: Add<Output = T>,
    T: Div<Output = T>,
    T: Neg<Output = T>,
    T: Clone,
{
    type Output = Self;

    fn not(self) -> Self::Output {
        let len_sq = self.real.clone().sq() + self.imag.clone().sq();
        (self.real / len_sq.clone(), -self.imag / len_sq).into()
    }
}