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

Originally committed as revision 8474 to svn://svn.ffmpeg.org/ffmpeg/trunk

diff --git a/libavcodec/resample2.c b/libavcodec/resample2.c
index ebda7e5..70d8a1d 100644 (file)
--- a/libavcodec/resample2.c
+++ b/libavcodec/resample2.c
@@ -71,9 +71,10 @@ static double bessel(double x){
     double t=1;
     int i;
 
+    x= x*x/4;
     for(i=1; i<50; i++){
-        t *= i;
-        v += pow(x*x/4, i)/(t*t);
+        t *= x/(i*i);
+        v += t;
     }
     return v;
 }