| 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697 |
- ---
- ieee-utils/fp-win.c | 70 +++++++++++++++++++++++++++++++++++++++++++++++++++++
- ieee-utils/fp.c | 2 ++
- 2 files changed, 72 insertions(+)
- create mode 100644 ieee-utils/fp-win.c
- diff --git a/ieee-utils/fp-win.c b/ieee-utils/fp-win.c
- new file mode 100644
- index 0000000..e024eae
- --- /dev/null
- +++ b/ieee-utils/fp-win.c
- @@ -0,0 +1,70 @@
- +/* fp-win.c
- + *
- + * Author: Brian Gladman
- + *
- + * This program is free software; you can redistribute it and/or modify
- + * it under the terms of the GNU General Public License as published by
- + * the Free Software Foundation; either version 2 of the License, or (at
- + * your option) any later version.
- + *
- + * This program is distributed in the hope that it will be useful, but
- + * WITHOUT ANY WARRANTY; without even the implied warranty of
- + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- + * General Public License for more details.
- + *
- + * You should have received a copy of the GNU General Public License
- + * along with this program; if not, write to the Free Software
- + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
- + */
- +
- +#include <float.h>
- +
- +#include <config.h>
- +#include <gsl/gsl_ieee_utils.h>
- +#include <gsl/gsl_errno.h>
- +
- +const char *fp_env_string = "round-to-nearest,double-precision,mask-all";
- +
- +int
- +gsl_ieee_set_mode (int precision, int rounding, int exception_mask)
- +{
- + unsigned int old, mode = _DN_SAVE, mask = _MCW_DN | _MCW_RC | _MCW_EM;
- +
- + switch(precision)
- + {
- + case GSL_IEEE_SINGLE_PRECISION: mode |= _PC_24; break;
- + case GSL_IEEE_EXTENDED_PRECISION: mode |= _PC_64; break;
- + case GSL_IEEE_DOUBLE_PRECISION:
- + default: mode |= _PC_53;
- + }
- +#ifndef _M_AMD64
- + mask |= _MCW_PC;
- +#endif
- +
- + switch(rounding)
- + {
- + case GSL_IEEE_ROUND_DOWN: mode |= _RC_DOWN; break;
- + case GSL_IEEE_ROUND_UP: mode |= _RC_UP; break;
- + case GSL_IEEE_ROUND_TO_ZERO: mode |= _RC_CHOP; break;
- + case GSL_IEEE_ROUND_TO_NEAREST:
- + default: mode |= _RC_NEAR;
- + }
- +
- + if(exception_mask & GSL_IEEE_MASK_INVALID)
- + mode |= _EM_INVALID;
- + if(exception_mask & GSL_IEEE_MASK_DENORMALIZED)
- + mode |= _EM_DENORMAL;
- + if(exception_mask & GSL_IEEE_MASK_DIVISION_BY_ZERO)
- + mode |= _EM_ZERODIVIDE;
- + if(exception_mask & GSL_IEEE_MASK_OVERFLOW)
- + mode |= _EM_OVERFLOW;
- + if(exception_mask & GSL_IEEE_MASK_UNDERFLOW)
- + mode |= _EM_UNDERFLOW;
- + if(exception_mask & GSL_IEEE_TRAP_INEXACT)
- + mode &= ~_EM_INEXACT;
- + else
- + mode |= _EM_INEXACT;
- +
- + _controlfp_s( &old, mode, mask);
- + return GSL_SUCCESS;
- +}
- diff --git a/ieee-utils/fp.c b/ieee-utils/fp.c
- index 445a14f..b6ae5af 100644
- --- a/ieee-utils/fp.c
- +++ b/ieee-utils/fp.c
- @@ -45,6 +45,8 @@
- #endif
- #elif HAVE_DECL_FEENABLEEXCEPT || HAVE_DECL_FESETTRAPENABLE
- #include "fp-gnuc99.c"
- +#elif _MSC_VER
- +#include "fp-win.c"
- #else
- #include "fp-unknown.c"
- #endif
- --
|
