44 lines
775 B
ArmAsm
44 lines
775 B
ArmAsm
|
/*
|
||
|
|
* ====================================================
|
||
|
|
* 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)
|
||
|
|
|
||
|
|
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
|