|
kusano |
7d535a |
#include <stdio.h></stdio.h>
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
main()
|
|
kusano |
7d535a |
{
|
|
kusano |
7d535a |
/* Local variables */
|
|
kusano |
7d535a |
double base, emin, prec, emax, rmin, rmax, t, sfmin;
|
|
kusano |
7d535a |
extern double dlamch_(char *);
|
|
kusano |
7d535a |
double rnd, eps;
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
eps = dlamch_("Epsilon");
|
|
kusano |
7d535a |
sfmin = dlamch_("Safe minimum");
|
|
kusano |
7d535a |
base = dlamch_("Base");
|
|
kusano |
7d535a |
prec = dlamch_("Precision");
|
|
kusano |
7d535a |
t = dlamch_("Number of digits in mantissa");
|
|
kusano |
7d535a |
rnd = dlamch_("Rounding mode");
|
|
kusano |
7d535a |
emin = dlamch_("Minnimum exponent");
|
|
kusano |
7d535a |
rmin = dlamch_("Underflow threshold");
|
|
kusano |
7d535a |
emax = dlamch_("Largest exponent");
|
|
kusano |
7d535a |
rmax = dlamch_("Overflow threshold");
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
printf(" Epsilon = %e\n", eps);
|
|
kusano |
7d535a |
printf(" Safe minimum = %e\n", sfmin);
|
|
kusano |
7d535a |
printf(" Base = %.0f\n", base);
|
|
kusano |
7d535a |
printf(" Precision = %e\n", prec);
|
|
kusano |
7d535a |
printf(" Number of digits in mantissa = %.0f\n", t);
|
|
kusano |
7d535a |
printf(" Rounding mode = %.0f\n", rnd);
|
|
kusano |
7d535a |
printf(" Minimum exponent = %.0f\n", emin);
|
|
kusano |
7d535a |
printf(" Underflow threshold = %e\n", rmin);
|
|
kusano |
7d535a |
printf(" Largest exponent = %.0f\n", emax);
|
|
kusano |
7d535a |
printf(" Overflow threshold = %e\n", rmax);
|
|
kusano |
7d535a |
printf(" Reciprocal of safe minimum = %e\n", 1./sfmin);
|
|
kusano |
7d535a |
return 0;
|
|
kusano |
7d535a |
}
|