contrib/gcc/sreal.c

26 	x = sig * 2 ^ exp
30 	exp = exponent
44 	-SREAL_MAX_EXP <= exp && exp <= SREAL_MAX_EXP
50 	sig == 0 && exp == -SREAL_MAX_EXP
71 	   x->sig_hi, x->sig_lo, x->exp);
73   fprintf (file, "(" HOST_WIDE_INT_PRINT_UNSIGNED " * 2^%d)", x->sig, x->exp);
88   r->exp = a->exp;
102   gcc_assert (x->exp + s <= SREAL_MAX_EXP);
104   x->exp += s;
143       x->exp = -SREAL_MAX_EXP;
152 	  x->exp -= SREAL_PART_BITS;
158 	  x->exp--;
162       if (x->exp < -SREAL_MAX_EXP)
164 	  x->exp = -SREAL_MAX_EXP;
195       x->exp += shift;
206 	      x->exp++;
211       if (x->exp > SREAL_MAX_EXP)
213 	  x->exp = SREAL_MAX_EXP;
221       x->exp = -SREAL_MAX_EXP;
228 	  x->exp--;
233       if (x->exp < -SREAL_MAX_EXP)
235 	  x->exp = -SREAL_MAX_EXP;
246 	  x->exp++;
255 	  x->exp++;
259       if (x->exp > SREAL_MAX_EXP)
261 	  x->exp = SREAL_MAX_EXP;
271 sreal_init (sreal *r, unsigned HOST_WIDE_INT sig, signed int exp)
276   r->exp = exp - 16;
279   r->exp = exp;
291   if (r->exp <= -SREAL_BITS)
293   if (r->exp >= 0)
295   return ((r->sig_hi << SREAL_PART_BITS) + r->sig_lo) >> -r->exp;
297   if (r->exp <= -SREAL_BITS)
299   if (r->exp >= SREAL_PART_BITS)
301   if (r->exp > 0)
302     return r->sig << r->exp;
303   if (r->exp < 0)
304     return r->sig >> -r->exp;
314   if (a->exp > b->exp)
316   if (a->exp < b->exp)
353   dexp = a->exp - b->exp;
354   r->exp = a->exp;
401   dexp = a->exp - b->exp;
402   r->exp = a->exp;
450       r->exp = -SREAL_MAX_EXP;
463       r->exp = a->exp + b->exp + SREAL_PART_BITS;
484       r->exp = -SREAL_MAX_EXP;
489       r->exp = a->exp + b->exp;
509       r->exp = -SREAL_MAX_EXP;
534       r->exp = a->exp - b->exp - SREAL_BITS - SREAL_PART_BITS / 2;
540   r->exp = a->exp - b->exp - SREAL_PART_BITS;