%% options copyright owner = Dirk Krause copyright year = 2015-xxxx SPDX-License-Identifier: BSD-3-Clause %% header /** @file Arithmetic operations on dk4_im_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 /** Absolute value. @param a Operand. @param erp Error report, may be NULL. @return Result on success, DK4_IM_MAX if a is DK4_IM_MIN (overflow is reported in erp). Error codes: DK4_E_OVERFLOW. */ dk4_im_t dk4ma_im_abs(dk4_im_t a, dk4_er_t *erp); /** 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_im_t dk4ma_im_add(dk4_im_t a, dk4_im_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_im_t dk4ma_im_sub(dk4_im_t a, dk4_im_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_im_t dk4ma_im_mul(dk4_im_t a, dk4_im_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_IM_MIN, DK4_IM_MAX, or 0 on division by zero or overflow (see the ec component of erp). Error codes: DK4_E_OVERFLOW, DK4_E_DIVZERO. */ dk4_im_t dk4ma_im_div(dk4_im_t a, dk4_im_t b, dk4_er_t *erp); /** Calculate greatest common divisor. @param a One number. @param b Other number. @param erp Error report, may be NULL. @return Greatest common divisor of a and b. Error codes: DK4_E_OVERFLOW, DK4_E_DIVZERO. */ dk4_im_t dk4ma_im_gcd(dk4_im_t a, dk4_im_t b, dk4_er_t *erp); /** Calculate least common multiple. @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_im_t dk4ma_im_lcm(dk4_im_t a, dk4_im_t b, dk4_er_t *erp); #ifdef __cplusplus } #endif %% module #include "dk4conf.h" #include "dk4types.h" #include "dk4error.h" #include "dk4numco.h" #if DK4_HAVE_INTMAX_T #include "dk4maami.h" #else #if DK4_HAVE_LONG_LONG #include "dk4maall.h" #else #include "dk4maal.h" #endif #endif #include "dk4maadi.h" #ifndef DK4MAADU_H_INCLUDED #include "dk4maadu.h" #endif dk4_im_t dk4ma_im_abs(dk4_im_t a, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_intmax_t_abs(a, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ll_abs(a, erp)); #else return (dk4ma_long_abs(a, erp)); #endif #endif } dk4_im_t dk4ma_im_add(dk4_im_t a, dk4_im_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_intmax_t_add(a, b, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ll_add(a, b, erp)); #else return (dk4ma_long_add(a, b, erp)); #endif #endif } dk4_im_t dk4ma_im_sub(dk4_im_t a, dk4_im_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_intmax_t_sub(a, b, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ll_sub(a, b, erp)); #else return (dk4ma_long_sub(a, b, erp)); #endif #endif } dk4_im_t dk4ma_im_mul(dk4_im_t a, dk4_im_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_intmax_t_mul(a, b, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ll_mul(a, b, erp)); #else return (dk4ma_long_mul(a, b, erp)); #endif #endif } dk4_im_t dk4ma_im_div(dk4_im_t a, dk4_im_t b, dk4_er_t *erp) { #if DK4_HAVE_INTMAX_T return (dk4ma_intmax_t_div(a, b, erp)); #else #if DK4_HAVE_LONG_LONG return (dk4ma_ll_div(a, b, erp)); #else return (dk4ma_long_div(a, b, erp)); #endif #endif } dk4_im_t dk4ma_im_gcd(dk4_im_t a, dk4_im_t b, dk4_er_t *erp) { dk4_im_t back; dk4_um_t v; v = dk4ma_um_gcd(dk4ma_im_abs(a, erp), dk4ma_im_abs(b, erp)); back = (dk4_im_t)v; if ((dk4_um_t)(DK4_IM_MAX) < v) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } return back; } dk4_im_t dk4ma_im_lcm(dk4_im_t a, dk4_im_t b, dk4_er_t *erp) { return ( dk4ma_im_mul( dk4ma_im_div(dk4ma_im_abs(a, erp), dk4ma_im_gcd(a, b, erp), erp), dk4ma_im_abs(b, erp), erp ) ); } /* vim: set ai sw=4 ts=4 : */