%% options copyright owner = Dirk Krause copyright year = 2015-xxxx license = bsd %% header /** @file Arithmetic operations on unsigned long long. CRT on Windows: Optional. */ #ifndef DK4CONF_H_INCLUDED #include "dk4conf.h" #endif #ifndef DK4TYPES_H_INCLUDED #include "dk4types.h" #endif #ifndef DK4ERROR_H_INCLUDED #include "dk4error.h" #endif #ifdef __cplusplus extern "C" { #endif #if DK4_HAVE_LONG_LONG /** Addition. @param a Left operand. @param b Right operand. @param erp Error report, may be NULL. @return Result on success, least significant result bits on overflow (see ec component of erp). Error codes: DK4_E_OVERFLOW. */ unsigned long long dk4ma_ull_add(unsigned long long a, unsigned long long b, dk4_er_t *erp); /** Subtraction. @param a Left operand. @param b Right operand. @param erp Error report, may be NULL. @return Result on success, least significant result bits on overflow (see ec component of erp). Error codes: DK4_E_OVERFLOW. */ unsigned long long dk4ma_ull_sub(unsigned long long a, unsigned long long b, dk4_er_t *erp); /** Multiplication. @param a Left factor. @param b right factor. @param erp Error report, may be NULL. @return Result on success, least significant bits on overflow (see ec component of erp). Error codes: DK4_E_OVERFLOW. */ unsigned long long dk4ma_ull_mul(unsigned long long a, unsigned long long b, dk4_er_t *erp); /** Division. @param a Counter (left operand). @param b Denominator (right operand). @param erp Error report, may be NULL. @return Division result on success, ULLONG_MAX on division by zero (see the ec component of erp). Error codes: DK4_E_DIVZERO. */ unsigned long long dk4ma_ull_div(unsigned long long a, unsigned long long b, dk4_er_t *erp); #endif #ifdef __cplusplus } #endif %% module #include "dk4conf.h" #include "dk4types.h" #include "dk4error.h" #include "dk4numco.h" #include "dk4maalu.h" $!trace-include #if DK4_HAVE_LONG_LONG unsigned long long dk4ma_ull_add(unsigned long long a, unsigned long long b, dk4_er_t *erp) { if ((ULLONG_MAX - a) < b) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } return (unsigned long long)(a + b); } unsigned long long dk4ma_ull_sub(unsigned long long a, unsigned long long b, dk4_er_t *erp) { if (b > a) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } return (unsigned long long)(a - b); } unsigned long long dk4ma_ull_mul(unsigned long long a, unsigned long long b, dk4_er_t *erp) { if (0ULL < a) { if ((ULLONG_MAX / a) < b) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } return (unsigned long long)(a * b); } unsigned long long dk4ma_ull_div(unsigned long long a, unsigned long long b, dk4_er_t *erp) { $? "+ dk4ma_ull_div %llu %llu", a, b if (0ULL < b) { $? "- dk4ma_ull_div (normal) %llu", (a / b) return (unsigned long long)(a / b); } else { dk4error_set_simple_error_code(erp, DK4_E_MATH_DIVZERO); $? "- dk4ma_ull_div (divzero) %llu", ULLONG_MAX return (unsigned long long)(ULLONG_MAX); } } #endif