
(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 (cbrt (pow (/ (+ 1.0 x) (- 1.0 x)) -1.5)))))
double code(double x) {
return 2.0 * atan(cbrt(pow(((1.0 + x) / (1.0 - x)), -1.5)));
}
public static double code(double x) {
return 2.0 * Math.atan(Math.cbrt(Math.pow(((1.0 + x) / (1.0 - x)), -1.5)));
}
function code(x) return Float64(2.0 * atan(cbrt((Float64(Float64(1.0 + x) / Float64(1.0 - x)) ^ -1.5)))) end
code[x_] := N[(2.0 * N[ArcTan[N[Power[N[Power[N[(N[(1.0 + x), $MachinePrecision] / N[(1.0 - x), $MachinePrecision]), $MachinePrecision], -1.5], $MachinePrecision], 1/3], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\sqrt[3]{{\left(\frac{1 + x}{1 - x}\right)}^{-1.5}}\right)
\end{array}
Initial program 99.9%
pow1/299.9%
metadata-eval99.9%
metadata-eval99.9%
pow-pow100.0%
clear-num99.9%
inv-pow99.9%
pow-pow99.9%
metadata-eval99.9%
metadata-eval99.9%
Applied egg-rr99.9%
unpow1/3100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* 2.0 (atan (/ (sqrt (- 1.0 (* x x))) (+ 1.0 x)))))
double code(double x) {
return 2.0 * atan((sqrt((1.0 - (x * x))) / (1.0 + x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((sqrt((1.0d0 - (x * x))) / (1.0d0 + x)))
end function
public static double code(double x) {
return 2.0 * Math.atan((Math.sqrt((1.0 - (x * x))) / (1.0 + x)));
}
def code(x): return 2.0 * math.atan((math.sqrt((1.0 - (x * x))) / (1.0 + x)))
function code(x) return Float64(2.0 * atan(Float64(sqrt(Float64(1.0 - Float64(x * x))) / Float64(1.0 + x)))) end
function tmp = code(x) tmp = 2.0 * atan((sqrt((1.0 - (x * x))) / (1.0 + x))); end
code[x_] := N[(2.0 * N[ArcTan[N[(N[Sqrt[N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{\sqrt{1 - x \cdot x}}{1 + x}\right)
\end{array}
Initial program 99.9%
clear-num99.9%
sqrt-div99.9%
metadata-eval99.9%
Applied egg-rr99.9%
flip--100.0%
associate-/r/100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-*l/100.0%
associate-/l*100.0%
Simplified100.0%
add-sqr-sqrt99.9%
sqrt-unprod100.0%
frac-times100.0%
metadata-eval100.0%
add-sqr-sqrt100.0%
Applied egg-rr100.0%
associate-/l*100.0%
*-lft-identity100.0%
associate-/l/99.9%
Simplified99.9%
Taylor expanded in x around 0 100.0%
associate-*l/100.0%
*-lft-identity100.0%
unpow2100.0%
Simplified100.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 99.9%
Final simplification99.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 99.9%
Taylor expanded in x around 0 99.3%
neg-mul-199.3%
unsub-neg99.3%
*-commutative99.3%
unpow299.3%
associate-*l*99.3%
Simplified99.3%
Final simplification99.3%
(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(Float64(1.0 + 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[(N[(1.0 + N[(x * N[(x * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\left(1 + x \cdot \left(x \cdot 0.5\right)\right) - x\right)
\end{array}
Initial program 99.9%
Taylor expanded in x around 0 99.3%
neg-mul-199.3%
unsub-neg99.3%
*-commutative99.3%
unpow299.3%
associate-*l*99.3%
Simplified99.3%
associate-+r-99.3%
Applied egg-rr99.3%
Final simplification99.3%
(FPCore (x) :precision binary64 (* 2.0 (atan (+ 1.0 (- (* x x) x)))))
double code(double x) {
return 2.0 * atan((1.0 + ((x * x) - x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 + ((x * x) - x)))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 + ((x * x) - x)));
}
def code(x): return 2.0 * math.atan((1.0 + ((x * x) - x)))
function code(x) return Float64(2.0 * atan(Float64(1.0 + Float64(Float64(x * x) - x)))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 + ((x * x) - x))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 + N[(N[(x * x), $MachinePrecision] - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(1 + \left(x \cdot x - x\right)\right)
\end{array}
Initial program 99.9%
clear-num99.9%
sqrt-div99.9%
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in x around 0 99.1%
Taylor expanded in x around 0 99.2%
neg-mul-199.2%
+-commutative99.2%
unsub-neg99.2%
unpow299.2%
Simplified99.2%
Final simplification99.2%
(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 99.9%
clear-num99.9%
sqrt-div99.9%
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in x around 0 99.1%
Final simplification99.1%
(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 99.9%
Taylor expanded in x around 0 99.1%
neg-mul-199.1%
sub-neg99.1%
Simplified99.1%
Final simplification99.1%
herbie shell --seed 2023230
(FPCore (x)
:name "arccos"
:precision binary64
(* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))