def cdfNormal(z):
if (abs(z) < ERF_CODY_LIMIT1):<EOL><INDENT>return <NUM_LIT:0.5> * (<NUM_LIT:1.0> + (z / M_SQRT2) * _erfRationalHelperR3(<NUM_LIT:0.5> * z * z))<EOL><DEDENT>elif (z < <NUM_LIT:0.0>):<EOL><INDENT>return np.exp(logPdfNormal(z)) * _erfRationalHelper(-z) / (-z)<EOL><DEDENT>else:<EOL><INDENT>return <NUM_LIT:1.0> - np.exp(logPdfNormal(z)) * _erfRationalHelper(z) / z<EOL><DEDENT>
Robust implementations of cdf of a standard normal.

 @see [[https://github.com/mseeger/apbsint/blob/master/src/eptools/potentials/SpecfunServices.h original implementation]]
 in C from Matthias Seeger.
*/