Re: The old round off problem?

2006-03-05 Thread David Treadwell
On Mar 5, 2006, at 1:01 AM, sam wrote:David Treadwell wrote: exp(x) is implemented by:1.  reducing x into the range |r| <=  0.5 * ln(2), such that x = k *ln(2) + r2.  approximating exp(r) with a fifth-order polynomial,3.  re-scaling by multiplying by 2^k: exp(x) = 2^k * exp(r)sinh(x) is mathematica

Re: The old round off problem?

2006-03-04 Thread David Treadwell
I wish I knew! So I asked Google. Here's what I learned: Most implementations are based on, or similar to the implementation in the fdlibm package. sinh(x) and cosh(x) are both based on exp(x). See http:// www.netlib.org/cgi-bin/netlibfiles.pl?filename=/fdlibm/e_sinh.c exp(x) is implemented

Re: Python advocacy in scientific computation

2006-03-04 Thread David Treadwell
On Mar 4, 2006, at 11:16 PM, Dennis Lee Bieber wrote: > On Sat, 4 Mar 2006 14:23:10 -0500, David Treadwell > <[EMAIL PROTECTED]> declaimed the following in > comp.lang.python: > >> needed programming, be it CS or chemical engineering, taught it in >> a [EMAIL PRO

Re: The old round off problem?

2006-03-04 Thread David Treadwell
On Mar 4, 2006, at 4:33 PM, Paul Rubin wrote: > "sam" <[EMAIL PROTECTED]> writes: >> Hello all, I am taking a class in scientific programming at the local >> college. My problem is that the following difference produces >> round off >> errors as the value of x increases. For x >= 19 the diferenc

Re: Python advocacy in scientific computation

2006-03-04 Thread David Treadwell
On Mar 4, 2006, at 5:55 AM, Dennis Lee Bieber wrote: > On Fri, 3 Mar 2006 22:05:19 -0500, David Treadwell > <[EMAIL PROTECTED]> declaimed the following in > comp.lang.python: > > >> My ability to think of data structures was stunted BECAUSE of >> Fortran and BAS

Re: Python advocacy in scientific computation

2006-03-03 Thread David Treadwell
o far) has made sense to me, even if it took a few days of thought to figure it out. > > Now please go ahead and tell me how Python can help me become a better > scientist. And try to steer clear of the computer science buzzwords > that don't mean anyting to me. > My pleasure.