kusano 7d535a
#include <stdio.h></stdio.h>
kusano 7d535a
kusano 7d535a
int main()
kusano 7d535a
{
kusano 7d535a
    /* Local variables */
kusano 7d535a
    float base, emin, prec, emax, rmin, rmax, t, sfmin;
kusano 7d535a
    extern float slamch_(char *);
kusano 7d535a
    float rnd, eps;
kusano 7d535a
kusano 7d535a
    eps = slamch_("Epsilon");
kusano 7d535a
    sfmin = slamch_("Safe minimum");
kusano 7d535a
    base = slamch_("Base");
kusano 7d535a
    prec = slamch_("Precision");
kusano 7d535a
    t = slamch_("Number of digits in mantissa");
kusano 7d535a
    rnd = slamch_("Rounding mode");
kusano 7d535a
    emin = slamch_("Minnimum exponent");
kusano 7d535a
    rmin = slamch_("Underflow threshold");
kusano 7d535a
    emax = slamch_("Largest exponent");
kusano 7d535a
    rmax = slamch_("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
}