00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015
00016
00017
00018 static integer c__9 = 9;
00019 static integer c__1 = 1;
00020 static integer c__5 = 5;
00021
00022 int MAIN__(void)
00023 {
00024
00025 doublereal d__1;
00026
00027
00028 integer s_wsle(cilist *), do_lio(integer *, integer *, char *, ftnlen),
00029 e_wsle(void);
00030
00031
00032 doublereal t, rnd, eps, base, emin, prec, emax, rmin, rmax, sfmin;
00033 extern doublereal dlamch_(char *);
00034
00035
00036 static cilist io___11 = { 0, 6, 0, 0, 0 };
00037 static cilist io___12 = { 0, 6, 0, 0, 0 };
00038 static cilist io___13 = { 0, 6, 0, 0, 0 };
00039 static cilist io___14 = { 0, 6, 0, 0, 0 };
00040 static cilist io___15 = { 0, 6, 0, 0, 0 };
00041 static cilist io___16 = { 0, 6, 0, 0, 0 };
00042 static cilist io___17 = { 0, 6, 0, 0, 0 };
00043 static cilist io___18 = { 0, 6, 0, 0, 0 };
00044 static cilist io___19 = { 0, 6, 0, 0, 0 };
00045 static cilist io___20 = { 0, 6, 0, 0, 0 };
00046 static cilist io___21 = { 0, 6, 0, 0, 0 };
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060 eps = dlamch_("Epsilon");
00061 sfmin = dlamch_("Safe minimum");
00062 base = dlamch_("Base");
00063 prec = dlamch_("Precision");
00064 t = dlamch_("Number of digits in mantissa");
00065 rnd = dlamch_("Rounding mode");
00066 emin = dlamch_("Minimum exponent");
00067 rmin = dlamch_("Underflow threshold");
00068 emax = dlamch_("Largest exponent");
00069 rmax = dlamch_("Overflow threshold");
00070
00071 s_wsle(&io___11);
00072 do_lio(&c__9, &c__1, " Epsilon = ", (ftnlen)32);
00073 do_lio(&c__5, &c__1, (char *)&eps, (ftnlen)sizeof(doublereal));
00074 e_wsle();
00075 s_wsle(&io___12);
00076 do_lio(&c__9, &c__1, " Safe minimum = ", (ftnlen)32);
00077 do_lio(&c__5, &c__1, (char *)&sfmin, (ftnlen)sizeof(doublereal));
00078 e_wsle();
00079 s_wsle(&io___13);
00080 do_lio(&c__9, &c__1, " Base = ", (ftnlen)32);
00081 do_lio(&c__5, &c__1, (char *)&base, (ftnlen)sizeof(doublereal));
00082 e_wsle();
00083 s_wsle(&io___14);
00084 do_lio(&c__9, &c__1, " Precision = ", (ftnlen)32);
00085 do_lio(&c__5, &c__1, (char *)&prec, (ftnlen)sizeof(doublereal));
00086 e_wsle();
00087 s_wsle(&io___15);
00088 do_lio(&c__9, &c__1, " Number of digits in mantissa = ", (ftnlen)32);
00089 do_lio(&c__5, &c__1, (char *)&t, (ftnlen)sizeof(doublereal));
00090 e_wsle();
00091 s_wsle(&io___16);
00092 do_lio(&c__9, &c__1, " Rounding mode = ", (ftnlen)32);
00093 do_lio(&c__5, &c__1, (char *)&rnd, (ftnlen)sizeof(doublereal));
00094 e_wsle();
00095 s_wsle(&io___17);
00096 do_lio(&c__9, &c__1, " Minimum exponent = ", (ftnlen)32);
00097 do_lio(&c__5, &c__1, (char *)&emin, (ftnlen)sizeof(doublereal));
00098 e_wsle();
00099 s_wsle(&io___18);
00100 do_lio(&c__9, &c__1, " Underflow threshold = ", (ftnlen)32);
00101 do_lio(&c__5, &c__1, (char *)&rmin, (ftnlen)sizeof(doublereal));
00102 e_wsle();
00103 s_wsle(&io___19);
00104 do_lio(&c__9, &c__1, " Largest exponent = ", (ftnlen)32);
00105 do_lio(&c__5, &c__1, (char *)&emax, (ftnlen)sizeof(doublereal));
00106 e_wsle();
00107 s_wsle(&io___20);
00108 do_lio(&c__9, &c__1, " Overflow threshold = ", (ftnlen)32);
00109 do_lio(&c__5, &c__1, (char *)&rmax, (ftnlen)sizeof(doublereal));
00110 e_wsle();
00111 s_wsle(&io___21);
00112 do_lio(&c__9, &c__1, " Reciprocal of safe minimum = ", (ftnlen)32);
00113 d__1 = 1 / sfmin;
00114 do_lio(&c__5, &c__1, (char *)&d__1, (ftnlen)sizeof(doublereal));
00115 e_wsle();
00116
00117 return 0;
00118 }
00119
00120 int test3_ () { MAIN__ (); return 0; }