The R Project SVN R

Rev

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

Rev 19500 Rev 24735
Line 172... Line 172...
172 
 
 172 
 
173 
	    qexpo = ac * ac;
 173 
	    qexpo = ac * ac;
174 
	    if (qexpo > C3)
 174 
	    if (qexpo > C3)
175 
		break;
 175 
		break;
176 
 
 176 
 
177 
#ifdef OLD
 - 
 
178 
	    if (ac > 0.)	pplus = 1.0 + erf(ac / M_SQRT2);
 - 
 
179 
	    else		pplus = erfc(-(ac / M_SQRT2));
 - 
 
180 
 
 - 
 
181 
	    if (ac > w) pminus = 1.0 + erf((ac / M_SQRT2) - (w / M_SQRT2));
 - 
 
182 
	    else	pminus = erfc((w / M_SQRT2) - (ac / M_SQRT2));
 - 
 
183 
#else
 - 
 
184 
	    pplus = 2 * pnorm(ac, 0., 1., 1,0);
 177 
	    pplus = 2 * pnorm(ac, 0., 1., 1,0);
185 
	    pminus= 2 * pnorm(ac, w,  1., 1,0);
 178 
	    pminus= 2 * pnorm(ac, w,  1., 1,0);
186 
#endif
 - 
 
187 
 
 179 
 
188 
	    /* if rinsum ^ (cc-1) < 9e-14, */
 180 
	    /* if rinsum ^ (cc-1) < 9e-14, */
189 
	    /* then doesn't contribute to integral */
 181 
	    /* then doesn't contribute to integral */
190 
 
 182 
 
191 
	    rinsum = (pplus * 0.5) - (pminus * 0.5);
 183 
	    rinsum = (pplus * 0.5) - (pminus * 0.5);