Files
mercury/library/int64.m
Zoltan Somogyi 0fb9165ac3 Fix typo.
2026-03-07 00:32:26 +11:00

827 lines
22 KiB
Mathematica

%---------------------------------------------------------------------------%
% vim: ts=4 sw=4 et ft=mercury
%---------------------------------------------------------------------------%
% Copyright (C) 2018-2021, 2025-2026 The Mercury team.
% This file is distributed under the terms specified in COPYING.LIB.
%---------------------------------------------------------------------------%
%
% File: int64.m
% Main author: juliensf
% Stability: high.
%
% Predicates and functions for dealing with signed 64-bit integer numbers.
%
%---------------------------------------------------------------------------%
:- module int64.
:- interface.
:- import_module pretty_printer.
%---------------------------------------------------------------------------%
%
% Conversion from int.
%
% from_int(I) = I64:
%
% Convert an int to an int64.
%
% Since an int can be only 32 or 64 bits, this is guaranteed to yield
% a result that is mathematically equal to the original.
%
:- func from_int(int) = int64.
% cast_from_int(I) = I64:
%
% Convert an int to an int64.
%
% While a cast from int to intN for N =< 32 may yield a result
% that is not mathematically equal to the original (because
% the original integer may be too big to be representable),
% casting an int to int64 *will* yield a result that is mathematically
% equal to the original. It is therefore a synonym for the from_int
% function. It is provided only for uniformity, to allow an int
% to be cast to intN for *all* of int8, int16, int32 and int64.
%
:- func cast_from_int(int) = int64.
%---------------------------------------------------------------------------%
%
% Conversion to int.
%
% to_int(I64, I):
%
% Convert an int64 to a mathematically equivalent int.
%
% If ints are 64 bits, always succeed, since the two types
% are effectively identical.
%
% If ints are 32 bits, succeed only if I64 is in the range
% [int.min_int, int.max_int], meaning [-2^31, 2^31 - 1].
%
:- pred to_int(int64::in, int::out) is semidet.
% det_to_int(I64) = I:
%
% Convert an int64 into an int. If the given value of I64
% does not fit into an int, throw an exception.
%
:- func det_to_int(int64) = int.
% cast_to_int(I64) = I:
%
% Convert an int64 to an int.
% Always succeeds.
%
% - If ints are 64 bits, I will always be mathematically equal to I64.
%
% - If ints are 32 bits, then I will be mathematically equal to I64
% only if I64 is in the range [-(2^31), 2^31 - 1].
%
:- func cast_to_int(int64) = int.
%---------------------------------------------------------------------------%
%
% Change of signedness.
%
% cast_from_uint64(U64) = I64:
%
% Convert a uint64 to an int64. This will yield a result that is
% mathematically equal to U64 only if U64 is in the range [0, 2^63 - 1].
%
:- func cast_from_uint64(uint64) = int64.
%---------------------------------------------------------------------------%
%
% Conversion from byte sequence.
%
% from_bytes_le(Byte0, Byte1, ..., Byte7) = I64:
%
% I64 is the int64 whose bytes are given in little-endian order by the
% arguments from left-to-right (i.e. Byte0 is the least significant byte
% and Byte7 is the most significant byte).
%
:- func from_bytes_le(uint8, uint8, uint8, uint8, uint8, uint8, uint8, uint8)
= int64.
% from_bytes_be(Byte0, Byte1, ..., Byte7) = I64:
%
% I64 is the int64 whose bytes are given in big-endian order by the
% arguments in left-to-right order (i.e. Byte0 is the most significant
% byte and Byte7 is the least significant byte).
%
:- func from_bytes_be(uint8, uint8, uint8, uint8, uint8, uint8, uint8, uint8)
= int64.
%---------------------------------------------------------------------------%
%
% Comparisons and related operations.
%
% Less than.
%
:- pred (int64::in) < (int64::in) is semidet.
% Greater than.
%
:- pred (int64::in) > (int64::in) is semidet.
% Less than or equal.
%
:- pred (int64::in) =< (int64::in) is semidet.
% Greater than or equal.
%
:- pred (int64::in) >= (int64::in) is semidet.
% Maximum.
%
:- func max(int64, int64) = int64.
% Minimum.
%
:- func min(int64, int64) = int64.
% clamp(Min, Max, N):
%
% Clamp N to the range [Min, Max] (inclusive).
% Returns Min if N < Min, Max if N > Max, and N otherwise.
% Throws an exception if Max < Min.
%
:- func clamp(int64, int64, int64) = int64.
%---------------------------------------------------------------------------%
%
% Absolute values.
%
% abs(X) returns the absolute value of X.
% Throws an exception if X = int64.min_int64.
%
:- func abs(int64) = int64.
% unchecked_abs(X) returns the absolute value of X, except that the result
% is undefined if X = int64.min_int64.
%
:- func unchecked_abs(int64) = int64.
% nabs(X) returns the negative of the absolute value of X.
% Unlike abs/1 this function is defined for X = int64.min_int64.
%
:- func nabs(int64) = int64.
%---------------------------------------------------------------------------%
%
% Arithmetic operations.
%
% Unary plus.
%
:- func + (int64::in) = (int64::uo) is det.
% Unary minus.
%
:- func - (int64::in) = (int64::uo) is det.
% Addition.
%
:- func int64 + int64 = int64.
:- mode in + in = uo is det.
:- mode uo + in = in is det.
:- mode in + uo = in is det.
:- func plus(int64, int64) = int64.
% Subtraction.
%
:- func int64 - int64 = int64.
:- mode in - in = uo is det.
:- mode uo - in = in is det.
:- mode in - uo = in is det.
:- func minus(int64, int64) = int64.
% Multiplication.
%
:- func (int64::in) * (int64::in) = (int64::uo) is det.
:- func times(int64, int64) = int64.
% Flooring integer division.
% Truncates towards minus infinity, e.g. -10_i64 div 3_i64 = -4_i64.
%
% Throws a `domain_error' exception if the right operand is zero.
%
:- func (int64::in) div (int64::in) = (int64::uo) is det.
% Truncating integer division.
% Truncates towards zero, e.g. -10_i64 // 3_i64 = -3_i64.
% `div' has nicer mathematical properties for negative operands,
% but `//' is typically more efficient.
%
% Throws a `domain_error' exception if the right operand is zero.
%
:- func (int64::in) // (int64::in) = (int64::uo) is det.
% (/)/2 is a synonym for (//)/2.
%
:- func (int64::in) / (int64::in) = (int64::uo) is det.
% unchecked_quotient(X, Y) is the same as X // Y, but the behaviour
% is undefined if the right operand is zero.
%
:- func unchecked_quotient(int64::in, int64::in) = (int64::uo) is det.
% Modulus.
% X mod Y = X - (X div Y) * Y
%
% Throws a `domain_error' exception if the right operand is zero.
%
:- func (int64::in) mod (int64::in) = (int64::uo) is det.
% Remainder.
% X rem Y = X - (X // Y) * Y.
%
% Throws a `domain_error' exception if the right operand is zero.
%
:- func (int64::in) rem (int64::in) = (int64::uo) is det.
% unchecked_rem(X, Y) is the same as X rem Y, but the behaviour is
% undefined if the right operand is zero.
%
:- func unchecked_rem(int64::in, int64::in) = (int64::uo) is det.
% even(X) is equivalent to (X mod 2 = 0).
%
:- pred even(int64::in) is semidet.
% odd(X) is equivalent to (not even(X)), i.e. (X mod 2 = 1).
%
:- pred odd(int64::in) is semidet.
%---------------------------------------------------------------------------%
%
% Shift operations.
%
% Left shift.
% X << Y returns X "left shifted" by Y bits.
% The bit positions vacated by the shift are filled by zeros.
% Throws an exception if Y is not in the range [0, 64).
%
:- func (int64::in) << (int::in) = (int64::uo) is det.
:- func (int64::in) <<u (uint::in) = (int64::uo) is det.
% unchecked_left_shift(X, Y) is the same as X << Y except that
% the behaviour is undefined if Y is not in the range [0, 64).
% It will typically be implemented more efficiently than X << Y.
%
:- func unchecked_left_shift(int64::in, int::in) = (int64::uo) is det.
:- func unchecked_left_ushift(int64::in, uint::in) = (int64::uo) is det.
% Right shift.
% X >> Y returns X "right shifted" by Y bits.
% The bit positions vacated by the shift are filled by the sign bit.
% Throws an exception if Y is not in the range [0, 64).
%
:- func (int64::in) >> (int::in) = (int64::uo) is det.
:- func (int64::in) >>u (uint::in) = (int64::uo) is det.
% unchecked_right_shift(X, Y) is the same as X >> Y except that
% the behaviour is undefined if Y is not in the range [0, 64).
% It will typically be implemented more efficiently than X >> Y.
%
:- func unchecked_right_shift(int64::in, int::in) = (int64::uo) is det.
:- func unchecked_right_ushift(int64::in, uint::in) = (int64::uo) is det.
%---------------------------------------------------------------------------%
%
% Logical operations.
%
% Bitwise and.
%
:- func (int64::in) /\ (int64::in) = (int64::uo) is det.
% Bitwise or.
%
:- func (int64::in) \/ (int64::in) = (int64::uo) is det.
% Bitwise exclusive or (xor).
%
:- func xor(int64, int64) = int64.
:- mode xor(in, in) = uo is det.
:- mode xor(in, uo) = in is det.
:- mode xor(uo, in) = in is det.
% Bitwise complement.
%
:- func \ (int64::in) = (int64::uo) is det.
%---------------------------------------------------------------------------%
%
% Operations on bits and bytes.
%
% num_zeros(I) = N:
%
% N is the number of zeros in the binary representation of I.
%
:- func num_zeros(int64) = int.
% num_ones(I) = N:
%
% N is the number of ones in the binary representation of I.
%
:- func num_ones(int64) = int.
% num_leading_zeros(I) = N:
%
% N is the number of leading zeros in the binary representation of I,
% starting at the most significant bit position.
% Note that num_leading_zeros(0i64) = 64.
%
:- func num_leading_zeros(int64) = int.
% num_trailing_zeros(I) = N:
%
% N is the number of trailing zeros in the binary representation of I,
% starting at the least significant bit position.
% Note that num_trailing_zeros(0i64) = 64.
%
:- func num_trailing_zeros(int64) = int.
% reverse_bytes(A) = B:
%
% B is the value that results from reversing the bytes
% in the binary representation of A.
%
:- func reverse_bytes(int64) = int64.
% reverse_bits(A) = B:
%
% B is the value that results from reversing the bits
% in the binary representation of A.
%
:- func reverse_bits(int64) = int64.
%---------------------------------------------------------------------------%
%
% Limits.
%
:- func min_int64 = int64.
:- func max_int64 = int64.
%---------------------------------------------------------------------------%
%
% Prettyprinting.
%
% Convert an int64 to a pretty_printer.doc for formatting.
%
:- func int64_to_doc(int64) = pretty_printer.doc.
:- pragma obsolete(func(int64_to_doc/1), [pretty_printer.int64_to_doc/1]).
%---------------------------------------------------------------------------%
%---------------------------------------------------------------------------%
:- implementation.
:- import_module exception.
:- import_module int.
:- import_module require.
:- import_module uint.
:- import_module uint64.
%---------------------------------------------------------------------------%
from_int(I) = cast_from_int(I).
:- pragma foreign_proc("C",
cast_from_int(I::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe, will_not_modify_trail],
"
I64 = (int64_t) I;
").
:- pragma foreign_proc("C#",
cast_from_int(I::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I64 = (long) I; // Mercury's 'int' type in the C# grade is 32-bits.
").
:- pragma foreign_proc("Java",
cast_from_int(I::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I64 = (long) I; // Mercury's 'int' type in the Java grade is 32-bits.
").
%---------------------------------------------------------------------------%
:- pragma foreign_proc("C",
to_int(I64::in, I::out),
[will_not_call_mercury, promise_pure, thread_safe, will_not_modify_trail,
does_not_affect_liveness],
"
#if MR_BYTES_PER_WORD == 8
// Every bit in I64 means the same in I.
I = (MR_Integer) I64;
SUCCESS_INDICATOR = MR_TRUE;
#else
// This is the code of int32.from_int, but with different types.
//
// In int32.from_int,
// - the input is an int, which may be 32 or 64 bits, and
// - the output is an int32.
//
// Here,
// - the input is an int64, and
// - the output is an int, which happens to be 32 bits.
//
// So while the types are different, the bit counts are the same.
if (I64 > (int64_t) INT32_MAX) {
SUCCESS_INDICATOR = MR_FALSE;
} else if (I64 < (int64_t) INT32_MIN) {
SUCCESS_INDICATOR = MR_FALSE;
} else {
I = (MR_Integer) I64;
SUCCESS_INDICATOR = MR_TRUE;
}
#endif
").
:- pragma foreign_proc("C#",
to_int(I64::in, I::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
// In C#, Mercury's int is always 32 bits (int), and Mercury's int64 is
// a C# long. Check that the value fits in the range of a C# int.
if (I64 > (long) System.Int32.MaxValue) {
I = 0;
SUCCESS_INDICATOR = false;
} else if (I64 < (long) System.Int32.MinValue) {
I = 0;
SUCCESS_INDICATOR = false;
} else {
I = (int) I64;
SUCCESS_INDICATOR = true;
}
").
:- pragma foreign_proc("Java",
to_int(I64::in, I::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
// In Java, Mercury's int is always 32 bits (int), and Mercury's int64 is
// a Java long. Check that the value fits in the range of a Java int.
if (I64 > (long) java.lang.Integer.MAX_VALUE) {
I = 0;
SUCCESS_INDICATOR = false;
} else if (I64 < (long) java.lang.Integer.MIN_VALUE) {
I = 0;
SUCCESS_INDICATOR = false;
} else {
I = (int) I64;
SUCCESS_INDICATOR = true;
}
").
%---------------------------------------------------------------------------%
det_to_int(I64) = I :-
( if to_int(I64, IPrime) then
I = IPrime
else
error($pred, "cannot convert int64 to int")
).
%---------------------------------------------------------------------------%
:- pragma foreign_proc("C",
cast_to_int(I64::in) = (I::out),
[will_not_call_mercury, promise_pure, thread_safe, will_not_modify_trail],
"
I = (MR_Integer) I64;
").
:- pragma foreign_proc("C#",
cast_to_int(I64::in) = (I::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I = (int) I64;
").
:- pragma foreign_proc("Java",
cast_to_int(I64::in) = (I::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I = (int) I64;
").
%---------------------------------------------------------------------------%
:- pragma foreign_proc("C",
cast_from_uint64(U64::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe, will_not_modify_trail],
"
I64 = (int64_t) U64;
").
:- pragma foreign_proc("C#",
cast_from_uint64(U64::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I64 = (long) U64;
").
:- pragma foreign_proc("Java",
cast_from_uint64(U64::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I64 = (long) U64;
").
%---------------------------------------------------------------------------%
:- pragma foreign_proc("C",
from_bytes_le(Byte0::in, Byte1::in, Byte2::in, Byte3::in,
Byte4::in, Byte5::in, Byte6::in, Byte7::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe, will_not_modify_trail],
"
unsigned char *int64_bytes = (unsigned char *) &I64;
#if defined(MR_BIG_ENDIAN)
int64_bytes[0] = Byte7;
int64_bytes[1] = Byte6;
int64_bytes[2] = Byte5;
int64_bytes[3] = Byte4;
int64_bytes[4] = Byte3;
int64_bytes[5] = Byte2;
int64_bytes[6] = Byte1;
int64_bytes[7] = Byte0;
#else
int64_bytes[0] = Byte0;
int64_bytes[1] = Byte1;
int64_bytes[2] = Byte2;
int64_bytes[3] = Byte3;
int64_bytes[4] = Byte4;
int64_bytes[5] = Byte5;
int64_bytes[6] = Byte6;
int64_bytes[7] = Byte7;
#endif
").
:- pragma foreign_proc("Java",
from_bytes_le(Byte0::in, Byte1::in, Byte2::in, Byte3::in,
Byte4::in, Byte5::in, Byte6::in, Byte7::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I64 =
(long) (Byte7 & 0xff) << 56 |
(long) (Byte6 & 0xff) << 48 |
(long) (Byte5 & 0xff) << 40 |
(long) (Byte4 & 0xff) << 32 |
(long) (Byte3 & 0xff) << 24 |
(long) (Byte2 & 0xff) << 16 |
(long) (Byte1 & 0xff) << 8 |
(long) (Byte0 & 0xff);
").
:- pragma foreign_proc("C#",
from_bytes_le(Byte0::in, Byte1::in, Byte2::in, Byte3::in,
Byte4::in, Byte5::in, Byte6::in, Byte7::in) = (I64::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
I64 = (long) (
(ulong) Byte7 << 56 |
(ulong) Byte6 << 48 |
(ulong) Byte5 << 40 |
(ulong) Byte4 << 32 |
(ulong) Byte3 << 24 |
(ulong) Byte2 << 16 |
(ulong) Byte1 << 8 |
(ulong) Byte0);
").
from_bytes_be(Byte7, Byte6, Byte5,Byte4, Byte3, Byte2, Byte1, Byte0) =
from_bytes_le(Byte0, Byte1, Byte2, Byte3, Byte4, Byte5, Byte6, Byte7).
%---------------------------------------------------------------------------%
% The comparison operations <, >, =< and >= are builtins.
max(X, Y) =
( if X > Y then X else Y ).
min(X, Y) =
( if X < Y then X else Y ).
clamp(Min, Max, N) =
( if Max >= Min then
( if N < Min then Min else if N > Max then Max else N )
else
func_error($pred, "Max < Min")
).
%---------------------------------------------------------------------------%
abs(Num) = Abs :-
( if Num = int64.min_int64 then
error($pred, "abs(min_int64) would overflow")
else
Abs = unchecked_abs(Num)
).
unchecked_abs(Num) =
( if Num < 0i64 then
0i64 - Num
else
Num
).
nabs(Num) =
( if Num > 0i64 then
-Num
else
Num
).
%---------------------------------------------------------------------------%
% The operations + and - (both unary and binary), plus, minus, *, and times
% are builtins.
X div Y = Div :-
Trunc = X // Y,
( if
( X >= 0i64, Y >= 0i64
; X < 0i64, Y < 0i64
; X rem Y = 0i64
)
then
Div = Trunc
else
Div = Trunc - 1i64
).
:- pragma inline(func('//'/2)).
X // Y = Div :-
( if Y = 0i64 then
throw(domain_error("int64.'//': division by zero"))
else
Div = unchecked_quotient(X, Y)
).
:- pragma inline(func('/'/2)).
X / Y = X // Y.
% The operations unchecked_quotient and unchecked_rem are builtins.
X mod Y = X - (X div Y) * Y.
:- pragma inline(func(rem/2)).
X rem Y = Rem :-
( if Y = 0i64 then
throw(domain_error("int64.rem: division by zero"))
else
Rem = unchecked_rem(X, Y)
).
:- pragma inline(pred(even/1)).
even(X) :-
(X /\ 1i64) = 0i64.
:- pragma inline(pred(odd/1)).
odd(X) :-
(X /\ 1i64) \= 0i64.
%---------------------------------------------------------------------------%
% The unchecked shift operations are builtins.
X << Y = Result :-
( if cast_from_int(Y) < 64u then
Result = unchecked_left_shift(X, Y)
else
Msg = "int64.(<<): second operand is out of range",
throw(domain_error(Msg))
).
X <<u Y = Result :-
( if Y < 64u then
Result = unchecked_left_ushift(X, Y)
else
Msg = "int64.(<<u): second operand is out of range",
throw(domain_error(Msg))
).
X >> Y = Result :-
( if cast_from_int(Y) < 64u then
Result = unchecked_right_shift(X, Y)
else
Msg = "int64.(>>): second operand is out of range",
throw(domain_error(Msg))
).
X >>u Y = Result :-
( if Y < 64u then
Result = unchecked_right_ushift(X, Y)
else
Msg = "int64.(>>u): second operand is out of range",
throw(domain_error(Msg))
).
% The operations /\, \/, xor and \ are builtins.
%---------------------------------------------------------------------------%
num_zeros(U) = 64 - num_ones(U).
:- pragma foreign_proc("Java",
num_ones(U::in) = (N::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
N = java.lang.Long.bitCount(U);
").
num_ones(I64) = N :-
U64 = uint64.cast_from_int64(I64),
N = uint64.num_ones(U64).
%---------------------%
:- pragma foreign_proc("Java",
num_leading_zeros(U::in) = (N::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
N = java.lang.Long.numberOfLeadingZeros(U);
").
num_leading_zeros(I64) = N :-
U64 = uint64.cast_from_int64(I64),
N = uint64.num_leading_zeros(U64).
:- pragma foreign_proc("Java",
num_trailing_zeros(U::in) = (N::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
N = java.lang.Long.numberOfTrailingZeros(U);
").
num_trailing_zeros(I64) = N :-
U64 = uint64.cast_from_int64(I64),
N = uint64.num_trailing_zeros(U64).
%---------------------%
:- pragma foreign_proc("C",
reverse_bytes(A::in) = (B::out),
[will_not_call_mercury, promise_pure, thread_safe, will_not_modify_trail],
"
B = (int64_t) MR_uint64_reverse_bytes((uint64_t) A);
").
:- pragma foreign_proc("Java",
reverse_bytes(A::in) = (B::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
B = java.lang.Long.reverseBytes(A);
").
reverse_bytes(I64) = Result :-
U64 = uint64.cast_from_int64(I64),
Result0 = uint64.reverse_bytes(U64),
Result = int64.cast_from_uint64(Result0).
:- pragma foreign_proc("Java",
reverse_bits(A::in) = (B::out),
[will_not_call_mercury, promise_pure, thread_safe],
"
B = java.lang.Long.reverse(A);
").
reverse_bits(I64) = RevI64 :-
U64 = uint64.cast_from_int64(I64),
RevU64 = uint64.reverse_bits(U64),
RevI64 = int64.cast_from_uint64(RevU64).
%---------------------------------------------------------------------------%
min_int64 = -9_223_372_036_854_775_808_i64.
max_int64 = 9_223_372_036_854_775_807_i64.
%---------------------------------------------------------------------------%
int64_to_doc(I) = pretty_printer.int64_to_doc(I).
%---------------------------------------------------------------------------%
:- end_module int64.
%---------------------------------------------------------------------------%