
(FPCore (x) :precision binary64 (* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))
double code(double x) {
return 2.0 * atan(sqrt(((1.0 - x) / (1.0 + x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(sqrt(((1.0d0 - x) / (1.0d0 + x))))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.sqrt(((1.0 - x) / (1.0 + x))));
}
def code(x): return 2.0 * math.atan(math.sqrt(((1.0 - x) / (1.0 + x))))
function code(x) return Float64(2.0 * atan(sqrt(Float64(Float64(1.0 - x) / Float64(1.0 + x))))) end
function tmp = code(x) tmp = 2.0 * atan(sqrt(((1.0 - x) / (1.0 + x)))); end
code[x_] := N[(2.0 * N[ArcTan[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\sqrt{\frac{1 - x}{1 + x}}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))
double code(double x) {
return 2.0 * atan(sqrt(((1.0 - x) / (1.0 + x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(sqrt(((1.0d0 - x) / (1.0d0 + x))))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.sqrt(((1.0 - x) / (1.0 + x))));
}
def code(x): return 2.0 * math.atan(math.sqrt(((1.0 - x) / (1.0 + x))))
function code(x) return Float64(2.0 * atan(sqrt(Float64(Float64(1.0 - x) / Float64(1.0 + x))))) end
function tmp = code(x) tmp = 2.0 * atan(sqrt(((1.0 - x) / (1.0 + x)))); end
code[x_] := N[(2.0 * N[ArcTan[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\sqrt{\frac{1 - x}{1 + x}}\right)
\end{array}
(FPCore (x) :precision binary64 (* 2.0 (atan (pow (/ (+ 1.0 x) (- 1.0 x)) -0.5))))
double code(double x) {
return 2.0 * atan(pow(((1.0 + x) / (1.0 - x)), -0.5));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((((1.0d0 + x) / (1.0d0 - x)) ** (-0.5d0)))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.pow(((1.0 + x) / (1.0 - x)), -0.5));
}
def code(x): return 2.0 * math.atan(math.pow(((1.0 + x) / (1.0 - x)), -0.5))
function code(x) return Float64(2.0 * atan((Float64(Float64(1.0 + x) / Float64(1.0 - x)) ^ -0.5))) end
function tmp = code(x) tmp = 2.0 * atan((((1.0 + x) / (1.0 - x)) ^ -0.5)); end
code[x_] := N[(2.0 * N[ArcTan[N[Power[N[(N[(1.0 + x), $MachinePrecision] / N[(1.0 - x), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left({\left(\frac{1 + x}{1 - x}\right)}^{-0.5}\right)
\end{array}
Initial program 100.0%
pow1/2100.0%
clear-num100.0%
inv-pow100.0%
pow-pow100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))
double code(double x) {
return 2.0 * atan(sqrt(((1.0 - x) / (1.0 + x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(sqrt(((1.0d0 - x) / (1.0d0 + x))))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.sqrt(((1.0 - x) / (1.0 + x))));
}
def code(x): return 2.0 * math.atan(math.sqrt(((1.0 - x) / (1.0 + x))))
function code(x) return Float64(2.0 * atan(sqrt(Float64(Float64(1.0 - x) / Float64(1.0 + x))))) end
function tmp = code(x) tmp = 2.0 * atan(sqrt(((1.0 - x) / (1.0 + x)))); end
code[x_] := N[(2.0 * N[ArcTan[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\sqrt{\frac{1 - x}{1 + x}}\right)
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (let* ((t_0 (+ x (* x (* x -0.5))))) (* 2.0 (atan (/ (- 1.0 (* t_0 t_0)) (+ 1.0 t_0))))))
double code(double x) {
double t_0 = x + (x * (x * -0.5));
return 2.0 * atan(((1.0 - (t_0 * t_0)) / (1.0 + t_0)));
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
t_0 = x + (x * (x * (-0.5d0)))
code = 2.0d0 * atan(((1.0d0 - (t_0 * t_0)) / (1.0d0 + t_0)))
end function
public static double code(double x) {
double t_0 = x + (x * (x * -0.5));
return 2.0 * Math.atan(((1.0 - (t_0 * t_0)) / (1.0 + t_0)));
}
def code(x): t_0 = x + (x * (x * -0.5)) return 2.0 * math.atan(((1.0 - (t_0 * t_0)) / (1.0 + t_0)))
function code(x) t_0 = Float64(x + Float64(x * Float64(x * -0.5))) return Float64(2.0 * atan(Float64(Float64(1.0 - Float64(t_0 * t_0)) / Float64(1.0 + t_0)))) end
function tmp = code(x) t_0 = x + (x * (x * -0.5)); tmp = 2.0 * atan(((1.0 - (t_0 * t_0)) / (1.0 + t_0))); end
code[x_] := Block[{t$95$0 = N[(x + N[(x * N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(2.0 * N[ArcTan[N[(N[(1.0 - N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + x \cdot \left(x \cdot -0.5\right)\\
2 \cdot \tan^{-1} \left(\frac{1 - t_0 \cdot t_0}{1 + t_0}\right)
\end{array}
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.9%
associate-+r+98.9%
neg-mul-198.9%
sub-neg98.9%
*-commutative98.9%
unpow298.9%
Simplified98.9%
associate-+l-98.9%
flip--98.9%
metadata-eval98.9%
*-commutative98.9%
cancel-sign-sub-inv98.9%
metadata-eval98.9%
*-commutative98.9%
cancel-sign-sub-inv98.9%
metadata-eval98.9%
*-commutative98.9%
cancel-sign-sub-inv98.9%
metadata-eval98.9%
Applied egg-rr98.9%
associate-*r*98.9%
associate-*r*98.9%
associate-*r*98.9%
Simplified98.9%
Final simplification98.9%
(FPCore (x) :precision binary64 (* 2.0 (atan (/ (- 1.0 (* x x)) (+ 1.0 (+ x (* x (* x -0.5))))))))
double code(double x) {
return 2.0 * atan(((1.0 - (x * x)) / (1.0 + (x + (x * (x * -0.5))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(((1.0d0 - (x * x)) / (1.0d0 + (x + (x * (x * (-0.5d0)))))))
end function
public static double code(double x) {
return 2.0 * Math.atan(((1.0 - (x * x)) / (1.0 + (x + (x * (x * -0.5))))));
}
def code(x): return 2.0 * math.atan(((1.0 - (x * x)) / (1.0 + (x + (x * (x * -0.5))))))
function code(x) return Float64(2.0 * atan(Float64(Float64(1.0 - Float64(x * x)) / Float64(1.0 + Float64(x + Float64(x * Float64(x * -0.5))))))) end
function tmp = code(x) tmp = 2.0 * atan(((1.0 - (x * x)) / (1.0 + (x + (x * (x * -0.5)))))); end
code[x_] := N[(2.0 * N[ArcTan[N[(N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(x + N[(x * N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{1 - x \cdot x}{1 + \left(x + x \cdot \left(x \cdot -0.5\right)\right)}\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.9%
associate-+r+98.9%
neg-mul-198.9%
sub-neg98.9%
*-commutative98.9%
unpow298.9%
Simplified98.9%
associate-+l-98.9%
flip--98.9%
metadata-eval98.9%
*-commutative98.9%
cancel-sign-sub-inv98.9%
metadata-eval98.9%
*-commutative98.9%
cancel-sign-sub-inv98.9%
metadata-eval98.9%
*-commutative98.9%
cancel-sign-sub-inv98.9%
metadata-eval98.9%
Applied egg-rr98.9%
associate-*r*98.9%
associate-*r*98.9%
associate-*r*98.9%
Simplified98.9%
Taylor expanded in x around 0 98.9%
unpow298.9%
Simplified98.9%
Final simplification98.9%
(FPCore (x) :precision binary64 (* 2.0 (atan (+ 1.0 (- (* x (* x 0.5)) x)))))
double code(double x) {
return 2.0 * atan((1.0 + ((x * (x * 0.5)) - x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 + ((x * (x * 0.5d0)) - x)))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 + ((x * (x * 0.5)) - x)));
}
def code(x): return 2.0 * math.atan((1.0 + ((x * (x * 0.5)) - x)))
function code(x) return Float64(2.0 * atan(Float64(1.0 + Float64(Float64(x * Float64(x * 0.5)) - x)))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 + ((x * (x * 0.5)) - x))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 + N[(N[(x * N[(x * 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(1 + \left(x \cdot \left(x \cdot 0.5\right) - x\right)\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.9%
+-commutative98.9%
neg-mul-198.9%
unsub-neg98.9%
*-commutative98.9%
unpow298.9%
associate-*l*98.9%
Simplified98.9%
Final simplification98.9%
(FPCore (x) :precision binary64 (* 2.0 (atan (+ (- 1.0 x) (* (* x x) 0.5)))))
double code(double x) {
return 2.0 * atan(((1.0 - x) + ((x * x) * 0.5)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(((1.0d0 - x) + ((x * x) * 0.5d0)))
end function
public static double code(double x) {
return 2.0 * Math.atan(((1.0 - x) + ((x * x) * 0.5)));
}
def code(x): return 2.0 * math.atan(((1.0 - x) + ((x * x) * 0.5)))
function code(x) return Float64(2.0 * atan(Float64(Float64(1.0 - x) + Float64(Float64(x * x) * 0.5)))) end
function tmp = code(x) tmp = 2.0 * atan(((1.0 - x) + ((x * x) * 0.5))); end
code[x_] := N[(2.0 * N[ArcTan[N[(N[(1.0 - x), $MachinePrecision] + N[(N[(x * x), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\left(1 - x\right) + \left(x \cdot x\right) \cdot 0.5\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.9%
associate-+r+98.9%
neg-mul-198.9%
sub-neg98.9%
*-commutative98.9%
unpow298.9%
Simplified98.9%
Final simplification98.9%
(FPCore (x) :precision binary64 (* 2.0 (atan (- 1.0 x))))
double code(double x) {
return 2.0 * atan((1.0 - x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 - x))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 - x));
}
def code(x): return 2.0 * math.atan((1.0 - x))
function code(x) return Float64(2.0 * atan(Float64(1.0 - x))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 - x)); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(1 - x\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.5%
neg-mul-198.5%
sub-neg98.5%
Simplified98.5%
Final simplification98.5%
(FPCore (x) :precision binary64 (* 2.0 (atan 1.0)))
double code(double x) {
return 2.0 * atan(1.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(1.0d0)
end function
public static double code(double x) {
return 2.0 * Math.atan(1.0);
}
def code(x): return 2.0 * math.atan(1.0)
function code(x) return Float64(2.0 * atan(1.0)) end
function tmp = code(x) tmp = 2.0 * atan(1.0); end
code[x_] := N[(2.0 * N[ArcTan[1.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} 1
\end{array}
Initial program 100.0%
div-inv100.0%
sqrt-prod100.0%
pow1/2100.0%
inv-pow100.0%
pow-pow100.0%
metadata-eval100.0%
Applied egg-rr100.0%
*-commutative100.0%
Simplified100.0%
Taylor expanded in x around 0 97.4%
Final simplification97.4%
herbie shell --seed 2023297
(FPCore (x)
:name "arccos"
:precision binary64
(* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))