
(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 5 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) (+ 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 (pow (- 1.0 (* x 0.6666666666666666)) 1.5))))
double code(double x) {
return 2.0 * atan(pow((1.0 - (x * 0.6666666666666666)), 1.5));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan(((1.0d0 - (x * 0.6666666666666666d0)) ** 1.5d0))
end function
public static double code(double x) {
return 2.0 * Math.atan(Math.pow((1.0 - (x * 0.6666666666666666)), 1.5));
}
def code(x): return 2.0 * math.atan(math.pow((1.0 - (x * 0.6666666666666666)), 1.5))
function code(x) return Float64(2.0 * atan((Float64(1.0 - Float64(x * 0.6666666666666666)) ^ 1.5))) end
function tmp = code(x) tmp = 2.0 * atan(((1.0 - (x * 0.6666666666666666)) ^ 1.5)); end
code[x_] := N[(2.0 * N[ArcTan[N[Power[N[(1.0 - N[(x * 0.6666666666666666), $MachinePrecision]), $MachinePrecision], 1.5], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left({\left(1 - x \cdot 0.6666666666666666\right)}^{1.5}\right)
\end{array}
Initial program 100.0%
add-cube-cbrt100.0%
pow3100.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%
+-commutative99.1%
*-commutative99.1%
fma-undefine99.1%
metadata-eval99.1%
pow-sqr99.1%
rem-sqrt-square99.1%
unpow199.1%
sqr-pow99.1%
fabs-sqr99.1%
sqr-pow99.1%
unpow199.1%
Simplified99.1%
Taylor expanded in x around inf 99.1%
Final simplification99.1%
(FPCore (x) :precision binary64 (* 2.0 (atan (+ 1.0 (* x (+ -1.0 (* x 0.16666666666666666)))))))
double code(double x) {
return 2.0 * atan((1.0 + (x * (-1.0 + (x * 0.16666666666666666)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * atan((1.0d0 + (x * ((-1.0d0) + (x * 0.16666666666666666d0)))))
end function
public static double code(double x) {
return 2.0 * Math.atan((1.0 + (x * (-1.0 + (x * 0.16666666666666666)))));
}
def code(x): return 2.0 * math.atan((1.0 + (x * (-1.0 + (x * 0.16666666666666666)))))
function code(x) return Float64(2.0 * atan(Float64(1.0 + Float64(x * Float64(-1.0 + Float64(x * 0.16666666666666666)))))) end
function tmp = code(x) tmp = 2.0 * atan((1.0 + (x * (-1.0 + (x * 0.16666666666666666))))); end
code[x_] := N[(2.0 * N[ArcTan[N[(1.0 + N[(x * N[(-1.0 + N[(x * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \tan^{-1} \left(1 + x \cdot \left(-1 + x \cdot 0.16666666666666666\right)\right)
\end{array}
Initial program 100.0%
add-cube-cbrt100.0%
pow3100.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%
unpow299.1%
associate-*r*99.1%
distribute-rgt-out99.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.0%
*-commutative99.0%
Simplified99.0%
add-cbrt-cube99.0%
pow1/399.0%
add-sqr-sqrt99.0%
pow199.0%
pow1/299.0%
pow-prod-up99.0%
+-commutative99.0%
fma-define99.0%
metadata-eval99.0%
Applied egg-rr99.0%
unpow1/399.0%
Simplified99.0%
Taylor expanded in x around 0 99.0%
mul-1-neg99.0%
unsub-neg99.0%
Simplified99.0%
Final simplification99.0%
(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%
Taylor expanded in x around 0 99.0%
*-commutative99.0%
Simplified99.0%
add-cbrt-cube99.0%
pow1/399.0%
add-sqr-sqrt99.0%
pow199.0%
pow1/299.0%
pow-prod-up99.0%
+-commutative99.0%
fma-define99.0%
metadata-eval99.0%
Applied egg-rr99.0%
unpow1/399.0%
Simplified99.0%
Taylor expanded in x around 0 98.2%
Final simplification98.2%
herbie shell --seed 2024041
(FPCore (x)
:name "arccos"
:precision binary64
(* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))