llml 0.2.3

Implementation of basic math data types with high level frontend and low level backend
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

use crate::others::{Complxd, Complxf};

macro_rules! impl_complx_type {
    ($($target:ident, $ty:ty),+) => {
        $(
            impl ComplxType<$ty> for $target {
                fn re (&self) -> $ty {
                    <$target>::re(self)
                }

                fn im (&self) -> $ty {
                    <$target>::im(self)
                }
            }
        )*
    }
}

macro_rules! impl_sqrtc {
    ($($target:ident, $ty:ty),+) => {
        $(
            impl ComplexSqrt for $ty {
                type Output = $target;

                #[inline(always)]
                fn sqrtc (self) -> Self::Output {
                    $target::sqrtc(self)
                }
            }
        )*
    };
}

macro_rules! impl_powc {
    ($($target:ident, $ty:ty),+) => {
        $(
            impl ComplexPow<$target> for $ty {
                type Output = $target;

                #[inline(always)]
                fn expi (self) -> Self::Output {
                    $target::expi(self)
                }

                #[inline(always)]
                fn powci (self, rhs: $ty) -> Self::Output {
                    $target::powci(self, rhs)
                }

                #[inline(always)]
                fn powc (self, rhs: $target) -> Self::Output {
                    self.powf(rhs.re()) * self.powci(rhs.im())
                }
            }
        )*
    };
}

macro_rules! impl_all {
    ($($target:ident, $ty:ty),+) => {
        $(
            impl_complx_type!($target, $ty);
            impl_sqrtc!($target, $ty);
            impl_powc!($target, $ty);
        )*
    }
}

pub trait ComplxType<T>: Copy {
    fn re (&self) -> T;
    fn im (&self) -> T;
}

pub trait ComplexSqrt<T = Self> {
    type Output: ComplxType<T>;

    /// Computes the square root of the value, returning the complex result
    fn sqrtc (self) -> Self::Output;
}

pub trait ComplexPow<Rhs: ComplxType<T>, T = Self> {
    type Output: ComplxType<T>;

    /// Computes ```exp(self * i)```
    fn expi (self) -> Self::Output;

    /// Computes ```pow(self, rhs * i)```
    fn powci (self, rhs: T) -> Self::Output;

    /// Computes ```pow(self, rhs)```, where ```rhs``` is a complex number
    fn powc (self, rhs: Rhs) -> Self::Output;
}

impl_all!(
    Complxf, f32,
    Complxd, f64
);