
(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 10 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 (/ (sqrt (- 1.0 x)) (sqrt (+ 1.0 x))))))
double code(double x) {
return 2.0 * atan((sqrt((1.0 - x)) / sqrt((1.0 + x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((sqrt((1.0d0 - x)) / sqrt((1.0d0 + x))))
end function
public static double code(double x) {
return 2.0 * Math.atan((Math.sqrt((1.0 - x)) / Math.sqrt((1.0 + x))));
}
def code(x): return 2.0 * math.atan((math.sqrt((1.0 - x)) / math.sqrt((1.0 + x))))
function code(x) return Float64(2.0 * atan(Float64(sqrt(Float64(1.0 - x)) / sqrt(Float64(1.0 + x))))) end
function tmp = code(x) tmp = 2.0 * atan((sqrt((1.0 - x)) / sqrt((1.0 + x)))); end
code[x_] := N[(2.0 * N[ArcTan[N[(N[Sqrt[N[(1.0 - x), $MachinePrecision]], $MachinePrecision] / N[Sqrt[N[(1.0 + x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{\sqrt{1 - x}}{\sqrt{1 + x}}\right)
\end{array}
Initial program 100.0%
sqrt-div100.0%
div-inv100.0%
add-sqr-sqrt50.4%
hypot-1-def50.4%
Applied egg-rr50.4%
associate-*r/50.4%
*-rgt-identity50.4%
hypot-undefine50.4%
metadata-eval50.4%
rem-square-sqrt100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* 2.0 (atan (/ 1.0 (sqrt (/ (+ 1.0 x) (- 1.0 x)))))))
double code(double x) {
return 2.0 * atan((1.0 / sqrt(((1.0 + x) / (1.0 - x)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 / sqrt(((1.0d0 + x) / (1.0d0 - x)))))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 / Math.sqrt(((1.0 + x) / (1.0 - x)))));
}
def code(x): return 2.0 * math.atan((1.0 / math.sqrt(((1.0 + x) / (1.0 - x)))))
function code(x) return Float64(2.0 * atan(Float64(1.0 / sqrt(Float64(Float64(1.0 + x) / Float64(1.0 - x)))))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 / sqrt(((1.0 + x) / (1.0 - x))))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 / N[Sqrt[N[(N[(1.0 + x), $MachinePrecision] / N[(1.0 - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{1}{\sqrt{\frac{1 + x}{1 - x}}}\right)
\end{array}
Initial program 100.0%
clear-num100.0%
sqrt-div100.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 (* 2.0 (atan (/ 1.0 (+ 1.0 (* x (+ 1.0 (* x (+ 0.5 (* x 0.5))))))))))
double code(double x) {
return 2.0 * atan((1.0 / (1.0 + (x * (1.0 + (x * (0.5 + (x * 0.5))))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 / (1.0d0 + (x * (1.0d0 + (x * (0.5d0 + (x * 0.5d0))))))))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 / (1.0 + (x * (1.0 + (x * (0.5 + (x * 0.5))))))));
}
def code(x): return 2.0 * math.atan((1.0 / (1.0 + (x * (1.0 + (x * (0.5 + (x * 0.5))))))))
function code(x) return Float64(2.0 * atan(Float64(1.0 / Float64(1.0 + Float64(x * Float64(1.0 + Float64(x * Float64(0.5 + Float64(x * 0.5))))))))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 / (1.0 + (x * (1.0 + (x * (0.5 + (x * 0.5)))))))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 / N[(1.0 + N[(x * N[(1.0 + N[(x * N[(0.5 + N[(x * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{1}{1 + x \cdot \left(1 + x \cdot \left(0.5 + x \cdot 0.5\right)\right)}\right)
\end{array}
Initial program 100.0%
clear-num100.0%
sqrt-div100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0 99.5%
*-commutative99.5%
Simplified99.5%
Final simplification99.5%
(FPCore (x) :precision binary64 (* 2.0 (atan (+ 1.0 (* x (+ (* x (+ 0.5 (* x -0.5))) -1.0))))))
double code(double x) {
return 2.0 * atan((1.0 + (x * ((x * (0.5 + (x * -0.5))) + -1.0))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 + (x * ((x * (0.5d0 + (x * (-0.5d0)))) + (-1.0d0)))))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 + (x * ((x * (0.5 + (x * -0.5))) + -1.0))));
}
def code(x): return 2.0 * math.atan((1.0 + (x * ((x * (0.5 + (x * -0.5))) + -1.0))))
function code(x) return Float64(2.0 * atan(Float64(1.0 + Float64(x * Float64(Float64(x * Float64(0.5 + Float64(x * -0.5))) + -1.0))))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 + (x * ((x * (0.5 + (x * -0.5))) + -1.0)))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 + N[(x * N[(N[(x * N[(0.5 + N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(1 + x \cdot \left(x \cdot \left(0.5 + x \cdot -0.5\right) + -1\right)\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 99.5%
Final simplification99.5%
(FPCore (x) :precision binary64 (* 2.0 (atan (/ 1.0 (+ 1.0 (* x (+ 1.0 (* x 0.5))))))))
double code(double x) {
return 2.0 * atan((1.0 / (1.0 + (x * (1.0 + (x * 0.5))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 / (1.0d0 + (x * (1.0d0 + (x * 0.5d0))))))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 / (1.0 + (x * (1.0 + (x * 0.5))))));
}
def code(x): return 2.0 * math.atan((1.0 / (1.0 + (x * (1.0 + (x * 0.5))))))
function code(x) return Float64(2.0 * atan(Float64(1.0 / Float64(1.0 + Float64(x * Float64(1.0 + Float64(x * 0.5))))))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 / (1.0 + (x * (1.0 + (x * 0.5)))))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 / N[(1.0 + N[(x * N[(1.0 + N[(x * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{1}{1 + x \cdot \left(1 + x \cdot 0.5\right)}\right)
\end{array}
Initial program 100.0%
clear-num100.0%
sqrt-div100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0 99.1%
Final simplification99.1%
(FPCore (x) :precision binary64 (* 2.0 (atan (+ 1.0 (* x (+ (* x 0.5) -1.0))))))
double code(double x) {
return 2.0 * atan((1.0 + (x * ((x * 0.5) + -1.0))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 + (x * ((x * 0.5d0) + (-1.0d0)))))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 + (x * ((x * 0.5) + -1.0))));
}
def code(x): return 2.0 * math.atan((1.0 + (x * ((x * 0.5) + -1.0))))
function code(x) return Float64(2.0 * atan(Float64(1.0 + Float64(x * Float64(Float64(x * 0.5) + -1.0))))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 + (x * ((x * 0.5) + -1.0)))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 + N[(x * N[(N[(x * 0.5), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(1 + x \cdot \left(x \cdot 0.5 + -1\right)\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 99.1%
Final simplification99.1%
(FPCore (x) :precision binary64 (* 2.0 (atan (/ 1.0 (+ 1.0 x)))))
double code(double x) {
return 2.0 * atan((1.0 / (1.0 + x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 / (1.0d0 + x)))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 / (1.0 + x)));
}
def code(x): return 2.0 * math.atan((1.0 / (1.0 + x)))
function code(x) return Float64(2.0 * atan(Float64(1.0 / Float64(1.0 + x)))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 / (1.0 + x))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{1}{1 + x}\right)
\end{array}
Initial program 100.0%
clear-num100.0%
sqrt-div100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0 98.6%
Final simplification98.6%
(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.6%
neg-mul-198.6%
sub-neg98.6%
Simplified98.6%
Final simplification98.6%
(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%
sqrt-div100.0%
div-inv100.0%
add-sqr-sqrt50.4%
hypot-1-def50.4%
Applied egg-rr50.4%
associate-*r/50.4%
*-rgt-identity50.4%
hypot-undefine50.4%
metadata-eval50.4%
rem-square-sqrt100.0%
Simplified100.0%
Taylor expanded in x around 0 97.3%
Final simplification97.3%
herbie shell --seed 2024130
(FPCore (x)
:name "arccos"
:precision binary64
(* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))