%% options copyright owner = Dirk Krause copyright year = 2015-xxxx SPDX-License-Identifier: BSD-3-Clause %% header /** @file Arithmetic operations on dk4_um_t. CRT on Windows: Optional. */ #ifndef DK4CONF_H_INCLUDED #if DK4_BUILDING_DKTOOLS4 #include "dk4conf.h" #else #include #endif #endif #ifndef DK4TYPES_H_INCLUDED #if DK4_BUILDING_DKTOOLS4 #include "dk4types.h" #else #include #endif #endif #ifndef DK4ERROR_H_INCLUDED #if DK4_BUILDING_DKTOOLS4 #include "dk4error.h" #else #include #endif #endif #ifdef __cplusplus extern "C" { #endif /** 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. */ dk4_um_t dk4ma_um_add(dk4_um_t a, dk4_um_t 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. */ dk4_um_t dk4ma_um_sub(dk4_um_t a, dk4_um_t 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. */ dk4_um_t dk4ma_um_mul(dk4_um_t a, dk4_um_t 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, DK4_UM_MAX on division by zero (see the ec component of erp). Error codes: DK4_E_DIVZERO. */ dk4_um_t dk4ma_um_div(dk4_um_t a, dk4_um_t b, dk4_er_t *erp); /** Calculate greatest common divisor of two numbers. @param a One number. @param b Other number. @return Greatest common divisor. */ dk4_um_t dk4ma_um_gcd(dk4_um_t a, dk4_um_t b); /** Calculate least common multiple of two numbers. @param a One number. @param b Other number. @param erp Error report, may be NULL. @return Least common multiple of a and b. Error codes: DK4_E_OVERFLOW, DK4_E_DIVZERO. */ dk4_um_t dk4ma_um_lcm(dk4_um_t a, dk4_um_t b, dk4_er_t *erp); /** Comparison: Is a greater than b. @param a Left comparison operand. @param b Right comparison operand. @return 1 if a>b, 0 otherwise. */ int dk4ma_um_isgr(dk4_um_t a, dk4_um_t b); /** Comparison: Is a greater than or equal to b. @param a Left comparison operand. @param b Right comparison operand. @return 1 if a>=b, 0 otherwise. */ int dk4ma_um_isgeq(dk4_um_t a, dk4_um_t b); #ifdef __cplusplus } #endif %% module #include "dk4conf.h" #include "dk4types.h" #include "dk4error.h" #include "dk4numco.h" #if DK4_HAVE_INTMAX_T #include "dk4maamu.h" #else #if DK4_HAVE_LONG_LONG #include "dk4maalu.h" #else #include "dk4maaul.h" #endif #endif #include "dk4maadu.h" $!trace-include dk4_um_t dk4ma_um_add(dk4_um_t a, dk4_um_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_uintmax_t_add(a, b, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ull_add(a, b, erp)); #else return (dk4ma_ulong_add(a, b, erp)); #endif #endif } dk4_um_t dk4ma_um_sub(dk4_um_t a, dk4_um_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_uintmax_t_sub(a, b, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ull_sub(a, b, erp)); #else return (dk4ma_ulong_sub(a, b, erp)); #endif #endif } dk4_um_t dk4ma_um_mul(dk4_um_t a, dk4_um_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_uintmax_t_mul(a, b, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ull_mul(a, b, erp)); #else return (dk4ma_ulong_mul(a, b, erp)); #endif #endif } dk4_um_t dk4ma_um_div(dk4_um_t a, dk4_um_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_uintmax_t_div(a, b,erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ull_div(a, b, erp)); #else return (dk4ma_ulong_div(a, b, erp)); #endif #endif } dk4_um_t dk4ma_um_gcd(dk4_um_t a, dk4_um_t b) { dk4_um_t h; $? "+ dk4ma_um_gcd %lu %lu", (unsigned long)a, (unsigned long)b while (0UL < b) { h = a % b; a = b; b = h; } if (0UL == a) { a = 1UL; } $? "- dk4ma_um_gcd %lu", (unsigned long)a return a; } dk4_um_t dk4ma_um_lcm(dk4_um_t a, dk4_um_t b, dk4_er_t *erp) { return( dk4ma_um_mul( dk4ma_um_div(a, dk4ma_um_gcd(a, b), erp), b, erp ) ); } int dk4ma_um_isgr(dk4_um_t a, dk4_um_t b) { return ((a > b) ? ((int)1) : ((int)0)); } int dk4ma_um_isgeq(dk4_um_t a, dk4_um_t b) { return ((a >= b) ? ((int)1) : ((int)0)); } /* vim: set ai sw=4 ts=4 : */