fast_reals/timing_basic_operations.cpp at master · psg-mit/fast_reals · GitHub

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
// Header file for input output functions
#include<iostream>
#include <mpfr.h>
#include <functional>
#include <math.h>
#include <list>

using namespace std;


int main()
{
    // register constants
    mpfr_t pi, one, e, result;

    function<double(double, double)> convtime = [](clock_t begin, clock_t end)
    {
        return double(end - begin) / CLOCKS_PER_SEC;
    };

    // Add
    int refinement_iterations = 17;
    int precision = 3;
    for (int i = 1; i < refinement_iterations; i++){

        mpfr_init2(pi, precision);
        mpfr_init2(one, precision);
        mpfr_init2(e, precision);
        mpfr_init2(result, precision);
        mpfr_free_cache ();

        clock_t begin = clock();
        mpfr_const_pi(pi, GMP_RNDN);
        mpfr_exp(e, one, GMP_RNDN);
        mpfr_add(result, e, pi, GMP_RNDN);
        clock_t end = clock();

        mpfr_clear(pi);
        mpfr_clear(one);
        mpfr_clear(e);
        mpfr_clear(result);

        printf("%.2g sec at precision %d\n", convtime(begin, end), precision);
        precision *= 2;
    }


    printf("\n***===***\n");

    // Sub
    precision = 3;
    for (int i = 1; i < refinement_iterations; i++){

        mpfr_init2(pi, precision);
        mpfr_init2(one, precision);
        mpfr_init2(e, precision);
        mpfr_init2(result, precision);
        mpfr_free_cache ();

        clock_t begin = clock();
        mpfr_const_pi(pi, GMP_RNDN);
        mpfr_exp(e, one, GMP_RNDN);
        mpfr_sub(result, e, pi, GMP_RNDN);
        clock_t end = clock();

        mpfr_clear(pi);
        mpfr_clear(one);
        mpfr_clear(e);
        mpfr_clear(result);

        printf("%.2g sec at precision %d\n", convtime(begin, end), precision);
        precision *= 2;
    }
    printf("\n***===***\n");

    // Mul
    precision = 3;

    for (int i = 1; i < refinement_iterations; i++){
        mpfr_init2(pi, precision);
        mpfr_init2(one, precision);
        mpfr_init2(e, precision);
        mpfr_init2(result, precision);
        mpfr_free_cache ();

        clock_t begin = clock();
        mpfr_const_pi(pi, GMP_RNDN);
        mpfr_exp(e, one, GMP_RNDN);
        mpfr_mul(result, e, pi, GMP_RNDN);
        clock_t end = clock();

        mpfr_clear(pi);
        mpfr_clear(one);
        mpfr_clear(e);
        mpfr_clear(result);

        printf("%.2g sec at precision %d\n", convtime(begin, end), precision);
        precision *= 2;
    }

    printf("\n***===***\n");

    // Div
    precision = 3;
    for (int i = 1; i < refinement_iterations; i++){

        mpfr_init2(pi, precision);
        mpfr_init2(one, precision);
        mpfr_init2(e, precision);
        mpfr_init2(result, precision);
        mpfr_free_cache ();

        clock_t begin = clock();
        mpfr_const_pi(pi, GMP_RNDN);
        mpfr_exp(e, one, GMP_RNDN);
        mpfr_div(result, e, pi, GMP_RNDN);
        clock_t end = clock();

        mpfr_clear(pi);
        mpfr_clear(one);
        mpfr_clear(e);
        mpfr_clear(result);

        printf("%.2g sec at precision %d\n", convtime(begin, end), precision);
        precision *= 2;
    }
}