27#include <ql/math/errorfunction.hpp>
198 Real R,
S,P,Q,s,y,z,r, ax;
200 if (!std::isfinite(x)) {
204 return ( x > 0 ? 1 : -1);
210 if(ax < 3.7252902984e-09) {
212 return 0.125*(8.0*x+
efx8*x);
225 if(x>=0)
return erx + P/Q;
else return -
erx - P/Q;
234 if(ax < 2.85714285714285) {
241 r = std::exp( -ax*ax-0.5625 +R/
S);
242 if(x>=0)
return one-r/ax;
else return r/ax-
one;
Real operator()(Real x) const