arpfloat 0.1.11

Arbitrary-precision floating point 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
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
119
120
121
122
123
124
125
126
127

use arpfloat::{BigInt, Float, RoundingMode, Semantics};

use RoundingMode::NearestTiesToEven as rme;

fn test_e() {
    let sem = Semantics::new(32, 2000, rme);
    black_box(Float::e(sem));
}

fn test_sqrt() {
    let sem = Semantics::new(32, 10000, rme);
    black_box(Float::one(sem, false).scale(1, rme).sqrt());
}

fn test_pi() {
    let sem = Semantics::new(32, 2000, rme);
    black_box(Float::pi(sem));
}

fn test_powi() {
    let a = BigInt::from_u64(1275563424);
    black_box(a.powi(11000));
}

fn test_bigint_as_dec() {
    let a = BigInt::from_u64(197123);
    black_box(a.powi(100).as_decimal());
}

fn test_bigint_div() {
    let a = BigInt::pseudorandom(1000, 12345);
    let b = BigInt::pseudorandom(500, 67890);
    black_box(a / b);
}

fn test_cos() {
    let sem = Semantics::new(32, 90, rme);
    for i in 0..100 {
        let a = Float::from_u64(sem, i).cos();
        black_box(a);
    }
}

fn test_sin() {
    let sem = Semantics::new(32, 90, rme);
    for i in 0..100 {
        let a = Float::from_u64(sem, i).sin();
        black_box(a);
    }
}

fn test_log() {
    let sem = Semantics::new(32, 100, rme);
    for i in 0..100 {
        let a = Float::from_u64(sem, i).log();
        black_box(a);
    }
}

fn test_exp() {
    let sem = Semantics::new(32, 100, rme);
    for i in 0..1000 {
        let a = Float::from_u64(sem, 100 - i).exp();
        let b = Float::from_u64(sem, i).exp();
        black_box(a + b);
    }
}

fn test_bigint_mul_1() {
    let a = BigInt::pseudorandom(1000, 98765);
    let b = BigInt::pseudorandom(1000, 43210);
    black_box(a * b);
}

fn test_bigint_mul_2() {
    let a = BigInt::pseudorandom(10, 98765);
    let b = BigInt::pseudorandom(10, 43210);
    black_box(a * b);
}

fn test_bigint_mul_3() {
    let a = BigInt::pseudorandom(100, 98765);
    let b = BigInt::pseudorandom(100, 43210);
    black_box(a * b);
}

fn test_bigint_mul_4() {
    let a = BigInt::pseudorandom(5000, 98765);
    let b = BigInt::pseudorandom(1, 43210);
    black_box(a * b);
}

fn test_bigint_div_1() {
    let a = BigInt::pseudorandom(1000, 98765);
    let b = BigInt::pseudorandom(1000, 43210);
    black_box(a / b);
}

fn test_bigint_div_2() {
    let a = BigInt::pseudorandom(1000, 98765);
    let b = BigInt::pseudorandom(1, 43210);
    black_box(a / b);
}

use criterion::{black_box, criterion_group, criterion_main, Criterion};

pub fn criterion_benchmark(c: &mut Criterion) {
    c.bench_function("pi", |b| b.iter(test_pi));
    c.bench_function("e", |b| b.iter(test_e));
    c.bench_function("sqrt", |b| b.iter(test_sqrt));
    c.bench_function("powi", |b| b.iter(test_powi));
    c.bench_function("bigint_as_dec", |b| b.iter(test_bigint_as_dec));
    c.bench_function("bigint_div", |b| b.iter(test_bigint_div));
    c.bench_function("cos", |b| b.iter(test_cos));
    c.bench_function("sin", |b| b.iter(test_sin));
    c.bench_function("exp", |b| b.iter(test_exp));
    c.bench_function("log", |b| b.iter(test_log));
    c.bench_function("bigint_mul_1", |b| b.iter(test_bigint_mul_1));
    c.bench_function("bigint_mul_2", |b| b.iter(test_bigint_mul_2));
    c.bench_function("bigint_mul_3", |b| b.iter(test_bigint_mul_3));
    c.bench_function("bigint_mul_4", |b| b.iter(test_bigint_mul_4));
    c.bench_function("bigint_div_1", |b| b.iter(test_bigint_div_1));
    c.bench_function("bigint_div_2", |b| b.iter(test_bigint_div_2));
}

criterion_group!(benches, criterion_benchmark);
criterion_main!(benches);