[ROOT] problem in TMath::BesselI

From: Peter Steinberg (peter.steinberg@bnl.gov)
Date: Wed May 07 2003 - 17:55:31 MEST


hello -

in TMath::BesselI(Int_t n, Double_t x)

there is a line at the beginning of the function:

if (TMath::Abs(x) < kBigPositive) return 0;

this essentially renders the function useless for
all practical applications, since in the n=0,1
cases is simply uses TMath::BesselI0 and I1.

is this simply a bug?  can it be fixed?

thanks,
peter

-- 
Dr. Peter Steinberg

Fulbright Visting Professor
Department of Physics, University of Cape Town, South Africa
(currently on leave from Brookhaven National Laboratory, USA)

Office: +21 650 5091, Mobile: +83 426 9743 (+ = 011 27 from the US)
Email: peter.steinberg@bnl.gov, Web: http://hep.phy.uct.ac.za/~steinber



This archive was generated by hypermail 2b29 : Thu Jan 01 2004 - 17:50:11 MET