%% options copyright owner = Dirk Krause copyright year = 2015-xxxx SPDX-License-Identifier: BSD-3-Clause %% header /** @file Arithmetic operations on short. 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, SHRT_MAX if a is SHRT_MIN (overflow is reported in erp). Error codes: DK4_E_OVERFLOW. */ short dk4ma_short_abs(short 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. */ short dk4ma_short_add(short a, short 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. */ short dk4ma_short_sub(short a, short 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. */ short dk4ma_short_mul(short a, short 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, SHRT_MIN, SHRT_MAX, or 0 on division by zero or overflow (see the ec component of erp). Error codes: DK4_E_OVERFLOW, DK4_E_DIVZERO. */ short dk4ma_short_div(short a, short b, dk4_er_t *erp); /** Convert dk4_im_t value to short. @param i Value to convert. @param erp Error report, may be NULL. @return Conversion result. Error codes: DK4_E_OVERFLOW. */ short dk4ma_short_from(dk4_im_t i, dk4_er_t *erp); /** Convert double to short. @param d Value to convert. @param erp Error report, may be NULL. @return Conversion result. Error codes: DK4_E_OVERFLOW. */ short dk4ma_short_from_double(double d, dk4_er_t *erp); #ifdef __cplusplus } #endif %% module #include "dk4conf.h" #include "dk4types.h" #include "dk4error.h" #include "dk4numco.h" #include "dk4maas.h" short dk4ma_short_abs(short a, dk4_er_t *erp) { if (SHRT_MIN == a) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } return ( ((short)0 <= a) ? (a) : ((SHRT_MIN == a) ? (SHRT_MAX) : ((short)0 - a)) ); } short dk4ma_short_add(short a, short b, dk4_er_t *erp) { if ((0 < a) && (0 < b)) { if ((SHRT_MAX - a) < b) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } else { if ((0 > a) && (0 > b)) { if (SHRT_MIN == a) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } else { if ((SHRT_MIN - a) > b) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } } } return (short)(a + b); } short dk4ma_short_sub(short a, short b, dk4_er_t *erp) { if ((0 < a) && (0 > b)) { if ((SHRT_MAX + b) < a) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } else { if((0 > a) && (0 < b)) { if ((SHRT_MIN + b) > a) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } } return (short)(a - b); } short dk4ma_short_mul(short a, short b, dk4_er_t *erp) { if ((0 != a) && (0 != b)) { if (SHRT_MIN == a) { if (1 != b) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } else { if (SHRT_MIN == b) { if (1 != a) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } else { if (0 < a) { if (0 < b) { if ((SHRT_MAX / a) < b) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } else { if ((SHRT_MIN / a) > b) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } } else { if (0 < b) { if ((SHRT_MIN / b) > a) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } else { if ((SHRT_MAX / (0 - a)) < (0 - b)) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } } } } } } return (short)(a * b); } short dk4ma_short_div(short a, short b, dk4_er_t *erp) { short back = 0; if (0 != b) { if ((SHRT_MIN == a) && (-1 == b)) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); back = SHRT_MAX; } else { back = a / b; } } else { dk4error_set_simple_error_code(erp, DK4_E_MATH_DIVZERO); if (0 < a) { back = SHRT_MAX; } else { if (0 > a) { back = SHRT_MIN; } } } return back; } short dk4ma_short_from(dk4_im_t i, dk4_er_t *erp) { if (((dk4_im_t)(SHRT_MAX) < i) || ((dk4_im_t)(SHRT_MIN) > i)) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); } return ((short)i); } short dk4ma_short_from_double(double d, dk4_er_t *erp) { short back = 0; if ((double)(SHRT_MIN) > d) { dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); back = SHRT_MIN; } else { if ((double)(SHRT_MAX) < d) { dk4error_set_simple_error_code(erp,DK4_E_MATH_OVERFLOW); back = SHRT_MAX; } else { back = (short)d; } } return back; }