Add BigInteger type

std/haxe/math/bigint/BigIntError.hx

@@ -0,0 +1,33 @@
+/*
+ * Copyright (C)2005-2023 Haxe Foundation
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+ * DEALINGS IN THE SOFTWARE.
+ */
+
+package haxe.math.bigint;
+
+enum abstract BigIntError(String) to String {
+	var INVALID_ARGUMENT = "Invalid argument";
+	var BUFFER_TOO_SMALL = "Buffer too small";
+	var DIVISION_BY_ZERO = "Division by zero";
+	var NEGATIVE_EXPONENT = "Negative exponent";
+	var INVALID_OPERATION = "Invalid operation";
+	var NEGATIVE_MODULUS = "Modulus should be positive";
+	var EVEN_VALUES = "Both values are even";
+}

std/haxe/math/bigint/BigIntException.hx

@@ -0,0 +1,29 @@
+/*
+ * Copyright (C)2005-2023 Haxe Foundation
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+ * DEALINGS IN THE SOFTWARE.
+ */
+
+package haxe.math.bigint;
+
+class BigIntException extends haxe.Exception {
+	public function new(error:BigIntError, ?previous:haxe.Exception) {
+		super( error, previous);
+	}
+}

std/haxe/math/bigint/BigIntHelper.hx

@@ -0,0 +1,156 @@
+/*
+ * Copyright (C)2005-2022 Haxe Foundation
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+ * DEALINGS IN THE SOFTWARE.
+ */
+
+package haxe.math.bigint;
+
+/* Original code courtesy Chuck Batson (github.com/cbatson) */
+class BigIntHelper {
+	/**
+		"Numbler of leading zeros" - return the number of leading
+		0-value bits in the binary representation of `x`.
+	**/
+	public static function nlz(x:Int):Int {
+		// From "Hacker's Delight", Second Edition; Henry S. Warren, Jr.; 2013. Figure 5-15, p. 102.
+		var y:Int, m:Int, n:Int;
+
+		y = -(x >>> 16);
+		m = (y >> 16) & 16;
+		n = 16 - m;
+		x = x >>> m;
+
+		y = x - 0x100;
+		m = (y >> 16) & 8;
+		n = n + m;
+		x = x << m;
+
+		y = x - 0x1000;
+		m = (y >> 16) & 4;
+		n = n + m;
+		x = x << m;
+
+		y = x - 0x4000;
+		m = (y >> 16) & 2;
+		n = n + m;
+		x = x << m;
+
+		y = x >> 14;
+		m = y & (~y >> 1);
+		return n + 2 - m;
+	}
+
+	public static function bitLen(x:Int):Int {
+		var sign:Int = (x<0)?-1:0;
+		return ( 32 - BigIntHelper.nlz(x^sign) );
+	}
+
+	/**
+		"Ceiling power of two" -- round up to the least power of two
+		greater than or equal to input `x`, which is interpreted as
+		unsigned.
+	**/
+	public static function clp2(x:Int32):Int {
+		// From "Hacker's Delight", Second Edition; Henry S. Warren, Jr.; 2013. Figure 3-3, p. 62.
+		x = x - 1;
+		x = x | (x >> 1);
+		x = x | (x >> 2);
+		x = x | (x >> 4);
+		x = x | (x >> 8);
+		x = x | (x >> 16);
+		return x + 1;
+	}
+
+	/**
+		Unsigned greater than comparison.
+
+		Returns `true` if `a > b` when both `a` and `b` are
+		interpreted as unsigned integers; `false` otherwise.
+	**/
+	public static inline function u32gtu32(a:Int, b:Int):Bool {
+		return (a ^ -2147483648) > (b ^ -2147483648); // unsigned comparison, see "Hacker's Delight" p. 25.
+	}
+
+	/**
+		Integer division of unsigned 32-bit integer by unsigned 16-bit integer.
+
+		Result is undefined when `divisor` <= 0 or `divisor` >= 2^16.
+	**/
+	public static function u32divu16(dividend:Int32, divisor:Int32):Int {
+		/*
+			Complicated because Haxe's division is always performed as
+			floating-point.  Here we rely on the ability to exactly represent
+			a 31-bit integer as a Float.  In other words, 64-bit floating
+			point is required.
+
+			TODO: Implement a method without this restriction.
+			TODO: Consider C++-specific optimization here.
+		 */
+
+		// From "Hacker's Delight", Second Edition; Henry S. Warren, Jr.; 2013. Section 9-3, p. 192.
+		var t:Int = divisor >> 31;
+		var nprime:Int = dividend & ~t;
+		var q:Int32 = Std.int((nprime >>> 1) / divisor) << 1;
+		var r:Int = dividend - q * divisor;
+		var c:Int = u32geu32(r, divisor) ? 1 : 0;
+		return q + c;
+	}
+
+	/**
+		Unsigned greater than or equal comparison.
+		Returns `true` if `a >= b` when both `a` and `b` are
+		interpreted as unsigned integers; `false` otherwise.
+	**/
+	public static inline function u32geu32(a:Int, b:Int):Bool {
+		return (a ^ -2147483648) >= (b ^ -2147483648); // unsigned comparison, see "Hacker's Delight" p. 25.
+	}
+
+	/**
+		Number of trailing zeros - return the number of trailing
+		0-value bits 
+	**/
+	public static function ntz(x:Int32):Int {
+		if (x == 0)
+			return 32;
+		var y:Int;
+		var n:Int = 31;
+		y = x << 16;
+		if (y != 0) {
+			n -= 16;
+			x = y;
+		}
+		y = x << 8;
+		if (y != 0) {
+			n -= 8;
+			x = y;
+		}
+		y = x << 4;
+		if (y != 0) {
+			n -= 4;
+			x = y;
+		}
+		y = x << 2;
+		if (y != 0) {
+			n -= 2;
+			x = y;
+		}
+		return (n - ((x << 1) >>> 31));
+	}
+}

