Mail Archives: djgpp-workers/2002/03/03/11:44:00
> From: Martin Str|mberg <ams AT ludd DOT luth DOT se>
> Date: Sun, 3 Mar 2002 15:19:10 +0100 (MET)
>
> How about this:
> Index: djgpp//src/libm/math/e_pow.c
> ===================================================================
> RCS file: /cvs/djgpp/djgpp/src/libm/math/e_pow.c,v
> retrieving revision 1.1
> diff -p -u -r1.1 e_pow.c
> --- djgpp//src/libm/math/e_pow.c 1998/10/06 10:04:10 1.1
> +++ djgpp//src/libm/math/e_pow.c 2002/03/03 14:17:46
> @@ -109,7 +109,7 @@ ivln2_l = 1.92596299112661746887e-08;
> double y_1,t1,t2,r,s,t,u,v,w;
> __int32_t i,j,k,yisint,n;
> __int32_t hx,hy,ix,iy;
> - __uint32_t lx,ly;
> + __uint32_t lx,ly,ju;
>
> EXTRACT_WORDS(hx,lx,x);
> EXTRACT_WORDS(hy,ly,y);
> @@ -134,8 +134,8 @@ ivln2_l = 1.92596299112661746887e-08;
> else if(iy>=0x3ff00000) {
> k = (iy>>20)-0x3ff; /* exponent */
> if(k>20) {
> - j = ly>>(52-k);
> - if((j<<(52-k))==ly) yisint = 2-(j&1);
> + ju = ly>>(52-k);
> + if((ju<<(52-k))==ly) yisint = 2-(ju&1);
This is fine, thanks.
- Raw text -