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/dk4fpe.c |
Initial commit
Diffstat (limited to 'support/dktools/dk4fpe.c')
-rw-r--r-- | support/dktools/dk4fpe.c | 224 |
1 files changed, 224 insertions, 0 deletions
diff --git a/support/dktools/dk4fpe.c b/support/dktools/dk4fpe.c new file mode 100644 index 0000000000..c1c392ac1f --- /dev/null +++ b/support/dktools/dk4fpe.c @@ -0,0 +1,224 @@ +/* + 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: dk4fpe.ctr +*/ + +/* +Copyright (C) 2016-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 dk4fpe.c The dk4fpe module. +*/ + + +#line 200 "dk4fpe.ctr" + +#include "dk4fpe.h" + + + +#ifdef _WIN32 + +/** Exception set for errors. +*/ +static const unsigned int dk4fpe_errors = +_SW_ZERODIVIDE | _SW_INVALID | _SW_OVERFLOW | _SW_UNDERFLOW; + +#else + +/** Exception set for errors. +*/ +static const int dk4fpe_errors = +FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW | FE_UNDERFLOW; + +#endif + + + + +int +dk4fpe_check_error(dk4_fpe_t fpe, dk4_er_t *erp) +{ + int back = 0; + if ((dk4_fpe_t)0 != (dk4fpe_errors & fpe)) { + back = 1; +#ifdef _WIN32 + if ((dk4_fpe_t)0 != (_SW_INVALID & fpe)) +#else + if ((dk4_fpe_t)0 != (FE_INVALID & fpe)) +#endif + { + dk4error_set_simple_error_code(erp, DK4_E_MATH_INVALID); + } + else + { +#ifdef _WIN32 + if ((dk4_fpe_t)0 != (_SW_ZERODIVIDE & fpe)) +#else + if ((dk4_fpe_t)0 != (FE_DIVBYZERO & fpe)) +#endif + { + dk4error_set_simple_error_code(erp, DK4_E_MATH_DIVZERO); + } + else + { +#ifdef _WIN32 + if ((dk4_fpe_t)0 != (_SW_OVERFLOW & fpe)) +#else + if ((dk4_fpe_t)0 != (FE_OVERFLOW & fpe)) +#endif + { + dk4error_set_simple_error_code(erp, DK4_E_MATH_OVERFLOW); + } + else + { +#ifdef _WIN32 + if ((dk4_fpe_t)0 != (_SW_UNDERFLOW & fpe)) +#else + if ((dk4_fpe_t)0 != (FE_UNDERFLOW & fpe)) +#endif + { + dk4error_set_simple_error_code(erp, DK4_E_MATH_UNDERFLOW); + } + } + } + } + } + return back; +} + + +int +dk4fpe_clear(dk4_er_t *erp) +{ +#ifdef _WIN32 + _clearfp(); + return 1; +#else +#if DK4_HAVE_FETESTEXCEPT + int back = 0; + if (0 == feclearexcept(FE_ALL_EXCEPT)) { + back = 1; + } else { + dk4error_set_simple_error_code(erp, DK4_E_SYSTEM); + } + return back; +#else + dk4error_set_simple_error_code(erp, DK4_E_NOT_SUPPORTED); + return 0; +#endif +#endif +} + + + +int +dk4fpe_get(dk4_fpe_t *dptr, dk4_er_t *erp) +{ + int back = 0; + if (NULL != dptr) { +#ifdef _WIN32 + *dptr = (dk4_fpe_t)_statusfp(); + back = 1; +#else +#if DK4_HAVE_FETESTEXCEPT + *dptr = (dk4_fpe_t)fetestexcept(FE_ALL_EXCEPT); + back = 1; +#else + dk4error_set_simple_error_code(erp, DK4_E_NOT_SUPPORTED); +#endif +#endif + } else { + dk4error_set_simple_error_code(erp, DK4_E_INVALID_ARGUMENTS); + } + return back; +} + + + +int +dk4fpe_get_clear(dk4_fpe_t *dptr, dk4_er_t *erp) +{ + int back = 0; + if (NULL != dptr) { +#ifdef _WIN32 + *dptr = _clearfp(); + back = 1; +#else +#if DK4_HAVE_FETESTEXCEPT + *dptr = (dk4_fpe_t)fetestexcept(FE_ALL_EXCEPT); + if (0 == feclearexcept(FE_ALL_EXCEPT)) { + back = 1; + } else { + dk4error_set_simple_error_code(erp, DK4_E_SYSTEM); + } +#else + dk4error_set_simple_error_code(erp, DK4_E_NOT_SUPPORTED); +#endif +#endif + } else { + dk4error_set_simple_error_code(erp, DK4_E_INVALID_ARGUMENTS); + } + return back; +} + + + +int +dk4fpe_is_finite(double d) +{ + int back = 0; +#ifdef _WIN32 + back = ((0 != _finite(d)) ? 1 : 0); +#else + back = ((0 != isfinite(d)) ? 1 : 0); +#endif + return back; +} + + + +int +dk4fpe_is_normal(double d) +{ + int back = 0; +#ifdef _WIN32 + switch (_fpclass(d)) { + case _FPCLASS_NN : case _FPCLASS_NZ : + case _FPCLASS_PZ : case _FPCLASS_PN : + { + back = 1; + } break; + } +#else + back = ((0 != isnormal(d)) ? 1 : 0); +#endif + return back; +} + |