
(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 14 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
(let* ((t_0 (cbrt (PI)))
(t_1 (/ 2.0 (PI)))
(t_2 (asin (- 1.0 x)))
(t_3 (fma 0.5 (PI) t_2)))
(/
(fma
(/ 2.0 (pow (pow (PI) 0.16666666666666666) 4.0))
(/ t_3 t_0)
(* (- t_3 (* (fma t_1 t_2 1.0) (acos (- 1.0 x)))) (/ -2.0 (PI))))
(*
(* t_1 t_3)
(/ 2.0 (* (pow (cbrt t_0) 3.0) (pow (cbrt (pow t_0 2.0)) 3.0)))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt[3]{\mathsf{PI}\left(\right)}\\
t_1 := \frac{2}{\mathsf{PI}\left(\right)}\\
t_2 := \sin^{-1} \left(1 - x\right)\\
t_3 := \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_2\right)\\
\frac{\mathsf{fma}\left(\frac{2}{{\left({\mathsf{PI}\left(\right)}^{0.16666666666666666}\right)}^{4}}, \frac{t\_3}{t\_0}, \left(t\_3 - \mathsf{fma}\left(t\_1, t\_2, 1\right) \cdot \cos^{-1} \left(1 - x\right)\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right)}{\left(t\_1 \cdot t\_3\right) \cdot \frac{2}{{\left(\sqrt[3]{t\_0}\right)}^{3} \cdot {\left(\sqrt[3]{{t\_0}^{2}}\right)}^{3}}}
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
clear-numN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites6.5%
Applied rewrites9.8%
lift-pow.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
metadata-evalN/A
metadata-evalN/A
sqrt-pow2N/A
lift-sqrt.f64N/A
pow-powN/A
pow1/3N/A
pow2N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f64N/A
pow1/3N/A
lift-sqrt.f64N/A
sqrt-pow2N/A
lower-pow.f64N/A
metadata-evalN/A
metadata-eval9.9
Applied rewrites9.9%
unpow1N/A
metadata-evalN/A
pow-powN/A
rem-3cbrt-lftN/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
unpow2N/A
lift-pow.f64N/A
unpow-prod-downN/A
metadata-evalN/A
pow-powN/A
lift-pow.f64N/A
cube-prodN/A
lower-*.f64N/A
Applied rewrites9.9%
Final simplification9.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (cbrt (PI)))
(t_1 (/ 2.0 (PI)))
(t_2 (asin (- 1.0 x)))
(t_3 (fma 0.5 (PI) t_2)))
(/
(fma
(/ 2.0 (pow (pow (PI) 0.16666666666666666) 4.0))
(/ t_3 t_0)
(* (- t_3 (* (fma t_1 t_2 1.0) (acos (- 1.0 x)))) (/ -2.0 (PI))))
(* (* (fma 0.5 (pow (pow t_0 2.0) 1.5) t_2) t_1) t_1))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt[3]{\mathsf{PI}\left(\right)}\\
t_1 := \frac{2}{\mathsf{PI}\left(\right)}\\
t_2 := \sin^{-1} \left(1 - x\right)\\
t_3 := \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_2\right)\\
\frac{\mathsf{fma}\left(\frac{2}{{\left({\mathsf{PI}\left(\right)}^{0.16666666666666666}\right)}^{4}}, \frac{t\_3}{t\_0}, \left(t\_3 - \mathsf{fma}\left(t\_1, t\_2, 1\right) \cdot \cos^{-1} \left(1 - x\right)\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right)}{\left(\mathsf{fma}\left(0.5, {\left({t\_0}^{2}\right)}^{1.5}, t\_2\right) \cdot t\_1\right) \cdot t\_1}
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
clear-numN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites6.5%
Applied rewrites9.8%
lift-pow.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
metadata-evalN/A
metadata-evalN/A
sqrt-pow2N/A
lift-sqrt.f64N/A
pow-powN/A
pow1/3N/A
pow2N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f64N/A
pow1/3N/A
lift-sqrt.f64N/A
sqrt-pow2N/A
lower-pow.f64N/A
metadata-evalN/A
metadata-eval9.9
Applied rewrites9.9%
unpow1N/A
metadata-evalN/A
pow-powN/A
pow1/3N/A
lift-cbrt.f64N/A
metadata-evalN/A
metadata-evalN/A
pow-powN/A
lift-pow.f64N/A
lower-pow.f64N/A
metadata-eval9.9
Applied rewrites9.9%
Final simplification9.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (asin (- 1.0 x)))
(t_1 (fma (/ 2.0 (PI)) t_0 1.0))
(t_2 (fma 0.5 (PI) t_0)))
(/
(fma
(/ 2.0 (pow (pow (PI) 0.16666666666666666) 4.0))
(/ t_2 (cbrt (PI)))
(* (- t_2 (* t_1 (acos (- 1.0 x)))) (/ -2.0 (PI))))
(/ t_1 (* 0.5 (PI))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(1 - x\right)\\
t_1 := \mathsf{fma}\left(\frac{2}{\mathsf{PI}\left(\right)}, t\_0, 1\right)\\
t_2 := \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_0\right)\\
\frac{\mathsf{fma}\left(\frac{2}{{\left({\mathsf{PI}\left(\right)}^{0.16666666666666666}\right)}^{4}}, \frac{t\_2}{\sqrt[3]{\mathsf{PI}\left(\right)}}, \left(t\_2 - t\_1 \cdot \cos^{-1} \left(1 - x\right)\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right)}{\frac{t\_1}{0.5 \cdot \mathsf{PI}\left(\right)}}
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
clear-numN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites6.5%
Applied rewrites9.8%
lift-pow.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
metadata-evalN/A
metadata-evalN/A
sqrt-pow2N/A
lift-sqrt.f64N/A
pow-powN/A
pow1/3N/A
pow2N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f64N/A
pow1/3N/A
lift-sqrt.f64N/A
sqrt-pow2N/A
lower-pow.f64N/A
metadata-evalN/A
metadata-eval9.9
Applied rewrites9.9%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
rem-3cbrt-lftN/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
lift-cbrt.f64N/A
unpow2N/A
lift-pow.f64N/A
frac-timesN/A
lift-/.f64N/A
lift-/.f64N/A
Applied rewrites9.9%
Final simplification9.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (asin (- 1.0 x))) (t_1 (fma 0.5 (PI) t_0)))
(/
(fma
(/ 2.0 (pow (pow (PI) 0.16666666666666666) 4.0))
(/ t_1 (cbrt (PI)))
(* (- t_1 (* (fma (/ 2.0 (PI)) t_0 1.0) (acos (- 1.0 x)))) (/ -2.0 (PI))))
(/ (* (* t_1 2.0) 2.0) (* (PI) (PI))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(1 - x\right)\\
t_1 := \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_0\right)\\
\frac{\mathsf{fma}\left(\frac{2}{{\left({\mathsf{PI}\left(\right)}^{0.16666666666666666}\right)}^{4}}, \frac{t\_1}{\sqrt[3]{\mathsf{PI}\left(\right)}}, \left(t\_1 - \mathsf{fma}\left(\frac{2}{\mathsf{PI}\left(\right)}, t\_0, 1\right) \cdot \cos^{-1} \left(1 - x\right)\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right)}{\frac{\left(t\_1 \cdot 2\right) \cdot 2}{\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)}}
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
clear-numN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites6.5%
Applied rewrites9.8%
lift-pow.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
metadata-evalN/A
metadata-evalN/A
sqrt-pow2N/A
lift-sqrt.f64N/A
pow-powN/A
pow1/3N/A
pow2N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f64N/A
pow1/3N/A
lift-sqrt.f64N/A
sqrt-pow2N/A
lower-pow.f64N/A
metadata-evalN/A
metadata-eval9.9
Applied rewrites9.9%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-/.f64N/A
associate-*l/N/A
lift-/.f64N/A
frac-timesN/A
*-commutativeN/A
lift-*.f64N/A
lower-/.f64N/A
Applied rewrites9.9%
Final simplification9.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (asin (- 1.0 x))) (t_1 (fma (PI) 0.5 t_0)) (t_2 (/ 2.0 (PI))))
(/
(/
(fma
(* (- t_1 (* (fma t_2 t_0 1.0) (acos (- 1.0 x)))) -2.0)
(PI)
(* (* t_1 2.0) (PI)))
(* (PI) (PI)))
(* (* t_2 (fma 0.5 (PI) t_0)) t_2))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(1 - x\right)\\
t_1 := \mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, t\_0\right)\\
t_2 := \frac{2}{\mathsf{PI}\left(\right)}\\
\frac{\frac{\mathsf{fma}\left(\left(t\_1 - \mathsf{fma}\left(t\_2, t\_0, 1\right) \cdot \cos^{-1} \left(1 - x\right)\right) \cdot -2, \mathsf{PI}\left(\right), \left(t\_1 \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right)}{\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)}}{\left(t\_2 \cdot \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_0\right)\right) \cdot t\_2}
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
clear-numN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites6.5%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
*-lft-identityN/A
lift-*.f64N/A
*-commutativeN/A
Applied rewrites9.9%
Applied rewrites9.9%
Final simplification9.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (asin (- 1.0 x))) (t_1 (fma 0.5 (PI) t_0)) (t_2 (/ 2.0 (PI))))
(/
(fma
t_1
t_2
(* (- t_1 (* (fma t_2 t_0 1.0) (acos (- 1.0 x)))) (/ -2.0 (PI))))
(/ (* (* (fma (PI) 0.5 t_0) 2.0) 2.0) (* (PI) (PI))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(1 - x\right)\\
t_1 := \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_0\right)\\
t_2 := \frac{2}{\mathsf{PI}\left(\right)}\\
\frac{\mathsf{fma}\left(t\_1, t\_2, \left(t\_1 - \mathsf{fma}\left(t\_2, t\_0, 1\right) \cdot \cos^{-1} \left(1 - x\right)\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right)}{\frac{\left(\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, t\_0\right) \cdot 2\right) \cdot 2}{\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)}}
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
clear-numN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites6.5%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
*-lft-identityN/A
lift-*.f64N/A
*-commutativeN/A
Applied rewrites9.9%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lift-/.f64N/A
associate-*l/N/A
lift-/.f64N/A
frac-timesN/A
lift-PI.f64N/A
lift-PI.f64N/A
lower-/.f64N/A
Applied rewrites9.9%
Final simplification9.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (* (PI) (PI))) (t_1 (asin (- 1.0 x))))
(/
(fma (* 0.25 t_0) (* 0.5 (PI)) (- (pow t_1 3.0)))
(fma 0.25 t_0 (* (fma 0.5 (PI) t_1) t_1)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\\
t_1 := \sin^{-1} \left(1 - x\right)\\
\frac{\mathsf{fma}\left(0.25 \cdot t\_0, 0.5 \cdot \mathsf{PI}\left(\right), -{t\_1}^{3}\right)}{\mathsf{fma}\left(0.25, t\_0, \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_1\right) \cdot t\_1\right)}
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
flip3--N/A
lower-/.f64N/A
Applied rewrites6.4%
lift-pow.f64N/A
rem-cube-cbrtN/A
lift-cbrt.f64N/A
pow-powN/A
lower-pow.f64N/A
metadata-eval4.6
Applied rewrites4.6%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
*-commutativeN/A
lift-pow.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
pow-powN/A
metadata-evalN/A
lift-pow.f64N/A
lift-pow.f64N/A
metadata-evalN/A
unpow-prod-downN/A
*-commutativeN/A
unpow3N/A
lower-fma.f64N/A
Applied rewrites9.9%
Final simplification9.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (sqrt (PI))))
(fma
(pow (PI) 0.6666666666666666)
(* (cbrt (* t_0 t_0)) 0.5)
(- (asin (- 1.0 x))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{PI}\left(\right)}\\
\mathsf{fma}\left({\mathsf{PI}\left(\right)}^{0.6666666666666666}, \sqrt[3]{t\_0 \cdot t\_0} \cdot 0.5, -\sin^{-1} \left(1 - x\right)\right)
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
sub-negN/A
div-invN/A
add-cube-cbrtN/A
associate-*l*N/A
lower-fma.f64N/A
pow2N/A
lower-pow.f64N/A
lower-cbrt.f64N/A
lower-PI.f64N/A
lower-*.f64N/A
lower-cbrt.f64N/A
lower-PI.f64N/A
metadata-evalN/A
lower-neg.f64N/A
lower-asin.f644.6
Applied rewrites4.6%
lift-PI.f64N/A
add-sqr-sqrtN/A
lower-*.f64N/A
lift-PI.f64N/A
lower-sqrt.f64N/A
lift-PI.f64N/A
lower-sqrt.f644.6
Applied rewrites4.6%
lift-pow.f64N/A
lift-cbrt.f64N/A
pow1/3N/A
pow-powN/A
lower-pow.f64N/A
metadata-eval9.8
Applied rewrites9.8%
(FPCore (x) :precision binary64 (let* ((t_0 (sqrt (PI)))) (fma (PI) 0.5 (- (fma (* t_0 0.5) t_0 (- (acos (- 1.0 x))))))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{PI}\left(\right)}\\
\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, -\mathsf{fma}\left(t\_0 \cdot 0.5, t\_0, -\cos^{-1} \left(1 - x\right)\right)\right)
\end{array}
\end{array}
Initial program 6.5%
lift-acos.f64N/A
acos-asinN/A
sub-negN/A
div-invN/A
lower-fma.f64N/A
lower-PI.f64N/A
metadata-evalN/A
lower-neg.f64N/A
lower-asin.f646.5
Applied rewrites6.5%
lift-asin.f64N/A
asin-acosN/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
*-commutativeN/A
lift-acos.f64N/A
sub-negN/A
rem-square-sqrtN/A
lift-sqrt.f64N/A
lift-sqrt.f64N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower-neg.f649.8
Applied rewrites9.8%
Final simplification9.8%
(FPCore (x)
:precision binary64
(let* ((t_0 (sqrt (PI))))
(if (<= x 5.6e-17)
(acos (- x))
(fma (* t_0 t_0) 0.5 (- (asin (- 1.0 x)))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{PI}\left(\right)}\\
\mathbf{if}\;x \leq 5.6 \cdot 10^{-17}:\\
\;\;\;\;\cos^{-1} \left(-x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t\_0 \cdot t\_0, 0.5, -\sin^{-1} \left(1 - x\right)\right)\\
\end{array}
\end{array}
if x < 5.5999999999999998e-17Initial program 3.9%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f646.5
Applied rewrites6.5%
if 5.5999999999999998e-17 < x Initial program 69.4%
lift-acos.f64N/A
acos-asinN/A
sub-negN/A
div-invN/A
lower-fma.f64N/A
lower-PI.f64N/A
metadata-evalN/A
lower-neg.f64N/A
lower-asin.f6469.3
Applied rewrites69.3%
rem-square-sqrtN/A
lift-sqrt.f64N/A
lift-sqrt.f64N/A
lift-*.f6469.6
Applied rewrites69.6%
(FPCore (x)
:precision binary64
(let* ((t_0 (sqrt (PI))))
(if (<= x 5.6e-17)
(acos (- x))
(fma t_0 (* t_0 0.5) (- (asin (- 1.0 x)))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{PI}\left(\right)}\\
\mathbf{if}\;x \leq 5.6 \cdot 10^{-17}:\\
\;\;\;\;\cos^{-1} \left(-x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t\_0, t\_0 \cdot 0.5, -\sin^{-1} \left(1 - x\right)\right)\\
\end{array}
\end{array}
if x < 5.5999999999999998e-17Initial program 3.9%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f646.5
Applied rewrites6.5%
if 5.5999999999999998e-17 < x Initial program 69.4%
lift-acos.f64N/A
acos-asinN/A
sub-negN/A
div-invN/A
add-sqr-sqrtN/A
associate-*l*N/A
lower-fma.f64N/A
lower-sqrt.f64N/A
lower-PI.f64N/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-PI.f64N/A
metadata-evalN/A
lower-neg.f64N/A
lower-asin.f6469.5
Applied rewrites69.5%
(FPCore (x) :precision binary64 (if (<= x 5.6e-17) (acos (- x)) (acos (- 1.0 x))))
double code(double x) {
double tmp;
if (x <= 5.6e-17) {
tmp = acos(-x);
} else {
tmp = acos((1.0 - x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 5.6d-17) then
tmp = acos(-x)
else
tmp = acos((1.0d0 - x))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 5.6e-17) {
tmp = Math.acos(-x);
} else {
tmp = Math.acos((1.0 - x));
}
return tmp;
}
def code(x): tmp = 0 if x <= 5.6e-17: tmp = math.acos(-x) else: tmp = math.acos((1.0 - x)) return tmp
function code(x) tmp = 0.0 if (x <= 5.6e-17) tmp = acos(Float64(-x)); else tmp = acos(Float64(1.0 - x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 5.6e-17) tmp = acos(-x); else tmp = acos((1.0 - x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 5.6e-17], N[ArcCos[(-x)], $MachinePrecision], N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.6 \cdot 10^{-17}:\\
\;\;\;\;\cos^{-1} \left(-x\right)\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(1 - x\right)\\
\end{array}
\end{array}
if x < 5.5999999999999998e-17Initial program 3.9%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f646.5
Applied rewrites6.5%
if 5.5999999999999998e-17 < x Initial program 69.4%
(FPCore (x) :precision binary64 (acos (- x)))
double code(double x) {
return acos(-x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = acos(-x)
end function
public static double code(double x) {
return Math.acos(-x);
}
def code(x): return math.acos(-x)
function code(x) return acos(Float64(-x)) end
function tmp = code(x) tmp = acos(-x); end
code[x_] := N[ArcCos[(-x)], $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} \left(-x\right)
\end{array}
Initial program 6.5%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f646.8
Applied rewrites6.8%
(FPCore (x) :precision binary64 (acos 1.0))
double code(double x) {
return acos(1.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = acos(1.0d0)
end function
public static double code(double x) {
return Math.acos(1.0);
}
def code(x): return math.acos(1.0)
function code(x) return acos(1.0) end
function tmp = code(x) tmp = acos(1.0); end
code[x_] := N[ArcCos[1.0], $MachinePrecision]
\begin{array}{l}
\\
\cos^{-1} 1
\end{array}
Initial program 6.5%
Taylor expanded in x around 0
Applied rewrites3.9%
(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 2024278
(FPCore (x)
:name "bug323 (missed optimization)"
:precision binary64
:pre (and (<= 0.0 x) (<= x 0.5))
:alt
(! :herbie-platform default (* 2 (asin (sqrt (/ x 2)))))
(acos (- 1.0 x)))