Lines Matching refs:r
36 unsigned r[4], *y, q, c;
72 r[0] = px->significand[0];
73 r[1] = px->significand[1];
74 r[2] = px->significand[2];
75 r[3] = px->significand[3];
78 if (fpu_cmpli(r, y, 4) >= 0)
86 if (fpu_cmpli(r, y, 4) >= 0) {
87 q += 1; /* if r>y do r-=y and q+=1 */
89 c = fpu_sub3wc(&r[3], r[3], y[3], c);
90 c = fpu_sub3wc(&r[2], r[2], y[2], c);
91 c = fpu_sub3wc(&r[1], r[1], y[1], c);
92 c = fpu_sub3wc(&r[0], r[0], y[0], c);
94 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
95 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
96 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
97 r[3] = (r[3]<<1);
104 if (fpu_cmpli(r, y, 4) >= 0) {
105 q += 1; /* if r>y do r-=y and q+=1 */
107 c = fpu_sub3wc(&r[3], r[3], y[3], c);
108 c = fpu_sub3wc(&r[2], r[2], y[2], c);
109 c = fpu_sub3wc(&r[1], r[1], y[1], c);
110 c = fpu_sub3wc(&r[0], r[0], y[0], c);
112 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
113 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
114 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
115 r[3] = (r[3]<<1);
122 if (fpu_cmpli(r, y, 4) >= 0) {
123 q += 1; /* if r>y do r-=y and q+=1 */
125 c = fpu_sub3wc(&r[3], r[3], y[3], c);
126 c = fpu_sub3wc(&r[2], r[2], y[2], c);
127 c = fpu_sub3wc(&r[1], r[1], y[1], c);
128 c = fpu_sub3wc(&r[0], r[0], y[0], c);
130 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
131 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
132 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
133 r[3] = (r[3]<<1);
140 if (fpu_cmpli(r, y, 4) >= 0) {
141 q += 1; /* if r>y do r-=y and q+=1 */
143 c = fpu_sub3wc(&r[3], r[3], y[3], c);
144 c = fpu_sub3wc(&r[2], r[2], y[2], c);
145 c = fpu_sub3wc(&r[1], r[1], y[1], c);
146 c = fpu_sub3wc(&r[0], r[0], y[0], c);
148 r[0] = (r[0]<<1)|((r[1]&0x80000000)>>31); /* r << 1 */
149 r[1] = (r[1]<<1)|((r[2]&0x80000000)>>31);
150 r[2] = (r[2]<<1)|((r[3]&0x80000000)>>31);
151 r[3] = (r[3]<<1);
154 if ((r[0]|r[1]|r[2]|r[3]) == 0) pz->sticky = pz->rounded = 0;
157 if (fpu_cmpli(r, y, 4) >= 0) pz->rounded = 1;
167 unsigned *x, r, c, q, t[4], s[4];
209 r = 0x00010000;
210 while (r != 0) { /* compute sqrt[0] */
211 t[0] = s[0] + r;
213 s[0] = t[0] + r;
215 q += r;
221 r >>= 1;
225 r = (unsigned)0x80000000;
226 while (r != 0) { /* compute sqrt[1] */
227 t[1] = s[1] + r; /* no carry */
231 c = fpu_add3wc(&s[1], t[1], r, c);
236 q += r;
242 r >>= 1;
246 r = (unsigned)0x80000000;
247 while (r != 0) { /* compute sqrt[2] */
248 t[2] = s[2] + r; /* no carry */
253 c = fpu_add3wc(&s[2], t[2], r, c);
260 q += r;
266 r >>= 1;
270 r = (unsigned)0x80000000;
271 while (r != 0) { /* compute sqrt[3] */
272 t[3] = s[3] + r; /* no carry */
278 c = fpu_add3wc(&s[3], t[3], r, c);
287 q += r;
293 r >>= 1;