
(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 (acos (- 1.0 x)))
(t_1 (asin (- 1.0 x)))
(t_2 (fma t_1 2.0 (PI)))
(t_3 (* t_2 t_0))
(t_4 (* (fma (PI) 0.5 t_1) (PI)))
(t_5 (fma t_3 (fma t_0 t_2 t_4) (pow t_4 2.0)))
(t_6 (fma t_1 (- t_1 (* 0.5 (PI))) (* 0.25 (* (PI) (PI))))))
(/
(/
(fma
(* 2.0 (fma 0.125 (pow (PI) 3.0) (pow t_1 3.0)))
t_5
(* (* (- (pow t_4 3.0) (pow t_3 3.0)) (/ -2.0 (PI))) t_6))
(* t_6 t_5))
(* (* (fma 0.5 (PI) t_1) 2.0) (/ 2.0 (PI))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(1 - x\right)\\
t_1 := \sin^{-1} \left(1 - x\right)\\
t_2 := \mathsf{fma}\left(t\_1, 2, \mathsf{PI}\left(\right)\right)\\
t_3 := t\_2 \cdot t\_0\\
t_4 := \mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, t\_1\right) \cdot \mathsf{PI}\left(\right)\\
t_5 := \mathsf{fma}\left(t\_3, \mathsf{fma}\left(t\_0, t\_2, t\_4\right), {t\_4}^{2}\right)\\
t_6 := \mathsf{fma}\left(t\_1, t\_1 - 0.5 \cdot \mathsf{PI}\left(\right), 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)\\
\frac{\frac{\mathsf{fma}\left(2 \cdot \mathsf{fma}\left(0.125, {\mathsf{PI}\left(\right)}^{3}, {t\_1}^{3}\right), t\_5, \left(\left({t\_4}^{3} - {t\_3}^{3}\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right) \cdot t\_6\right)}{t\_6 \cdot t\_5}}{\left(\mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_1\right) \cdot 2\right) \cdot \frac{2}{\mathsf{PI}\left(\right)}}
\end{array}
\end{array}
Initial program 7.1%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites7.1%
Applied rewrites10.6%
Final simplification10.6%
(FPCore (x)
:precision binary64
(let* ((t_0 (asin (- 1.0 x)))
(t_1 (fma t_0 2.0 (PI)))
(t_2 (* (fma (PI) 0.5 t_0) (PI)))
(t_3 (acos (- 1.0 x)))
(t_4 (* t_1 t_3))
(t_5 (fma t_4 (fma t_3 t_1 t_2) (pow t_2 2.0)))
(t_6 (* 0.25 (* (PI) (PI)))))
(/
(/
(fma
(* 2.0 (fma 0.125 (pow (PI) 3.0) (pow t_0 3.0)))
t_5
(*
(* (- (pow t_2 3.0) (pow t_4 3.0)) (/ -2.0 (PI)))
(fma t_0 (- t_0 (* 0.5 (PI))) t_6)))
(* (fma t_0 (- t_3) t_6) t_5))
(* (* (fma 0.5 (PI) t_0) 2.0) (/ 2.0 (PI))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(1 - x\right)\\
t_1 := \mathsf{fma}\left(t\_0, 2, \mathsf{PI}\left(\right)\right)\\
t_2 := \mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, t\_0\right) \cdot \mathsf{PI}\left(\right)\\
t_3 := \cos^{-1} \left(1 - x\right)\\
t_4 := t\_1 \cdot t\_3\\
t_5 := \mathsf{fma}\left(t\_4, \mathsf{fma}\left(t\_3, t\_1, t\_2\right), {t\_2}^{2}\right)\\
t_6 := 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\
\frac{\frac{\mathsf{fma}\left(2 \cdot \mathsf{fma}\left(0.125, {\mathsf{PI}\left(\right)}^{3}, {t\_0}^{3}\right), t\_5, \left(\left({t\_2}^{3} - {t\_4}^{3}\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right) \cdot \mathsf{fma}\left(t\_0, t\_0 - 0.5 \cdot \mathsf{PI}\left(\right), t\_6\right)\right)}{\mathsf{fma}\left(t\_0, -t\_3, t\_6\right) \cdot t\_5}}{\left(\mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_0\right) \cdot 2\right) \cdot \frac{2}{\mathsf{PI}\left(\right)}}
\end{array}
\end{array}
Initial program 7.1%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites7.1%
Applied rewrites10.6%
lift--.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate--r-N/A
lift-*.f64N/A
metadata-evalN/A
div-invN/A
lift-PI.f64N/A
lift-asin.f64N/A
acos-asinN/A
lift-acos.f64N/A
neg-sub0N/A
lift-neg.f6410.6
Applied rewrites10.6%
Final simplification10.6%
(FPCore (x)
:precision binary64
(let* ((t_0 (acos (- 1.0 x)))
(t_1 (- t_0))
(t_2 (* 0.25 (* (PI) (PI))))
(t_3 (asin (- 1.0 x)))
(t_4 (fma 0.5 (PI) t_3))
(t_5 (* t_4 (PI)))
(t_6 (* (fma (PI) 0.5 t_3) (PI)))
(t_7 (* (fma 2.0 t_3 (PI)) t_0))
(t_8 (fma t_3 2.0 (PI))))
(/
(/
(fma
(fma (fma t_4 (PI) t_7) t_7 (pow t_5 2.0))
(* (fma (pow (PI) 3.0) 0.125 (pow t_3 3.0)) 2.0)
(* (- (pow t_5 3.0) (pow t_7 3.0)) (* (fma t_1 t_3 t_2) (/ -2.0 (PI)))))
(* (fma t_3 t_1 t_2) (fma (* t_8 t_0) (fma t_0 t_8 t_6) (pow t_6 2.0))))
(* (* t_4 2.0) (/ 2.0 (PI))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(1 - x\right)\\
t_1 := -t\_0\\
t_2 := 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\\
t_3 := \sin^{-1} \left(1 - x\right)\\
t_4 := \mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_3\right)\\
t_5 := t\_4 \cdot \mathsf{PI}\left(\right)\\
t_6 := \mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, t\_3\right) \cdot \mathsf{PI}\left(\right)\\
t_7 := \mathsf{fma}\left(2, t\_3, \mathsf{PI}\left(\right)\right) \cdot t\_0\\
t_8 := \mathsf{fma}\left(t\_3, 2, \mathsf{PI}\left(\right)\right)\\
\frac{\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(t\_4, \mathsf{PI}\left(\right), t\_7\right), t\_7, {t\_5}^{2}\right), \mathsf{fma}\left({\mathsf{PI}\left(\right)}^{3}, 0.125, {t\_3}^{3}\right) \cdot 2, \left({t\_5}^{3} - {t\_7}^{3}\right) \cdot \left(\mathsf{fma}\left(t\_1, t\_3, t\_2\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right)\right)}{\mathsf{fma}\left(t\_3, t\_1, t\_2\right) \cdot \mathsf{fma}\left(t\_8 \cdot t\_0, \mathsf{fma}\left(t\_0, t\_8, t\_6\right), {t\_6}^{2}\right)}}{\left(t\_4 \cdot 2\right) \cdot \frac{2}{\mathsf{PI}\left(\right)}}
\end{array}
\end{array}
Initial program 7.1%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites7.1%
Applied rewrites10.6%
lift--.f64N/A
sub-negN/A
+-commutativeN/A
neg-sub0N/A
associate--r-N/A
lift-*.f64N/A
metadata-evalN/A
div-invN/A
lift-PI.f64N/A
lift-asin.f64N/A
acos-asinN/A
lift-acos.f64N/A
neg-sub0N/A
lift-neg.f6410.6
Applied rewrites10.6%
Applied rewrites10.6%
Final simplification10.6%
(FPCore (x)
:precision binary64
(let* ((t_0 (asin (- 1.0 x))))
(/
(fma
(pow (PI) 0.25)
(sqrt (pow (PI) 1.5))
(fma
t_0
2.0
(*
(- (* (fma (PI) 0.5 t_0) (PI)) (* (fma t_0 2.0 (PI)) (acos (- 1.0 x))))
(/ -2.0 (PI)))))
(* (* (fma 0.5 (PI) t_0) 2.0) (/ 2.0 (PI))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(1 - x\right)\\
\frac{\mathsf{fma}\left({\mathsf{PI}\left(\right)}^{0.25}, \sqrt{{\mathsf{PI}\left(\right)}^{1.5}}, \mathsf{fma}\left(t\_0, 2, \left(\mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, t\_0\right) \cdot \mathsf{PI}\left(\right) - \mathsf{fma}\left(t\_0, 2, \mathsf{PI}\left(\right)\right) \cdot \cos^{-1} \left(1 - x\right)\right) \cdot \frac{-2}{\mathsf{PI}\left(\right)}\right)\right)}{\left(\mathsf{fma}\left(0.5, \mathsf{PI}\left(\right), t\_0\right) \cdot 2\right) \cdot \frac{2}{\mathsf{PI}\left(\right)}}
\end{array}
\end{array}
Initial program 7.1%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites7.1%
Applied rewrites10.6%
Final simplification10.6%
(FPCore (x) :precision binary64 (let* ((t_0 (asin (- 1.0 x))) (t_1 (/ 2.0 (PI))) (t_2 (pow t_0 2.0))) (/ (fma (- (pow t_0 3.0)) t_1 t_2) (* (+ (* 0.0 t_0) t_2) t_1))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin^{-1} \left(1 - x\right)\\
t_1 := \frac{2}{\mathsf{PI}\left(\right)}\\
t_2 := {t\_0}^{2}\\
\frac{\mathsf{fma}\left(-{t\_0}^{3}, t\_1, t\_2\right)}{\left(0 \cdot t\_0 + t\_2\right) \cdot t\_1}
\end{array}
\end{array}
Initial program 7.1%
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.f647.1
Applied rewrites7.1%
Applied rewrites10.5%
lift-*.f64N/A
lift-+.f64N/A
+-lft-identityN/A
*-rgt-identity10.5
lift-+.f64N/A
lift-*.f64N/A
mul0-lftN/A
+-commutativeN/A
+-lft-identity10.5
Applied rewrites10.5%
Final simplification10.5%
(FPCore (x)
:precision binary64
(let* ((t_0 (acos (- 1.0 x))))
(fma
(fma -0.125 (pow (PI) 3.0) (pow t_0 3.0))
(/ 1.0 (fma t_0 (- t_0 (* -0.5 (PI))) (* 0.25 (* (PI) (PI)))))
(* 0.5 (PI)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos^{-1} \left(1 - x\right)\\
\mathsf{fma}\left(\mathsf{fma}\left(-0.125, {\mathsf{PI}\left(\right)}^{3}, {t\_0}^{3}\right), \frac{1}{\mathsf{fma}\left(t\_0, t\_0 - -0.5 \cdot \mathsf{PI}\left(\right), 0.25 \cdot \left(\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)\right)\right)}, 0.5 \cdot \mathsf{PI}\left(\right)\right)
\end{array}
\end{array}
Initial program 7.1%
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.f647.1
Applied rewrites7.1%
lift-asin.f64N/A
asin-acosN/A
lift-acos.f64N/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
rem-cube-cbrtN/A
lift-cbrt.f64N/A
pow3N/A
unpow2N/A
lift-pow.f64N/A
associate-*r*N/A
lift-*.f64N/A
sub-negN/A
lower-fma.f64N/A
lower-neg.f6410.5
Applied rewrites10.5%
lift-neg.f64N/A
lift-fma.f64N/A
distribute-neg-inN/A
lift-*.f64N/A
associate-*r*N/A
lift-pow.f64N/A
unpow2N/A
pow3N/A
lift-cbrt.f64N/A
rem-cube-cbrtN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
lift-neg.f64N/A
remove-double-negN/A
Applied rewrites7.1%
Applied rewrites10.5%
(FPCore (x)
:precision binary64
(let* ((t_0 (sqrt (PI))))
(fma
(pow (cbrt (* t_0 t_0)) 2.0)
(* (cbrt (PI)) 0.5)
(- (asin (- 1.0 x))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{PI}\left(\right)}\\
\mathsf{fma}\left({\left(\sqrt[3]{t\_0 \cdot t\_0}\right)}^{2}, \sqrt[3]{\mathsf{PI}\left(\right)} \cdot 0.5, -\sin^{-1} \left(1 - x\right)\right)
\end{array}
\end{array}
Initial program 7.1%
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.f645.3
Applied rewrites5.3%
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.f6410.5
Applied rewrites10.5%
(FPCore (x) :precision binary64 (let* ((t_0 (fma (PI) 0.5 (acos (- 1.0 x))))) (fma (* t_0 (PI)) (/ 0.5 t_0) (- (asin (- 1.0 x))))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{PI}\left(\right), 0.5, \cos^{-1} \left(1 - x\right)\right)\\
\mathsf{fma}\left(t\_0 \cdot \mathsf{PI}\left(\right), \frac{0.5}{t\_0}, -\sin^{-1} \left(1 - x\right)\right)
\end{array}
\end{array}
Initial program 7.1%
lift-acos.f64N/A
acos-asinN/A
asin-acosN/A
lift-acos.f64N/A
flip--N/A
frac-subN/A
lower-/.f64N/A
Applied rewrites7.1%
Applied rewrites10.5%
Final simplification10.5%
(FPCore (x) :precision binary64 (if (<= (- 1.0 x) 0.9999999999999999) (* (* (fma -2.0 (/ (fma -2.0 (acos (- 1.0 x)) (PI)) (PI)) 2.0) 0.25) (PI)) (acos (- x))))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;1 - x \leq 0.9999999999999999:\\
\;\;\;\;\left(\mathsf{fma}\left(-2, \frac{\mathsf{fma}\left(-2, \cos^{-1} \left(1 - x\right), \mathsf{PI}\left(\right)\right)}{\mathsf{PI}\left(\right)}, 2\right) \cdot 0.25\right) \cdot \mathsf{PI}\left(\right)\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(-x\right)\\
\end{array}
\end{array}
if (-.f64 #s(literal 1 binary64) x) < 0.999999999999999889Initial program 63.2%
lift-acos.f64N/A
acos-asinN/A
clear-numN/A
asin-acosN/A
acos-asinN/A
flip--N/A
frac-subN/A
frac-subN/A
lower-/.f64N/A
Applied rewrites63.6%
Taylor expanded in x around 0
Applied rewrites63.2%
if 0.999999999999999889 < (-.f64 #s(literal 1 binary64) x) Initial program 3.9%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f646.6
Applied rewrites6.6%
(FPCore (x) :precision binary64 (let* ((t_0 (sqrt (PI)))) (fma (* t_0 -0.5) t_0 (+ (* 0.5 (PI)) (acos (- 1.0 x))))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{\mathsf{PI}\left(\right)}\\
\mathsf{fma}\left(t\_0 \cdot -0.5, t\_0, 0.5 \cdot \mathsf{PI}\left(\right) + \cos^{-1} \left(1 - x\right)\right)
\end{array}
\end{array}
Initial program 7.1%
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.f647.1
Applied rewrites7.1%
lift-asin.f64N/A
asin-acosN/A
lift-acos.f64N/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
rem-cube-cbrtN/A
lift-cbrt.f64N/A
pow3N/A
unpow2N/A
lift-pow.f64N/A
associate-*r*N/A
lift-*.f64N/A
sub-negN/A
lower-fma.f64N/A
lower-neg.f6410.5
Applied rewrites10.5%
lift-neg.f64N/A
lift-fma.f64N/A
distribute-neg-inN/A
lift-*.f64N/A
associate-*r*N/A
lift-pow.f64N/A
unpow2N/A
pow3N/A
lift-cbrt.f64N/A
rem-cube-cbrtN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
lift-neg.f64N/A
remove-double-negN/A
Applied rewrites7.1%
lift-fma.f64N/A
lift-*.f64N/A
lift-fma.f64N/A
associate-+r+N/A
+-commutativeN/A
metadata-evalN/A
div-invN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
associate-+l+N/A
div-invN/A
metadata-evalN/A
*-commutativeN/A
Applied rewrites10.5%
Final simplification10.5%
(FPCore (x) :precision binary64 (let* ((t_0 (sqrt (PI)))) (fma (PI) 0.5 (fma (* t_0 t_0) -0.5 (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 t\_0, -0.5, \cos^{-1} \left(1 - x\right)\right)\right)
\end{array}
\end{array}
Initial program 7.1%
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.f647.1
Applied rewrites7.1%
lift-asin.f64N/A
asin-acosN/A
lift-acos.f64N/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
rem-cube-cbrtN/A
lift-cbrt.f64N/A
pow3N/A
unpow2N/A
lift-pow.f64N/A
associate-*r*N/A
lift-*.f64N/A
sub-negN/A
lower-fma.f64N/A
lower-neg.f6410.5
Applied rewrites10.5%
lift-neg.f64N/A
lift-fma.f64N/A
distribute-neg-inN/A
lift-*.f64N/A
associate-*r*N/A
lift-pow.f64N/A
unpow2N/A
pow3N/A
lift-cbrt.f64N/A
rem-cube-cbrtN/A
distribute-rgt-neg-inN/A
metadata-evalN/A
metadata-evalN/A
lift-neg.f64N/A
remove-double-negN/A
Applied rewrites7.1%
rem-square-sqrtN/A
lift-sqrt.f64N/A
lift-sqrt.f64N/A
lower-*.f6410.5
Applied rewrites10.5%
(FPCore (x) :precision binary64 (if (<= (- 1.0 x) 0.9999999999999999) (acos (- 1.0 x)) (acos (- x))))
double code(double x) {
double tmp;
if ((1.0 - x) <= 0.9999999999999999) {
tmp = acos((1.0 - x));
} else {
tmp = acos(-x);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if ((1.0d0 - x) <= 0.9999999999999999d0) then
tmp = acos((1.0d0 - x))
else
tmp = acos(-x)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if ((1.0 - x) <= 0.9999999999999999) {
tmp = Math.acos((1.0 - x));
} else {
tmp = Math.acos(-x);
}
return tmp;
}
def code(x): tmp = 0 if (1.0 - x) <= 0.9999999999999999: tmp = math.acos((1.0 - x)) else: tmp = math.acos(-x) return tmp
function code(x) tmp = 0.0 if (Float64(1.0 - x) <= 0.9999999999999999) tmp = acos(Float64(1.0 - x)); else tmp = acos(Float64(-x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if ((1.0 - x) <= 0.9999999999999999) tmp = acos((1.0 - x)); else tmp = acos(-x); end tmp_2 = tmp; end
code[x_] := If[LessEqual[N[(1.0 - x), $MachinePrecision], 0.9999999999999999], N[ArcCos[N[(1.0 - x), $MachinePrecision]], $MachinePrecision], N[ArcCos[(-x)], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;1 - x \leq 0.9999999999999999:\\
\;\;\;\;\cos^{-1} \left(1 - x\right)\\
\mathbf{else}:\\
\;\;\;\;\cos^{-1} \left(-x\right)\\
\end{array}
\end{array}
if (-.f64 #s(literal 1 binary64) x) < 0.999999999999999889Initial program 63.2%
if 0.999999999999999889 < (-.f64 #s(literal 1 binary64) x) Initial program 3.9%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f646.6
Applied rewrites6.6%
(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 7.1%
Taylor expanded in x around inf
mul-1-negN/A
lower-neg.f646.9
Applied rewrites6.9%
(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 7.1%
Taylor expanded in x around 0
Applied rewrites3.8%
(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 2024308
(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)))