47 const double oneIfXXZero = (xx==0.);
48 double t=yy/(xx+oneIfXXZero);
53 t = (t-1.0) / (t+1.0);
66 const double t2 = t * t;
73 const double poq=px / qx;
75 double ret = t * t2 * poq + t;
81 ret*= (1.-oneIfXXZero);
106 const float oneIfXXZero = (xx==0.f);
110 if( t > 0.4142135623730950
f )
111 z = (t-1.0f)/(t+1.0
f);
116 float ret =(((( 8.05374449538e-2
f * z2
117 - 1.38776856032E-1
f) * z2
118 + 1.99777106478
E-1
f) * z2
119 - 3.33329491539E-1
f) * z2 * z
123 ret*= (1.f - oneIfXXZero);
130 if (y<0.
f) ret = -ret;
double get_atan_px(const double x2)
VecExpr< UnaryOp< Fabs< T >, VecExpr< A, T, D >, T >, T, D > fabs(const VecExpr< A, T, D > &rhs)
you should not use this method at all Int_t Int_t z
double get_atan_qx(const double x2)
double fast_atan2(double y, double x)
float fast_atan2f(float y, float x)