1 /* divmod.c: bcmath library file. */
2 /*
3 Copyright (C) 1991, 1992, 1993, 1994, 1997 Free Software Foundation, Inc.
4 Copyright (C) 2000 Philip A. Nelson
5
6 This library is free software; you can redistribute it and/or
7 modify it under the terms of the GNU Lesser General Public
8 License as published by the Free Software Foundation; either
9 version 2 of the License, or (at your option) any later version.
10
11 This library is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 Lesser General Public License for more details. (LICENSE)
15
16 You should have received a copy of the GNU Lesser General Public
17 License along with this library; if not, write to:
18
19 The Free Software Foundation, Inc.
20 59 Temple Place, Suite 330
21 Boston, MA 02111-1307 USA.
22
23 You may contact the author by:
24 e-mail: philnelson@acm.org
25 us-mail: Philip A. Nelson
26 Computer Science Department, 9062
27 Western Washington University
28 Bellingham, WA 98226-9062
29
30 *************************************************************************/
31
32 #include <config.h>
33 #include <stdio.h>
34 #include <stdlib.h>
35 #include <ctype.h>
36 #include <stdarg.h>
37 #include "bcmath.h"
38 #include "private.h"
39
40
41 /* Division *and* modulo for numbers. This computes both NUM1 / NUM2 and
42 NUM1 % NUM2 and puts the results in QUOT and REM, except that if QUOT
43 is NULL then that store will be omitted.
44 */
45
46 int
bc_divmod(bc_num num1,bc_num num2,bc_num * quot,bc_num * rem,int scale)47 bc_divmod (bc_num num1, bc_num num2, bc_num *quot, bc_num *rem, int scale)
48 {
49 bc_num quotient = NULL;
50 bc_num temp;
51 int rscale;
52
53 /* Check for correct numbers. */
54 if (bc_is_zero (num2)) return -1;
55
56 /* Calculate final scale. */
57 rscale = MAX (num1->n_scale, num2->n_scale+scale);
58 bc_init_num(&temp);
59
60 /* Calculate it. */
61 bc_divide (num1, num2, &temp, 0);
62 if (quot)
63 quotient = bc_copy_num (temp);
64 bc_multiply (temp, num2, &temp, rscale);
65 bc_sub (num1, temp, rem, rscale);
66 bc_free_num (&temp);
67
68 if (quot)
69 {
70 bc_free_num (quot);
71 *quot = quotient;
72 }
73
74 return 0; /* Everything is OK. */
75 }
76
77
78 /* Modulo for numbers. This computes NUM1 % NUM2 and puts the
79 result in RESULT. */
80
81 int
bc_modulo(bc_num num1,bc_num num2,bc_num * result,int scale)82 bc_modulo (bc_num num1, bc_num num2, bc_num *result, int scale)
83 {
84 return bc_divmod (num1, num2, NULL, result, scale);
85 }
86
