annotate src/decimal/bcd.d @ 7:b304232c476c

Coded and tested div and mod. Changed setNumDigits.
author Paul (paul.d.anderson@comcast.net)
date Tue, 23 Mar 2010 20:50:07 -0700
parents 14badf9104b9
children c991b9fde45c
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1 /**
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
2 * A D programming language implementation of the
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
3 * General Decimal Arithmetic Specification,
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
4 * Version 1.70, (25 March 2009).
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
5 *
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
6 * by Paul D. Anderson
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
7 *
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
8 * Boost Software License - Version 1.0 - August 17th, 2003
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
9 *
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
10 * Permission is hereby granted, free of charge, to any person or organization
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
11 * obtaining a copy of the software and accompanying documentation covered by
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
12 * this license (the "Software") to use, reproduce, display, distribute,
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
13 * execute, and transmit the Software, and to prepare derivative works of the
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
14 * Software, and to permit third-parties to whom the Software is furnished to
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
15 * do so, all subject to the following:
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
16 *
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
17 * The copyright notices in the Software and this entire statement, including
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
18 * the above license grant, this restriction and the following disclaimer,
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
19 * must be included in all copies of the Software, in whole or in part, and
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
20 * all derivative works of the Software, unless such copies or derivative
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
21 * works are solely in the form of machine-executable object code generated by
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
22 * a source language processor.
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
23 *
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
24 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
25 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
26 * FITNESS FOR A PARTICULAR PURPOSE, TITLE AND NON-INFRINGEMENT. IN NO EVENT
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
27 * SHALL THE COPYRIGHT HOLDERS OR ANYONE DISTRIBUTING THE SOFTWARE BE LIABLE
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
28 * FOR ANY DAMAGES OR OTHER LIABILITY, WHETHER IN CONTRACT, TORT OR OTHERWISE,
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
29 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
30 * DEALINGS IN THE SOFTWARE.
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
31 **/
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
32
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
33 module decimal.bcd;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
34
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
35 //import decimal.context: RoundingMode;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
36 import std.algorithm: max;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
37 import std.conv: ConvError;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
38 import std.math;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
39 import std.stdio: write, writeln;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
40 import std.string: strip;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
41
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
42 alias ubyte Digit;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
43
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
44 public immutable Bcd ZERO = { digits:[0] };
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
45 public immutable Bcd ONE = { digits:[1] };
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
46 public immutable Bcd NEG_ONE = { sign:true, digits:[1] };
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
47
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
48 /**
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
49 * Enumeration of available rounding modes.
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
50 */
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
51 public enum Rounding {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
52 HALF_EVEN,
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
53 HALF_DOWN,
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
54 HALF_UP,
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
55 DOWN,
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
56 UP,
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
57 FLOOR,
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
58 CEILING,
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
59 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
60
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
61 /**
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
62 * Provides BCD-encoded integral values of arbitrary length.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
63 *
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
64 * Advantages of BCD:
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
65 * To/from string is easy.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
66 * Easy to pull values of individual digits.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
67 * Easy to multiply and divide by powers of 10.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
68 *
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
69 * Disadvantages of BCD:
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
70 * Significantly less compact representation.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
71 * Operations are carried out bytewise, requiring more iterations.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
72 * Extra steps required for addition.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
73 *
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
74 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
75 public struct Bcd {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
76
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
77 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
78 // members
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
79 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
80
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
81 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
82 * The sign of the BCD integer. Sign is kept explicitly
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
83 * and not encoded in the number.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
84 **/
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
85 private bool sign = false;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
86
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
87 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
88 * An array of decimal digits in reverse (right-to-left) order
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
89 * representing an integral value.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
90 * Trailing zeros have no effect on the value of the number;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
91 * they correspond to leading zeros in a left-to-right representation.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
92 **/
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
93 private Digit[] digits = [ 0 ];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
94
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
95 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
96 // construction
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
97 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
98
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
99 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
100 // copy construction
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
101 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
102
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
103
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
104 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
105 * Copy constructor. Returns a (non-const) copy of the argument.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
106 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
107 public this(const Bcd bcd) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
108 sign = bcd.sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
109 digits = bcd.digits.dup;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
110 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
111
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
112 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
113 * Copies a BCD integer and adjusts the number of digits,
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
114 * padding or turncating, if necessary.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
115 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
116 public this(const Bcd bcd, uint numDigits) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
117 this(bcd);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
118 setNumDigits(numDigits);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
119 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
120
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
121 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
122 // construction from integer types.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
123 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
124
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
125 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
126 * Constructs a BCD integer from a (signed) byte, short, int or long value.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
127 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
128 public this(const long n) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
129 ulong m;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
130 if (n < 0) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
131 sign = true;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
132 m = std.math.abs(n);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
133 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
134 else {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
135 sign = false;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
136 m = n;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
137 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
138 Digit[20] dig;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
139 int i = 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
140 do {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
141 ulong q = m/10;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
142 ulong r = m%10;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
143 dig[i] = cast(Digit) r;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
144 m = q;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
145 i++;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
146 } while (m > 0);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
147 digits = dig[0..i].dup;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
148 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
149
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
150 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
151 * Constructs a BCD integer from a (signed) byte, short,
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
152 * int or long value, and then adjusts the number of digits,
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
153 * padding or truncating as needed.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
154 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
155 public this(const long n, uint numDigits) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
156 this = n;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
157 setNumDigits(numDigits);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
158 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
159
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
160 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
161 write("this(long)..");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
162 Bcd p;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
163 p = -1234;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
164 assert(p == Bcd("-1234"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
165 p = 12345678L;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
166 assert(p == Bcd("12345678"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
167 p = Bcd(-19166586400L);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
168 assert(p == Bcd("-19166586400"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
169 p = Bcd(139676498390L);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
170 assert(p == Bcd("139676498390"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
171 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
172 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
173
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
174 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
175 // construction from a string.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
176 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
177
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
178 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
179 * Constructs a BCD integer from a string representation.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
180 * If the string has leading zeros they are retained internally.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
181 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
182 public this(const string str) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
183 this = parse(str);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
184 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
185
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
186 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
187 // member access functions
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
188 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
189
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
190 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
191 // access to the sign of the number
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
192 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
193
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
194 const bool isSigned() {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
195 return sign;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
196 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
197
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
198 void setSign(bool sign) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
199 this.sign = sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
200 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
201
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
202 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
203 // access to the digits
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
204 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
205
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
206 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
207 * Returns the number of digits in this BCD integer.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
208 * Includes leading zero digits.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
209 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
210 const uint numDigits() {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
211 return digits.length;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
212 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
213
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
214 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
215 * Adjusts the number of digits in this BCD integer,
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
216 * padding or truncating if necessary.
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
217 * If truncating, first leading zeros are stripped away, then
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
218 * the remainder is clipped.
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
219 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
220 void setNumDigits(uint n) {
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
221 this = stripLeadingZeros(this);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
222 if (n > digits.length) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
223 digits.length = n;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
224 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
225 else {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
226 this = truncate(this, n);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
227 }
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
228 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
229
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
230 unittest {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
231 write("setWidth...");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
232 Bcd a;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
233 a = 1234567;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
234 a.setNumDigits(5);
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
235 assert(a == Bcd(12345));
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
236 a.setNumDigits(9);
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
237 assert(a.toString == "000012345");
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
238 writeln("passed");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
239 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
240
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
241 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
242 * Returns the first digit of this BCD integer. If the
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
243 * integer has leading zeros this function will return
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
244 * zero, but the value of the number is not necessarily zero.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
245 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
246 const uint firstDigit() {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
247 return digits[$-1];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
248 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
249
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
250 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
251 * Returns the last digit of this BCD integer.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
252 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
253 const uint lastDigit() {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
254 return digits[0];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
255 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
256
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
257 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
258 * Returns the specified digit of this BCD integer. The index is
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
259 * zero based, i.e., getDigit(0) returns the first digit.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
260 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
261 const int getDigit(int n) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
262 return digits[$-n-1];
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
263 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
264
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
265 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
266 * Sets the specified digit of this BCD integer to the specified value.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
267 * The index is zero based, i.e., getDigit(0) returns the first digit.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
268 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
269 void setDigit(uint n, Digit value) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
270 assert(value < 10);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
271 digits[$-n-1] = value;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
272 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
273
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
274 unittest {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
275 write("digits...");
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
276 Bcd bcd;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
277 bcd = 12345678;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
278 assert(bcd.numDigits() == 8);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
279 assert(bcd.firstDigit() == 1);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
280 assert(bcd.lastDigit() == 8);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
281 assert(bcd.getDigit(3) == 4);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
282 bcd.setDigit(3, 7);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
283 assert(bcd.getDigit(3) == 7);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
284 bcd.setNumDigits(10);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
285 assert(bcd.numDigits == 10);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
286 assert(bcd.firstDigit() == 0);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
287 assert(bcd.getDigit(5) == 7);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
288 bcd.setNumDigits(5);
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
289 assert(bcd.getDigit(2) == 3);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
290 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
291 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
292
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
293 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
294 // Common object functions.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
295 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
296
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
297 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
298 * Returns a string representation of this BCD integer.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
299 * Leading zeros are displayed. The sign is displayed if
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
300 * the number is negative.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
301 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
302 const string toString() {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
303 return format(this);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
304 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
305
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
306 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
307 write("toString...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
308 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
309 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
310
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
311 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
312 * Returns a hash code for this BCD integer.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
313 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
314 const hash_t toHash() {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
315 hash_t hash = 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
316 foreach(Digit digit; digits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
317 hash += digit;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
318 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
319 return hash;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
320 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
321
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
322 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
323 write("toHash...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
324 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
325 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
326
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
327 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
328 * Returns a mutable copy of this BCD integer.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
329 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
330 const Bcd dup() {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
331 Bcd bcd;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
332 bcd.sign = this.sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
333 bcd.digits = digits.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
334 return bcd;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
335 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
336
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
337 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
338 // Utility functions.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
339 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
340
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
341 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
342 * Returns true if the internal representation of
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
343 * this BCD integer has leading zeros.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
344 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
345 const bool hasLeadingZeros() {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
346 return numDigits > 0 && firstDigit == 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
347 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
348
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
349 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
350 * Returns true if the value of this BCD integer is zero.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
351 * Ignores and leading zeros and the sign of the number.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
352 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
353 const bool isZero() {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
354 Bcd temp = stripLeadingZeros(this);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
355 if (temp.numDigits > 1) return false;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
356 return temp.lastDigit == 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
357 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
358
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
359 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
360 // Assignment operators.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
361 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
362
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
363 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
364 * Assignment operator for BCD integer to BCD integer operations.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
365 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
366 Bcd opAssign(const Bcd that) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
367 this.sign = that.sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
368 this.digits = that.digits.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
369 return this;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
370 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
371
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
372 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
373 * Converts the string argument into a BCD integer and assigns it
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
374 * to this BCD integer.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
375 */
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
376 void opAssign(const string str) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
377 this = Bcd(str);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
378 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
379
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
380 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
381 * Converts the integer argument into a BCD integer and assigns it
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
382 * to this BCD integer.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
383 */
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
384 void opAssign(const long n) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
385 this = Bcd(n);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
386 /+ uint m;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
387 if (n < 0) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
388 sign = true;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
389 m = std.math.abs(n);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
390 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
391 else {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
392 sign = false;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
393 m = n;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
394 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
395 Digit[20] dig;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
396 int i = 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
397 do {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
398 uint q = m/10;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
399 uint r = m%10;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
400 dig[i] = cast(Digit) r;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
401 m = q;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
402 i++;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
403 } while (m > 0);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
404 digits = dig[0..i].dup;+/
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
405 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
406
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
407
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
408 // index operator
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
409 Digit opIndex(uint index) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
410 return this.digits[index];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
411 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
412
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
413 // index assign operator
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
414 void opIndexAssign(Digit value, uint index) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
415 this.digits[index] = value;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
416 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
417
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
418 /+ int opApply(int delegate (ref Digit) dg) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
419 int result = 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
420 for (int i = 0; i < numDigits; i++) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
421 result = dg(digits[i]);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
422 if (result) break;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
423 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
424 return result;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
425 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
426
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
427 int opApply(int delegate (ref Digit) dg) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
428 int result = 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
429 for (int i = 0; i < numDigits; i++) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
430 result = dg(digits[i]);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
431 if (result) break;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
432 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
433 return result;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
434 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
435
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
436 int opApplyReverse(int delegate (ref Digit) dg) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
437 int result = 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
438 for (int i = 0; i < numDigits; i++) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
439 result = dg(digits[i]);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
440 if (result) break;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
441 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
442 return result;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
443 }+/
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
444
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
445 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
446 // unary operators
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
447 //--------------------------------
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
448
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
449 const Bcd opPos() {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
450 return this.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
451 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
452
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
453 const Bcd opNeg() {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
454 Bcd copy = this.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
455 copy.sign = !this.sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
456 return copy;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
457 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
458
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
459 const Bcd opCom() {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
460 return twosComp(this);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
461 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
462
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
463 const Bcd opNot() {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
464 return not(this);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
465 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
466
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
467 Bcd opPostInc() {
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
468 this = this + ONE;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
469 return this;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
470 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
471
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
472 Bcd opPostDec() {
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
473 this = this - ONE;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
474 return this;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
475 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
476
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
477 const Bcd opShl(int n) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
478 return shift(this, n);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
479 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
480
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
481 Bcd opShlAssign(int n) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
482 this = this << n;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
483 return this;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
484 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
485
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
486 const Bcd opShr(int n) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
487 return shift(this, -n);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
488 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
489
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
490 Bcd opShrAssign(int n) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
491 this = this >> n;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
492 return this;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
493 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
494
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
495 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
496 // binary operators
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
497 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
498
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
499 const Bcd opAdd(T:Bcd)(const T addend) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
500 return add(this, addend);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
501 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
502
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
503 const Bcd opAdd(T)(const T addend) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
504 return add(this, Bcd(addend));
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
505 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
506
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
507 Bcd opAddAssign(T)(const T addend) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
508 this = this + addend;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
509 return this;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
510 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
511
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
512 const Bcd opSub(T:Bcd)(const T addend) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
513 return subtract(this, addend);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
514 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
515
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
516 const Bcd opSub(T)(const T addend) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
517 return subtract(this, Bcd(addend));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
518 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
519
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
520 Bcd opSubAssign(T)(const T addend) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
521 this = this - addend;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
522 return this;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
523 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
524
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
525 const Bcd opMul(T:Bcd)(const T b) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
526 return multiply(this, b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
527 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
528
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
529 const Bcd opMul(T)(const T b) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
530 return multiply(this, Bcd(b));
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
531 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
532
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
533 Bcd opMulAssign(T)(const T b) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
534 this = this * b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
535 return this;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
536 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
537
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
538 unittest {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
539 write("multiply...");
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
540 Bcd a, b, c;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
541 a = 105, b = -22;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
542 c = a * b;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
543 assert(c == -2310);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
544 c = 45;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
545 c *= 1205;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
546 // writeln("c = ", c);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
547 assert(c == 54225);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
548 writeln("passed");
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
549 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
550
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
551 const Bcd opDiv(T:Bcd)(const T b) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
552 Bcd rem;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
553 return divide(this, b, rem);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
554 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
555
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
556 unittest {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
557 write("opDiv...");
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
558 Bcd a, b;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
559 a = 7;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
560 b = 3;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
561 assert(a/b == 2);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
562 writeln("passed");
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
563 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
564
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
565 const Bcd opDiv(T)(const T b) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
566 Bcd rem;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
567 return divide(this, Bcd(b), rem);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
568 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
569
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
570 Bcd opDivAssign(T)(const T b) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
571 this = this / b;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
572 return this;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
573 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
574
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
575 const Bcd opMod(T:Bcd)(const T b) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
576 Bcd rem;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
577 divide(this, b, rem);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
578 return rem;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
579 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
580
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
581 unittest {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
582 write("opMod...");
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
583 Bcd a, b;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
584 a = 7;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
585 b = 3;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
586 assert(a % b == 1);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
587 writeln("passed");
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
588 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
589
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
590 const Bcd opMod(T)(const T b) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
591 Bcd rem;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
592 divide(this, Bcd(b), rem);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
593 return rem;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
594 }
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
595
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
596 Bcd opModAssign(T)(const T b) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
597 this = this % b;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
598 return this;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
599 }
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
600
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
601 const Bcd opAnd(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
602 return and(this, a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
603 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
604
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
605 Bcd opAndAssign(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
606 this = this & a;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
607 return this;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
608 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
609
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
610 const Bcd opOr(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
611 return or(this, a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
612 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
613
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
614 Bcd opOrAssign(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
615 this = this | a;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
616 return this;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
617 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
618
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
619 const Bcd opXor(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
620 return xor(this, a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
621 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
622
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
623 Bcd opXorAssign(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
624 this = this ^ a;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
625 return this;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
626 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
627
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
628 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
629 // comparison, equality operators
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
630 //--------------------------------
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
631
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
632 const int opCmp(T:Bcd)(const T that) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
633 return compare(this, that);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
634 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
635
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
636 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
637 write("compare...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
638 Bcd a,b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
639 a = 100;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
640 b = 5;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
641 assert(a > b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
642 a = 5;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
643 b = 100;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
644 assert(a < b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
645 assert(b > a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
646 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
647 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
648
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
649 const int opCmp(T)(const T that) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
650 return opCmp(Bcd(that));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
651 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
652
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
653 const bool opEquals(T:Bcd)(const T that) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
654 Bcd a = stripLeadingZeros(this);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
655 Bcd b = stripLeadingZeros(that);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
656 if (this.sign != that.sign) return false;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
657 if (a.digits.length != b.digits.length) return false;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
658 foreach_reverse(int i, Digit digit; a.digits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
659 if (digit != b.digits[i]) return false;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
660 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
661 return true;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
662 // return compare(this, that) == 0;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
663 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
664
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
665 const bool opEquals(T)(const T that) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
666 return opEquals(Bcd(that));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
667 }
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
668
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
669 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
670 write("equals...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
671 assert(Bcd(0) == Bcd(0));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
672 assert(Bcd(4) == Bcd(4));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
673 assert(Bcd(40) == Bcd(40));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
674 assert(Bcd(-400) == Bcd(-400));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
675 assert(Bcd(12345678) == Bcd(+12345678));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
676 assert(Bcd(40) != Bcd(53));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
677 assert(Bcd(402) != Bcd(531));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
678 assert(Bcd(402) != Bcd(-402));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
679 assert(Bcd(65432) != Bcd(65431));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
680 assert(Bcd(2) != Bcd(1));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
681 assert(Bcd(1) != Bcd(2));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
682 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
683 }
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
684
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
685 } // end struct Bcd
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
686
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
687 //public string format(const Bcd bcd, bool showPlus = false, bool showLeadingZeros = false) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
688 public string format(const Bcd bcd) {
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
689 int len = bcd.digits.length;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
690 if (bcd.isSigned) len++;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
691 // if (bcd.isSigned || bcd.hasLeadingZeros) len++;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
692 char[] str = new char[len];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
693
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
694 int index = 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
695 if (bcd.sign) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
696 str[index] = '-';
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
697 index++;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
698 }
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
699 /+ if (bcd.hasLeadingZeros) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
700 str[index] = '+';
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
701 index++;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
702 }+/
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
703 foreach_reverse(Digit digit; bcd.digits) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
704 str[index] = cast(char) digit + '0';
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
705 index++;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
706 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
707 return str.idup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
708 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
709
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
710 public Bcd parse(string str) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
711 Bcd bcd;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
712 // are leading zeros retained?
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
713 bool lz = false;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
714
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
715 // strip whitespace, convert to char array
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
716 char[] str1 = str.strip.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
717
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
718 // check for leading '-'
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
719 if (str1[0] == '-') {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
720 bcd.sign = true;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
721 str1 = str1[1..$];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
722 lz = true;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
723 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
724
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
725 // check for leading '+'
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
726 if (str1[0] == '+') {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
727 bcd.sign = false;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
728 str1 = str1[1..$];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
729 lz = true;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
730 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
731
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
732 // verify chars are digits, strip out underscores
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
733 int index = 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
734 char[] str2 = new char[str1.length];
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
735 foreach_reverse(char ch; str1) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
736 if (isDigit(ch)) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
737 str2[index] = ch;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
738 index++;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
739 continue;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
740 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
741 if (ch == '_') {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
742 continue;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
743 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
744 throw (new ConvError("Invalid character: " ~ ch));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
745 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
746 str2.length = index;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
747
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
748 bcd.digits.length = index;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
749 foreach(int i, char ch; str2) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
750 bcd.digits[i] = cast(Digit)(ch - '0');
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
751 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
752
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
753 if (!lz) return stripLeadingZeros(bcd);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
754 return bcd;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
755 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
756
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
757 private bool isDigit(const char ch) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
758 return (ch >= '0' && ch <= '9');
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
759 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
760
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
761 public Bcd negate(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
762 return -a;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
763 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
764
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
765 public Bcd abs(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
766 return a.isSigned ? -a: +a;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
767 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
768
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
769 public int compare(const Bcd m, const Bcd n) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
770 Bcd a = stripLeadingZeros(m);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
771 Bcd b = stripLeadingZeros(n);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
772 if (!a.isSigned && b.isSigned) return 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
773 if (a.isSigned && !b.isSigned) return -1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
774 if (a.digits.length > b.digits.length) return 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
775 if (a.digits.length < b.digits.length) return -1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
776 foreach_reverse(int i, Digit digit; a.digits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
777 if (digit > b.digits[i]) return 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
778 if (digit < b.digits[i]) return -1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
779 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
780 return 0;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
781 }
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
782
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
783 // TODO: does this need a const and a dup?
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
784 // should it be public?? private would be okay
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
785 public Bcd stripLeadingZeros(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
786 Bcd d = a.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
787 if (!a.hasLeadingZeros) return d;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
788 int len = a.numDigits;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
789 int i = 0;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
790 while(i < len-1 && a.getDigit(i) == 0) {
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
791 i++;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
792 }
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
793 d.digits.length = (len-i); // - 1);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
794 if (d.numDigits == 1 && d.firstDigit == 0) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
795 d.sign = false;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
796 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
797 return d;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
798 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
799
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
800 unittest {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
801 write("strip...");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
802 Bcd bcd;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
803 bcd = "+00123";
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
804 assert(bcd.toString == "00123");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
805 bcd = stripLeadingZeros(bcd);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
806 // writeln("bcd = ", bcd);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
807 assert(bcd.toString == "123");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
808 bcd = "+000";
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
809 // writeln("bcd = ", bcd);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
810 assert(bcd.toString == "000");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
811 bcd = stripLeadingZeros(bcd);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
812 // writeln("bcd = ", bcd);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
813 writeln("passed");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
814 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
815
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
816 public bool sameLength(const Bcd a, const Bcd b) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
817 return a.numDigits == b.numDigits;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
818 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
819
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
820 public int setSameLength(ref Bcd a, ref Bcd b) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
821 if (sameLength(a, b)) return a.numDigits;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
822 uint alen = a.numDigits;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
823 uint blen = b.numDigits;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
824 if (alen > blen) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
825 b.setNumDigits(alen);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
826 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
827 else {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
828 a.setNumDigits(blen);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
829 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
830 return a.numDigits();
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
831 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
832
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
833 public Bcd twosComp(const Bcd a) {
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
834 Bcd b = Bcd(0, a.numDigits);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
835 foreach(int i, Digit digit; a.digits) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
836 b.digits[i] = digit == 0 ? 1 : 0 ;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
837 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
838 return b;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
839 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
840
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
841 unittest {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
842 write("com...");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
843 Bcd bcd;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
844 bcd = 123;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
845 assert(twosComp(bcd).toString == "000");
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
846 bcd = 40509;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
847 assert(twosComp(bcd).toString == "01010");
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
848 // writeln("bcd = ", bcd);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
849 // writeln("~bcd = ", twosComp(bcd));
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
850 writeln("passed");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
851 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
852
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
853 public int sgn(const Bcd bcd) {
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
854 if (bcd.isZero) return 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
855 return bcd.isSigned ? -1 : 1;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
856 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
857
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
858 public Bcd not(const Bcd x) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
859 Bcd result = cast(Bcd) x;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
860 for (int i = 0; i < x.numDigits; i++) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
861 result[i] = not(result[i]);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
862 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
863 return result;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
864 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
865
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
866 private Digit not(const Digit a) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
867 return a != 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
868 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
869
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
870 public Bcd and(const Bcd x, const Bcd y) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
871 Bcd a = x.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
872 Bcd b = y.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
873 int len = setSameLength(a,b);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
874 Bcd result = Bcd(0, len);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
875 for (int i = 0; i < len; i++) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
876 result[i] = and(a[i], b[i]);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
877 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
878 return result;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
879 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
880
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
881 private Digit and(const Digit a, const Digit b) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
882 return a != 0 && b != 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
883 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
884
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
885 unittest {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
886 write("and...");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
887 Bcd a;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
888 Bcd b;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
889 a = 123000;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
890 b = 123000;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
891 assert(and(a,b).toString == "111000");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
892 b = 12300;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
893 assert(and(a,b).toString == "011000");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
894 a = 1234567;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
895 b = 7654321;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
896 assert(and(a,b).toString == "1111111");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
897 a = 1010;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
898 b = 101;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
899 // writeln("a = ", a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
900 // writeln("b = ", b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
901 // writeln("and = ", and(a,b));
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
902 assert(and(a,b).isZero);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
903 writeln("passed");
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
904 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
905
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
906 // TODO: looks like there's some room for templates or mixins here.
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
907 public Bcd or(const Bcd x, const Bcd y) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
908 Bcd a = x.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
909 Bcd b = y.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
910 int len = setSameLength(a,b);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
911 Bcd result = Bcd(0, len);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
912 for (int i = 0; i < len; i++) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
913 result[i] = or(a[i], b[i]);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
914 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
915 return result;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
916 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
917
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
918 private Digit or(const Digit a, const Digit b) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
919 return a != 0 || b != 0;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
920 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
921
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
922 public Bcd xor(const Bcd x, const Bcd y) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
923 Bcd a = x.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
924 Bcd b = y.dup;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
925 int len = setSameLength(a,b);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
926 Bcd result = Bcd(0, len);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
927 for (int i = 0; i < len; i++) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
928 result[i] = xor(a[i], b[i]);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
929 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
930 return result;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
931 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
932
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
933 private Digit xor(const Digit a, const Digit b) {
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
934 return (a == 0 && b != 0) || (a != 0 && b == 0);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
935 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
936
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
937 // TODO: watch the leading zeros
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
938 public Bcd add(const Bcd x, const Bcd y) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
939 Bcd a = x.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
940 Bcd b = y.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
941 Bcd sum;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
942 if (a.sign == b.sign) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
943 sum = simpleAdd(a, b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
944 sum = stripLeadingZeros(sum);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
945 sum.sign = a.sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
946 return sum;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
947 }
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
948 else { // signs differ
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
949 bool sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
950 switch (compare(abs(a), abs(b))) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
951 case 0:
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
952 return ZERO.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
953 case 1:
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
954 sign = a.isSigned;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
955 break;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
956 case -1:
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
957 sign = b.isSigned;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
958 break;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
959 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
960 setSameLength(a, b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
961 sum = simpleAdd(a, tensComp(b));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
962 if (sum.firstDigit == 1) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
963 sum.digits = sum.digits[0..$-1];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
964 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
965 else {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
966 sum = tensComp(sum);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
967 sum = stripLeadingZeros(sum);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
968 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
969 sum.sign = sign;
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
970 }
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
971 return stripLeadingZeros(sum);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
972 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
973
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
974 unittest {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
975 write("add....");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
976 Bcd a, b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
977 a = 1234;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
978 b = 4321;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
979 assert(a + b == 5555);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
980 a = 2000;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
981 b = 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
982 assert(a + b == 2001);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
983 a = -2000;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
984 b = -1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
985 assert(a + b == -2001);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
986 a = 2000;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
987 b = -1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
988 assert(a + b == 1999);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
989 a = -123;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
990 b = 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
991 assert(a + b == -122);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
992 a = 123;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
993 b = -1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
994 assert(a + b == 122);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
995 a = -1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
996 b = 123;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
997 // writeln();
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
998 // writeln("a + b = ", a + b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
999 assert(a + b == 122);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1000 a = -123;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1001 b = 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1002 assert(a + b == -122);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1003 a = 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1004 b = -123;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1005 assert(a + b == -122);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1006 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1007 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1008
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1009 Bcd subtract(const Bcd a, const Bcd b) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1010 return add(a, negate(b));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1011 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1012
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1013 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1014 * Adds two BCD integers without regard to sign
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1015 **/
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1016 private Bcd simpleAdd(const Bcd a, const Bcd b) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1017 Bcd x = a.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1018 Bcd y = b.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1019 int len = setSameLength(x, y);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1020 Bcd sum = Bcd(0, len+1);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1021 Digit carry = 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1022 uint i;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1023 // TODO: this is the place to stop adding zeros and
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1024 // just copy the rest of the digits.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1025 for (i = 0; i < len; i++) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1026 sum[i] = add(x[i], y[i], carry);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1027 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1028 if (carry) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1029 sum[i] = 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1030 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1031 else {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1032 sum = stripLeadingZeros(sum);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1033 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1034 return sum;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1035 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1036
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1037 /**
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1038 * Adds two digits and a carry digit.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1039 * Returns the (single-digit) sum and sets or resets the carry digit.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1040 **/
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1041 private Digit add(const Digit a, const Digit b, ref Digit carry) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1042 Digit sum = a + b + carry;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1043 carry = sum > 9;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1044 if (carry) sum -= 10;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1045 return sum;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1046 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1047
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1048 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1049 write("simpleAdd...");
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1050 Bcd a, b;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1051 a = 123;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1052 b = 222;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1053 Bcd sum = simpleAdd(a,b);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1054 assert(sum == 345);
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1055 a = 2;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1056 b = 102000;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1057 sum = simpleAdd(a,b);
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1058 assert(sum == 102002);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1059 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1060 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1061
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1062 private Bcd tensComp(const Bcd a) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1063 /+ foreach(Digit digit; a.digits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1064 digit = 9 - digit;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1065 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1066 a = simpleAdd(a, ONE);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1067 return a;+/
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1068
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1069 Bcd x = Bcd(0, a.numDigits);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1070 foreach (int i, digit; a.digits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1071 x[i] = 9 - digit;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1072 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1073 return simpleAdd(x, ONE);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1074 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1075
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1076 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1077 write("tensComp...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1078 Bcd a, b, c;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1079 a = 123456;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1080 b = tensComp(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1081 // writeln("a = ", a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1082 // writeln("b = ", b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1083 c = 876544;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1084 assert(b == c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1085 a = Bcd(0, 4);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1086 b = tensComp(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1087 c = 10000;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1088 assert(b == c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1089 a = Bcd(1, 4);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1090 b = tensComp(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1091 c = 9999;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1092 assert(b == c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1093 a = 9999;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1094 b = tensComp(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1095 c = 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1096 assert(b == c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1097 a = 10000;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1098 b = tensComp(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1099 c = 90000;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1100 assert(b == c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1101 a = 1;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1102 b = tensComp(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1103 c = 9;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1104 // writeln("a = ", a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1105 // writeln("b = ", b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1106 assert(b == c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1107 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1108 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1109
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1110 // TODO: there are lots of ways to speed this up.
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1111 public Bcd multiply(const Bcd a, const Bcd b) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1112 Bcd n, m;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1113 Bcd p = Bcd(0, a.numDigits + b.numDigits);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1114 if (a.numDigits > b.numDigits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1115 n = stripLeadingZeros(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1116 m = stripLeadingZeros(b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1117 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1118 else {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1119 n = stripLeadingZeros(b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1120 m = stripLeadingZeros(a);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1121 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1122 foreach(uint i, Digit d; m.digits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1123 if (d == 0) continue;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1124 // writeln("mul(n, m[i]) = ", mul(n, m[i]));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1125 // writeln("pow10(mul(n, m[i]), i) = ", pow10(mul(n, m[i]), i));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1126 p = simpleAdd(p, pow10(mul(n, m[i]), i));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1127 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1128 p.sign = a.sign ^ b.sign;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1129 return p;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1130 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1131
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1132 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1133 write("mul...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1134 Bcd a, b, p;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1135 a = 2105, b = 301;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1136 p = multiply(a, b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1137 assert(p == Bcd(633605));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1138 b = 800, a = 23958233;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1139 p = multiply(a, b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1140 assert(p == Bcd("19166586400"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1141 b = 5830, a = 23958233;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1142 p = multiply(a, b);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1143 assert(p == Bcd("139676498390"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1144 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1145 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1146
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1147 private Bcd mul(const Bcd a, Digit m) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1148 Bcd p = Bcd(0, a.numDigits+1);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1149 Digit c = 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1150 int i = 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1151 foreach (Digit d; a.digits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1152 p[i] = mul(m, d, c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1153 i++;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1154 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1155 if (c) p[i] = c;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1156 return p;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1157 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1158
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1159 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1160 write("mul...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1161 Bcd a, p;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1162 Digit m;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1163 a = 2105, m = 3;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1164 p = mul(a, m);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1165 assert(p == Bcd(6315));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1166 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1167 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1168
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1169 private Digit mul(Digit a, Digit b, ref Digit c) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1170 Digit p = a * b + c;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1171 c = p / 10;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1172 return p % 10;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1173 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1174
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1175 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1176 write("mul...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1177 Digit a, b, c, d;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1178 a = 2, b = 3, c = 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1179 d = mul(a,b,c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1180 assert(d == 6);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1181 assert(c == 0);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1182 a = 9, b = 8, c = 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1183 d = mul(a,b,c);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1184 assert(d == 2);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1185 assert(c == 7);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1186 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1187 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1188
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1189 public Bcd truncate(const Bcd a, int n) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1190 /+ if (n <= 0) return ZERO.dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1191 if (n >= a.numDigits) return a.dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1192 Bcd b;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1193 b.digits = a.digits[$-n..$].dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1194 b.sign = a.sign;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1195 return b;+/
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1196 Bcd dummy;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1197 return truncate(a, n, dummy);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1198 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1199
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1200 unittest {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1201 write("truncate....");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1202 Bcd a;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1203 a = 1234567;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1204 assert(truncate(a,3) == 123);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1205 a = 10256;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1206 assert(truncate(a,5) == 10256);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1207 a = 8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1208 assert(truncate(a,1) == 8);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1209 a = -8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1210 assert(truncate(a,1) == -8);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1211 a = -8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1212 assert(truncate(a,-1) == 0);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1213 writeln("passed");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1214 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1215
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1216 /**
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1217 * Truncate to the specified number of digits
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1218 */
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1219 public Bcd truncate(const Bcd a, int n, out Bcd r) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1220 if (n <= 0) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1221 r = stripLeadingZeros(a);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1222 return ZERO.dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1223 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1224 if (n >= a.numDigits) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1225 r = ZERO.dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1226 return stripLeadingZeros(a);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1227 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1228 Bcd b;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1229 b.digits = a.digits[$-n..$].dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1230 b.sign = a.sign;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1231 r.digits = a.digits[0..$-n].dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1232 return b;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1233 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1234
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1235 unittest {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1236 write("truncRem...");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1237 Bcd a;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1238 Bcd r;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1239 a = 1234567;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1240 assert(truncate(a, 3, r) == 123);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1241 assert(r == 4567);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1242 a = 10256;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1243 assert(truncate(a,5,r) == 10256);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1244 assert(r == 0);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1245 a = 8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1246 assert(truncate(a,1,r) == 8);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1247 assert(r == 500);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1248 a = -8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1249 assert(truncate(a,1,r) == -8);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1250 assert(r == 500);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1251 a = -8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1252 // writeln("truncate(a, 6, r) = ", truncate(a, 6, r));
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1253 // writeln("r = ", r);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1254 assert(truncate(a,6,r) == -8500);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1255 assert(r == 0);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1256 writeln("passed");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1257 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1258
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1259 public Bcd pad(Bcd a, int n) {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1260 if (n <= a.numDigits) return a.dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1261 return Bcd(a,n);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1262 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1263
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1264 unittest {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1265 write("pad.........");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1266 Bcd a;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1267 a = 1234567;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1268 assert(pad(a,3) == 1234567);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1269 a = 10256;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1270 // writeln("pad(a,9) = ", pad(a,9));
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1271 assert(pad(a,9).toString == "000010256");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1272 a = 8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1273 assert(pad(a,6).toString == "008500");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1274 a = -8500;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1275 assert(pad(a,6).toString == "-008500");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1276 writeln("passed");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1277 }
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1278 public Bcd divide(const Bcd a, const Bcd b, out Bcd r) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1279
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1280 if (b == ZERO) throw new Exception("Divide by zero");
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1281 if (a == ZERO) return ZERO.dup;
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1282
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1283 Bcd d = stripLeadingZeros(abs(b));
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1284 r = stripLeadingZeros(abs(a));
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1285 Bcd q;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1286 Bcd p = d;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1287 Bcd m;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1288
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1289 // shift the divisor
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1290 do {
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1291 m = p;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1292 p <<= 1;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1293 } while (p <= r && p > m);
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1294
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1295 while (m >= d) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1296 q <<= 1;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1297 // subtract repeatedly
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1298 while (m <= r) {
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1299 r -= m;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1300 q++;
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1301 }
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1302 m >>= 1;
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1303 }
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1304
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1305 // adjust the sign
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1306 r.sign = a.sign;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1307 q.sign = a.sign != b.sign;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1308 return q;
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1309 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1310
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1311 unittest {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1312 write("divide...");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1313 Bcd a, b, q, r;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1314 a = 144;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1315 b = 12;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1316 q = divide(a, b, r);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1317 assert(q == 12);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1318 assert(r == 0);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1319 a = 144;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1320 b = 14;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1321 q = divide(a, b, r);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1322 assert(q == 10);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1323 assert(r == 4);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1324 a = 142351;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1325 b = 12;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1326 q = divide(a, b, r);
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1327 // writeln("q = ", q);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1328 // writeln("r = ", r);
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1329 assert(q == 11862);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1330 assert(r == 7);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1331 a = 2;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1332 b = 12;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1333 q = divide(a, b, r);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1334 assert(q == 0);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1335 assert(r == 2);
7
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1336 a = 7;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1337 b = 3;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1338 q = divide(a, b, r);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1339 assert(q == 2);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1340 assert(r == 1);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1341 a = 7;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1342 b = -3;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1343 q = divide(a, b, r);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1344 assert(q == -2);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1345 assert(r == 1);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1346 a = -7;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1347 b = 3;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1348 q = divide(a, b, r);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1349 assert(q == -2);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1350 assert(r == -1);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1351 a = -7;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1352 b = -3;
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1353 q = divide(a, b, r);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1354 assert(q == 2);
b304232c476c Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents: 6
diff changeset
1355 assert(r == -1);
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1356 writeln("passed");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1357 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1358
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1359
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1360 public Bcd pow10(const Bcd a, uint n) {
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1361 Bcd b = stripLeadingZeros(a);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1362 if (n == 0 || b == 0) return b;
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1363 Digit[] zeros = new Digit[n];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1364 b.digits = zeros ~ b.digits;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1365 return b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1366 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1367
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1368 public Bcd pow10(uint n) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1369 return pow10(ONE, n);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1370 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1371
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1372 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1373 write("pow10...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1374 Bcd a, b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1375 a = 21345;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1376 a = pow10(a, 12);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1377 assert(a == Bcd("21345000000000000"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1378 a = pow10(a, 0);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1379 assert(a == "21345000000000000");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1380 a = pow10(9);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1381 assert(a == Bcd("1000000000"));
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1382 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1383 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1384
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1385 public Bcd shift(const Bcd a, int n) {
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1386 if (n == 0) return stripLeadingZeros(a);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1387 if (n > 0) { // shift left
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1388 return pow10(a, n);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1389 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1390 else { // shift right
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1391 if (-n >= a.numDigits) return ZERO.dup;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1392 return truncate(a, a.numDigits + n);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1393 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1394 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1395
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1396 unittest {
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1397 write("shift.......");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1398 Bcd a, b;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1399 a = 1234;
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1400 b = shift(a, 1);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1401 assert(b.toString == "12340");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1402 b = shift(a, 2);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1403 assert(b.toString == "123400");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1404 b = shift(a, -1);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1405 // writeln("b = ", b);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1406 assert(b.toString == "123");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1407 b = shift(a, 0);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1408 assert(b.toString == "1234");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1409 b = shift(a, -12);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1410 /+ writeln("b = ", b);+/
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1411 assert(b.toString == "0");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1412 b = shift(a, -4);
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1413 assert(b.toString == "0");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1414 writeln("passed");
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1415 }
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1416
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1417 public Bcd shiftFixed(const Bcd a, int n) {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1418 if (n == 0) return a.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1419 Bcd b = Bcd(0, a.numDigits);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1420 if (std.math.abs(n) >= a.numDigits) return b;
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1421 if (n > 0) { // shiftFixed left
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1422 b.digits[0..$-n] = a.digits[n..$];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1423 }
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1424 else { // shiftFixed right
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1425 n = -n;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1426 b.digits[n..$] = a.digits[0..$-n];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1427 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1428 return b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1429 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1430
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1431 unittest {
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1432 write("shiftFixed...");
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1433 Bcd a, b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1434 a = 1234;
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1435 b = shiftFixed(a, 1);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1436 assert(b.toString == "0123");
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1437 b = shiftFixed(a, 2);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1438 assert(b.toString == "0012");
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1439 b = shiftFixed(a, -1);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1440 assert(b.toString == "2340");
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1441 b = shiftFixed(a, 0);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1442 assert(b.toString == "1234");
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1443 b = shiftFixed(a, 12);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1444 assert(b.toString == "0000");
6
14badf9104b9 shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents: 5
diff changeset
1445 b = shiftFixed(a, -4);
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1446 assert(b.toString == "0000");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1447 writeln("passed");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1448 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1449
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1450 public Bcd rotate(const Bcd a, int n) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1451 if (n == 0) return a.dup;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1452 Bcd b = Bcd(0, a.numDigits);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1453 if (std.math.abs(n) >= a.numDigits) {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1454 bool sign = n < 0;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1455 n = std.math.abs(n) % a.numDigits;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1456 if (sign) n = -n;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1457 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1458 if (n > 0) { // rotate left
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1459 b.digits[0..$-n] = a.digits[n..$];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1460 b.digits[$-n..$] = a.digits[0..n];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1461 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1462 else { // rotate right
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1463 n = -n;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1464 b.digits[n..$] = a.digits[0..$-n];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1465 b.digits[0..n] = a.digits[$-n..$];
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1466 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1467 return b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1468 }
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1469
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1470 unittest {
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1471 write("rotate...");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1472 Bcd a, b;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1473 a = 1234;
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1474 b = rotate(a, 1);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1475 assert(b.toString == "4123");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1476 b = rotate(a, 2);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1477 assert(b.toString == "3412");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1478 b = rotate(a, -1);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1479 assert(b.toString == "2341");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1480 b = rotate(a, 0);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1481 assert(b.toString == "1234");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1482 b = rotate(a, 12);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1483 assert(b.toString == "1234");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1484 b = rotate(a, -4);
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1485 assert(b.toString == "1234");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1486 writeln("passed");
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1487 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1488
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1489 //==========================================
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1490
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1491 public void main() {
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1492 writeln();
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1493 writeln("If you got this far all the unit tests were successful.");
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1494 writeln();
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1495 writeln(" Congratulations!");
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1496 Bcd bcd;
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1497 writeln("bcd = ", format(bcd));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1498 bcd = parse("12");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1499 writeln("bcd = ", format(bcd));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1500 bcd = parse("012");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1501 writeln("bcd = ", format(bcd));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1502 bcd = parse("-012");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1503 writeln("bcd = ", format(bcd));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1504 bcd = parse("+012");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1505 writeln("bcd = ", format(bcd));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1506 bcd = parse("-0");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1507 writeln("bcd = ", format(bcd));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1508 bcd = parse("012_345_678");
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1509 writeln("bcd = ", format(bcd));
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1510
5
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1511 /+
c021ed211f89 bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents: 4
diff changeset
1512 writeln("OK!");+/
4
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1513 }
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1514
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1515
b37c218c1442 Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff changeset
1516