Distribution::Normal.cdf(x) was returning 0.5 for absolute values over 38, probably due to some overflow or numeric error. These cases should be covered by the changed lines, except that these lines weren't actually doing anything other than evaluating and carrying on.
Distribution::Normal.cdf(x) was returning 0.5 for absolute values over 38, probably due to some overflow or numeric error. These cases should be covered by the changed lines, except that these lines weren't actually doing anything other than evaluating and carrying on.