The R Project SVN R

Rev

Rev 13117 | Rev 36986 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 13117 Rev 13434
Line 50... Line 50...
50 
   ----------------------------------------------------------------------*/
 50 
   ----------------------------------------------------------------------*/
51 
 
 51 
 
52 
/* ----------------------------------------------------------------------
 52 
/* ----------------------------------------------------------------------
53 
   Mathematical constants
 53 
   Mathematical constants
54 
   ----------------------------------------------------------------------*/
 54 
   ----------------------------------------------------------------------*/
55 
    static double sqrtpi = .9189385332046727417803297; /* == ??? */
 55 
    const double sqrtpi = .9189385332046727417803297; /* == ??? */
56 
 
 56 
 
57 
/* *******************************************************************
 57 
/* *******************************************************************
58 
 
 58 
 
59 
   Explanation of machine-dependent constants
 59 
   Explanation of machine-dependent constants
60 
 
 60 
 
Line 103... Line 103...
103 
   Machine dependent parameters
 103 
   Machine dependent parameters
104 
   ----------------------------------------------------------------------
 104 
   ----------------------------------------------------------------------
105 
   */
 105 
   */
106 
 
 106 
 
107 
 
 107 
 
108 
    static double xbig = 171.624;
 108 
    const double xbig = 171.624;
109 
    /* ML_POSINF ==   static double xinf = 1.79e308;*/
 109 
    /* ML_POSINF ==   const double xinf = 1.79e308;*/
110 
    /* DBL_EPSILON = static double eps = 2.22e-16;*/
 110 
    /* DBL_EPSILON = const double eps = 2.22e-16;*/
111 
    /* DBL_MIN ==   static double xminin = 2.23e-308;*/
 111 
    /* DBL_MIN ==   const double xminin = 2.23e-308;*/
112 
 
 112 
 
113 
    /*----------------------------------------------------------------------
 113 
    /*----------------------------------------------------------------------
114 
      Numerator and denominator coefficients for rational minimax
 114 
      Numerator and denominator coefficients for rational minimax
115 
      approximation over (1,2).
 115 
      approximation over (1,2).
116 
      ----------------------------------------------------------------------*/
 116 
      ----------------------------------------------------------------------*/
117 
    static double p[8] = {
 117 
    const double p[8] = {
118 
	-1.71618513886549492533811,
 118 
	-1.71618513886549492533811,
119 
	24.7656508055759199108314,-379.804256470945635097577,
 119 
	24.7656508055759199108314,-379.804256470945635097577,
120 
	629.331155312818442661052,866.966202790413211295064,
 120 
	629.331155312818442661052,866.966202790413211295064,
121 
	-31451.2729688483675254357,-36144.4134186911729807069,
 121 
	-31451.2729688483675254357,-36144.4134186911729807069,
122 
	66456.1438202405440627855 };
 122 
	66456.1438202405440627855 };
123 
    static double q[8] = {
 123 
    const double q[8] = {
124 
	-30.8402300119738975254353,
 124 
	-30.8402300119738975254353,
125 
	315.350626979604161529144,-1015.15636749021914166146,
 125 
	315.350626979604161529144,-1015.15636749021914166146,
126 
	-3107.77167157231109440444,22538.1184209801510330112,
 126 
	-3107.77167157231109440444,22538.1184209801510330112,
127 
	4755.84627752788110767815,-134659.959864969306392456,
 127 
	4755.84627752788110767815,-134659.959864969306392456,
128 
	-115132.259675553483497211 };
 128 
	-115132.259675553483497211 };
129 
    /*----------------------------------------------------------------------
 129 
    /*----------------------------------------------------------------------
130 
      Coefficients for minimax approximation over (12, INF).
 130 
      Coefficients for minimax approximation over (12, INF).
131 
      ----------------------------------------------------------------------*/
 131 
      ----------------------------------------------------------------------*/
132 
    static double c[7] = {
 132 
    const double c[7] = {
133 
	-.001910444077728,8.4171387781295e-4,
 133 
	-.001910444077728,8.4171387781295e-4,
134 
	-5.952379913043012e-4,7.93650793500350248e-4,
 134 
	-5.952379913043012e-4,7.93650793500350248e-4,
135 
	-.002777777777777681622553,.08333333333333333331554247,
 135 
	-.002777777777777681622553,.08333333333333333331554247,
136 
	.0057083835261 };
 136 
	.0057083835261 };
137 
 
 137