22 static double cos_sign_tbl[] = {1,-1,-1,1};
25 static double sin_sign_tbl[] = {1,-1,-1,1};
39     quadrant = (quadrant - 1) & 0x3;  in sin()
49      * result = 1 - x^2/2! + x^4/4! - x^6/6! + x^8/8!  in sin()
52      * result = 1 - x^2 * (1/2! - x^2 * (1/4! - x^2 * (1/6! - x^2 * (1/8!))))  in sin()
59     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18*19*20);  in sin()
63     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18);  in sin()
67     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16);  in sin()
71     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14);  in sin()
75     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12);  in sin()
79     result += 1./(1.*2*3*4*5*6*7*8*9*10);  in sin()
82     result += 1./(1.*2*3*4*5*6*7*8);  in sin()
85     result += 1./(1.*2*3*4*5*6);  in sin()
88     result += 1./(1.*2*3*4);  in sin()
91     result += 1./(1.*2);  in sin()
94     result += 1;  in sin()
124      * result = 1 - x^2/2! + x^4/4! - x^6/6! + x^8/8!  in cos()
127      * result = 1 - x^2 * (1/2! - x^2 * (1/4! - x^2 * (1/6! - x^2 * (1/8!))))  in cos()
134     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18*19*20);  in cos()
138     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18);  in cos()
142     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16);  in cos()
146     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12*13*14);  in cos()
150     result += 1./(1.*2*3*4*5*6*7*8*9*10*11*12);  in cos()
154     result += 1./(1.*2*3*4*5*6*7*8*9*10);  in cos()
157     result += 1./(1.*2*3*4*5*6*7*8);  in cos()
160     result += 1./(1.*2*3*4*5*6);  in cos()
163     result += 1./(1.*2*3*4);  in cos()
166     result += 1./(1.*2);  in cos()
169     result += 1;  in cos()