
(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 (exp (* (- (log1p (- x)) (log1p x)) 0.5)))))
double code(double x) {
return 2.0 * atan(exp(((log1p(-x) - log1p(x)) * 0.5)));
}
public static double code(double x) {
return 2.0 * Math.atan(Math.exp(((Math.log1p(-x) - Math.log1p(x)) * 0.5)));
}
def code(x): return 2.0 * math.atan(math.exp(((math.log1p(-x) - math.log1p(x)) * 0.5)))
function code(x) return Float64(2.0 * atan(exp(Float64(Float64(log1p(Float64(-x)) - log1p(x)) * 0.5)))) end
code[x_] := N[(2.0 * N[ArcTan[N[Exp[N[(N[(N[Log[1 + (-x)], $MachinePrecision] - N[Log[1 + x], $MachinePrecision]), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(e^{\left(\mathsf{log1p}\left(-x\right) - \mathsf{log1p}\left(x\right)\right) \cdot 0.5}\right)
\end{array}
Initial program 100.0%
pow1/2100.0%
pow-to-exp100.0%
log-div100.0%
sub-neg100.0%
log1p-def100.0%
log1p-udef100.0%
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* 2.0 (atan (sqrt (* (- 1.0 x) (/ 1.0 (+ x 1.0)))))))
double code(double x) {
return 2.0 * atan(sqrt(((1.0 - x) * (1.0 / (x + 1.0)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(sqrt(((1.0d0 - x) * (1.0d0 / (x + 1.0d0)))))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.sqrt(((1.0 - x) * (1.0 / (x + 1.0)))));
}
def code(x): return 2.0 * math.atan(math.sqrt(((1.0 - x) * (1.0 / (x + 1.0)))))
function code(x) return Float64(2.0 * atan(sqrt(Float64(Float64(1.0 - x) * Float64(1.0 / Float64(x + 1.0)))))) end
function tmp = code(x) tmp = 2.0 * atan(sqrt(((1.0 - x) * (1.0 / (x + 1.0))))); end
code[x_] := N[(2.0 * N[ArcTan[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] * N[(1.0 / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\sqrt{\left(1 - x\right) \cdot \frac{1}{x + 1}}\right)
\end{array}
Initial program 100.0%
clear-num100.0%
associate-/r/100.0%
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* 2.0 (atan (/ 1.0 (sqrt (/ (+ x 1.0) (- 1.0 x)))))))
double code(double x) {
return 2.0 * atan((1.0 / sqrt(((x + 1.0) / (1.0 - x)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 / sqrt(((x + 1.0d0) / (1.0d0 - x)))))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 / Math.sqrt(((x + 1.0) / (1.0 - x)))));
}
def code(x): return 2.0 * math.atan((1.0 / math.sqrt(((x + 1.0) / (1.0 - x)))))
function code(x) return Float64(2.0 * atan(Float64(1.0 / sqrt(Float64(Float64(x + 1.0) / Float64(1.0 - x)))))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 / sqrt(((x + 1.0) / (1.0 - x))))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 / N[Sqrt[N[(N[(x + 1.0), $MachinePrecision] / N[(1.0 - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{1}{\sqrt{\frac{x + 1}{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) (+ x 1.0))))))
double code(double x) {
return 2.0 * atan(sqrt(((1.0 - x) / (x + 1.0))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(sqrt(((1.0d0 - x) / (x + 1.0d0))))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.sqrt(((1.0 - x) / (x + 1.0))));
}
def code(x): return 2.0 * math.atan(math.sqrt(((1.0 - x) / (x + 1.0))))
function code(x) return Float64(2.0 * atan(sqrt(Float64(Float64(1.0 - x) / Float64(x + 1.0))))) end
function tmp = code(x) tmp = 2.0 * atan(sqrt(((1.0 - x) / (x + 1.0)))); end
code[x_] := N[(2.0 * N[ArcTan[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\sqrt{\frac{1 - x}{x + 1}}\right)
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* 2.0 (atan (exp (- x)))))
double code(double x) {
return 2.0 * atan(exp(-x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(exp(-x))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.exp(-x));
}
def code(x): return 2.0 * math.atan(math.exp(-x))
function code(x) return Float64(2.0 * atan(exp(Float64(-x)))) end
function tmp = code(x) tmp = 2.0 * atan(exp(-x)); end
code[x_] := N[(2.0 * N[ArcTan[N[Exp[(-x)], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(e^{-x}\right)
\end{array}
Initial program 100.0%
pow1/2100.0%
pow-to-exp100.0%
log-div100.0%
sub-neg100.0%
log1p-def100.0%
log1p-udef100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0 99.3%
*-commutative99.3%
Simplified99.3%
Taylor expanded in x around 0 99.3%
mul-1-neg99.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(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 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 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(Float64(1.0 - x) + Float64(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[(N[(1.0 - x), $MachinePrecision] + N[(x * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\left(1 - x\right) + x \cdot 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 99.1%
+-commutative99.1%
Simplified99.1%
Taylor expanded in x around 0 99.1%
associate-+r+99.1%
mul-1-neg99.1%
sub-neg99.1%
unpow299.1%
Simplified99.1%
Final simplification99.1%
(FPCore (x) :precision binary64 (* 2.0 (atan (/ 1.0 (+ x 1.0)))))
double code(double x) {
return 2.0 * atan((1.0 / (x + 1.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 / (x + 1.0d0)))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 / (x + 1.0)));
}
def code(x): return 2.0 * math.atan((1.0 / (x + 1.0)))
function code(x) return Float64(2.0 * atan(Float64(1.0 / Float64(x + 1.0)))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 / (x + 1.0))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 / N[(x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(\frac{1}{x + 1}\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%
+-commutative99.1%
Simplified99.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 100.0%
Taylor expanded in x around 0 99.1%
neg-mul-199.1%
sub-neg99.1%
Simplified99.1%
Final simplification99.1%
(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%
pow1/2100.0%
pow-to-exp100.0%
log-div100.0%
sub-neg100.0%
log1p-def100.0%
log1p-udef100.0%
Applied egg-rr100.0%
Applied egg-rr97.9%
expm1-def96.4%
expm1-log1p97.9%
*-inverses97.9%
metadata-eval97.9%
Simplified97.9%
Final simplification97.9%
herbie shell --seed 2023238
(FPCore (x)
:name "arccos"
:precision binary64
(* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))