Use expansion for exp(z^2) * erfc(-z)

[jlapeyre]

Recently, we have avoided the particular form exp(z^2) * erfc(-z) for z < 0 and z large in magnitude, which suffers from overflow, by using the generic algorithm instead. But, this is much slower. Using a term or two in the asymptotic expansion is probably sufficient to avoid the overflow.