快速sin cos

在优化算法时时间项费非常关键，如果需要大量使用sin cos，那么对sin cos的优化将大大提高运算速度

最近发现一个快速计算sin cos的方法，

#include"cmath"
#include<iostream>
#include<cstdio>
#include<ctime>
using namespace std;
#define MILLION 1000000


#define CV_PI 3.14159265359
#define halfPi ((float)(CV_PI*0.5))
#define Pi     ((float)CV_PI)
#define a0  0 /*-4.172325e-7f*/   /*(-(float)0x7)/((float)0x1000000); */
#define a1 1.000025f        /*((float)0x1922253)/((float)0x1000000)*2/Pi; */
#define a2 -2.652905e-4f    /*(-(float)0x2ae6)/((float)0x1000000)*4/(Pi*Pi); */
#define a3 -0.165624f       /*(-(float)0xa45511)/((float)0x1000000)*8/(Pi*Pi*Pi); */
#define a4 -1.964532e-3f    /*(-(float)0x30fd3)/((float)0x1000000)*16/(Pi*Pi*Pi*Pi); */
#define a5 1.02575e-2f      /*((float)0x191cac)/((float)0x1000000)*32/(Pi*Pi*Pi*Pi*Pi); */
#define a6 -9.580378e-4f    /*(-(float)0x3af27)/((float)0x1000000)*64/(Pi*Pi*Pi*Pi*Pi*Pi); */

#define _sin(x) ((((((a6*(x) + a5)*(x) + a4)*(x) + a3)*(x) + a2)*(x) + a1)*(x) + a0)
#define _cos(x) _sin(halfPi - (x))

int main(){
    struct timespec tpstart;
    struct timespec tpend;
    long timedif;


	//fastsin cos
	{
		
		float sine, cose;
		float dsin, dcos;
		clock_gettime(CLOCK_MONOTONIC, &tpstart);
		for(float ang = 0.0; ang < CV_PI; ang += 0.0001)
		{
			sine = _sin(ang);
			cose = _cos(ang);		
		}
		clock_gettime(CLOCK_MONOTONIC, &tpend);
        timedif = MILLION*(tpend.tv_sec-tpstart.tv_sec)+(tpend.tv_nsec-tpstart.tv_nsec)/1000;
        printf("it took %ld microseconds\n", timedif);
	}	
	//sin cos
	{
		double sine,cose;
		clock_gettime(CLOCK_MONOTONIC, &tpstart);
		for(double ang = 0.0; ang < CV_PI; ang += 0.0001)
		{
			sine = sin(ang);
			cose = cos(ang);
		}
		clock_gettime(CLOCK_MONOTONIC, &tpend);
        timedif = MILLION*(tpend.tv_sec-tpstart.tv_sec)+(tpend.tv_nsec-tpstart.tv_nsec)/1000;
        printf("it took %ld microseconds\n", timedif);
	}	
	return 0;
}

Makefile文件

fastsin : fastsin.cpp
	g++ -o fastsin fastsin.cpp -lrt
clean:
	rm -f fastsin

计算结果：

it took 1300 microseconds
it took 6297 microseconds

系统自带的sin cos计算（double），大概是快速sin cos计算（float）时间的5倍。