Author mark.dickinson
Recipients mark.dickinson, ned.deily, nirinA, rhettinger, steven.daprano, stutzbach, terry.reedy, tim.peters
Date 2009-12-13.12:53:43
SpamBayes Score 4.46264e-06
Marked as misclassified No
Message-id <1260708828.34.0.963757970672.issue3366@psf.upfronthosting.co.za>
In-reply-to
Content
nirinA: thanks for the erf patch.  It's fine as far as it goes;  the main 
technical issue is that the series development only converges at a 
sensible rate for smallish x;  for larger x something more is needed.

Here's a prototype patch for the erf and erfc functions that uses a power 
series (not quite the same as the series nirinA used;  by pulling out a 
factor of exp(-x*2) you can make all coefficients in the remaining series 
positive, which is good for the numerics) for small x, and a continued 
fraction expansion for larger x.   The erf and erfc functions are 
currently written in Python.  I'll do some more testing (and possibly 
tweaking) of this patch and then translate the Python to C.
History
Date User Action Args
2009-12-13 12:53:50mark.dickinsonsetrecipients: + mark.dickinson, tim.peters, rhettinger, terry.reedy, ned.deily, stutzbach, nirinA, steven.daprano
2009-12-13 12:53:48mark.dickinsonsetmessageid: <1260708828.34.0.963757970672.issue3366@psf.upfronthosting.co.za>
2009-12-13 12:53:46mark.dickinsonlinkissue3366 messages
2009-12-13 12:53:45mark.dickinsoncreate