gcc/fortran/arith.cc

136   gfc_real_info *real_info;
180   for (real_info = gfc_real_kinds; real_info->kind != 0; real_info++)
182       gfc_set_model_kind (real_info->kind);
189       mpfr_init (real_info->huge);
190       mpfr_set_ui (real_info->huge, 1, GFC_RND_MODE);
191       mpfr_set_ui (a, real_info->radix, GFC_RND_MODE);
192       mpfr_pow_si (a, a, -real_info->digits, GFC_RND_MODE);
193       mpfr_sub (real_info->huge, real_info->huge, a, GFC_RND_MODE);
196       mpfr_set_ui (a, real_info->radix, GFC_RND_MODE);
197       mpfr_pow_ui (a, a, real_info->max_exponent - 1, GFC_RND_MODE);
200       mpfr_mul (real_info->huge, real_info->huge, a, GFC_RND_MODE);
203       mpfr_mul_ui (real_info->huge, real_info->huge, real_info->radix,
207       mpfr_init (real_info->tiny);
208       mpfr_set_ui (real_info->tiny, real_info->radix, GFC_RND_MODE);
209       mpfr_pow_si (real_info->tiny, real_info->tiny,
210 		   real_info->min_exponent - 1, GFC_RND_MODE);
213       mpfr_init (real_info->subnormal);
214       mpfr_set_ui (real_info->subnormal, real_info->radix, GFC_RND_MODE);
215       mpfr_pow_si (real_info->subnormal, real_info->subnormal,
216 		   real_info->min_exponent - real_info->digits, GFC_RND_MODE);
219       mpfr_init (real_info->epsilon);
220       mpfr_set_ui (real_info->epsilon, real_info->radix, GFC_RND_MODE);
221       mpfr_pow_si (real_info->epsilon, real_info->epsilon,
222 		   1 - real_info->digits, GFC_RND_MODE);
225       mpfr_log10 (a, real_info->huge, GFC_RND_MODE);
226       mpfr_log10 (b, real_info->tiny, GFC_RND_MODE);
232       real_info->range = (int) mpfr_get_si (a, GFC_RND_MODE);
235       mpfr_set_ui (a, real_info->radix, GFC_RND_MODE);
237       mpfr_mul_ui (a, a, real_info->digits - 1, GFC_RND_MODE);
239       real_info->precision = (int) mpfr_get_si (a, GFC_RND_MODE);
242       for (i = 10; i <= real_info->radix; i *= 10)
243 	if (i == real_info->radix)
244 	  real_info->precision++;