decimal-scaled 0.3.3

Const-generic base-10 fixed-point decimals (D9/D18/D38/D76/D153/D307) with integer-only transcendentals correctly rounded to within 0.5 ULP — exact at the type's last representable place. Deterministic across every platform; no_std-friendly.
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
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

//! Micro-bench: wide-tier `sqrt_strict` per tier and midpoint scale.
//!
//! Run with:
//!     cargo run --release --example sqrtbench --features wide,x-wide

use std::hint::black_box;
use std::time::Instant;

use decimal_scaled::{D153, D307, D38, D76};

const ITERS: u32 = 200_000;

fn time<F: FnMut()>(label: &str, mut f: F) {
    for _ in 0..2_000 {
        f();
    }
    let start = Instant::now();
    for _ in 0..ITERS {
        f();
    }
    let elapsed = start.elapsed();
    let per = elapsed.as_nanos() as f64 / ITERS as f64;
    println!("  {:<40} {:>10.1} ns/op", label, per);
}

fn main() {
    println!("== D38<19>::sqrt(2) — baseline ==");
    {
        let a = D38::<19>::from_int(2);
        time("D38<19>::sqrt_strict(2)", || {
            black_box(black_box(a).sqrt_strict());
        });
    }

    println!();
    println!("== D76<35>::sqrt(2) ==");
    {
        let a = D76::<35>::from_int(2);
        time("D76<35>::sqrt_strict(2)", || {
            black_box(black_box(a).sqrt_strict());
        });
    }

    println!();
    println!("== D76<76>::sqrt(2) ==");
    {
        let a = D76::<76>::from_int(2);
        time("D76<76>::sqrt_strict(2)", || {
            black_box(black_box(a).sqrt_strict());
        });
    }

    println!();
    println!("== D153<75>::sqrt(2) ==");
    {
        let a = D153::<75>::from_int(2);
        time("D153<75>::sqrt_strict(2)", || {
            black_box(black_box(a).sqrt_strict());
        });
    }

    println!();
    println!("== D153<153>::sqrt(2) ==");
    {
        let a = D153::<153>::from_int(2);
        time("D153<153>::sqrt_strict(2)", || {
            black_box(black_box(a).sqrt_strict());
        });
    }

    println!();
    println!("== D307<150>::sqrt(2) ==");
    {
        let a = D307::<150>::from_int(2);
        time("D307<150>::sqrt_strict(2)", || {
            black_box(black_box(a).sqrt_strict());
        });
    }

    println!();
    println!("== D307<307>::sqrt(2) ==");
    {
        let a = D307::<307>::from_int(2);
        time("D307<307>::sqrt_strict(2)", || {
            black_box(black_box(a).sqrt_strict());
        });
    }
}