Hi all, I'm looking for a C or Fortran routine that computes log(erf(x)) without ever computing erf(x) directly, for use in PyMC. Does anyone have one laying around? It looks like GSL has such a thing, but it's GPL and we're using the MIT license. Thanks, Anand