def f(x,returntype=LINALGMAXFLOAT): ''' evaluates f(x)=1/(1+exp(-x)) as accurately as possible ''' return returntype(sexp(-np.logaddexp(ZERO128,-np.float128(x))))