• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 /*
2  * Copyright (c) 2008-2020 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  *
11  * 2. Redistributions in binary form must reproduce the above copyright
12  *    notice, this list of conditions and the following disclaimer in the
13  *    documentation and/or other materials provided with the distribution.
14  *
15  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS" AND
16  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
19  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
20  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
21  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
22  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
23  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
24  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
25  * SUCH DAMAGE.
26  */
27 
28 
29 #include <stdio.h>
30 #include <stdlib.h>
31 #include <time.h>
32 #include <mpdecimal.h>
33 
34 
35 int
main(int argc,char ** argv)36 main(int argc, char **argv)
37 {
38 	mpd_context_t ctx;
39 	mpd_t *a, *b;
40 	mpd_t *q, *r;
41 	char *qs, *rs;
42 	char status_str[MPD_MAX_FLAG_STRING];
43 	clock_t start_clock, end_clock;
44 
45 	if (argc != 3) {
46 		fprintf(stderr, "divmod: usage: ./divmod x y\n");
47 		exit(1);
48 	}
49 
50 	mpd_init(&ctx, 38);
51 	ctx.traps = 0;
52 
53 	q = mpd_new(&ctx);
54 	r = mpd_new(&ctx);
55 	a = mpd_new(&ctx);
56 	b = mpd_new(&ctx);
57 	mpd_set_string(a, argv[1], &ctx);
58 	mpd_set_string(b, argv[2], &ctx);
59 
60 	start_clock = clock();
61 	mpd_divmod(q, r, a, b, &ctx);
62 	end_clock = clock();
63 	fprintf(stderr, "time: %f\n\n",
64 	           (double)(end_clock-start_clock)/(double)CLOCKS_PER_SEC);
65 
66 	qs = mpd_to_sci(q, 1);
67 	rs = mpd_to_sci(r, 1);
68 
69 	mpd_snprint_flags(status_str, MPD_MAX_FLAG_STRING, ctx.status);
70 	printf("%s  %s  %s\n", qs, rs, status_str);
71 
72 	mpd_del(q);
73 	mpd_del(r);
74 	mpd_del(a);
75 	mpd_del(b);
76 	mpd_free(qs);
77 	mpd_free(rs);
78 
79 	return 0;
80 }
81 
82 
83