29 	int64_t n,hx,hy,hz,ix,iy,sx,i;  in __ieee754_fmodl()  local
32 	GET_LDOUBLE_WORDS64(hx,lx,x);  in __ieee754_fmodl()
34 	sx = hx&0x8000000000000000ULL;		/* sign of x */  in __ieee754_fmodl()
35 	hx ^=sx;				/* |x| */  in __ieee754_fmodl()
39 	if((hy|ly)==0||(hx>=0x7fff000000000000LL)|| /* y=0,or x not finite */  in __ieee754_fmodl()
42 	if(hx<=hy) {  in __ieee754_fmodl()
43 	    if((hx<hy)||(lx<ly)) return x;	/* |x|<|y| return x */  in __ieee754_fmodl()
49 	if(hx<0x0001000000000000LL) {	/* subnormal x */  in __ieee754_fmodl()
50 	    if(hx==0) {  in __ieee754_fmodl()
53 		for (ix = -16382, i=hx<<15; i>0; i<<=1) ix -=1;  in __ieee754_fmodl()
55 	} else ix = (hx>>48)-0x3fff;  in __ieee754_fmodl()
68 	    hx = 0x0001000000000000LL|(0x0000ffffffffffffLL&hx);  in __ieee754_fmodl()
72 		hx = (hx<<n)|(lx>>(64-n));  in __ieee754_fmodl()
75 		hx = lx<<(n-64);  in __ieee754_fmodl()
95 	    hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;  in __ieee754_fmodl()
96 	    if(hz<0){hx = hx+hx+(lx>>63); lx = lx+lx;}  in __ieee754_fmodl()
100 		hx = hz+hz+(lz>>63); lx = lz+lz;  in __ieee754_fmodl()
103 	hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;  in __ieee754_fmodl()
104 	if(hz>=0) {hx=hz;lx=lz;}  in __ieee754_fmodl()
107 	if((hx|lx)==0)			/* return sign(x)*0 */  in __ieee754_fmodl()
109 	while(hx<0x0001000000000000LL) {	/* normalize x */  in __ieee754_fmodl()
110 	    hx = hx+hx+(lx>>63); lx = lx+lx;  in __ieee754_fmodl()
114 	    hx = ((hx-0x0001000000000000LL)|((iy+16383)<<48));  in __ieee754_fmodl()
115 	    SET_LDOUBLE_WORDS64(x,hx|sx,lx);  in __ieee754_fmodl()
119 		lx = (lx>>n)|((uint64_t)hx<<(64-n));  in __ieee754_fmodl()
120 		hx >>= n;  in __ieee754_fmodl()
122 		lx = (hx<<(64-n))|(lx>>n); hx = sx;  in __ieee754_fmodl()
124 		lx = hx>>(n-64); hx = sx;  in __ieee754_fmodl()
126 	    SET_LDOUBLE_WORDS64(x,hx|sx,lx);  in __ieee754_fmodl()