Mercurial > projects > decimal
annotate src/decimal/bcd.d @ 11:8e5f172ec55c
added comparable
author | Paul (paul.d.anderson@comcast.net) |
---|---|
date | Thu, 25 Mar 2010 21:37:14 -0700 |
parents | f925fe996255 |
children |
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 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
35 import std.algorithm: max; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
36 import std.conv: ConvError; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
37 import std.math; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
38 import std.stdio: write, writeln; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
39 import std.string: strip; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
40 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
41 alias ubyte Digit; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
42 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
43 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
|
44 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
|
45 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
|
46 |
10 | 47 /+static this() { |
48 immutable(Bcd) MAX_LONG = cast(immutable) Bcd(long.max); | |
49 // MAX_LONG = cast(immutable) Bcd(long.max); | |
50 }+/ | |
51 | |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
52 /** |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
53 * Enumeration of available rounding modes. |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
54 */ |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
55 public enum Rounding { |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
56 HALF_EVEN, |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
57 HALF_DOWN, |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
58 HALF_UP, |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
59 DOWN, |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
60 UP, |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
61 FLOOR, |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
62 CEILING, |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
63 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
64 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
65 /** |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
66 * 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
|
67 * |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
68 * Advantages of BCD: |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
69 * To/from string is easy. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
70 * 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
|
71 * 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
|
72 * |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
73 * Disadvantages of BCD: |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
74 * Significantly less compact representation. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
75 * 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
|
76 * Extra steps required for addition. |
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 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
79 public struct Bcd { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
80 |
5
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 // members |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
83 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
84 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
85 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
86 * 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
|
87 * 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
|
88 **/ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
89 private bool sign = false; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
90 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
91 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
92 * 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
|
93 * representing an integral value. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
94 * 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
|
95 * 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
|
96 **/ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
97 private Digit[] digits = [ 0 ]; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
98 |
5
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 // 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 // copy construction |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
105 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
106 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
107 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
108 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
109 * 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
|
110 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
111 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
|
112 sign = bcd.sign; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
113 digits = bcd.digits.dup; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
114 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
115 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
116 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
117 * 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
|
118 * padding or turncating, if necessary. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
119 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
120 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
|
121 this(bcd); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
122 setNumDigits(numDigits); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
123 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
124 |
5
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 // construction from integer types. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
127 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
128 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
129 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
130 * 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
|
131 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
132 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
|
133 ulong m; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
134 if (n < 0) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
135 sign = true; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
136 m = std.math.abs(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 else { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
139 sign = false; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
140 m = n; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
141 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
142 Digit[20] dig; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
143 int i = 0; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
144 do { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
145 ulong q = m/10; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
146 ulong r = m%10; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
147 dig[i] = cast(Digit) r; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
148 m = q; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
149 i++; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
150 } while (m > 0); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
151 digits = dig[0..i].dup; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
152 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
153 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
154 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
155 * 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
|
156 * 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
|
157 * padding or truncating as needed. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
158 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
159 public this(const long n, uint numDigits) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
160 this = n; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
161 setNumDigits(numDigits); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
162 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
163 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
164 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
165 write("this(long).."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
166 Bcd p; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
167 p = -1234; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
168 assert(p == Bcd("-1234")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
169 p = 12345678L; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
170 assert(p == Bcd("12345678")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
171 p = Bcd(-19166586400L); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
172 assert(p == Bcd("-19166586400")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
173 p = Bcd(139676498390L); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
174 assert(p == Bcd("139676498390")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
175 writeln("passed"); |
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 // construction from a string. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
180 //-------------------------------- |
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 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
183 * 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
|
184 * 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
|
185 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
186 public this(const string str) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
187 this = parse(str); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
188 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
189 |
5
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 // member access functions |
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 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
195 // 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
|
196 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
197 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
198 const bool isSigned() { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
199 return sign; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
200 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
201 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
202 void setSign(bool sign) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
203 this.sign = sign; |
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 // access to the digits |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
208 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
209 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
210 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
211 * 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
|
212 * Includes leading zero digits. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
213 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
214 const uint numDigits() { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
215 return digits.length; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
216 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
217 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
218 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
219 * 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
|
220 * padding or truncating if necessary. |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
221 * 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
|
222 * the remainder is clipped. |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
223 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
224 void setNumDigits(uint n) { |
10 | 225 if (n == 0) n++; |
226 digits.length = n; | |
227 } | |
228 | |
229 /+ /** | |
230 * Adjusts the number of digits in this BCD integer, | |
231 * padding or truncating if necessary. | |
232 * If truncating, first leading zeros are stripped away, then | |
233 * the remainder is clipped. | |
234 */ | |
235 void setNumDigits(uint n) { | |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
236 stripl(this); |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
237 if (n > digits.length) { |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
238 digits.length = n; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
239 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
240 else { |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
241 this = truncate(this, n); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
242 } |
10 | 243 }+/ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
244 |
10 | 245 /+unittest { |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
246 write("setWidth..."); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
247 Bcd a; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
248 a = 1234567; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
249 a.setNumDigits(5); |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
250 assert(a == Bcd(12345)); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
251 a.setNumDigits(9); |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
252 assert(a.toString == "000012345"); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
253 writeln("passed"); |
10 | 254 }+/ |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
255 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
256 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
257 * 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
|
258 * 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
|
259 * 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
|
260 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
261 const uint firstDigit() { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
262 return digits[$-1]; |
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 * 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
|
267 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
268 const uint lastDigit() { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
269 return digits[0]; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
270 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
271 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
272 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
273 * 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
|
274 * 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
|
275 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
276 const int getDigit(int n) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
277 return digits[$-n-1]; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
278 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
279 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
280 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
281 * 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
|
282 * 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
|
283 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
284 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
|
285 assert(value < 10); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
286 digits[$-n-1] = value; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
287 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
288 |
10 | 289 /+ unittest { |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
290 write("digits..."); |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
291 Bcd bcd; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
292 bcd = 12345678; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
293 assert(bcd.numDigits() == 8); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
294 assert(bcd.firstDigit() == 1); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
295 assert(bcd.lastDigit() == 8); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
296 assert(bcd.getDigit(3) == 4); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
297 bcd.setDigit(3, 7); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
298 assert(bcd.getDigit(3) == 7); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
299 bcd.setNumDigits(10); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
300 assert(bcd.numDigits == 10); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
301 assert(bcd.firstDigit() == 0); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
302 assert(bcd.getDigit(5) == 7); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
303 bcd.setNumDigits(5); |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
304 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
|
305 writeln("passed"); |
10 | 306 }+/ |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
307 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
308 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
309 // Common object functions. |
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 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
313 * 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
|
314 * 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
|
315 * the number is negative. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
316 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
317 const string toString() { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
318 return format(this); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
319 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
320 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
321 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
322 write("toString..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
323 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
324 } |
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 * 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
|
328 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
329 const hash_t toHash() { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
330 hash_t hash = 0; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
331 foreach(Digit digit; digits) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
332 hash += digit; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
333 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
334 return hash; |
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 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
338 write("toHash..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
339 writeln("passed"); |
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 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
343 * 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
|
344 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
345 const Bcd dup() { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
346 Bcd bcd; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
347 bcd.sign = this.sign; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
348 bcd.digits = digits.dup; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
349 return bcd; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
350 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
351 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
352 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
353 // Utility functions. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
354 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
355 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
356 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
357 * 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
|
358 * 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
|
359 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
360 const bool hasLeadingZeros() { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
361 return numDigits > 1 && firstDigit == 0; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
362 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
363 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
364 unittest { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
365 write("hasLeadingZeros..."); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
366 assert(!ZERO.hasLeadingZeros); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
367 Bcd b = Bcd(0); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
368 assert(!b.hasLeadingZeros); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
369 b = "-00000001"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
370 assert(b.hasLeadingZeros); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
371 b = "-00000000"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
372 assert(b.hasLeadingZeros); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
373 writeln("passed"); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
374 } |
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 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
377 * Returns true if the value of this BCD integer is zero. |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
378 * Ignores leading zeros and sign. |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
379 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
380 const bool isZero() { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
381 Bcd b = canonical(this); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
382 if (b.numDigits > 1) return false; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
383 return b.lastDigit == 0; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
384 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
385 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
386 unittest { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
387 write("isZero........."); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
388 assert(ZERO.isZero); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
389 Bcd b = Bcd(0); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
390 assert(b.isZero); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
391 b = "-00000001"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
392 assert(!b.isZero); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
393 b = "-00000000"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
394 assert(b.isZero); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
395 writeln("passed"); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
396 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
397 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
398 //-------------------------------- |
10 | 399 // casting operators. |
400 //-------------------------------- | |
401 | |
402 const long opCast() { | |
403 Bcd MAX_LONG = Bcd(long.max); | |
404 Bcd MIN_LONG = Bcd(long.min); | |
405 | |
406 if (this > MAX_LONG || this < MIN_LONG) | |
407 throw new Exception("Can't cast -- out of range"); | |
408 | |
409 long n = 0; | |
410 foreach_reverse(Digit digit; digits) { | |
411 n = 10*n + digit; | |
412 } | |
413 return this.sign ? -n : n; | |
414 } | |
415 | |
416 unittest { | |
417 write("cast(long)..."); | |
418 Bcd a; | |
419 a = "12345678901234567890123456"; | |
420 // long n = cast(long) a; | |
421 a = 12345678L; | |
422 long n = cast(long) a; | |
423 assert(n == 12345678L); | |
424 int m = cast(int) a; | |
425 assert(m == 12345678L); | |
426 a = -a; | |
427 m = cast(int) a; | |
428 assert(m == -12345678L); | |
429 | |
430 writeln("passed"); | |
431 } | |
432 | |
433 | |
434 //-------------------------------- | |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
435 // Assignment operators. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
436 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
437 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
438 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
439 * 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
|
440 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
441 Bcd opAssign(const Bcd that) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
442 this.sign = that.sign; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
443 this.digits = that.digits.dup; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
444 return this; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
445 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
446 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
447 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
448 * 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
|
449 * to this BCD integer. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
450 */ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
451 void opAssign(const string str) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
452 this = Bcd(str); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
453 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
454 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
455 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
456 * 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
|
457 * to this BCD integer. |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
458 */ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
459 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
|
460 this = Bcd(n); |
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 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
464 // index operator |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
465 Digit opIndex(uint index) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
466 return this.digits[index]; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
467 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
468 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
469 // index assign operator |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
470 void opIndexAssign(Digit value, uint index) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
471 this.digits[index] = value; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
472 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
473 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
474 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
475 // unary operators |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
476 //-------------------------------- |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
477 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
478 const Bcd opPos() { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
479 return this.dup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
480 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
481 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
482 const Bcd opNeg() { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
483 return negate(this); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
484 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
485 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
486 const Bcd opCom() { |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
487 return twosComp(this); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
488 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
489 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
490 const Bcd opNot() { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
491 return not(this); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
492 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
493 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
494 Bcd opPostInc() { |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
495 this = this + ONE; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
496 return this; |
5
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 Bcd opPostDec() { |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
500 this = this - ONE; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
501 return this; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
502 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
503 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
504 const Bcd opShl(int n) { |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
505 return shift(this, n); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
506 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
507 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
508 Bcd opShlAssign(int n) { |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
509 this = this << n; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
510 return this; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
511 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
512 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
513 const Bcd opShr(int n) { |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
514 return shift(this, -n); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
515 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
516 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
517 Bcd opShrAssign(int n) { |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
518 this = this >> n; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
519 return this; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
520 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
521 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
522 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
523 // binary operators |
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 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
526 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
|
527 return add(this, addend); |
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 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
530 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
|
531 return add(this, Bcd(addend)); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
532 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
533 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
534 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
|
535 this = this + addend; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
536 return this; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
537 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
538 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
539 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
|
540 return subtract(this, addend); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
541 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
542 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
543 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
|
544 return subtract(this, Bcd(addend)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
545 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
546 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
547 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
|
548 this = this - addend; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
549 return this; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
550 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
551 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
552 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
|
553 return multiply(this, b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
554 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
555 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
556 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
|
557 return multiply(this, Bcd(b)); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
558 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
559 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
560 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
|
561 this = this * b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
562 return this; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
563 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
564 |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
565 unittest { |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
566 write("multiply..."); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
567 Bcd a, b, c; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
568 a = 105, b = -22; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
569 c = a * b; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
570 assert(c == -2310); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
571 c = 45; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
572 c *= 1205; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
573 // writeln("c = ", c); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
574 assert(c == 54225); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
575 writeln("passed"); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
576 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
577 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
578 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
|
579 Bcd rem; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
580 return divide(this, b, rem); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
581 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
582 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
583 unittest { |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
584 write("opDiv..."); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
585 Bcd a, b; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
586 a = 7; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
587 b = 3; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
588 assert(a/b == 2); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
589 writeln("passed"); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
590 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
591 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
592 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
|
593 Bcd rem; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
594 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
|
595 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
596 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
597 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
|
598 this = this / b; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
599 return this; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
600 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
601 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
602 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
|
603 Bcd rem; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
604 divide(this, b, rem); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
605 return rem; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
606 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
607 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
608 unittest { |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
609 write("opMod..."); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
610 Bcd a, b; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
611 a = 7; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
612 b = 3; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
613 assert(a % b == 1); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
614 writeln("passed"); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
615 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
616 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
617 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
|
618 Bcd rem; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
619 divide(this, Bcd(b), rem); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
620 return rem; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
621 } |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
622 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
623 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
|
624 this = this % b; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
625 return this; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
626 } |
5
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 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
|
629 return and(this, a); |
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 Bcd opAndAssign(const Bcd a) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
633 this = this & a; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
634 return this; |
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 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
637 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
|
638 return or(this, a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
639 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
640 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
641 Bcd opOrAssign(const Bcd a) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
642 this = this | a; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
643 return this; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
644 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
645 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
646 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
|
647 return xor(this, a); |
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 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
650 Bcd opXorAssign(const Bcd a) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
651 this = this ^ a; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
652 return this; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
653 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
654 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
655 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
656 // comparison, equality operators |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
657 //-------------------------------- |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
658 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
659 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
|
660 return compare(this, that); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
661 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
662 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
663 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
664 write("compare..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
665 Bcd a,b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
666 a = 100; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
667 b = 5; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
668 assert(a > b); |
10 | 669 a = -100; |
670 b = 5; | |
671 assert(a < b); | |
672 a = -100; | |
673 b = -5; | |
674 assert(a < b); | |
675 a = 100; | |
676 b = -5; | |
677 assert(a > b); | |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
678 a = 5; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
679 b = 100; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
680 assert(a < b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
681 assert(b > a); |
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 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
684 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
685 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
|
686 return opCmp(Bcd(that)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
687 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
688 |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
689 const bool opEquals(T:Bcd)(const T that) { |
11 | 690 Bcd a = comparable(this); |
691 Bcd b = comparable(that); | |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
692 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
|
693 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
|
694 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
|
695 if (digit != b.digits[i]) return false; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
696 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
697 return true; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
698 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
699 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
700 const bool opEquals(T)(const T that) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
701 return opEquals(Bcd(that)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
702 } |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
703 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
704 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
705 write("equals..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
706 assert(Bcd(0) == Bcd(0)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
707 assert(Bcd(4) == Bcd(4)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
708 assert(Bcd(40) == Bcd(40)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
709 assert(Bcd(-400) == Bcd(-400)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
710 assert(Bcd(12345678) == Bcd(+12345678)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
711 assert(Bcd(40) != Bcd(53)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
712 assert(Bcd(402) != Bcd(531)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
713 assert(Bcd(402) != Bcd(-402)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
714 assert(Bcd(65432) != Bcd(65431)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
715 assert(Bcd(2) != Bcd(1)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
716 assert(Bcd(1) != Bcd(2)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
717 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
718 } |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
719 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
720 } // end struct Bcd |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
721 |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
722 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
723 public bool isCanonical(const Bcd a) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
724 // no leading zeros |
11 | 725 return !a.hasLeadingZeros; |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
726 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
727 |
11 | 728 // Strips leading zeros |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
729 public Bcd canonical(const Bcd a) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
730 Bcd d = a.dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
731 if (isCanonical(a)) return d; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
732 stripl(d); |
11 | 733 return d; |
734 } | |
735 | |
736 // Strips leading zeros and changes sign if == -0 | |
737 private Bcd comparable(const Bcd a) { | |
738 Bcd d = canonical(a); | |
739 if (d.numDigits == 1 && d.firstDigit == 0) { | |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
740 d.sign = false; |
11 | 741 } |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
742 return d; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
743 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
744 |
10 | 745 |
746 private void clipFirst(ref Bcd a, uint n = 1) { | |
747 if (n == 0) return; | |
748 if (n >= a.digits.length) a = ZERO; | |
749 else a.digits = a.digits[0..$-n]; | |
750 } | |
751 | |
752 private void clipLast(ref Bcd a, uint n = 1) { | |
753 if (n == 0) return; | |
754 if (n >= a.digits.length) a = ZERO; | |
755 else a.digits = a.digits[n..$]; | |
756 } | |
757 | |
758 unittest { | |
759 write("clipping...."); | |
760 Bcd a; | |
761 a = 12345; | |
762 clipFirst(a); | |
763 assert(a == 2345); | |
764 clipLast(a); | |
765 assert(a == 234); | |
766 a = 1234567890; | |
767 clipFirst(a, 3); | |
768 assert(a == 4567890); | |
769 clipLast(a, 5); | |
770 assert(a == 45); | |
771 clipLast(a, 5); | |
772 assert(a == 0); | |
773 writeln("passed"); | |
774 } | |
775 | |
776 | |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
777 /** |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
778 * Strips leading zeros from the specified decint. |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
779 */ |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
780 private void stripl(ref Bcd a) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
781 if (a.hasLeadingZeros) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
782 foreach_reverse(int i, Digit digit; a.digits) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
783 if (i == 0 || digit != 0) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
784 a.digits.length = i+1; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
785 break; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
786 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
787 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
788 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
789 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
790 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
791 unittest { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
792 write("stripl......"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
793 Bcd b = Bcd("+00023"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
794 assert(b.toString == "00023"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
795 stripl(b); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
796 assert(b.toString == "23"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
797 b = "-00050432"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
798 stripl(b); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
799 assert(b.toString == "-50432"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
800 b = "-000"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
801 stripl(b); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
802 assert(b.toString == "-0"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
803 writeln("passed"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
804 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
805 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
806 /** |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
807 * Strips trailing zeros from the specified decint. |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
808 */ |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
809 private void stripr(ref Bcd a) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
810 int len = a.digits.length; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
811 foreach(int i, Digit digit; a.digits) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
812 if (i == len-1 || digit != 0) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
813 a.digits = a.digits[i..$]; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
814 break; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
815 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
816 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
817 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
818 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
819 unittest { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
820 write("stripl......"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
821 Bcd b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
822 b = 23; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
823 stripr(b); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
824 assert(b == 23); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
825 b = "-5000"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
826 stripr(b); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
827 assert(b.toString == "-5"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
828 b = "-000"; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
829 stripr(b); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
830 assert(b.toString == "-0"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
831 writeln("passed"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
832 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
833 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
834 //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
|
835 public string format(const Bcd bcd) { |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
836 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
|
837 if (bcd.isSigned) len++; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
838 // if (bcd.isSigned || bcd.hasLeadingZeros) len++; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
839 char[] str = new char[len]; |
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 int index = 0; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
842 if (bcd.sign) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
843 str[index] = '-'; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
844 index++; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
845 } |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
846 /+ if (bcd.hasLeadingZeros) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
847 str[index] = '+'; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
848 index++; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
849 }+/ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
850 foreach_reverse(Digit digit; bcd.digits) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
851 str[index] = cast(char) digit + '0'; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
852 index++; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
853 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
854 return str.idup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
855 } |
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 public Bcd parse(string str) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
858 Bcd bcd; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
859 // are leading zeros retained? |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
860 bool lz = false; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
861 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
862 // strip whitespace, convert to char array |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
863 char[] str1 = str.strip.dup; |
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 // check for leading '-' |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
866 if (str1[0] == '-') { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
867 bcd.sign = true; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
868 str1 = str1[1..$]; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
869 lz = true; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
870 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
871 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
872 // check for leading '+' |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
873 if (str1[0] == '+') { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
874 bcd.sign = false; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
875 str1 = str1[1..$]; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
876 lz = true; |
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 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
879 // verify chars are digits, strip out underscores |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
880 int index = 0; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
881 char[] str2 = new char[str1.length]; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
882 foreach_reverse(char ch; str1) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
883 if (isDigit(ch)) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
884 str2[index] = ch; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
885 index++; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
886 continue; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
887 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
888 if (ch == '_') { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
889 continue; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
890 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
891 throw (new ConvError("Invalid character: " ~ ch)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
892 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
893 str2.length = index; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
894 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
895 bcd.digits.length = index; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
896 foreach(int i, char ch; str2) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
897 bcd.digits[i] = cast(Digit)(ch - '0'); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
898 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
899 |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
900 if (!lz) return canonical(bcd); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
901 return bcd; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
902 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
903 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
904 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
|
905 return (ch >= '0' && ch <= '9'); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
906 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
907 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
908 public Bcd negate(const Bcd a) { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
909 Bcd b = a.dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
910 b.sign = !a.sign; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
911 return b; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
912 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
913 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
914 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
|
915 return a.isSigned ? -a: +a; |
4
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 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
918 public int compare(const Bcd m, const Bcd n) { |
11 | 919 Bcd a = comparable(m); |
920 Bcd b = comparable(n); | |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
921 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
|
922 if (a.isSigned && !b.isSigned) return -1; |
10 | 923 if (!a.isSigned) { |
924 if (a.digits.length > b.digits.length) return 1; | |
925 if (a.digits.length < b.digits.length) return -1; | |
926 foreach_reverse(int i, Digit digit; a.digits) { | |
927 if (digit > b.digits[i]) return 1; | |
928 if (digit < b.digits[i]) return -1; | |
929 } | |
930 } | |
931 else { | |
932 if (a.digits.length > b.digits.length) return -1; | |
933 if (a.digits.length < b.digits.length) return 11; | |
934 foreach_reverse(int i, Digit digit; a.digits) { | |
935 if (digit > b.digits[i]) return -1; | |
936 if (digit < b.digits[i]) return 1; | |
937 } | |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
938 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
939 return 0; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
940 } |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
941 |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
942 unittest { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
943 write("strip..."); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
944 Bcd bcd; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
945 bcd = "+00123"; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
946 assert(bcd.toString == "00123"); |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
947 bcd = canonical(bcd); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
948 // writeln("bcd = ", bcd); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
949 assert(bcd.toString == "123"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
950 bcd = "+000"; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
951 // writeln("bcd = ", bcd); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
952 assert(bcd.toString == "000"); |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
953 bcd = canonical(bcd); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
954 // writeln("bcd = ", bcd); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
955 writeln("passed"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
956 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
957 |
10 | 958 private bool sameLength(const Bcd a, const Bcd b) { |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
959 return a.numDigits == b.numDigits; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
960 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
961 |
10 | 962 private 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
|
963 if (sameLength(a, b)) return a.numDigits; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
964 uint alen = a.numDigits; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
965 uint blen = b.numDigits; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
966 if (alen > blen) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
967 b.setNumDigits(alen); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
968 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
969 else { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
970 a.setNumDigits(blen); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
971 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
972 return a.numDigits(); |
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 |
10 | 975 // TODO: really? private? |
976 private Bcd twosComp(const Bcd a) { | |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
977 Bcd b = Bcd(0, a.numDigits); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
978 foreach(int i, Digit digit; a.digits) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
979 b.digits[i] = digit == 0 ? 1 : 0 ; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
980 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
981 return b; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
982 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
983 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
984 unittest { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
985 write("com..."); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
986 Bcd bcd; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
987 bcd = 123; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
988 assert(twosComp(bcd).toString == "000"); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
989 bcd = 40509; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
990 assert(twosComp(bcd).toString == "01010"); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
991 // writeln("bcd = ", bcd); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
992 // writeln("~bcd = ", twosComp(bcd)); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
993 writeln("passed"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
994 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
995 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
996 public int sgn(const Bcd bcd) { |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
997 if (bcd.isZero) return 0; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
998 return bcd.isSigned ? -1 : 1; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
999 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1000 |
10 | 1001 public Bcd not(const Bcd a) { |
1002 Bcd b = a.dup; | |
1003 foreach(Digit digit; b.digits) { | |
1004 digit = not(digit); | |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1005 } |
10 | 1006 return b; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1007 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1008 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1009 private Digit not(const Digit a) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1010 return a != 0; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1011 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1012 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1013 public Bcd and(const Bcd x, const Bcd y) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1014 Bcd a = x.dup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1015 Bcd b = y.dup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1016 int len = setSameLength(a,b); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1017 Bcd result = Bcd(0, len); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1018 for (int i = 0; i < len; i++) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1019 result[i] = and(a[i], b[i]); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1020 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1021 return result; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1022 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1023 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1024 private Digit and(const Digit a, const Digit b) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1025 return a != 0 && b != 0; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1026 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1027 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1028 unittest { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1029 write("and..."); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1030 Bcd a; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1031 Bcd b; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1032 a = 123000; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1033 b = 123000; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1034 assert(and(a,b).toString == "111000"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1035 b = 12300; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1036 assert(and(a,b).toString == "011000"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1037 a = 1234567; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1038 b = 7654321; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1039 assert(and(a,b).toString == "1111111"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1040 a = 1010; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1041 b = 101; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1042 // writeln("a = ", a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1043 // writeln("b = ", b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1044 // writeln("and = ", and(a,b)); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1045 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
|
1046 writeln("passed"); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1047 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1048 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1049 // 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
|
1050 public Bcd or(const Bcd x, const Bcd y) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1051 Bcd a = x.dup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1052 Bcd b = y.dup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1053 int len = setSameLength(a,b); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1054 Bcd result = Bcd(0, len); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1055 for (int i = 0; i < len; i++) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1056 result[i] = or(a[i], b[i]); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1057 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1058 return result; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1059 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1060 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1061 private Digit or(const Digit a, const Digit b) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1062 return a != 0 || b != 0; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1063 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1064 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1065 public Bcd xor(const Bcd x, const Bcd y) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1066 Bcd a = x.dup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1067 Bcd b = y.dup; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1068 int len = setSameLength(a,b); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1069 Bcd result = Bcd(0, len); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1070 for (int i = 0; i < len; i++) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1071 result[i] = xor(a[i], b[i]); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1072 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1073 return result; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1074 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1075 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1076 private Digit xor(const Digit a, const Digit b) { |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1077 return (a == 0 && b != 0) || (a != 0 && b == 0); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1078 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1079 |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1080 // 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
|
1081 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
|
1082 Bcd a = x.dup; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1083 Bcd b = y.dup; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1084 Bcd sum; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1085 if (a.sign == b.sign) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1086 sum = simpleAdd(a, b); |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1087 stripl(sum); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1088 sum.sign = a.sign; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1089 return sum; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1090 } |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1091 else { // signs differ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1092 bool sign; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1093 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
|
1094 case 0: |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1095 return ZERO.dup; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1096 case 1: |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1097 sign = a.isSigned; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1098 break; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1099 case -1: |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1100 sign = b.isSigned; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1101 break; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1102 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1103 setSameLength(a, b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1104 sum = simpleAdd(a, tensComp(b)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1105 if (sum.firstDigit == 1) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1106 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
|
1107 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1108 else { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1109 sum = tensComp(sum); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1110 } |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1111 stripl(sum); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1112 sum.sign = sign; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1113 } |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1114 return sum; |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1115 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1116 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1117 unittest { |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1118 write("add...."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1119 Bcd a, b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1120 a = 1234; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1121 b = 4321; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1122 assert(a + b == 5555); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1123 a = 2000; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1124 b = 1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1125 assert(a + b == 2001); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1126 a = -2000; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1127 b = -1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1128 assert(a + b == -2001); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1129 a = 2000; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1130 b = -1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1131 assert(a + b == 1999); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1132 a = -123; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1133 b = 1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1134 assert(a + b == -122); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1135 a = 123; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1136 b = -1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1137 assert(a + b == 122); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1138 a = -1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1139 b = 123; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1140 // writeln(); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1141 // writeln("a + b = ", a + b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1142 assert(a + b == 122); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1143 a = -123; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1144 b = 1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1145 assert(a + b == -122); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1146 a = 1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1147 b = -123; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1148 assert(a + b == -122); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1149 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1150 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1151 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1152 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
|
1153 return add(a, negate(b)); |
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 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1156 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1157 * 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
|
1158 **/ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1159 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
|
1160 Bcd x = a.dup; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1161 Bcd y = b.dup; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1162 int len = setSameLength(x, y); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1163 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
|
1164 Digit carry = 0; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1165 uint i; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1166 // 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
|
1167 // 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
|
1168 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
|
1169 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
|
1170 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1171 if (carry) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1172 sum[i] = 1; |
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 else { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1175 stripl(sum); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1176 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1177 return sum; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1178 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1179 |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1180 // TODO: need a clipFirst, clipLast method |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1181 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1182 /** |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1183 * 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
|
1184 * 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
|
1185 **/ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1186 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
|
1187 Digit sum = a + b + carry; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1188 carry = sum > 9; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1189 if (carry) sum -= 10; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1190 return sum; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1191 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1192 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1193 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1194 write("simpleAdd..."); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1195 Bcd a, b; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1196 a = 123; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1197 b = 222; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1198 Bcd sum = simpleAdd(a,b); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1199 assert(sum == 345); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1200 a = 2; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1201 b = 102000; |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1202 sum = simpleAdd(a,b); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1203 assert(sum == 102002); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1204 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1205 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1206 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1207 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
|
1208 /+ foreach(Digit digit; a.digits) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1209 digit = 9 - digit; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1210 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1211 a = simpleAdd(a, ONE); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1212 return a;+/ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1213 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1214 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
|
1215 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
|
1216 x[i] = 9 - digit; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1217 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1218 return simpleAdd(x, ONE); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1219 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1220 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1221 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1222 write("tensComp..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1223 Bcd a, b, c; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1224 a = 123456; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1225 b = tensComp(a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1226 // writeln("a = ", a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1227 // writeln("b = ", b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1228 c = 876544; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1229 assert(b == c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1230 a = Bcd(0, 4); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1231 b = tensComp(a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1232 c = 10000; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1233 assert(b == c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1234 a = Bcd(1, 4); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1235 b = tensComp(a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1236 c = 9999; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1237 assert(b == c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1238 a = 9999; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1239 b = tensComp(a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1240 c = 1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1241 assert(b == c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1242 a = 10000; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1243 b = tensComp(a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1244 c = 90000; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1245 assert(b == c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1246 a = 1; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1247 b = tensComp(a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1248 c = 9; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1249 // writeln("a = ", a); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1250 // writeln("b = ", b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1251 assert(b == c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1252 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1253 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1254 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1255 // 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
|
1256 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
|
1257 Bcd n, m; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1258 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
|
1259 if (a.numDigits > b.numDigits) { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1260 n = canonical(a); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1261 m = canonical(b); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1262 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1263 else { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1264 n = canonical(b); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1265 m = canonical(a); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1266 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1267 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
|
1268 if (d == 0) continue; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1269 // 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
|
1270 // 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
|
1271 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
|
1272 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1273 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
|
1274 return p; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1275 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1276 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1277 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1278 write("mul..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1279 Bcd a, b, p; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1280 a = 2105, b = 301; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1281 p = multiply(a, b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1282 assert(p == Bcd(633605)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1283 b = 800, a = 23958233; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1284 p = multiply(a, b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1285 assert(p == Bcd("19166586400")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1286 b = 5830, a = 23958233; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1287 p = multiply(a, b); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1288 assert(p == Bcd("139676498390")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1289 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1290 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1291 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1292 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
|
1293 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
|
1294 Digit c = 0; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1295 int i = 0; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1296 foreach (Digit d; a.digits) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1297 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
|
1298 i++; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1299 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1300 if (c) p[i] = c; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1301 return p; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1302 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1303 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1304 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1305 write("mul..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1306 Bcd a, p; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1307 Digit m; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1308 a = 2105, m = 3; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1309 p = mul(a, m); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1310 assert(p == Bcd(6315)); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1311 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1312 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1313 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1314 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
|
1315 Digit p = a * b + c; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1316 c = p / 10; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1317 return p % 10; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1318 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1319 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1320 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1321 write("mul..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1322 Digit a, b, c, d; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1323 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
|
1324 d = mul(a,b,c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1325 assert(d == 6); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1326 assert(c == 0); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1327 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
|
1328 d = mul(a,b,c); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1329 assert(d == 2); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1330 assert(c == 7); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1331 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1332 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1333 |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1334 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
|
1335 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1336 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
|
1337 if (a == ZERO) return ZERO.dup; |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1338 |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1339 Bcd d = canonical(abs(b)); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1340 r = canonical(abs(a)); |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1341 Bcd q; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1342 Bcd p = d; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1343 Bcd m; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1344 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1345 // shift the divisor |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1346 do { |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1347 m = p; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1348 p <<= 1; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1349 } while (p <= r && p > m); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1350 |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1351 while (m >= d) { |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1352 q <<= 1; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1353 // subtract repeatedly |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1354 while (m <= r) { |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1355 r -= m; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1356 q++; |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1357 } |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1358 m >>= 1; |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1359 } |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1360 |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1361 // adjust the sign |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1362 r.sign = a.sign; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1363 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
|
1364 return q; |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1365 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1366 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1367 unittest { |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1368 write("divide..."); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1369 Bcd a, b, q, r; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1370 a = 144; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1371 b = 12; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1372 q = divide(a, b, r); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1373 assert(q == 12); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1374 assert(r == 0); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1375 a = 144; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1376 b = 14; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1377 q = divide(a, b, r); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1378 assert(q == 10); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1379 assert(r == 4); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1380 a = 142351; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1381 b = 12; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1382 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
|
1383 // writeln("q = ", q); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1384 // writeln("r = ", r); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1385 assert(q == 11862); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1386 assert(r == 7); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1387 a = 2; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1388 b = 12; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1389 q = divide(a, b, r); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1390 assert(q == 0); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1391 assert(r == 2); |
7
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1392 a = 7; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1393 b = 3; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1394 q = divide(a, b, r); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1395 assert(q == 2); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1396 assert(r == 1); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1397 a = 7; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1398 b = -3; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1399 q = divide(a, b, r); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1400 assert(q == -2); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1401 assert(r == 1); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1402 a = -7; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1403 b = 3; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1404 q = divide(a, b, r); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1405 assert(q == -2); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1406 assert(r == -1); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1407 a = -7; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1408 b = -3; |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1409 q = divide(a, b, r); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1410 assert(q == 2); |
b304232c476c
Coded and tested div and mod. Changed setNumDigits.
Paul (paul.d.anderson@comcast.net)
parents:
6
diff
changeset
|
1411 assert(r == -1); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1412 writeln("passed"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1413 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1414 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1415 |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1416 public Bcd truncate(const Bcd a, int n) { |
10 | 1417 return truncate(a, n, Bcd()); |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1418 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1419 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1420 unittest { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1421 write("truncate...."); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1422 Bcd a; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1423 a = 1234567; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1424 assert(truncate(a,3) == 123); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1425 a = 10256; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1426 assert(truncate(a,5) == 10256); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1427 a = 8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1428 assert(truncate(a,1) == 8); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1429 a = -8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1430 assert(truncate(a,1) == -8); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1431 a = -8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1432 assert(truncate(a,-1) == 0); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1433 writeln("passed"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1434 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1435 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1436 /** |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1437 * Truncate to the specified number of digits |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1438 */ |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1439 public Bcd truncate(const Bcd a, int n, out Bcd r) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1440 if (n <= 0) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1441 r = canonical(a); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1442 return ZERO.dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1443 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1444 if (n >= a.numDigits) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1445 r = ZERO.dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1446 return canonical(a); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1447 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1448 Bcd b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1449 b.digits = a.digits[$-n..$].dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1450 b.sign = a.sign; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1451 r.digits = a.digits[0..$-n].dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1452 return b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1453 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1454 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1455 unittest { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1456 write("truncRem..."); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1457 Bcd a; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1458 Bcd r; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1459 a = 1234567; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1460 assert(truncate(a, 3, r) == 123); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1461 assert(r == 4567); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1462 a = 10256; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1463 assert(truncate(a,5,r) == 10256); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1464 assert(r == 0); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1465 a = 8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1466 assert(truncate(a,1,r) == 8); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1467 assert(r == 500); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1468 a = -8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1469 assert(truncate(a,1,r) == -8); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1470 assert(r == 500); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1471 a = -8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1472 // writeln("truncate(a, 6, r) = ", truncate(a, 6, r)); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1473 // writeln("r = ", r); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1474 assert(truncate(a,6,r) == -8500); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1475 assert(r == 0); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1476 writeln("passed"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1477 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1478 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1479 public Bcd pad(Bcd a, int n) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1480 if (n <= a.numDigits) return a.dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1481 return Bcd(a,n); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1482 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1483 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1484 unittest { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1485 write("pad........."); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1486 Bcd a; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1487 a = 1234567; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1488 assert(pad(a,3) == 1234567); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1489 a = 10256; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1490 // writeln("pad(a,9) = ", pad(a,9)); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1491 assert(pad(a,9).toString == "000010256"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1492 a = 8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1493 assert(pad(a,6).toString == "008500"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1494 a = -8500; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1495 assert(pad(a,6).toString == "-008500"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1496 writeln("passed"); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1497 } |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1498 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1499 public Bcd pow10(const Bcd a, uint n) { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1500 Bcd b = canonical(a); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1501 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
|
1502 Digit[] zeros = new Digit[n]; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1503 b.digits = zeros ~ b.digits; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1504 return b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1505 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1506 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1507 public Bcd pow10(uint n) { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1508 return pow10(ONE, n); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1509 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1510 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1511 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1512 write("pow10..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1513 Bcd a, b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1514 a = 21345; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1515 a = pow10(a, 12); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1516 assert(a == Bcd("21345000000000000")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1517 a = pow10(a, 0); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1518 assert(a == "21345000000000000"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1519 a = pow10(9); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1520 assert(a == Bcd("1000000000")); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1521 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1522 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1523 |
10 | 1524 public Bcd power(const Bcd a, uint n) { |
1525 Bcd b = canonical(a); | |
1526 Bcd p = ONE.dup; | |
1527 while (n > 0) { | |
1528 if (n & 1) { | |
1529 p *= b; | |
1530 if (n == 1) return p; | |
1531 } | |
1532 b *= b; | |
1533 n /= 2; | |
1534 } | |
1535 return p; | |
1536 } | |
1537 | |
1538 unittest { | |
1539 write("power......."); | |
1540 Bcd a; | |
1541 a = 2; | |
1542 int n = 3; | |
1543 assert(power(a,n) == 8); | |
1544 writeln("passed"); | |
1545 } | |
1546 | |
1547 | |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1548 public Bcd shift(const Bcd a, int n) { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1549 if (n == 0) return canonical(a); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1550 if (n > 0) { // shift left |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1551 return pow10(a, n); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1552 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1553 else { // shift right |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1554 if (-n >= a.numDigits) return ZERO.dup; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1555 return truncate(a, a.numDigits + n); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1556 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1557 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1558 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1559 unittest { |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1560 write("shift......."); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1561 Bcd a, b; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1562 a = 1234; |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1563 b = shift(a, 1); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1564 assert(b.toString == "12340"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1565 b = shift(a, 2); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1566 assert(b.toString == "123400"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1567 b = shift(a, -1); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1568 // writeln("b = ", b); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1569 assert(b.toString == "123"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1570 b = shift(a, 0); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1571 assert(b.toString == "1234"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1572 b = shift(a, -12); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1573 /+ writeln("b = ", b);+/ |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1574 assert(b.toString == "0"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1575 b = shift(a, -4); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1576 assert(b.toString == "0"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1577 writeln("passed"); |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1578 } |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1579 |
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1580 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
|
1581 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
|
1582 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
|
1583 if (std.math.abs(n) >= a.numDigits) return b; |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1584 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
|
1585 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
|
1586 } |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1587 else { // shiftFixed right |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1588 n = -n; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1589 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
|
1590 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1591 return b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1592 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1593 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1594 unittest { |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1595 write("shiftFixed..."); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1596 Bcd a, b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1597 a = 1234; |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1598 b = shiftFixed(a, 1); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1599 assert(b.toString == "0123"); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1600 b = shiftFixed(a, 2); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1601 assert(b.toString == "0012"); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1602 b = shiftFixed(a, -1); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1603 assert(b.toString == "2340"); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1604 b = shiftFixed(a, 0); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1605 assert(b.toString == "1234"); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1606 b = shiftFixed(a, 12); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1607 assert(b.toString == "0000"); |
6
14badf9104b9
shifts, rotations, division
Paul (paul.d.anderson@comcast.net)
parents:
5
diff
changeset
|
1608 b = shiftFixed(a, -4); |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1609 assert(b.toString == "0000"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1610 writeln("passed"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1611 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1612 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1613 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
|
1614 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
|
1615 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
|
1616 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
|
1617 bool sign = n < 0; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1618 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
|
1619 if (sign) n = -n; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1620 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1621 if (n > 0) { // rotate left |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1622 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
|
1623 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
|
1624 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1625 else { // rotate right |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1626 n = -n; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1627 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
|
1628 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
|
1629 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1630 return b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1631 } |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1632 |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1633 unittest { |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1634 write("rotate..."); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1635 Bcd a, b; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1636 a = 1234; |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1637 b = rotate(a, 1); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1638 assert(b.toString == "4123"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1639 b = rotate(a, 2); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1640 assert(b.toString == "3412"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1641 b = rotate(a, -1); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1642 assert(b.toString == "2341"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1643 b = rotate(a, 0); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1644 assert(b.toString == "1234"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1645 b = rotate(a, 12); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1646 assert(b.toString == "1234"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1647 b = rotate(a, -4); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1648 assert(b.toString == "1234"); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1649 writeln("passed"); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1650 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1651 |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1652 /** |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1653 * Rounds the specified integer to the specified number of digits. |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1654 * The rounding mode governs the method of rounding. |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1655 */ |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1656 public Bcd round(const Bcd a, int n, Rounding mode) { |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1657 Bcd r; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1658 Bcd b = truncate(a, n, r); |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1659 return round(b, r, mode); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1660 // return b; |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1661 } |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1662 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1663 /** |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1664 * Rounds the (truncated) integer based on its remainder and the rounding |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1665 * mode. |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1666 */ |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1667 public Bcd round(const Bcd a, const Bcd r, const Rounding mode) { |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1668 Bcd b = canonical(a); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1669 ubyte rfd = r.firstDigit; |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1670 switch (mode) { |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1671 case Rounding.DOWN: |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1672 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1673 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1674 case Rounding.HALF_UP: |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1675 if (rfd >= 5) b++; |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1676 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1677 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1678 case Rounding.HALF_EVEN: |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1679 if (rfd < 5) return b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1680 if (rfd > 5) return ++b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1681 Bcd p = r.dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1682 stripr(p); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1683 if (p.numDigits > 1) return ++b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1684 // to reach this point the remainder must be exactly 5 |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1685 // if odd, increment |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1686 if (b.lastDigit & 1 != 0) b++; |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1687 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1688 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1689 case Rounding.CEILING: |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1690 if (!b.isSigned && !r.isZero) b++; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1691 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1692 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1693 case Rounding.FLOOR: |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1694 if (b.isSigned && !r.isZero) b++; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1695 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1696 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1697 case Rounding.HALF_DOWN: |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1698 if (rfd < 5) return b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1699 if (rfd > 5) return ++b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1700 Bcd p = r.dup; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1701 stripr(p); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1702 if (p.numDigits > 1) return ++b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1703 // to reach this point the remainder must be exactly 5 |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1704 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1705 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1706 case Rounding.UP: |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1707 if (!r.isZero) b++; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1708 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1709 } |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1710 return b; |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1711 } |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1712 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1713 unittest { |
9
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1714 write("rounding...."); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1715 Bcd a, b; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1716 a = 12345; |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1717 b = round(a, 4, Rounding.DOWN); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1718 assert(b == 1234); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1719 b = round(a, 4, Rounding.HALF_UP); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1720 assert(b == 1235); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1721 b = round(a, 4, Rounding.CEILING); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1722 assert(b == 1235); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1723 b = round(a, 4, Rounding.HALF_EVEN); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1724 assert(b == 1234); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1725 b = round(a, 4, Rounding.FLOOR); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1726 assert(b == 1234); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1727 b = round(a, 4, Rounding.HALF_DOWN); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1728 assert(b == 1234); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1729 b = round(a, 4, Rounding.UP); |
48d564218e05
Added rounding, additional unit tests, strip leading and trailing zeros.
Paul (paul.d.anderson@comcast.net)
parents:
8
diff
changeset
|
1730 assert(b == 1235); |
8
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1731 writeln("passed"); |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1732 } |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1733 |
c991b9fde45c
added rounding (partially implemented)
Paul (paul.d.anderson@comcast.net)
parents:
7
diff
changeset
|
1734 |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1735 //========================================== |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1736 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1737 public void main() { |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1738 writeln(); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1739 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
|
1740 writeln(); |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1741 writeln(" Congratulations!"); |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1742 Bcd bcd; |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1743 writeln("bcd = ", format(bcd)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1744 bcd = parse("12"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1745 writeln("bcd = ", format(bcd)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1746 bcd = parse("012"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1747 writeln("bcd = ", format(bcd)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1748 bcd = parse("-012"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1749 writeln("bcd = ", format(bcd)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1750 bcd = parse("+012"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1751 writeln("bcd = ", format(bcd)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1752 bcd = parse("-0"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1753 writeln("bcd = ", format(bcd)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1754 bcd = parse("012_345_678"); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1755 writeln("bcd = ", format(bcd)); |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1756 |
5
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1757 /+ |
c021ed211f89
bcd add, sub, multiply, logic and shift operations
Paul (paul.d.anderson@comcast.net)
parents:
4
diff
changeset
|
1758 writeln("OK!");+/ |
4
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1759 } |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1760 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1761 |
b37c218c1442
Initial development of BCD integers
Paul (paul.d.anderson@comcast.net)
parents:
diff
changeset
|
1762 |