1 #include <stdio.h>
2
3
main()4 main()
5 {
6 /* Local variables */
7 double base, emin, prec, emax, rmin, rmax, t, sfmin;
8 extern double dlamch_(char *);
9 double rnd, eps;
10
11 eps = dlamch_("Epsilon");
12 sfmin = dlamch_("Safe minimum");
13 base = dlamch_("Base");
14 prec = dlamch_("Precision");
15 t = dlamch_("Number of digits in mantissa");
16 rnd = dlamch_("Rounding mode");
17 emin = dlamch_("Minnimum exponent");
18 rmin = dlamch_("Underflow threshold");
19 emax = dlamch_("Largest exponent");
20 rmax = dlamch_("Overflow threshold");
21
22 printf(" Epsilon = %e\n", eps);
23 printf(" Safe minimum = %e\n", sfmin);
24 printf(" Base = %.0f\n", base);
25 printf(" Precision = %e\n", prec);
26 printf(" Number of digits in mantissa = %.0f\n", t);
27 printf(" Rounding mode = %.0f\n", rnd);
28 printf(" Minimum exponent = %.0f\n", emin);
29 printf(" Underflow threshold = %e\n", rmin);
30 printf(" Largest exponent = %.0f\n", emax);
31 printf(" Overflow threshold = %e\n", rmax);
32 printf(" Reciprocal of safe minimum = %e\n", 1./sfmin);
33 return 0;
34 }
35