1 /*
2 * Copyright (c) 2020-2024 Stefan Krah. All rights reserved.
3 *
4 * Redistribution and use in source and binary forms, with or without
5 * modification, are permitted provided that the following conditions
6 * are met:
7 *
8 * 1. Redistributions of source code must retain the above copyright
9 * notice, this list of conditions and the following disclaimer.
10 * 2. Redistributions in binary form must reproduce the above copyright
11 * notice, this list of conditions and the following disclaimer in the
12 * documentation and/or other materials provided with the distribution.
13 *
14 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
15 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
16 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
17 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
18 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
19 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
20 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
21 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
22 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
23 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
24 * SUCH DAMAGE.
25 */
26
27
28 #include <cstdio>
29 #include <cstdlib>
30 #include <cstdint>
31 #include <cassert>
32 #include <ctime>
33
34 #include "mpdecimal.h"
35 #include "decimal.hh"
36
37
38 using decimal::Decimal;
39 using decimal::Context;
40 using decimal::context;
41
42
43 /* Nonsense version of escape-time algorithm for calculating a Mandelbrot
44 * set. Just for benchmarking. */
45 void
color_point(Decimal & x0,const Decimal & y0,long maxiter)46 color_point(Decimal& x0, const Decimal& y0, long maxiter)
47 {
48 Decimal x = 0;
49 Decimal y = 0;
50
51 Decimal sq_x = 0;
52 Decimal sq_y = 0;
53
54 Decimal two{2};
55
56 for (long i = 0; i < maxiter; i++) {
57 y = x * y;
58 y = y * two;
59 y = y + y0;
60
61 x = sq_x - sq_y;
62 x = x + x0;
63
64 sq_x = x * x;
65 sq_y = y * y;
66 }
67
68 x0 = x;
69 }
70
71 int
main(int argc,char ** argv)72 main(int argc, char **argv)
73 {
74 const double clocks_per_sec = CLOCKS_PER_SEC;
75 clock_t start_clock, end_clock;
76 uint32_t prec;
77 long iter;
78
79 assert(MPD_MINALLOC == 4);
80
81 if (argc != 3) {
82 fprintf(stderr, "usage: bench prec iter\n");
83 exit(1);
84 }
85 prec = strtoul(argv[1], NULL, 10);
86 iter = strtol(argv[2], NULL, 10);
87
88 context.prec(prec);
89
90 Decimal x0{"0.222"};
91 Decimal y0{"0.333"};
92
93 start_clock = clock();
94 color_point(x0, y0, iter);
95 end_clock = clock();
96
97 std::cout << x0 << std::endl;
98
99 fprintf(stderr, "time: %f\n\n", (end_clock-start_clock)/clocks_per_sec);
100
101 return 0;
102 }
103
104