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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118

use llml::vec::{EucVecd2, EucVecf2};
use rand::random;

macro_rules! test_arith {
    ($sy:tt) => {
        let (alpha, beta) = get_vecs();
        assert_eq!(alpha $sy beta, EucVecd2::new([alpha.x() $sy beta.x(), alpha.y() $sy beta.y()]))
    }
}

#[cfg(feature = "llml_serde")]
#[test]
fn serde () {
    let alpha : EucVecd2 = random();
    let json = serde_json::to_string(&alpha).unwrap();
    let beta : EucVecd2 = serde_json::from_str(json.as_str()).unwrap();

    assert!((alpha - beta).abs().sum() <= f64::EPSILON * 2.);
}

#[test]
fn eq () {
    assert_eq!(EucVecd2::new([1., 2.]), EucVecd2::new([1., 2.]));
    assert_ne!(EucVecd2::new([1., 2.]), EucVecd2::new([3., 3.]))
}

#[test]
fn into () {
    let alpha = EucVecd2::new([1., 2.]);
    assert_eq!(Into::<EucVecf2>::into(EucVecd2::new([1., 2.])), EucVecf2::new([1., 2.]));
    assert_eq!(Into::<[f64;2]>::into(alpha), [1., 2.])
}

#[test]
fn from_scalar () {
    let alpha : f64 = random();
    assert_eq!(EucVecd2::from_scal(alpha), EucVecd2::new([alpha, alpha]))
}

#[test]
fn add () {
    test_arith!(+);
}

#[test]
fn sub () {
    test_arith!(-);
}

#[test]
fn mul () {
    test_arith!(*);
}

#[test]
fn div () {
    test_arith!(/);
}

#[test]
fn neg () {
    let alpha : EucVecd2 = random();
    assert_eq!(-alpha, EucVecd2::new([-alpha.x(), -alpha.y()]))
}

#[test]
fn sum () {
    let alpha : EucVecd2 = random();
    assert_eq!(alpha.sum(), alpha.x() + alpha.y())
}


#[test]
fn dot () {
    let (alpha, beta) = get_vecs();
    assert_eq!(alpha.dot(beta), alpha.x() * beta.x() + alpha.y() * beta.y())
}

#[test]
fn norm () {
    let alpha : EucVecd2 = random();
    assert_eq!(alpha.norm(), alpha.x().hypot(alpha.y()))
}

#[test]
fn unit () {
    let alpha : EucVecd2 = random();

    let norm = alpha.x().hypot(alpha.y());
    assert_eq!(alpha.unit(), EucVecd2::new([alpha.x() / norm, alpha.y() / norm]))
}

#[test]
fn abs () {
    let alpha : EucVecd2 = random();
    assert_eq!(alpha.abs(), EucVecd2::new([alpha.x().abs(), alpha.y().abs()]))
}

#[test]
fn sqrt () {
    let alpha : EucVecd2 = random();
    assert_eq!(alpha.sqrt(), EucVecd2::new([alpha.x().sqrt(), alpha.y().sqrt()]))
}

const RSQRT_EPSILON : f64 = 0.0003662109375 + f64::EPSILON;

#[test]
fn sqrt_fast () {
    let alpha : EucVecd2 = random();
    let fast = alpha.sqrt_fast();

    assert!((fast.x() - alpha.x().sqrt()).abs() <= RSQRT_EPSILON);
    assert!((fast.y() - alpha.y().sqrt()).abs() <= RSQRT_EPSILON);
}

fn get_vecs () -> (EucVecd2, EucVecd2) {
    (random(), random())
}