Mercurial > mplayer.hg
view drivers/generic_math.h @ 23518:11b0be5cea0a
capital char and relative dot at end of phrase removed
author | ptt |
---|---|
date | Sun, 10 Jun 2007 17:58:53 +0000 |
parents | 544914dadb52 |
children | 3baf6a2283da |
line wrap: on
line source
/* generic implementation of sin(x) and cos(x) functions. Specially for Linux. Licence: GPL Copyright (C) 2002 Nick Kurshev */ #ifndef GENERIC_MATH_H #define GENERIC_MATH_H typedef struct _gen_sincos { double x; double sinx; double cosx; }gen_sincos_t; static gen_sincos_t g_sincos[201] = { { -3.141600e+00, 7.346410e-06, -1.000000e-00 }, { -3.110184e+00, -3.140349e-02, -9.995068e-01 }, { -3.078768e+00, -6.278333e-02, -9.980272e-01 }, { -3.047352e+00, -9.410122e-02, -9.955626e-01 }, { -3.015936e+00, -1.253262e-01, -9.921156e-01 }, { -2.984520e+00, -1.564276e-01, -9.876894e-01 }, { -2.953104e+00, -1.873745e-01, -9.822885e-01 }, { -2.921688e+00, -2.181366e-01, -9.759183e-01 }, { -2.890272e+00, -2.486833e-01, -9.685848e-01 }, { -2.858856e+00, -2.789847e-01, -9.602956e-01 }, { -2.827440e+00, -3.090107e-01, -9.510586e-01 }, { -2.796024e+00, -3.387318e-01, -9.408830e-01 }, { -2.764608e+00, -3.681185e-01, -9.297789e-01 }, { -2.733192e+00, -3.971420e-01, -9.177572e-01 }, { -2.701776e+00, -4.257736e-01, -9.048297e-01 }, { -2.670360e+00, -4.539849e-01, -8.910094e-01 }, { -2.638944e+00, -4.817483e-01, -8.763097e-01 }, { -2.607528e+00, -5.090362e-01, -8.607451e-01 }, { -2.576112e+00, -5.358217e-01, -8.443312e-01 }, { -2.544696e+00, -5.620785e-01, -8.270839e-01 }, { -2.513280e+00, -5.877805e-01, -8.090204e-01 }, { -2.481864e+00, -6.129025e-01, -7.901586e-01 }, { -2.450448e+00, -6.374196e-01, -7.705169e-01 }, { -2.419032e+00, -6.613076e-01, -7.501148e-01 }, { -2.387616e+00, -6.845430e-01, -7.289724e-01 }, { -2.356200e+00, -7.071029e-01, -7.071107e-01 }, { -2.324784e+00, -7.289649e-01, -6.845511e-01 }, { -2.293368e+00, -7.501075e-01, -6.613159e-01 }, { -2.261952e+00, -7.705099e-01, -6.374281e-01 }, { -2.230536e+00, -7.901518e-01, -6.129112e-01 }, { -2.199120e+00, -8.090140e-01, -5.877894e-01 }, { -2.167704e+00, -8.270777e-01, -5.620876e-01 }, { -2.136288e+00, -8.443252e-01, -5.358310e-01 }, { -2.104872e+00, -8.607395e-01, -5.090457e-01 }, { -2.073456e+00, -8.763043e-01, -4.817579e-01 }, { -2.042040e+00, -8.910044e-01, -4.539948e-01 }, { -2.010624e+00, -9.048251e-01, -4.257835e-01 }, { -1.979208e+00, -9.177528e-01, -3.971521e-01 }, { -1.947792e+00, -9.297748e-01, -3.681288e-01 }, { -1.916376e+00, -9.408793e-01, -3.387421e-01 }, { -1.884960e+00, -9.510552e-01, -3.090212e-01 }, { -1.853544e+00, -9.602925e-01, -2.789953e-01 }, { -1.822128e+00, -9.685821e-01, -2.486940e-01 }, { -1.790712e+00, -9.759158e-01, -2.181473e-01 }, { -1.759296e+00, -9.822865e-01, -1.873854e-01 }, { -1.727880e+00, -9.876877e-01, -1.564385e-01 }, { -1.696464e+00, -9.921142e-01, -1.253372e-01 }, { -1.665048e+00, -9.955616e-01, -9.411219e-02 }, { -1.633632e+00, -9.980265e-01, -6.279433e-02 }, { -1.602216e+00, -9.995064e-01, -3.141450e-02 }, { -1.570800e+00, -1.000000e-00, -3.673205e-06 }, { -1.539384e+00, -9.995067e-01, 3.140716e-02 }, { -1.507968e+00, -9.980269e-01, 6.278700e-02 }, { -1.476552e+00, -9.955623e-01, 9.410488e-02 }, { -1.445136e+00, -9.921151e-01, 1.253299e-01 }, { -1.413720e+00, -9.876889e-01, 1.564312e-01 }, { -1.382304e+00, -9.822879e-01, 1.873781e-01 }, { -1.350888e+00, -9.759175e-01, 2.181402e-01 }, { -1.319472e+00, -9.685839e-01, 2.486869e-01 }, { -1.288056e+00, -9.602945e-01, 2.789882e-01 }, { -1.256640e+00, -9.510574e-01, 3.090142e-01 }, { -1.225224e+00, -9.408817e-01, 3.387352e-01 }, { -1.193808e+00, -9.297775e-01, 3.681220e-01 }, { -1.162392e+00, -9.177557e-01, 3.971454e-01 }, { -1.130976e+00, -9.048282e-01, 4.257769e-01 }, { -1.099560e+00, -8.910077e-01, 4.539882e-01 }, { -1.068144e+00, -8.763079e-01, 4.817515e-01 }, { -1.036728e+00, -8.607433e-01, 5.090393e-01 }, { -1.005312e+00, -8.443292e-01, 5.358248e-01 }, { -9.738960e-01, -8.270819e-01, 5.620815e-01 }, { -9.424800e-01, -8.090183e-01, 5.877835e-01 }, { -9.110640e-01, -7.901563e-01, 6.129054e-01 }, { -8.796480e-01, -7.705146e-01, 6.374224e-01 }, { -8.482320e-01, -7.501124e-01, 6.613104e-01 }, { -8.168160e-01, -7.289699e-01, 6.845457e-01 }, { -7.854000e-01, -7.071081e-01, 7.071055e-01 }, { -7.539840e-01, -6.845484e-01, 7.289674e-01 }, { -7.225680e-01, -6.613131e-01, 7.501100e-01 }, { -6.911520e-01, -6.374252e-01, 7.705122e-01 }, { -6.597360e-01, -6.129083e-01, 7.901541e-01 }, { -6.283200e-01, -5.877864e-01, 8.090161e-01 }, { -5.969040e-01, -5.620845e-01, 8.270798e-01 }, { -5.654880e-01, -5.358279e-01, 8.443272e-01 }, { -5.340720e-01, -5.090425e-01, 8.607414e-01 }, { -5.026560e-01, -4.817547e-01, 8.763061e-01 }, { -4.712400e-01, -4.539915e-01, 8.910060e-01 }, { -4.398240e-01, -4.257802e-01, 9.048266e-01 }, { -4.084080e-01, -3.971488e-01, 9.177542e-01 }, { -3.769920e-01, -3.681254e-01, 9.297762e-01 }, { -3.455760e-01, -3.387387e-01, 9.408805e-01 }, { -3.141600e-01, -3.090177e-01, 9.510563e-01 }, { -2.827440e-01, -2.789917e-01, 9.602935e-01 }, { -2.513280e-01, -2.486905e-01, 9.685830e-01 }, { -2.199120e-01, -2.181437e-01, 9.759166e-01 }, { -1.884960e-01, -1.873817e-01, 9.822872e-01 }, { -1.570800e-01, -1.564348e-01, 9.876883e-01 }, { -1.256640e-01, -1.253335e-01, 9.921147e-01 }, { -9.424800e-02, -9.410853e-02, 9.955619e-01 }, { -6.283200e-02, -6.279067e-02, 9.980267e-01 }, { -3.141600e-02, -3.141083e-02, 9.995066e-01 }, { 0.000000e+00, 0.000000e+00, 1.000000e+00 }, { 3.141600e-02, 3.141083e-02, 9.995066e-01 }, { 6.283200e-02, 6.279067e-02, 9.980267e-01 }, { 9.424800e-02, 9.410853e-02, 9.955619e-01 }, { 1.256640e-01, 1.253335e-01, 9.921147e-01 }, { 1.570800e-01, 1.564348e-01, 9.876883e-01 }, { 1.884960e-01, 1.873817e-01, 9.822872e-01 }, { 2.199120e-01, 2.181437e-01, 9.759166e-01 }, { 2.513280e-01, 2.486905e-01, 9.685830e-01 }, { 2.827440e-01, 2.789917e-01, 9.602935e-01 }, { 3.141600e-01, 3.090177e-01, 9.510563e-01 }, { 3.455760e-01, 3.387387e-01, 9.408805e-01 }, { 3.769920e-01, 3.681254e-01, 9.297762e-01 }, { 4.084080e-01, 3.971488e-01, 9.177542e-01 }, { 4.398240e-01, 4.257802e-01, 9.048266e-01 }, { 4.712400e-01, 4.539915e-01, 8.910060e-01 }, { 5.026560e-01, 4.817547e-01, 8.763061e-01 }, { 5.340720e-01, 5.090425e-01, 8.607414e-01 }, { 5.654880e-01, 5.358279e-01, 8.443272e-01 }, { 5.969040e-01, 5.620845e-01, 8.270798e-01 }, { 6.283200e-01, 5.877864e-01, 8.090161e-01 }, { 6.597360e-01, 6.129083e-01, 7.901541e-01 }, { 6.911520e-01, 6.374252e-01, 7.705122e-01 }, { 7.225680e-01, 6.613131e-01, 7.501100e-01 }, { 7.539840e-01, 6.845484e-01, 7.289674e-01 }, { 7.854000e-01, 7.071081e-01, 7.071055e-01 }, { 8.168160e-01, 7.289699e-01, 6.845457e-01 }, { 8.482320e-01, 7.501124e-01, 6.613104e-01 }, { 8.796480e-01, 7.705146e-01, 6.374224e-01 }, { 9.110640e-01, 7.901563e-01, 6.129054e-01 }, { 9.424800e-01, 8.090183e-01, 5.877835e-01 }, { 9.738960e-01, 8.270819e-01, 5.620815e-01 }, { 1.005312e+00, 8.443292e-01, 5.358248e-01 }, { 1.036728e+00, 8.607433e-01, 5.090393e-01 }, { 1.068144e+00, 8.763079e-01, 4.817515e-01 }, { 1.099560e+00, 8.910077e-01, 4.539882e-01 }, { 1.130976e+00, 9.048282e-01, 4.257769e-01 }, { 1.162392e+00, 9.177557e-01, 3.971454e-01 }, { 1.193808e+00, 9.297775e-01, 3.681220e-01 }, { 1.225224e+00, 9.408817e-01, 3.387352e-01 }, { 1.256640e+00, 9.510574e-01, 3.090142e-01 }, { 1.288056e+00, 9.602945e-01, 2.789882e-01 }, { 1.319472e+00, 9.685839e-01, 2.486869e-01 }, { 1.350888e+00, 9.759175e-01, 2.181402e-01 }, { 1.382304e+00, 9.822879e-01, 1.873781e-01 }, { 1.413720e+00, 9.876889e-01, 1.564312e-01 }, { 1.445136e+00, 9.921151e-01, 1.253299e-01 }, { 1.476552e+00, 9.955623e-01, 9.410488e-02 }, { 1.507968e+00, 9.980269e-01, 6.278700e-02 }, { 1.539384e+00, 9.995067e-01, 3.140716e-02 }, { 1.570800e+00, 1.000000e-00, -3.673205e-06 }, { 1.602216e+00, 9.995064e-01, -3.141450e-02 }, { 1.633632e+00, 9.980265e-01, -6.279433e-02 }, { 1.665048e+00, 9.955616e-01, -9.411219e-02 }, { 1.696464e+00, 9.921142e-01, -1.253372e-01 }, { 1.727880e+00, 9.876877e-01, -1.564385e-01 }, { 1.759296e+00, 9.822865e-01, -1.873854e-01 }, { 1.790712e+00, 9.759158e-01, -2.181473e-01 }, { 1.822128e+00, 9.685821e-01, -2.486940e-01 }, { 1.853544e+00, 9.602925e-01, -2.789953e-01 }, { 1.884960e+00, 9.510552e-01, -3.090212e-01 }, { 1.916376e+00, 9.408793e-01, -3.387421e-01 }, { 1.947792e+00, 9.297748e-01, -3.681288e-01 }, { 1.979208e+00, 9.177528e-01, -3.971521e-01 }, { 2.010624e+00, 9.048251e-01, -4.257835e-01 }, { 2.042040e+00, 8.910044e-01, -4.539948e-01 }, { 2.073456e+00, 8.763043e-01, -4.817579e-01 }, { 2.104872e+00, 8.607395e-01, -5.090457e-01 }, { 2.136288e+00, 8.443252e-01, -5.358310e-01 }, { 2.167704e+00, 8.270777e-01, -5.620876e-01 }, { 2.199120e+00, 8.090140e-01, -5.877894e-01 }, { 2.230536e+00, 7.901518e-01, -6.129112e-01 }, { 2.261952e+00, 7.705099e-01, -6.374281e-01 }, { 2.293368e+00, 7.501075e-01, -6.613159e-01 }, { 2.324784e+00, 7.289649e-01, -6.845511e-01 }, { 2.356200e+00, 7.071029e-01, -7.071107e-01 }, { 2.387616e+00, 6.845430e-01, -7.289724e-01 }, { 2.419032e+00, 6.613076e-01, -7.501148e-01 }, { 2.450448e+00, 6.374196e-01, -7.705169e-01 }, { 2.481864e+00, 6.129025e-01, -7.901586e-01 }, { 2.513280e+00, 5.877805e-01, -8.090204e-01 }, { 2.544696e+00, 5.620785e-01, -8.270839e-01 }, { 2.576112e+00, 5.358217e-01, -8.443312e-01 }, { 2.607528e+00, 5.090362e-01, -8.607451e-01 }, { 2.638944e+00, 4.817483e-01, -8.763097e-01 }, { 2.670360e+00, 4.539849e-01, -8.910094e-01 }, { 2.701776e+00, 4.257736e-01, -9.048297e-01 }, { 2.733192e+00, 3.971420e-01, -9.177572e-01 }, { 2.764608e+00, 3.681185e-01, -9.297789e-01 }, { 2.796024e+00, 3.387318e-01, -9.408830e-01 }, { 2.827440e+00, 3.090107e-01, -9.510586e-01 }, { 2.858856e+00, 2.789847e-01, -9.602956e-01 }, { 2.890272e+00, 2.486833e-01, -9.685848e-01 }, { 2.921688e+00, 2.181366e-01, -9.759183e-01 }, { 2.953104e+00, 1.873745e-01, -9.822885e-01 }, { 2.984520e+00, 1.564276e-01, -9.876894e-01 }, { 3.015936e+00, 1.253262e-01, -9.921156e-01 }, { 3.047352e+00, 9.410122e-02, -9.955626e-01 }, { 3.078768e+00, 6.278333e-02, -9.980272e-01 }, { 3.110184e+00, 3.140349e-02, -9.995068e-01 }, { 3.141600e+00, -7.346410e-06, -1.000000e-00 } }; # define M_PI 3.14159265358979323846 /* pi */ static double inline gen_sin(double x) { int i; if(x < 0) while(x < -M_PI) x+= M_PI; else while(x > M_PI) x-= M_PI; for(i=0;i<sizeof(g_sincos)/sizeof(gen_sincos_t)-1;i++) { if(x>=g_sincos[i].x && x <= g_sincos[i+1].x) { return (g_sincos[i+1].sinx-g_sincos[i].sinx)*(x-g_sincos[i].x)/(g_sincos[i+1].x-g_sincos[i].x)+g_sincos[i].sinx; } } return x<0?1:-1; } #undef sin #define sin(x) gen_sin(x) static double inline gen_cos(double x) { int i; if(x < 0) while(x < -M_PI) x+= M_PI; else while(x > M_PI) x-= M_PI; for(i=0;i<sizeof(g_sincos)/sizeof(gen_sincos_t)-1;i++) { if(x>=g_sincos[i].x && x <= g_sincos[i+1].x) { return (g_sincos[i+1].cosx-g_sincos[i].cosx)*(x-g_sincos[i].x)/(g_sincos[i+1].x-g_sincos[i].x)+g_sincos[i].cosx; } } return x<0?1:-1; } #undef cos #define cos(x) gen_cos(x) #endif