boostvoronoi 0.12.1

Boost voronoi ported to 100% rust
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

#![allow(unused_imports)]
use crate::extended_scalar::robust_fpt::RobustDif;
use crate::extended_scalar::robust_fpt::RobustFpt;

#[test]
fn robustdif_sum_1() {
    let a_: f64 = 5.;
    let b_: f64 = 4.;

    let a = RobustDif::from(a_);
    let b = RobustDif::from(b_);
    let s = a + b;

    assert_eq!(s.dif().fpv(), a_ + b_);
    assert_eq!(s.dif().fpv(), 9.0)
}

#[test]
fn robustdif_sum_2() {
    type F2 = f64;

    let newf = RobustFpt::new;
    let c_x = RobustDif::from((newf(12.0, 2.0), newf(0.5, 13.0)));
    let mut r = RobustDif::from((newf(24.0, 2.0), newf(26.0, 2.0)));
    let line_a: F2 = -2.0;
    let line_a_2 = RobustFpt::from(line_a);

    assert_eq!(line_a_2.fpv(), line_a);
    dbg!(&c_x);
    dbg!(&line_a_2);
    let v = c_x * line_a_2;
    dbg!(&v);
    assert_eq!(v.dif().fpv(), (12.0 - 0.5) * -2.0);

    dbg!(&c_x, &line_a, &r);
    r += c_x * RobustFpt::from(line_a);
    dbg!(&r);

    assert_eq!(r.positive().fpv(), 25.0);
    //assert_eq!(r.positive().re(), 15.0);

    assert_eq!(r.negative().fpv(), 50.0);
    //assert_eq!(r.negative().re(), 4.0);
}

#[test]
fn robustdif_prod_1() {
    let a = RobustDif::from((RobustFpt::new(35058881.0, 7.0), RobustFpt::new(0.0, 0.0)));
    let b = RobustFpt::new(0.0, 3.0);
    let p = a * b;
    println!("p:{p:?}");
    assert_eq!(p.positive().fpv(), 0.0);
    assert_eq!(p.positive().re(), 11.0);
    assert_eq!(p.negative().fpv(), 0.0);
    assert_eq!(p.negative().re(), 4.0);
}