arch/generic/s_sin.c

74 static double sloww(double x, double dx, double orig);
75 static double sloww1(double x, double dx, double orig);
76 static double sloww2(double x, double dx, double orig, int n);
77 static double bsloww(double x, double dx, double orig, int n);
78 static double bsloww1(double x, double dx, double orig, int n);
79 static double bsloww2(double x, double dx, double orig, int n);
82 static double csloww(double x, double dx, double orig);
83 static double csloww1(double x, double dx, double orig);
84 static double csloww2(double x, double dx, double orig, int n);
697 static double sloww(double x,double dx, double orig) {
711   cor = (cor>0)? 1.0005*cor+ABS(orig)*3.1e-30 : 1.0005*cor-ABS(orig)*3.1e-30;
715     cor = (w[1]>0)? 1.000000001*w[1] + ABS(orig)*1.1e-30 : 1.000000001*w[1] - ABS(orig)*1.1e-30;
718       t = (orig*hpinv.x + toint.x);
721       y = (orig - xn*mp1.x) - xn*mp2.x;
731       cor = (w[1]>0)? 1.000000001*w[1] + ABS(orig)*1.1e-40 : 1.000000001*w[1] - ABS(orig)*1.1e-40;
733       else return __mpsin1(orig);
744 static double sloww1(double x, double dx, double orig) {
770   cor = (cor>0)? 1.0005*cor+3.1e-30*ABS(orig) : 1.0005*cor-3.1e-30*ABS(orig);
774     cor = (w[1]>0)? 1.000000005*w[1]+1.1e-30*ABS(orig) : 1.000000005*w[1]-1.1e-30*ABS(orig);
776   else  return __mpsin1(orig);
786 static double sloww2(double x, double dx, double orig, int n) {
813   cor = (cor>0)? 1.0005*cor+3.1e-30*ABS(orig) : 1.0005*cor-3.1e-30*ABS(orig);
817    cor = (w[1]>0)? 1.000000005*w[1]+1.1e-30*ABS(orig) : 1.000000005*w[1]-1.1e-30*ABS(orig);
819    else  return __mpsin1(orig);
830 static double bsloww(double x,double dx, double orig,int n) {
852     else return (n&1)?__mpcos1(orig):__mpsin1(orig);
863 static double bsloww1(double x, double dx, double orig,int n) {
895    else  return (n&1)?__mpcos1(orig):__mpsin1(orig);
906 static double bsloww2(double x, double dx, double orig, int n) {
939    else  return (n&1)?__mpsin1(orig):__mpcos1(orig);
991 static double csloww(double x,double dx, double orig) {
1006   cor = (cor>0)? 1.0005*cor+ABS(orig)*3.1e-30 : 1.0005*cor-ABS(orig)*3.1e-30;
1010     cor = (w[1]>0)? 1.000000001*w[1] + ABS(orig)*1.1e-30 : 1.000000001*w[1] - ABS(orig)*1.1e-30;
1013       t = (orig*hpinv.x + toint.x);
1016       y = (orig - xn*mp1.x) - xn*mp2.x;
1026       cor = (w[1]>0)? 1.000000001*w[1] + ABS(orig)*1.1e-40 : 1.000000001*w[1] - ABS(orig)*1.1e-40;
1028       else return __mpcos1(orig);
1040 static double csloww1(double x, double dx, double orig) {
1066   cor = (cor>0)? 1.0005*cor+3.1e-30*ABS(orig) : 1.0005*cor-3.1e-30*ABS(orig);
1070     cor = (w[1]>0)? 1.000000005*w[1]+1.1e-30*ABS(orig) : 1.000000005*w[1]-1.1e-30*ABS(orig);
1072     else  return __mpcos1(orig);
1084 static double csloww2(double x, double dx, double orig, int n) {
1111   cor = (cor>0)? 1.0005*cor+3.1e-30*ABS(orig) : 1.0005*cor-3.1e-30*ABS(orig);
1115     cor = (w[1]>0)? 1.000000005*w[1]+1.1e-30*ABS(orig) : 1.000000005*w[1]-1.1e-30*ABS(orig);
1117     else  return __mpcos1(orig);