🚀 go-pugleaf

In article <2324@l.cc.purdue.edu> cik@l.cc.purdue.edu (Herman Rubin) writes:
> The function is frexp(x,&n).  The argument x is a double precision floating
> point number.  The result is to be that number scaled down by a power of 
> 2 so that the absolute value is between .5 and 1, with the power stored in
> the location &n.

I suspect that something like this will work:

	double frexp(x,n)
	double x;
	int *n;
	{
		x = NORMALISE(x);
		*n = ((long *)&x)[EXPONENT_WORD] & EXPONENT_BITS;
		((long *)&x)[EXPONENT_WORD] &= ~EXPONENT_BITS;
		return x;
	}

You didn't say it had to be portable C. For most machines this will come
down to 2 or 3 instructions, given a decent compiler.
-- 
Peter da Silva.   `-_-'
+1 713 274 5180.
<peter@ficc.ferranti.com>