
(FPCore (x) :precision binary64 (acos (- 1.0 x)))
double code(double x) {
return acos((1.0 - x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = acos((1.0d0 - x))
end function
public static double code(double x) {
return Math.acos((1.0 - x));
}
def code(x): return math.acos((1.0 - x))
function code(x) return acos(Float64(1.0 - x)) end
function tmp = code(x) tmp = acos((1.0 - x)); end
code[x_] := N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(1 - x\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (acos (- 1.0 x)))
double code(double x) {
return acos((1.0 - x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = acos((1.0d0 - x))
end function
public static double code(double x) {
return Math.acos((1.0 - x));
}
def code(x): return math.acos((1.0 - x))
function code(x) return acos(Float64(1.0 - x)) end
function tmp = code(x) tmp = acos((1.0 - x)); end
code[x_] := N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(1 - x\right)
\end{array}
(FPCore (x) :precision binary64 (fma (pow (pow (* PI 0.5) 0.3333333333333333) 2.0) (cbrt (* PI 0.5)) (- (asin (- 1.0 x)))))
double code(double x) {
return fma(pow(pow((((double) M_PI) * 0.5), 0.3333333333333333), 2.0), cbrt((((double) M_PI) * 0.5)), -asin((1.0 - x)));
}
function code(x) return fma(((Float64(pi * 0.5) ^ 0.3333333333333333) ^ 2.0), cbrt(Float64(pi * 0.5)), Float64(-asin(Float64(1.0 - x)))) end
code[x_] := N[(N[Power[N[Power[N[(Pi * 0.5), $MachinePrecision], 0.3333333333333333], $MachinePrecision], 2.0], $MachinePrecision] * N[Power[N[(Pi * 0.5), $MachinePrecision], 1/3], $MachinePrecision] + (-N[ArcSin[N[(1.0 - x), $MachinePrecision]], $MachinePrecision])), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left({\left({\left(\pi \cdot 0.5\right)}^{0.3333333333333333}\right)}^{2}, \sqrt[3]{\pi \cdot 0.5}, -\sin^{-1} \left(1 - x\right)\right)
\end{array}
Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
add-cube-cbrt5.4%
fma-neg5.4%
pow25.4%
Applied egg-rr5.4%
pow1/310.6%
Applied egg-rr10.6%
Final simplification10.6%
(FPCore (x) :precision binary64 (- (* PI 0.5) (pow (cbrt (cbrt (asin (- 1.0 x)))) 9.0)))
double code(double x) {
return (((double) M_PI) * 0.5) - pow(cbrt(cbrt(asin((1.0 - x)))), 9.0);
}
public static double code(double x) {
return (Math.PI * 0.5) - Math.pow(Math.cbrt(Math.cbrt(Math.asin((1.0 - x)))), 9.0);
}
function code(x) return Float64(Float64(pi * 0.5) - (cbrt(cbrt(asin(Float64(1.0 - x)))) ^ 9.0)) end
code[x_] := N[(N[(Pi * 0.5), $MachinePrecision] - N[Power[N[Power[N[Power[N[ArcSin[N[(1.0 - x), $MachinePrecision]], $MachinePrecision], 1/3], $MachinePrecision], 1/3], $MachinePrecision], 9.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\pi \cdot 0.5 - {\left(\sqrt[3]{\sqrt[3]{\sin^{-1} \left(1 - x\right)}}\right)}^{9}
\end{array}
Initial program 7.2%
acos-asin7.2%
sub-neg7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
sub-neg7.2%
Simplified7.2%
add-cube-cbrt10.5%
pow310.5%
Applied egg-rr10.5%
add-cube-cbrt10.5%
pow310.5%
Applied egg-rr10.5%
pow-pow10.6%
pow-to-exp10.6%
metadata-eval10.6%
Applied egg-rr10.6%
exp-to-pow10.6%
Simplified10.6%
Final simplification10.6%
(FPCore (x) :precision binary64 (if (<= (- 1.0 x) 1.0) (- (* PI 0.5) (cbrt (pow (asin (- 1.0 x)) 3.0))) (- PI (acos (- 1.0 x)))))
double code(double x) {
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = (((double) M_PI) * 0.5) - cbrt(pow(asin((1.0 - x)), 3.0));
} else {
tmp = ((double) M_PI) - acos((1.0 - x));
}
return tmp;
}
public static double code(double x) {
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = (Math.PI * 0.5) - Math.cbrt(Math.pow(Math.asin((1.0 - x)), 3.0));
} else {
tmp = Math.PI - Math.acos((1.0 - x));
}
return tmp;
}
function code(x) tmp = 0.0 if (Float64(1.0 - x) <= 1.0) tmp = Float64(Float64(pi * 0.5) - cbrt((asin(Float64(1.0 - x)) ^ 3.0))); else tmp = Float64(pi - acos(Float64(1.0 - x))); end return tmp end
code[x_] := If[LessEqual[N[(1.0 - x), $MachinePrecision], 1.0], N[(N[(Pi * 0.5), $MachinePrecision] - N[Power[N[Power[N[ArcSin[N[(1.0 - x), $MachinePrecision]], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], N[(Pi - N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;1 - x \leq 1:\\
\;\;\;\;\pi \cdot 0.5 - \sqrt[3]{{\sin^{-1} \left(1 - x\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\pi - \cos^{-1} \left(1 - x\right)\\
\end{array}
\end{array}
if (-.f64 1 x) < 1Initial program 7.2%
acos-asin7.2%
sub-neg7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
sub-neg7.2%
Simplified7.2%
add-cbrt-cube5.4%
pow35.4%
Applied egg-rr5.4%
if 1 < (-.f64 1 x) Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
metadata-eval7.2%
div-inv7.2%
acos-asin7.2%
expm1-log1p-u7.2%
expm1-udef7.2%
log1p-udef7.2%
add-exp-log7.2%
associate--l+7.2%
+-commutative7.2%
sub-neg7.2%
metadata-eval7.2%
Applied egg-rr7.2%
associate-+l+7.2%
metadata-eval7.2%
+-rgt-identity7.2%
acos-asin7.2%
div-inv7.2%
metadata-eval7.2%
add-cube-cbrt10.5%
unpow310.5%
sub-neg10.5%
unpow310.5%
add-cube-cbrt7.2%
add-sqr-sqrt0.0%
sqrt-unprod6.8%
add-cube-cbrt6.8%
unpow36.8%
add-cube-cbrt6.8%
unpow36.8%
Applied egg-rr6.8%
distribute-lft-out6.8%
metadata-eval6.8%
*-rgt-identity6.8%
Simplified6.8%
Final simplification5.4%
(FPCore (x) :precision binary64 (- (* PI 0.5) (pow (cbrt (asin (- 1.0 x))) 3.0)))
double code(double x) {
return (((double) M_PI) * 0.5) - pow(cbrt(asin((1.0 - x))), 3.0);
}
public static double code(double x) {
return (Math.PI * 0.5) - Math.pow(Math.cbrt(Math.asin((1.0 - x))), 3.0);
}
function code(x) return Float64(Float64(pi * 0.5) - (cbrt(asin(Float64(1.0 - x))) ^ 3.0)) end
code[x_] := N[(N[(Pi * 0.5), $MachinePrecision] - N[Power[N[Power[N[ArcSin[N[(1.0 - x), $MachinePrecision]], $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\pi \cdot 0.5 - {\left(\sqrt[3]{\sin^{-1} \left(1 - x\right)}\right)}^{3}
\end{array}
Initial program 7.2%
acos-asin7.2%
sub-neg7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
sub-neg7.2%
Simplified7.2%
add-cube-cbrt10.5%
pow310.5%
Applied egg-rr10.5%
Final simplification10.5%
(FPCore (x) :precision binary64 (- (* PI 0.5) (pow (sqrt (asin (- 1.0 x))) 2.0)))
double code(double x) {
return (((double) M_PI) * 0.5) - pow(sqrt(asin((1.0 - x))), 2.0);
}
public static double code(double x) {
return (Math.PI * 0.5) - Math.pow(Math.sqrt(Math.asin((1.0 - x))), 2.0);
}
def code(x): return (math.pi * 0.5) - math.pow(math.sqrt(math.asin((1.0 - x))), 2.0)
function code(x) return Float64(Float64(pi * 0.5) - (sqrt(asin(Float64(1.0 - x))) ^ 2.0)) end
function tmp = code(x) tmp = (pi * 0.5) - (sqrt(asin((1.0 - x))) ^ 2.0); end
code[x_] := N[(N[(Pi * 0.5), $MachinePrecision] - N[Power[N[Sqrt[N[ArcSin[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\pi \cdot 0.5 - {\left(\sqrt{\sin^{-1} \left(1 - x\right)}\right)}^{2}
\end{array}
Initial program 7.2%
acos-asin7.2%
sub-neg7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
sub-neg7.2%
Simplified7.2%
add-sqr-sqrt10.6%
pow210.6%
Applied egg-rr10.6%
Final simplification10.6%
(FPCore (x) :precision binary64 (let* ((t_0 (acos (- 1.0 x)))) (if (<= (- 1.0 x) 1.0) (* 2.0 (log (exp (* 0.5 t_0)))) (- PI t_0))))
double code(double x) {
double t_0 = acos((1.0 - x));
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = 2.0 * log(exp((0.5 * t_0)));
} else {
tmp = ((double) M_PI) - t_0;
}
return tmp;
}
public static double code(double x) {
double t_0 = Math.acos((1.0 - x));
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = 2.0 * Math.log(Math.exp((0.5 * t_0)));
} else {
tmp = Math.PI - t_0;
}
return tmp;
}
def code(x): t_0 = math.acos((1.0 - x)) tmp = 0 if (1.0 - x) <= 1.0: tmp = 2.0 * math.log(math.exp((0.5 * t_0))) else: tmp = math.pi - t_0 return tmp
function code(x) t_0 = acos(Float64(1.0 - x)) tmp = 0.0 if (Float64(1.0 - x) <= 1.0) tmp = Float64(2.0 * log(exp(Float64(0.5 * t_0)))); else tmp = Float64(pi - t_0); end return tmp end
function tmp_2 = code(x) t_0 = acos((1.0 - x)); tmp = 0.0; if ((1.0 - x) <= 1.0) tmp = 2.0 * log(exp((0.5 * t_0))); else tmp = pi - t_0; end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(1.0 - x), $MachinePrecision], 1.0], N[(2.0 * N[Log[N[Exp[N[(0.5 * t$95$0), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(Pi - t$95$0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(1 - x\right)\\
\mathbf{if}\;1 - x \leq 1:\\
\;\;\;\;2 \cdot \log \left(e^{0.5 \cdot t_0}\right)\\
\mathbf{else}:\\
\;\;\;\;\pi - t_0\\
\end{array}
\end{array}
if (-.f64 1 x) < 1Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
metadata-eval7.2%
div-inv7.2%
acos-asin7.2%
add-log-exp7.2%
add-sqr-sqrt7.2%
log-prod7.2%
Applied egg-rr7.2%
count-27.2%
Simplified7.2%
pow1/27.2%
pow-exp7.2%
Applied egg-rr7.2%
if 1 < (-.f64 1 x) Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
metadata-eval7.2%
div-inv7.2%
acos-asin7.2%
expm1-log1p-u7.2%
expm1-udef7.2%
log1p-udef7.2%
add-exp-log7.2%
associate--l+7.2%
+-commutative7.2%
sub-neg7.2%
metadata-eval7.2%
Applied egg-rr7.2%
associate-+l+7.2%
metadata-eval7.2%
+-rgt-identity7.2%
acos-asin7.2%
div-inv7.2%
metadata-eval7.2%
add-cube-cbrt10.5%
unpow310.5%
sub-neg10.5%
unpow310.5%
add-cube-cbrt7.2%
add-sqr-sqrt0.0%
sqrt-unprod6.8%
add-cube-cbrt6.8%
unpow36.8%
add-cube-cbrt6.8%
unpow36.8%
Applied egg-rr6.8%
distribute-lft-out6.8%
metadata-eval6.8%
*-rgt-identity6.8%
Simplified6.8%
Final simplification7.2%
(FPCore (x) :precision binary64 (if (<= (- 1.0 x) 1.0) (- (* PI 0.5) (asin (- 1.0 x))) (- PI (acos (- 1.0 x)))))
double code(double x) {
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = (((double) M_PI) * 0.5) - asin((1.0 - x));
} else {
tmp = ((double) M_PI) - acos((1.0 - x));
}
return tmp;
}
public static double code(double x) {
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = (Math.PI * 0.5) - Math.asin((1.0 - x));
} else {
tmp = Math.PI - Math.acos((1.0 - x));
}
return tmp;
}
def code(x): tmp = 0 if (1.0 - x) <= 1.0: tmp = (math.pi * 0.5) - math.asin((1.0 - x)) else: tmp = math.pi - math.acos((1.0 - x)) return tmp
function code(x) tmp = 0.0 if (Float64(1.0 - x) <= 1.0) tmp = Float64(Float64(pi * 0.5) - asin(Float64(1.0 - x))); else tmp = Float64(pi - acos(Float64(1.0 - x))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if ((1.0 - x) <= 1.0) tmp = (pi * 0.5) - asin((1.0 - x)); else tmp = pi - acos((1.0 - x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[N[(1.0 - x), $MachinePrecision], 1.0], N[(N[(Pi * 0.5), $MachinePrecision] - N[ArcSin[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(Pi - N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;1 - x \leq 1:\\
\;\;\;\;\pi \cdot 0.5 - \sin^{-1} \left(1 - x\right)\\
\mathbf{else}:\\
\;\;\;\;\pi - \cos^{-1} \left(1 - x\right)\\
\end{array}
\end{array}
if (-.f64 1 x) < 1Initial program 7.2%
acos-asin7.2%
sub-neg7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
sub-neg7.2%
Simplified7.2%
if 1 < (-.f64 1 x) Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
metadata-eval7.2%
div-inv7.2%
acos-asin7.2%
expm1-log1p-u7.2%
expm1-udef7.2%
log1p-udef7.2%
add-exp-log7.2%
associate--l+7.2%
+-commutative7.2%
sub-neg7.2%
metadata-eval7.2%
Applied egg-rr7.2%
associate-+l+7.2%
metadata-eval7.2%
+-rgt-identity7.2%
acos-asin7.2%
div-inv7.2%
metadata-eval7.2%
add-cube-cbrt10.5%
unpow310.5%
sub-neg10.5%
unpow310.5%
add-cube-cbrt7.2%
add-sqr-sqrt0.0%
sqrt-unprod6.8%
add-cube-cbrt6.8%
unpow36.8%
add-cube-cbrt6.8%
unpow36.8%
Applied egg-rr6.8%
distribute-lft-out6.8%
metadata-eval6.8%
*-rgt-identity6.8%
Simplified6.8%
Final simplification7.2%
(FPCore (x) :precision binary64 (let* ((t_0 (acos (- 1.0 x)))) (if (<= (- 1.0 x) 1.0) (+ 1.0 (+ t_0 -1.0)) (- PI t_0))))
double code(double x) {
double t_0 = acos((1.0 - x));
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = 1.0 + (t_0 + -1.0);
} else {
tmp = ((double) M_PI) - t_0;
}
return tmp;
}
public static double code(double x) {
double t_0 = Math.acos((1.0 - x));
double tmp;
if ((1.0 - x) <= 1.0) {
tmp = 1.0 + (t_0 + -1.0);
} else {
tmp = Math.PI - t_0;
}
return tmp;
}
def code(x): t_0 = math.acos((1.0 - x)) tmp = 0 if (1.0 - x) <= 1.0: tmp = 1.0 + (t_0 + -1.0) else: tmp = math.pi - t_0 return tmp
function code(x) t_0 = acos(Float64(1.0 - x)) tmp = 0.0 if (Float64(1.0 - x) <= 1.0) tmp = Float64(1.0 + Float64(t_0 + -1.0)); else tmp = Float64(pi - t_0); end return tmp end
function tmp_2 = code(x) t_0 = acos((1.0 - x)); tmp = 0.0; if ((1.0 - x) <= 1.0) tmp = 1.0 + (t_0 + -1.0); else tmp = pi - t_0; end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(1.0 - x), $MachinePrecision], 1.0], N[(1.0 + N[(t$95$0 + -1.0), $MachinePrecision]), $MachinePrecision], N[(Pi - t$95$0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(1 - x\right)\\
\mathbf{if}\;1 - x \leq 1:\\
\;\;\;\;1 + \left(t_0 + -1\right)\\
\mathbf{else}:\\
\;\;\;\;\pi - t_0\\
\end{array}
\end{array}
if (-.f64 1 x) < 1Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
metadata-eval7.2%
div-inv7.2%
acos-asin7.2%
expm1-log1p-u7.2%
expm1-udef7.2%
log1p-udef7.2%
add-exp-log7.2%
associate--l+7.2%
+-commutative7.2%
sub-neg7.2%
metadata-eval7.2%
Applied egg-rr7.2%
if 1 < (-.f64 1 x) Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
metadata-eval7.2%
div-inv7.2%
acos-asin7.2%
expm1-log1p-u7.2%
expm1-udef7.2%
log1p-udef7.2%
add-exp-log7.2%
associate--l+7.2%
+-commutative7.2%
sub-neg7.2%
metadata-eval7.2%
Applied egg-rr7.2%
associate-+l+7.2%
metadata-eval7.2%
+-rgt-identity7.2%
acos-asin7.2%
div-inv7.2%
metadata-eval7.2%
add-cube-cbrt10.5%
unpow310.5%
sub-neg10.5%
unpow310.5%
add-cube-cbrt7.2%
add-sqr-sqrt0.0%
sqrt-unprod6.8%
add-cube-cbrt6.8%
unpow36.8%
add-cube-cbrt6.8%
unpow36.8%
Applied egg-rr6.8%
distribute-lft-out6.8%
metadata-eval6.8%
*-rgt-identity6.8%
Simplified6.8%
Final simplification7.2%
(FPCore (x) :precision binary64 (+ 1.0 (+ (acos (- 1.0 x)) -1.0)))
double code(double x) {
return 1.0 + (acos((1.0 - x)) + -1.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 + (acos((1.0d0 - x)) + (-1.0d0))
end function
public static double code(double x) {
return 1.0 + (Math.acos((1.0 - x)) + -1.0);
}
def code(x): return 1.0 + (math.acos((1.0 - x)) + -1.0)
function code(x) return Float64(1.0 + Float64(acos(Float64(1.0 - x)) + -1.0)) end
function tmp = code(x) tmp = 1.0 + (acos((1.0 - x)) + -1.0); end
code[x_] := N[(1.0 + N[(N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \left(\cos^{-1} \left(1 - x\right) + -1\right)
\end{array}
Initial program 7.2%
acos-asin7.2%
flip--7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
div-inv7.2%
metadata-eval7.2%
Applied egg-rr7.2%
flip--7.2%
metadata-eval7.2%
div-inv7.2%
acos-asin7.2%
expm1-log1p-u7.2%
expm1-udef7.2%
log1p-udef7.2%
add-exp-log7.2%
associate--l+7.2%
+-commutative7.2%
sub-neg7.2%
metadata-eval7.2%
Applied egg-rr7.2%
Final simplification7.2%
(FPCore (x) :precision binary64 (acos (- 1.0 x)))
double code(double x) {
return acos((1.0 - x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = acos((1.0d0 - x))
end function
public static double code(double x) {
return Math.acos((1.0 - x));
}
def code(x): return math.acos((1.0 - x))
function code(x) return acos(Float64(1.0 - x)) end
function tmp = code(x) tmp = acos((1.0 - x)); end
code[x_] := N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(1 - x\right)
\end{array}
Initial program 7.2%
Final simplification7.2%
(FPCore (x) :precision binary64 (* 2.0 (asin (sqrt (/ x 2.0)))))
double code(double x) {
return 2.0 * asin(sqrt((x / 2.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0 * asin(sqrt((x / 2.0d0)))
end function
public static double code(double x) {
return 2.0 * Math.asin(Math.sqrt((x / 2.0)));
}
def code(x): return 2.0 * math.asin(math.sqrt((x / 2.0)))
function code(x) return Float64(2.0 * asin(sqrt(Float64(x / 2.0)))) end
function tmp = code(x) tmp = 2.0 * asin(sqrt((x / 2.0))); end
code[x_] := N[(2.0 * N[ArcSin[N[Sqrt[N[(x / 2.0), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \sin^{-1} \left(\sqrt{\frac{x}{2}}\right)
\end{array}
herbie shell --seed 2023274
(FPCore (x)
:name "bug323 (missed optimization)"
:precision binary64
:pre (and (<= 0.0 x) (<= x 0.5))
:herbie-target
(* 2.0 (asin (sqrt (/ x 2.0))))
(acos (- 1.0 x)))