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
ext::shared_ptr< YieldTermStructure > r