diff options
Diffstat (limited to 'newlib/libm/machine/i386/f_frexp.S')
| -rw-r--r-- | newlib/libm/machine/i386/f_frexp.S | 44 |
1 files changed, 44 insertions, 0 deletions
diff --git a/newlib/libm/machine/i386/f_frexp.S b/newlib/libm/machine/i386/f_frexp.S new file mode 100644 index 000000000..58397b79c --- /dev/null +++ b/newlib/libm/machine/i386/f_frexp.S @@ -0,0 +1,44 @@ +/* + * ==================================================== + * Copyright (C) 1998 by Cygnus Solutions. All rights reserved. + * + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +/* +Fast version of frexp using Intel float instructions. + + double _f_frexp (double x, int *exp); + +Function splits x into y * 2 ** z. It then +returns the value of y and updates *exp with z. +There is no error checking or setting of errno. +*/ + + #include "i386mach.h" + + .global SYM (_f_frexp) + SOTYPE_FUNCTION(_f_frexp) + +SYM (_f_frexp): + pushl ebp + movl esp,ebp + fldl 8(ebp) + movl 16(ebp),eax + + fxtract + fld1 + fchs + fxch + fscale + fstp st1 + fxch + fld1 + faddp + fistpl 0(eax) + + leave + ret |