diff options
author | Norbert Preining <norbert@preining.info> | 2019-09-02 13:46:59 +0900 |
---|---|---|
committer | Norbert Preining <norbert@preining.info> | 2019-09-02 13:46:59 +0900 |
commit | e0c6872cf40896c7be36b11dcc744620f10adf1d (patch) | |
tree | 60335e10d2f4354b0674ec22d7b53f0f8abee672 /support/dktools/dk3mas.c |
Initial commit
Diffstat (limited to 'support/dktools/dk3mas.c')
-rw-r--r-- | support/dktools/dk3mas.c | 232 |
1 files changed, 232 insertions, 0 deletions
diff --git a/support/dktools/dk3mas.c b/support/dktools/dk3mas.c new file mode 100644 index 0000000000..4482aa85bf --- /dev/null +++ b/support/dktools/dk3mas.c @@ -0,0 +1,232 @@ +/* + WARNING: This file was generated by dkct. + 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 +*/ + +/* +Copyright (C) 2014-2017, Dirk Krause + +Redistribution and use in source and binary forms, with or without +modification, are permitted provided that the following conditions are met: + +* Redistributions of source code must retain the above copyright notice, + this list of conditions and the following disclaimer. +* Redistributions in binary form must reproduce the above opyright notice, + this list of conditions and the following disclaimer in the documentation + and/or other materials provided with the distribution. +* Neither the name of the author nor the names of contributors may be used + to endorse or promote products derived from this software without specific + prior written permission. + +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS ``AS IS'' +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE +ARE DISCLAIMED. +IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY +DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES +(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND +ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT +(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF +THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*/ + +/** @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; +} + |