Lines Matching refs:r
36 unsigned r[4],*y,q,c; local
71 r[0] = px->significand[0];
72 r[1] = px->significand[1];
73 r[2] = px->significand[2];
74 r[3] = px->significand[3];
77 if(fpu_cmpli(r,y,4)>=0)
85 if(fpu_cmpli(r,y,4)>=0) {
88 c = fpu_sub3wc(&r[3],r[3],y[3],c);
89 c = fpu_sub3wc(&r[2],r[2],y[2],c);
90 c = fpu_sub3wc(&r[1],r[1],y[1],c);
91 c = fpu_sub3wc(&r[0],r[0],y[0],c);
93 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
94 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
95 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
96 r[3] = (r[3]<<1);
103 if(fpu_cmpli(r,y,4)>=0) {
106 c = fpu_sub3wc(&r[3],r[3],y[3],c);
107 c = fpu_sub3wc(&r[2],r[2],y[2],c);
108 c = fpu_sub3wc(&r[1],r[1],y[1],c);
109 c = fpu_sub3wc(&r[0],r[0],y[0],c);
111 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
112 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
113 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
114 r[3] = (r[3]<<1);
121 if(fpu_cmpli(r,y,4)>=0) {
124 c = fpu_sub3wc(&r[3],r[3],y[3],c);
125 c = fpu_sub3wc(&r[2],r[2],y[2],c);
126 c = fpu_sub3wc(&r[1],r[1],y[1],c);
127 c = fpu_sub3wc(&r[0],r[0],y[0],c);
129 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
130 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
131 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
132 r[3] = (r[3]<<1);
139 if(fpu_cmpli(r,y,4)>=0) {
142 c = fpu_sub3wc(&r[3],r[3],y[3],c);
143 c = fpu_sub3wc(&r[2],r[2],y[2],c);
144 c = fpu_sub3wc(&r[1],r[1],y[1],c);
145 c = fpu_sub3wc(&r[0],r[0],y[0],c);
147 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
148 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
149 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
150 r[3] = (r[3]<<1);
153 if((r[0]|r[1]|r[2]|r[3])==0) pz->sticky = pz->rounded = 0;
156 if(fpu_cmpli(r,y,4)>=0) pz->rounded = 1;
166 unsigned *x,r,c,q,t[4],s[4]; local
202 r = 0x00010000;
203 while(r!=0) { /* compute sqrt[0] */
204 t[0] = s[0]+r;
206 s[0] = t[0]+r;
208 q += r;
214 r>>=1;
218 r = 0x80000000;
219 while(r!=0) { /* compute sqrt[1] */
220 t[1] = s[1]+r; /* no carry */
224 c = fpu_add3wc(&s[1],t[1],r,c);
229 q += r;
235 r>>=1;
239 r = 0x80000000;
240 while(r!=0) { /* compute sqrt[2] */
241 t[2] = s[2]+r; /* no carry */
246 c = fpu_add3wc(&s[2],t[2],r,c);
253 q += r;
259 r>>=1;
263 r = 0x80000000;
264 while(r!=0) { /* compute sqrt[3] */
265 t[3] = s[3]+r; /* no carry */
271 c = fpu_add3wc(&s[3],t[3],r,c);
280 q += r;
286 r>>=1;