/* Copyright (C) 2014-2020, Dirk Krause SPDX-License-Identifier: BSD-3-Clause */ /* WARNING: This file was generated by the dkct program (see http://dktools.sourceforge.net/ for details). Changes you make here will be lost if dkct is run again! You should modify the original source and run dkct on it. Original source: dk3mas.ctr */ /** @file dk3mas.c The dk3mas module. */ #line 108 "dk3mas.ctr" #include "dk3ma.h" #line 116 "dk3mas.ctr" short dk3ma_s_add_ok(short a, short b, int *ec) { if (NULL != ec) { if ((0 < a) && (0 < b)) { if((DK3_S_MAX - a) < b) { *ec = DK3_ERROR_MATH_OVERFLOW; } } else { if ((0 > a) && (0 > b)) { if (DK3_S_MIN == a) { *ec = DK3_ERROR_MATH_OVERFLOW; } else { if ((DK3_S_MIN - a) > b) { *ec = DK3_ERROR_MATH_OVERFLOW; } } } } } return (a + b); } short dk3ma_s_sub_ok(short a, short b, int *ec) { if (NULL != ec) { if ((0 < a) && (0 > b)) { if((DK3_S_MAX + b) < a) { *ec = DK3_ERROR_MATH_OVERFLOW; } } else { if ((0 > a) && (0 < b)) { if((DK3_S_MIN + b) > a) { *ec = DK3_ERROR_MATH_OVERFLOW; } } } } return (a - b); } short dk3ma_s_abs_ok(short a, int *ec) { if (NULL != ec) { if (DK3_S_MIN == a) { *ec = DK3_ERROR_MATH_OVERFLOW; } } return((0 <= a) ? a : ((short)(0 - a))); } #if VERSION_BEFORE_20141018 short dk3ma_s_mul_ok(short a, short b, int *ec) { if (NULL != ec) { if ((0 != a) && (0 != b)) { if ((DK3_S_MIN == a) && (1 != b)) { *ec = DK3_ERROR_MATH_OVERFLOW; } else { if ((DK3_S_MIN == b) && (1 != a)) { *ec = DK3_ERROR_MATH_OVERFLOW; } else { if ((DK3_S_MAX / dk3ma_s_abs_ok(a, ec)) < dk3ma_s_abs_ok(b, ec)) { *ec = DK3_ERROR_MATH_OVERFLOW; } } } } } return (a * b); } #else short dk3ma_s_mul_ok(short a, short b, int *ec) { if ((NULL != ec) && (0 != a) && (0 != b)) { if (DK3_S_MIN == a) { if (1 != b) { *ec = DK3_ERROR_MATH_OVERFLOW; } goto finished; } if (DK3_S_MIN == b) { if (1 != a) { *ec = DK3_ERROR_MATH_OVERFLOW; } goto finished; } if (0 < a) { if (0 < b) { if ((DK3_S_MAX / a) < b) { *ec = DK3_ERROR_MATH_OVERFLOW; } } else { if ((DK3_S_MIN / (0 - a)) > b) { *ec = DK3_ERROR_MATH_OVERFLOW; } } } else { if (0 < b) { if ((DK3_S_MIN / (0 - b)) > a) { *ec = DK3_ERROR_MATH_OVERFLOW; } } else { if ((DK3_S_MAX / (0 - a)) < (0 - b)) { *ec = DK3_ERROR_MATH_OVERFLOW; } } } } finished: return (a * b); } #endif short dk3ma_s_div_ok(short a, short b, int *ec) { short back = (short)0; if (0 != b) { if ((DK3_S_MIN == a) && (-1 == b)) { back = DK3_S_MAX; *ec = DK3_ERROR_MATH_OVERFLOW; } else { back = a / b; } } else { if (NULL != ec) { *ec = DK3_ERROR_MATH_DIVZERO; } if (0 <= a) { back = DK3_S_MAX; } else { back = DK3_S_MIN; } } return back; } short dk3ma_s_gcd_ok(short a, short b, int *ec) { short h; /* Check a */ if (0 > a) { if (DK3_S_MIN == a) { a = (short)0; if (NULL != ec) { *ec = DK3_ERROR_MATH_OVERFLOW; } goto finished; } else { a = (short)(0 - a); } } /* Check b */ if (0 > b) { if (DK3_S_MIN == b) { a = (short)0; if (NULL != ec) { *ec = DK3_ERROR_MATH_OVERFLOW; } goto finished; } else { b = (short)(0 - b); } } /* Calculate gcd. */ while (0 < b) { h = a % b; a = b; b = h; } if(0 == a) { a = (short)1; } /* Return result. */ finished: return a; }