Mercurial > mxe-octave
view src/build-msvctools/math/modff.c @ 7207:9ed6500e56d3 default tip @
maint: Merge release to default.
author | Markus Mützel <markus.muetzel@gmx.de> |
---|---|
date | Fri, 17 May 2024 20:16:41 +0200 |
parents | f8299bb6c872 |
children |
line wrap: on
line source
#include <fenv.h> #include <math.h> #include <errno.h> #define FE_ROUNDING_MASK \ (FE_TONEAREST | FE_DOWNWARD | FE_UPWARD | FE_TOWARDZERO) float modff (float value, float* iptr) { float int_part; unsigned short saved_cw; unsigned short tmp_cw; /* truncate */ asm ("fnstcw %0;" : "=m" (saved_cw)); /* save control word */ tmp_cw = (saved_cw & ~FE_ROUNDING_MASK) | FE_TOWARDZERO; asm ("fldcw %0;" : : "m" (tmp_cw)); asm ("frndint;" : "=t" (int_part) : "0" (value)); /* round */ asm ("fldcw %0;" : : "m" (saved_cw)); /* restore saved cw */ if (iptr) *iptr = int_part; return (isinf (value) ? 0.0F : value - int_part); }