std/haxe/math/bigint/BigIntTools.hx

@@ -0,0 +1,69 @@
+/*
+ * Copyright (C)2005-2023 Haxe Foundation
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+ * DEALINGS IN THE SOFTWARE.
+ */
+
+package haxe.math.bigint;
+
+import haxe.math.bigint.BigIntException;
+import haxe.math.bigint.BigIntError;
+
+/* Original code courtesy Chuck Batson (github.com/cbatson) */
+class BigIntTools {
+	public static inline function isNull(value:BigInt):Bool {
+		var a:BigInt_ = value;
+		return a == null;
+	}
+
+	public static inline function isBigInt(value:Dynamic):Bool {
+		return Std.isOfType(value, BigInt_);
+	}
+
+	public static inline function castFrom(value:Dynamic):BigInt {
+		return new BigInt(Std.downcast(value, BigInt_));
+	}
+
+	public static function parseValueUnsigned(value:Dynamic):BigInt {
+		var bi:BigInt;
+		if (Std.isOfType(value, String)) {
+			bi = parseStringUnsigned(cast(value, String));
+		} else if (isBigInt(value)) {
+			var t = new MutableBigInt_();
+			t.copyFrom(castFrom(value));
+			return new BigInt(t);
+		} else if (Std.isOfType(value, Int)) {
+			bi = BigInt.fromInt(cast(value, Int));
+		} else {
+			throw new BigIntException(BigIntError.INVALID_ARGUMENT);
+		}
+		return bi;
+	}
+
+	private static function parseStringUnsigned(value:String):BigInt {
+		var result = new MutableBigInt_();
+		if (StringTools.startsWith(value, "0x")) {
+			result.setFromHexUnsigned(value.substr(2));
+		} else {
+			result.setFromString(value);
+		}
+		var result2:MutableBigInt = result;
+		return result2;
+	}
+}