Commit 50df67d3 authored by Michael Niedermayer's avatar Michael Niedermayer
Browse files

optimize bessel function instead of trusting gcc to do trivial optimizations (as gcc doesnt ...)

Originally committed as revision 8474 to svn://
parent 14f887ef
...@@ -71,9 +71,10 @@ static double bessel(double x){ ...@@ -71,9 +71,10 @@ static double bessel(double x){
double t=1; double t=1;
int i; int i;
x= x*x/4;
for(i=1; i<50; i++){ for(i=1; i<50; i++){
t *= i; t *= x/(i*i);
v += pow(x*x/4, i)/(t*t); v += t;
} }
return v; return v;
} }
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment