hcomplex 0.2.0

Lightweight complex and hypercomplex algebra 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

use rand::{prelude::*};
use rand_xorshift::XorShiftRng;
use approx::*;
use crate::{transform::*, Complex, Quaternion, Octonion};


const TRANSFORM_ATTEMPTS: usize = 64;
const POINT_ATTEMPTS: usize = 16;

#[test]
fn moebius2() {
    let mut rng = XorShiftRng::seed_from_u64(0xDEAD0);
    for _ in 0..TRANSFORM_ATTEMPTS {
        let a: Moebius<Complex<f64>> = rng.sample(&Normalized);
        let b: Moebius<Complex<f64>> = rng.sample(&Normalized);
        let c = a.chain(b);
        for _ in 0..POINT_ATTEMPTS {
            let x: Complex<f64> = rng.sample(&StandardNormal);
            let y = a.apply(b.apply(x));
            let z = c.apply(x);
            assert_abs_diff_eq!(y, z, epsilon=1e-12);
        }
    }
}

#[test]
fn moebius4() {
    let mut rng = XorShiftRng::seed_from_u64(0xDEAD1);
    for _ in 0..TRANSFORM_ATTEMPTS {
        let a: Moebius<Quaternion<f64>> = rng.sample(&Normalized);
        let b: Moebius<Quaternion<f64>> = rng.sample(&Normalized);
        let c = a.chain(b);
        for _ in 0..POINT_ATTEMPTS {
            let x: Quaternion<f64> = rng.sample(&StandardNormal);
            let y = a.apply(b.apply(x));
            let z = c.apply(x);
            assert_abs_diff_eq!(y, z, epsilon=1e-12);
        }
    }
}

/// Moebuis transform over octonions isn't chainable and therefore should fail
#[test]
#[should_panic]
fn moebius8() {
    let mut rng = XorShiftRng::seed_from_u64(0xDEAD2);
    for _ in 0..TRANSFORM_ATTEMPTS {
        let a: Moebius<Octonion<f64>> = rng.sample(&Normalized);
        let b: Moebius<Octonion<f64>> = rng.sample(&Normalized);
        let c = a.chain(b);
        for _ in 0..POINT_ATTEMPTS {
            let x: Octonion<f64> = rng.sample(&StandardNormal);
            let y = a.apply(b.apply(x));
            let z = c.apply(x);
            assert_abs_diff_eq!(y, z, epsilon=1e-12);
        }
    }
}