projects/decimal: src/decimal/bcd.d annotate

annotate src/decimal/bcd.d @ 9:48d564218e05

Added rounding, additional unit tests, strip leading and trailing zeros.

author	Paul (paul.d.anderson@comcast.net)
date	Wed, 24 Mar 2010 22:21:38 -0700
parents	c991b9fde45c
children	f925fe996255

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

Mercurial > projects > decimal

annotate src/decimal/bcd.d @ 9:48d564218e05