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
|
use super::{Complex, Float};
impl Float for f32 {
fn sqrt(self) -> Self {
libm::sqrtf(self)
}
fn atan2(y: Self, x: Self) -> Self {
libm::atan2f(y, x)
}
}
impl Float for f64 {
fn sqrt(self) -> Self {
libm::sqrt(self)
}
fn atan2(y: Self, x: Self) -> Self {
libm::atan2(y, x)
}
}
impl Float for Complex<f32> {
fn sqrt(self) -> Self {
let r = self.abs();
let theta = self.phase();
let sqrt_r = libm::sqrtf(r);
let sqrt_theta = theta / 2.0;
Complex::new(
sqrt_r * libm::cosf(sqrt_theta),
sqrt_r * libm::sinf(sqrt_theta),
)
}
fn atan2(y: Self, x: Self) -> Self {
Complex::new(
libm::atan2f(y.real, x.real),
libm::atan2f(y.imaginary, x.imaginary),
)
}
}
impl Float for Complex<f64> {
fn sqrt(self) -> Self {
let r = self.abs();
let theta = self.phase();
let sqrt_r = libm::sqrt(r);
let sqrt_theta = theta / 2.0;
Complex::new(
sqrt_r * libm::cos(sqrt_theta),
sqrt_r * libm::sin(sqrt_theta),
)
}
fn atan2(y: Self, x: Self) -> Self {
Complex::new(
libm::atan2(y.real, x.real),
libm::atan2(y.imaginary, x.imaginary),
)
}
}
|
