aboutsummaryrefslogtreecommitdiff
path: root/core/math/big
diff options
context:
space:
mode:
authorJeroen van Rijn <Kelimion@users.noreply.github.com>2021-07-20 19:51:39 +0200
committerJeroen van Rijn <Kelimion@users.noreply.github.com>2021-08-11 20:59:50 +0200
commit4eadd0867d6529340ad1da3ab2af56d62b40f3d2 (patch)
treea63a0ec6d425e5fb8f5857b019afdf2ff750bfb4 /core/math/big
parent9dba17cf87a208f41f2d53bc9a2be14ed70a396a (diff)
big: Continuing to refactor.
Diffstat (limited to 'core/math/big')
-rw-r--r--core/math/big/api.odin155
-rw-r--r--core/math/big/basic.odin26
-rw-r--r--core/math/big/bigint.odin2
-rw-r--r--core/math/big/compare.odin75
-rw-r--r--core/math/big/helpers.odin64
-rw-r--r--core/math/big/log.odin2
-rw-r--r--core/math/big/logical.odin2
-rw-r--r--core/math/big/radix.odin2
8 files changed, 234 insertions, 94 deletions
diff --git a/core/math/big/api.odin b/core/math/big/api.odin
new file mode 100644
index 000000000..d2cbd8565
--- /dev/null
+++ b/core/math/big/api.odin
@@ -0,0 +1,155 @@
+package big
+
+/*
+ Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
+ Made available under Odin's BSD-2 license.
+
+ An arbitrary precision mathematics implementation in Odin.
+ For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
+ The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.
+
+ This file collects public procs into proc maps, as well as any of their aliases.
+/*
+
+*/
+ === === === === === === === === === === === === === === === === === === === === === === === ===
+ Basic arithmetic.
+ See `basic.odin`.
+ === === === === === === === === === === === === === === === === === === === === === === === ===
+*/
+
+/*
+ err = add(dest, a, b);
+*/
+add :: proc {
+ /*
+ int_add :: proc(dest, a, b: ^Int) -> (err: Error)
+ */
+ int_add,
+ /*
+ int_add_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error)
+ */
+ int_add_digit,
+};
+
+/*
+ err = sub(dest, a, b);
+*/
+sub :: proc {
+ /*
+ int_sub :: proc(dest, a, b: ^Int) -> (err: Error)
+ */
+ int_sub,
+ /*
+ int_sub_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error)
+ */
+ int_sub_digit,
+};
+
+/*
+ === === === === === === === === === === === === === === === === === === === === === === === ===
+ Comparisons.
+ See `compare.odin`.
+ === === === === === === === === === === === === === === === === === === === === === === === ===
+*/
+
+is_initialized :: proc {
+ /*
+ int_is_initialized :: proc(a: ^Int) -> bool
+ */
+ int_is_initialized,
+};
+
+is_zero :: proc {
+ /*
+ int_is_zero :: proc(a: ^Int) -> bool
+ */
+ int_is_zero,
+};
+
+is_positive :: proc {
+ /*
+ int_is_positive :: proc(a: ^Int) -> bool
+ */
+ int_is_positive,
+};
+is_pos :: is_positive;
+
+is_negative :: proc {
+ /*
+ int_is_negative :: proc(a: ^Int) -> bool
+ */
+ int_is_negative,
+};
+is_neg :: is_negative;
+
+is_even :: proc {
+ /*
+ int_is_even :: proc(a: ^Int) -> bool
+ */
+ int_is_even,
+};
+
+is_odd :: proc {
+ /*
+ int_is_odd :: proc(a: ^Int) -> bool
+ */
+ int_is_odd,
+};
+
+is_power_of_two :: proc {
+ /*
+ platform_int_is_power_of_two :: proc(a: int) -> bool
+ */
+ platform_int_is_power_of_two,
+ /*
+ int_is_power_of_two :: proc(a: ^Int) -> (res: bool)
+ */
+ int_is_power_of_two,
+};
+
+compare :: proc {
+ /*
+ Compare two `Int`s, signed.
+
+ int_compare :: proc(a, b: ^Int) -> Comparison_Flag
+ */
+ int_compare,
+ /*
+ Compare an `Int` to an unsigned number upto the size of the backing type.
+
+ int_compare_digit :: proc(a: ^Int, u: DIGIT) -> Comparison_Flag
+ */
+ int_compare_digit,
+};
+cmp :: compare;
+
+compare_magnitude :: proc {
+ /*
+ Compare the magnitude of two `Int`s, unsigned.
+ */
+ int_compare_magnitude,
+};
+cmp_mag :: compare_magnitude;
+
+/*
+ === === === === === === === === === === === === === === === === === === === === === === === ===
+ Initialization and other helpers.
+ See `helpers.odin`.
+ === === === === === === === === === === === === === === === === === === === === === === === ===
+*/
+
+destroy :: proc {
+ /*
+ Deallocates the backing memory of an `Int` and resets it.
+
+ int_destroy :: proc(a: ^Int, allocator_zeroes := false, allocator := context.allocator)
+ */
+ int_destroy,
+};
+
+init :: proc{
+ int_init_sized,
+ int_init_from_integer,
+ int_init_copy,
+};
diff --git a/core/math/big/basic.odin b/core/math/big/basic.odin
index 98f14687c..43f5a3715 100644
--- a/core/math/big/basic.odin
+++ b/core/math/big/basic.odin
@@ -4,7 +4,7 @@ package big
Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
Made available under Odin's BSD-2 license.
- A BigInt implementation in Odin.
+ An arbitrary precision mathematics implementation in Odin.
For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.
@@ -23,7 +23,7 @@ import "core:intrinsics"
/*
High-level addition. Handles sign.
*/
-add_two_ints :: proc(dest, a, b: ^Int) -> (err: Error) {
+int_add :: proc(dest, a, b: ^Int) -> (err: Error) {
dest := dest; x := a; y := b;
assert_initialized(dest); assert_initialized(a); assert_initialized(b);
@@ -32,7 +32,7 @@ add_two_ints :: proc(dest, a, b: ^Int) -> (err: Error) {
*/
if x.sign == y.sign {
dest.sign = x.sign;
- return _add(dest, x, y);
+ return _int_add(dest, x, y);
}
/*
@@ -45,7 +45,7 @@ add_two_ints :: proc(dest, a, b: ^Int) -> (err: Error) {
}
dest.sign = x.sign;
- return _sub(dest, x, y);
+ return _int_sub(dest, x, y);
}
/*
@@ -54,7 +54,7 @@ add_two_ints :: proc(dest, a, b: ^Int) -> (err: Error) {
dest = a + digit;
*/
-add_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error) {
+int_add_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error) {
dest := dest; digit := digit;
assert_initialized(dest); assert_initialized(a);
@@ -165,12 +165,10 @@ add_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error) {
return .OK;
}
-add :: proc{add_two_ints, add_digit};
-
/*
High-level subtraction, dest = number - decrease. Handles signs.
*/
-sub_two_ints :: proc(dest, number, decrease: ^Int) -> (err: Error) {
+int_sub :: proc(dest, number, decrease: ^Int) -> (err: Error) {
dest := dest; x := number; y := decrease;
assert_initialized(number); assert_initialized(decrease); assert_initialized(dest);
@@ -180,7 +178,7 @@ sub_two_ints :: proc(dest, number, decrease: ^Int) -> (err: Error) {
In either case, ADD their magnitudes and use the sign of the first number.
*/
dest.sign = x.sign;
- return _add(dest, x, y);
+ return _int_add(dest, x, y);
}
/*
@@ -201,7 +199,7 @@ sub_two_ints :: proc(dest, number, decrease: ^Int) -> (err: Error) {
*/
dest.sign = x.sign;
}
- return _sub(dest, x, y);
+ return _int_sub(dest, x, y);
}
/*
@@ -210,7 +208,7 @@ sub_two_ints :: proc(dest, number, decrease: ^Int) -> (err: Error) {
dest = a - digit;
*/
-sub_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error) {
+int_sub_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error) {
dest := dest; digit := digit;
assert_initialized(dest); assert_initialized(a);
@@ -296,8 +294,6 @@ sub_digit :: proc(dest, a: ^Int, digit: DIGIT) -> (err: Error) {
return .OK;
}
-sub :: proc{sub_two_ints, sub_digit};
-
/*
==========================
Low-level routines
@@ -308,7 +304,7 @@ sub :: proc{sub_two_ints, sub_digit};
Low-level addition, unsigned.
Handbook of Applied Cryptography, algorithm 14.7.
*/
-_add :: proc(dest, a, b: ^Int) -> (err: Error) {
+_int_add :: proc(dest, a, b: ^Int) -> (err: Error) {
dest := dest; x := a; y := b;
assert_initialized(a); assert_initialized(b); assert_initialized(dest);
@@ -389,7 +385,7 @@ _add :: proc(dest, a, b: ^Int) -> (err: Error) {
Low-level subtraction, dest = number - decrease. Assumes |number| > |decrease|.
Handbook of Applied Cryptography, algorithm 14.9.
*/
-_sub :: proc(dest, number, decrease: ^Int) -> (err: Error) {
+_int_sub :: proc(dest, number, decrease: ^Int) -> (err: Error) {
dest := dest; x := number; y := decrease;
assert_initialized(number); assert_initialized(decrease); assert_initialized(dest);
diff --git a/core/math/big/bigint.odin b/core/math/big/bigint.odin
index 82ddaeb79..0e5b11381 100644
--- a/core/math/big/bigint.odin
+++ b/core/math/big/bigint.odin
@@ -4,7 +4,7 @@ package big
Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
Made available under Odin's BSD-2 license.
- A BigInt implementation in Odin.
+ An arbitrary precision mathematics implementation in Odin.
For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.
*/
diff --git a/core/math/big/compare.odin b/core/math/big/compare.odin
index facfca2e8..b9c6dace2 100644
--- a/core/math/big/compare.odin
+++ b/core/math/big/compare.odin
@@ -4,32 +4,32 @@ package big
Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
Made available under Odin's BSD-2 license.
- A BigInt implementation in Odin.
+ An arbitrary precision mathematics implementation in Odin.
For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.
+
+ This file contains various comparison routines.
*/
import "core:intrinsics"
-is_initialized :: proc(a: ^Int) -> bool {
+int_is_initialized :: proc(a: ^Int) -> bool {
return a != rawptr(uintptr(0));
}
-is_zero :: proc(a: ^Int) -> bool {
+int_is_zero :: proc(a: ^Int) -> bool {
return is_initialized(a) && a.used == 0;
}
-is_positive :: proc(a: ^Int) -> bool {
+int_is_positive :: proc(a: ^Int) -> bool {
return is_initialized(a) && a.sign == .Zero_or_Positive;
}
-is_pos :: is_positive;;
-is_negative :: proc(a: ^Int) -> bool {
+int_is_negative :: proc(a: ^Int) -> bool {
return is_initialized(a) && a.sign == .Negative;
}
-is_neg :: is_negative;
-is_even :: proc(a: ^Int) -> bool {
+int_is_even :: proc(a: ^Int) -> bool {
if is_initialized(a) {
if is_zero(a) {
return true;
@@ -41,18 +41,18 @@ is_even :: proc(a: ^Int) -> bool {
return false;
}
-is_odd :: proc(a: ^Int) -> bool {
+int_is_odd :: proc(a: ^Int) -> bool {
if is_initialized(a) {
return !is_even(a);
}
return false;
}
-is_power_of_two_small :: proc(a: int) -> bool {
+platform_int_is_power_of_two :: proc(a: int) -> bool {
return ((a) != 0) && (((a) & ((a) - 1)) == 0);
}
-is_power_of_two_large :: proc(a: ^Int) -> (res: bool) {
+int_is_power_of_two :: proc(a: ^Int) -> (res: bool) {
/*
Early out for Int == 0.
*/
@@ -63,7 +63,7 @@ is_power_of_two_large :: proc(a: ^Int) -> (res: bool) {
/*
For an `Int` to be a power of two, its top limb has to be a power of two.
*/
- if !is_power_of_two_small(int(a.digit[a.used - 1])) {
+ if !platform_int_is_power_of_two(int(a.digit[a.used - 1])) {
return false;
}
@@ -84,12 +84,11 @@ is_power_of_two_large :: proc(a: ^Int) -> (res: bool) {
}
return true;
}
-is_power_of_two :: proc{is_power_of_two_small, is_power_of_two_large};
/*
Compare two `Int`s, signed.
*/
-compare :: proc(a, b: ^Int) -> Comparison_Flag {
+int_compare :: proc(a, b: ^Int) -> Comparison_Flag {
if !is_initialized(a) { return .Uninitialized; }
if !is_initialized(b) { return .Uninitialized; }
@@ -105,35 +104,11 @@ compare :: proc(a, b: ^Int) -> Comparison_Flag {
}
return cmp_mag(x, y);
}
-cmp :: compare;
-
-/*
- Compare the magnitude of two `Int`s, unsigned.
-*/
-compare_magnitude :: proc(a, b: ^Int) -> Comparison_Flag {
- if !is_initialized(a) { return .Uninitialized; }
- if !is_initialized(b) { return .Uninitialized; }
-
- /* Compare based on used digits */
- if a.used != b.used {
- return .Greater_Than if a.used > b.used else .Less_Than;
- }
-
- /* Same number of used digits, compare based on their value */
- for n := a.used - 1; n >= 0; n -= 1 {
- if a.digit[n] != b.digit[n] {
- return .Greater_Than if a.digit[n] > b.digit[n] else .Less_Than;
- }
- }
-
- return .Equal;
-}
-cmp_mag :: compare_magnitude;
/*
Compare an `Int` to an unsigned number upto the size of the backing type.
*/
-compare_digit :: proc(a: ^Int, u: DIGIT) -> Comparison_Flag {
+int_compare_digit :: proc(a: ^Int, u: DIGIT) -> Comparison_Flag {
if !is_initialized(a) { return .Uninitialized; }
/* Compare based on sign */
@@ -152,4 +127,26 @@ compare_digit :: proc(a: ^Int, u: DIGIT) -> Comparison_Flag {
}
return .Equal;
+}
+
+/*
+ Compare the magnitude of two `Int`s, unsigned.
+*/
+int_compare_magnitude :: proc(a, b: ^Int) -> Comparison_Flag {
+ if !is_initialized(a) { return .Uninitialized; }
+ if !is_initialized(b) { return .Uninitialized; }
+
+ /* Compare based on used digits */
+ if a.used != b.used {
+ return .Greater_Than if a.used > b.used else .Less_Than;
+ }
+
+ /* Same number of used digits, compare based on their value */
+ for n := a.used - 1; n >= 0; n -= 1 {
+ if a.digit[n] != b.digit[n] {
+ return .Greater_Than if a.digit[n] > b.digit[n] else .Less_Than;
+ }
+ }
+
+ return .Equal;
} \ No newline at end of file
diff --git a/core/math/big/helpers.odin b/core/math/big/helpers.odin
index 9d21295df..933b330a5 100644
--- a/core/math/big/helpers.odin
+++ b/core/math/big/helpers.odin
@@ -4,7 +4,7 @@ package big
Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
Made available under Odin's BSD-2 license.
- A BigInt implementation in Odin.
+ An arbitrary precision mathematics implementation in Odin.
For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.
*/
@@ -13,13 +13,10 @@ import "core:mem"
import "core:intrinsics"
/*
- Deallocates the backing memory of an Int.
+ Deallocates the backing memory of an `Int`.
*/
-destroy :: proc(a: ^Int, allocator_zeroes := false, free_int := true, loc := #caller_location) {
- if !is_initialized(a) {
- // Nothing to do.
- return;
- }
+int_destroy :: proc(a: ^Int, allocator_zeroes := false, allocator := context.allocator) {
+ if !is_initialized(a) { return; }
if !allocator_zeroes {
mem.zero_slice(a.digit[:]);
@@ -27,15 +24,15 @@ destroy :: proc(a: ^Int, allocator_zeroes := false, free_int := true, loc := #ca
free(&a.digit[0]);
a.used = 0;
a.allocated = 0;
- if free_int {
- free(a);
- }
+ free(a);
}
/*
Creates and returns a new `Int`.
+ Size is the last parameter so it doesn't conflict with `int_init_from_integer`,
+ and we can say `init(512)` to init & set to 512.
*/
-init_new :: proc(allocator_zeroes := true, allocator := context.allocator, size := _DEFAULT_DIGIT_COUNT) -> (a: ^Int, err: Error) {
+int_init_sized :: proc(allocator_zeroes := true, allocator := context.allocator, size := _DEFAULT_DIGIT_COUNT) -> (a: ^Int, err: Error) {
/*
Allocating a new variable.
*/
@@ -46,10 +43,10 @@ init_new :: proc(allocator_zeroes := true, allocator := context.allocator, size
a.used = 0;
a.sign = .Zero_or_Positive;
- if len(a.digit) != size {
+ if len(a.digit) != int(size) {
return a, .Out_of_Memory;
}
- a.allocated = size;
+ a.allocated = int(size);
if !allocator_zeroes {
_zero_unused(a);
@@ -58,17 +55,17 @@ init_new :: proc(allocator_zeroes := true, allocator := context.allocator, size
}
/*
- Initialize from a signed or unsigned integer.
+ Initialize from a signed or unsigned platform integer.
Inits a new `Int` and then calls the appropriate `set` routine.
*/
-init_from_integer :: proc(src: $T, minimize := false, allocator_zeroes := true, allocator := context.allocator) -> (a: ^Int, err: Error) where intrinsics.type_is_integer(T) {
+int_init_from_integer :: proc(src: $T, minimize := false, allocator_zeroes := true, allocator := context.allocator) -> (a: ^Int, err: Error) where intrinsics.type_is_integer(T) {
n := _DEFAULT_DIGIT_COUNT;
if minimize {
n = _MIN_DIGIT_COUNT;
}
- a, err = init_new(allocator_zeroes, allocator, n);
+ a, err = init(allocator_zeroes, allocator, n);
if err == .OK {
set(a, src, minimize);
}
@@ -78,43 +75,41 @@ init_from_integer :: proc(src: $T, minimize := false, allocator_zeroes := true,
/*
Initialize an `Int` as a copy from another `Int`.
*/
-init_copy :: proc(src: ^Int, minimize := false, allocator_zeroes := true, allocator := context.allocator) -> (a: ^Int, err: Error) {
+int_init_copy :: proc(src: ^Int, minimize := false, allocator_zeroes := true, allocator := context.allocator) -> (a: ^Int, err: Error) {
if !is_initialized(src) {
return nil, .Invalid_Input;
}
- a, err = init_new(allocator_zeroes, allocator, src.used);
+ a, err = init(allocator_zeroes, allocator, src.used);
if err == .OK {
copy(a, src);
}
return;
}
-init :: proc{init_new, init_from_integer, init_copy};
-
/*
Helpers to set an `Int` to a specific value.
*/
-set_integer :: proc(a: ^Int, n: $T, minimize := false, loc := #caller_location) where intrinsics.type_is_integer(T) {
- n := n;
- assert_initialized(a, loc);
+int_set_from_integer :: proc(dest: ^Int, src: $T, minimize := false, loc := #caller_location) where intrinsics.type_is_integer(T) {
+ src := src;
+ assert_initialized(dest, loc);
- a.used = 0;
- a.sign = .Zero_or_Positive if n >= 0 else .Negative;
- n = abs(n);
+ dest.used = 0;
+ dest.sign = .Zero_or_Positive if src >= 0 else .Negative;
+ src = abs(src);
- for n != 0 {
- a.digit[a.used] = DIGIT(n) & _MASK;
- a.used += 1;
- n >>= _DIGIT_BITS;
+ for src != 0 {
+ dest.digit[dest.used] = DIGIT(src) & _MASK;
+ dest.used += 1;
+ src >>= _DIGIT_BITS;
}
if minimize {
- shrink(a);
+ shrink(dest);
}
- _zero_unused(a);
+ _zero_unused(dest);
}
-set :: proc{set_integer};
+set :: proc{int_set_from_integer};
/*
Copy one `Int` to another.
@@ -375,9 +370,6 @@ minus_one :: proc(a: ^Int, minimize := false) -> (err: Error) {
power_of_two :: proc(a: ^Int, power: int) -> (err: Error) {
assert_initialized(a);
- /*
-
- */
if power < 0 || power > _MAX_BIT_COUNT {
return .Invalid_Input;
}
diff --git a/core/math/big/log.odin b/core/math/big/log.odin
index d4b794344..1ce2437bf 100644
--- a/core/math/big/log.odin
+++ b/core/math/big/log.odin
@@ -4,7 +4,7 @@ package big
Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
Made available under Odin's BSD-2 license.
- A BigInt implementation in Odin.
+ An arbitrary precision mathematics implementation in Odin.
For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.
*/
diff --git a/core/math/big/logical.odin b/core/math/big/logical.odin
index 9ecb4fba3..6c857a2f2 100644
--- a/core/math/big/logical.odin
+++ b/core/math/big/logical.odin
@@ -4,7 +4,7 @@ package big
Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
Made available under Odin's BSD-2 license.
- A BigInt implementation in Odin.
+ An arbitrary precision mathematics implementation in Odin.
For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.
diff --git a/core/math/big/radix.odin b/core/math/big/radix.odin
index be1be7d10..86080d041 100644
--- a/core/math/big/radix.odin
+++ b/core/math/big/radix.odin
@@ -4,7 +4,7 @@ package big
Copyright 2021 Jeroen van Rijn <nom@duclavier.com>.
Made available under Odin's BSD-2 license.
- A BigInt implementation in Odin.
+ An arbitrary precision mathematics implementation in Odin.
For the theoretical underpinnings, see Knuth's The Art of Computer Programming, Volume 2, section 4.3.
The code started out as an idiomatic source port of libTomMath, which is in the public domain, with thanks.