21 R Y(float_property)(
const float_property p)
23 const R base = FLT_RADIX;
24 static R eps = K(1.0);
26 const R emin = MIN_EXP;
27 const R emax = MAX_EXP;
28 const R prec = eps * base;
29 static R rmin = K(1.0);
30 static R rmax = K(1.0);
31 const R rnd = FLTROUND;
32 static R sfmin = K(-1.0);
33 static short first = TRUE;
39 for (INT i=0; i<MANT_DIG-1; i++)
44 const INT n = 1 - MIN_EXP;
46 for (i = 0; i < n; i++)
54 for (i = 0; i < emax; i++)
60 R small = K(1.0) / rmax;
63 sfmin = small * (eps + K(1.0));
69 if (p == NFFT_EPSILON)
71 else if (p == NFFT_SAFE__MIN)
73 else if (p == NFFT_BASE)
75 else if (p == NFFT_PRECISION)
77 else if (p == NFFT_MANT_DIG)
79 else if (p == NFFT_FLTROUND)
81 else if (p == NFFT_E_MIN)
83 else if (p == NFFT_R_MIN)
85 else if (p == NFFT_E_MAX)
87 else if (p == NFFT_R_MAX)
96 R Y(prod_real)(R *vec, INT d)
102 for (t = 0; t < d; t++)
Internal header file for auxiliary definitions and functions.