Good point, so we need a better solution that fixes all cases

>> I'll file a ticket.
>>
>> Incidentally, if tanh(z) is simply programmed as
>>
>> (1.0 - exp(-2.0*z)) / (1.0 + exp(-2.0*z))

>This will overflow as z -> -\infty. The solution is probably to use a
>different expression for Re(z) < 0, and to check how other libraries do
>this in case the above still misses something.
>
>       Pauli
_______________________________________________
NumPy-Discussion mailing list
[email protected]
http://mail.scipy.org/mailman/listinfo/numpy-discussion

Reply via email to