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
}