
(FPCore (x y z t a b) :precision binary64 (* (* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0))) (cos (/ (* (* (+ (* a 2.0) 1.0) b) t) 16.0))))
double code(double x, double y, double z, double t, double a, double b) {
return (x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos((((((a * 2.0) + 1.0) * b) * t) / 16.0));
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (x * cos((((((y * 2.0d0) + 1.0d0) * z) * t) / 16.0d0))) * cos((((((a * 2.0d0) + 1.0d0) * b) * t) / 16.0d0))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * Math.cos((((((a * 2.0) + 1.0) * b) * t) / 16.0));
}
def code(x, y, z, t, a, b): return (x * math.cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * math.cos((((((a * 2.0) + 1.0) * b) * t) / 16.0))
function code(x, y, z, t, a, b) return Float64(Float64(x * cos(Float64(Float64(Float64(Float64(Float64(y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(Float64(Float64(Float64(Float64(Float64(a * 2.0) + 1.0) * b) * t) / 16.0))) end
function tmp = code(x, y, z, t, a, b) tmp = (x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos((((((a * 2.0) + 1.0) * b) * t) / 16.0)); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Cos[N[(N[(N[(N[(N[(y * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(N[(N[(N[(a * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * b), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{\left(\left(a \cdot 2 + 1\right) \cdot b\right) \cdot t}{16}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (* (* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0))) (cos (/ (* (* (+ (* a 2.0) 1.0) b) t) 16.0))))
double code(double x, double y, double z, double t, double a, double b) {
return (x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos((((((a * 2.0) + 1.0) * b) * t) / 16.0));
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (x * cos((((((y * 2.0d0) + 1.0d0) * z) * t) / 16.0d0))) * cos((((((a * 2.0d0) + 1.0d0) * b) * t) / 16.0d0))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * Math.cos((((((a * 2.0) + 1.0) * b) * t) / 16.0));
}
def code(x, y, z, t, a, b): return (x * math.cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * math.cos((((((a * 2.0) + 1.0) * b) * t) / 16.0))
function code(x, y, z, t, a, b) return Float64(Float64(x * cos(Float64(Float64(Float64(Float64(Float64(y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(Float64(Float64(Float64(Float64(Float64(a * 2.0) + 1.0) * b) * t) / 16.0))) end
function tmp = code(x, y, z, t, a, b) tmp = (x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos((((((a * 2.0) + 1.0) * b) * t) / 16.0)); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Cos[N[(N[(N[(N[(N[(y * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(N[(N[(N[(a * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * b), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{\left(\left(a \cdot 2 + 1\right) \cdot b\right) \cdot t}{16}\right)
\end{array}
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (cbrt (* (* z (fma y 2.0 1.0)) (* t 0.0625)))))
(if (<=
(*
(* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0)))
(cos (/ (* t (* (+ 1.0 (* 2.0 a)) b)) 16.0)))
5e+301)
(*
x
(*
(cos (pow (* (pow (pow t_1 2.0) 0.3333333333333333) (cbrt t_1)) 3.0))
(cos (* t (/ (* b (fma 2.0 a 1.0)) 16.0)))))
x)))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = cbrt(((z * fma(y, 2.0, 1.0)) * (t * 0.0625)));
double tmp;
if (((x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(((t * ((1.0 + (2.0 * a)) * b)) / 16.0))) <= 5e+301) {
tmp = x * (cos(pow((pow(pow(t_1, 2.0), 0.3333333333333333) * cbrt(t_1)), 3.0)) * cos((t * ((b * fma(2.0, a, 1.0)) / 16.0))));
} else {
tmp = x;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = cbrt(Float64(Float64(z * fma(y, 2.0, 1.0)) * Float64(t * 0.0625))) tmp = 0.0 if (Float64(Float64(x * cos(Float64(Float64(Float64(Float64(Float64(y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(Float64(Float64(t * Float64(Float64(1.0 + Float64(2.0 * a)) * b)) / 16.0))) <= 5e+301) tmp = Float64(x * Float64(cos((Float64(((t_1 ^ 2.0) ^ 0.3333333333333333) * cbrt(t_1)) ^ 3.0)) * cos(Float64(t * Float64(Float64(b * fma(2.0, a, 1.0)) / 16.0))))); else tmp = x; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[Power[N[(N[(z * N[(y * 2.0 + 1.0), $MachinePrecision]), $MachinePrecision] * N[(t * 0.0625), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, If[LessEqual[N[(N[(x * N[Cos[N[(N[(N[(N[(N[(y * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(t * N[(N[(1.0 + N[(2.0 * a), $MachinePrecision]), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 5e+301], N[(x * N[(N[Cos[N[Power[N[(N[Power[N[Power[t$95$1, 2.0], $MachinePrecision], 0.3333333333333333], $MachinePrecision] * N[Power[t$95$1, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]], $MachinePrecision] * N[Cos[N[(t * N[(N[(b * N[(2.0 * a + 1.0), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sqrt[3]{\left(z \cdot \mathsf{fma}\left(y, 2, 1\right)\right) \cdot \left(t \cdot 0.0625\right)}\\
\mathbf{if}\;\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{t \cdot \left(\left(1 + 2 \cdot a\right) \cdot b\right)}{16}\right) \leq 5 \cdot 10^{+301}:\\
\;\;\;\;x \cdot \left(\cos \left({\left({\left({t\_1}^{2}\right)}^{0.3333333333333333} \cdot \sqrt[3]{t\_1}\right)}^{3}\right) \cdot \cos \left(t \cdot \frac{b \cdot \mathsf{fma}\left(2, a, 1\right)}{16}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) < 5.0000000000000004e301Initial program 43.3%
associate-*l*43.3%
*-commutative43.3%
associate-*r*42.7%
*-commutative42.7%
associate-*l*43.5%
associate-/l*43.5%
fma-define43.5%
*-commutative43.5%
associate-*r/43.5%
*-commutative43.5%
associate-*r*42.7%
*-commutative42.7%
Simplified43.5%
add-cube-cbrt42.7%
pow342.5%
associate-*r*42.7%
*-commutative42.7%
associate-*l*42.0%
div-inv42.0%
metadata-eval42.0%
Applied egg-rr42.0%
pow1/327.3%
add-cube-cbrt27.7%
unpow-prod-down27.1%
pow227.1%
associate-*r*27.3%
pow1/343.9%
associate-*r*44.7%
Applied egg-rr44.7%
if 5.0000000000000004e301 < (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) Initial program 0.0%
associate-*l*0.0%
*-commutative0.0%
associate-*r*1.0%
*-commutative1.0%
associate-*l*1.1%
associate-/l*1.3%
fma-define1.3%
*-commutative1.3%
associate-*r/1.3%
*-commutative1.3%
associate-*r*2.1%
*-commutative2.1%
Simplified1.3%
Taylor expanded in z around 0 4.8%
Taylor expanded in t around 0 11.1%
Final simplification30.1%
(FPCore (x y z t a b)
:precision binary64
(if (<=
(*
(* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0)))
(cos (/ (* t (* (+ 1.0 (* 2.0 a)) b)) 16.0)))
5e+301)
(*
x
(*
(log1p (expm1 (cos (* (* z (fma y 2.0 1.0)) (* t 0.0625)))))
(cos (pow (cbrt (* t (* (fma 2.0 a 1.0) (* b 0.0625)))) 3.0))))
x))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (((x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(((t * ((1.0 + (2.0 * a)) * b)) / 16.0))) <= 5e+301) {
tmp = x * (log1p(expm1(cos(((z * fma(y, 2.0, 1.0)) * (t * 0.0625))))) * cos(pow(cbrt((t * (fma(2.0, a, 1.0) * (b * 0.0625)))), 3.0)));
} else {
tmp = x;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(x * cos(Float64(Float64(Float64(Float64(Float64(y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(Float64(Float64(t * Float64(Float64(1.0 + Float64(2.0 * a)) * b)) / 16.0))) <= 5e+301) tmp = Float64(x * Float64(log1p(expm1(cos(Float64(Float64(z * fma(y, 2.0, 1.0)) * Float64(t * 0.0625))))) * cos((cbrt(Float64(t * Float64(fma(2.0, a, 1.0) * Float64(b * 0.0625)))) ^ 3.0)))); else tmp = x; end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(x * N[Cos[N[(N[(N[(N[(N[(y * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(t * N[(N[(1.0 + N[(2.0 * a), $MachinePrecision]), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 5e+301], N[(x * N[(N[Log[1 + N[(Exp[N[Cos[N[(N[(z * N[(y * 2.0 + 1.0), $MachinePrecision]), $MachinePrecision] * N[(t * 0.0625), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]], $MachinePrecision] * N[Cos[N[Power[N[Power[N[(t * N[(N[(2.0 * a + 1.0), $MachinePrecision] * N[(b * 0.0625), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{t \cdot \left(\left(1 + 2 \cdot a\right) \cdot b\right)}{16}\right) \leq 5 \cdot 10^{+301}:\\
\;\;\;\;x \cdot \left(\mathsf{log1p}\left(\mathsf{expm1}\left(\cos \left(\left(z \cdot \mathsf{fma}\left(y, 2, 1\right)\right) \cdot \left(t \cdot 0.0625\right)\right)\right)\right) \cdot \cos \left({\left(\sqrt[3]{t \cdot \left(\mathsf{fma}\left(2, a, 1\right) \cdot \left(b \cdot 0.0625\right)\right)}\right)}^{3}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) < 5.0000000000000004e301Initial program 43.3%
associate-*l*43.3%
*-commutative43.3%
associate-*r*42.7%
*-commutative42.7%
associate-*l*43.5%
associate-/l*43.5%
fma-define43.5%
*-commutative43.5%
associate-*r/43.5%
*-commutative43.5%
associate-*r*42.7%
*-commutative42.7%
Simplified43.5%
add-cube-cbrt42.7%
pow342.5%
associate-*r*42.7%
*-commutative42.7%
associate-*l*42.0%
div-inv42.0%
metadata-eval42.0%
Applied egg-rr42.0%
log1p-expm1-u42.0%
rem-cube-cbrt42.7%
associate-*r*43.3%
Applied egg-rr43.3%
associate-*r/43.5%
fma-undefine43.5%
*-commutative43.5%
*-commutative43.5%
add-cube-cbrt43.5%
pow343.9%
*-commutative43.9%
*-commutative43.9%
fma-undefine43.9%
associate-*r/43.9%
associate-/l*43.9%
div-inv43.9%
metadata-eval43.9%
Applied egg-rr43.9%
if 5.0000000000000004e301 < (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) Initial program 0.0%
associate-*l*0.0%
*-commutative0.0%
associate-*r*1.0%
*-commutative1.0%
associate-*l*1.1%
associate-/l*1.3%
fma-define1.3%
*-commutative1.3%
associate-*r/1.3%
*-commutative1.3%
associate-*r*2.1%
*-commutative2.1%
Simplified1.3%
Taylor expanded in z around 0 4.8%
Taylor expanded in t around 0 11.1%
Final simplification29.7%
(FPCore (x y z t a b)
:precision binary64
(if (<=
(*
(* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0)))
(cos (/ (* t (* (+ 1.0 (* 2.0 a)) b)) 16.0)))
5e+301)
(*
x
(*
(cos (pow (cbrt (* t (* (fma 2.0 a 1.0) (* b 0.0625)))) 3.0))
(cos (* (fma y 2.0 1.0) (* z (/ t 16.0))))))
x))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (((x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(((t * ((1.0 + (2.0 * a)) * b)) / 16.0))) <= 5e+301) {
tmp = x * (cos(pow(cbrt((t * (fma(2.0, a, 1.0) * (b * 0.0625)))), 3.0)) * cos((fma(y, 2.0, 1.0) * (z * (t / 16.0)))));
} else {
tmp = x;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(x * cos(Float64(Float64(Float64(Float64(Float64(y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(Float64(Float64(t * Float64(Float64(1.0 + Float64(2.0 * a)) * b)) / 16.0))) <= 5e+301) tmp = Float64(x * Float64(cos((cbrt(Float64(t * Float64(fma(2.0, a, 1.0) * Float64(b * 0.0625)))) ^ 3.0)) * cos(Float64(fma(y, 2.0, 1.0) * Float64(z * Float64(t / 16.0)))))); else tmp = x; end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(x * N[Cos[N[(N[(N[(N[(N[(y * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(t * N[(N[(1.0 + N[(2.0 * a), $MachinePrecision]), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 5e+301], N[(x * N[(N[Cos[N[Power[N[Power[N[(t * N[(N[(2.0 * a + 1.0), $MachinePrecision] * N[(b * 0.0625), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision]], $MachinePrecision] * N[Cos[N[(N[(y * 2.0 + 1.0), $MachinePrecision] * N[(z * N[(t / 16.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{t \cdot \left(\left(1 + 2 \cdot a\right) \cdot b\right)}{16}\right) \leq 5 \cdot 10^{+301}:\\
\;\;\;\;x \cdot \left(\cos \left({\left(\sqrt[3]{t \cdot \left(\mathsf{fma}\left(2, a, 1\right) \cdot \left(b \cdot 0.0625\right)\right)}\right)}^{3}\right) \cdot \cos \left(\mathsf{fma}\left(y, 2, 1\right) \cdot \left(z \cdot \frac{t}{16}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) < 5.0000000000000004e301Initial program 43.3%
associate-*l*43.3%
*-commutative43.3%
associate-*r*42.7%
*-commutative42.7%
associate-*l*43.5%
associate-/l*43.5%
fma-define43.5%
*-commutative43.5%
associate-*r/43.5%
*-commutative43.5%
associate-*r*42.7%
*-commutative42.7%
Simplified43.5%
associate-*r/43.5%
fma-undefine43.5%
*-commutative43.5%
*-commutative43.5%
add-cube-cbrt43.5%
pow343.9%
*-commutative43.9%
*-commutative43.9%
fma-undefine43.9%
associate-*r/43.9%
associate-/l*43.9%
div-inv43.9%
metadata-eval43.9%
Applied egg-rr43.9%
if 5.0000000000000004e301 < (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) Initial program 0.0%
associate-*l*0.0%
*-commutative0.0%
associate-*r*1.0%
*-commutative1.0%
associate-*l*1.1%
associate-/l*1.3%
fma-define1.3%
*-commutative1.3%
associate-*r/1.3%
*-commutative1.3%
associate-*r*2.1%
*-commutative2.1%
Simplified1.3%
Taylor expanded in z around 0 4.8%
Taylor expanded in t around 0 11.1%
Final simplification29.6%
(FPCore (x y z t a b)
:precision binary64
(if (<=
(*
(* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0)))
(cos (/ (* t (* (+ 1.0 (* 2.0 a)) b)) 16.0)))
5e+301)
(*
x
(*
(cos (* t (/ (* b (fma 2.0 a 1.0)) 16.0)))
(cos (* (fma y 2.0 1.0) (* z (/ t 16.0))))))
x))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (((x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(((t * ((1.0 + (2.0 * a)) * b)) / 16.0))) <= 5e+301) {
tmp = x * (cos((t * ((b * fma(2.0, a, 1.0)) / 16.0))) * cos((fma(y, 2.0, 1.0) * (z * (t / 16.0)))));
} else {
tmp = x;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(x * cos(Float64(Float64(Float64(Float64(Float64(y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(Float64(Float64(t * Float64(Float64(1.0 + Float64(2.0 * a)) * b)) / 16.0))) <= 5e+301) tmp = Float64(x * Float64(cos(Float64(t * Float64(Float64(b * fma(2.0, a, 1.0)) / 16.0))) * cos(Float64(fma(y, 2.0, 1.0) * Float64(z * Float64(t / 16.0)))))); else tmp = x; end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(x * N[Cos[N[(N[(N[(N[(N[(y * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(t * N[(N[(1.0 + N[(2.0 * a), $MachinePrecision]), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 5e+301], N[(x * N[(N[Cos[N[(t * N[(N[(b * N[(2.0 * a + 1.0), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Cos[N[(N[(y * 2.0 + 1.0), $MachinePrecision] * N[(z * N[(t / 16.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{t \cdot \left(\left(1 + 2 \cdot a\right) \cdot b\right)}{16}\right) \leq 5 \cdot 10^{+301}:\\
\;\;\;\;x \cdot \left(\cos \left(t \cdot \frac{b \cdot \mathsf{fma}\left(2, a, 1\right)}{16}\right) \cdot \cos \left(\mathsf{fma}\left(y, 2, 1\right) \cdot \left(z \cdot \frac{t}{16}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) < 5.0000000000000004e301Initial program 43.3%
associate-*l*43.3%
*-commutative43.3%
associate-*r*42.7%
*-commutative42.7%
associate-*l*43.5%
associate-/l*43.5%
fma-define43.5%
*-commutative43.5%
associate-*r/43.5%
*-commutative43.5%
associate-*r*42.7%
*-commutative42.7%
Simplified43.5%
if 5.0000000000000004e301 < (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) Initial program 0.0%
associate-*l*0.0%
*-commutative0.0%
associate-*r*1.0%
*-commutative1.0%
associate-*l*1.1%
associate-/l*1.3%
fma-define1.3%
*-commutative1.3%
associate-*r/1.3%
*-commutative1.3%
associate-*r*2.1%
*-commutative2.1%
Simplified1.3%
Taylor expanded in z around 0 4.8%
Taylor expanded in t around 0 11.1%
Final simplification29.4%
(FPCore (x y z t a b)
:precision binary64
(if (<=
(*
(* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0)))
(cos (/ (* t (* (+ 1.0 (* 2.0 a)) b)) 16.0)))
5e+230)
(*
x
(*
(cos (* t (/ (* b (fma 2.0 a 1.0)) 16.0)))
(cos (/ (* z (fma y 2.0 1.0)) (/ 16.0 t)))))
x))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (((x * cos((((((y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(((t * ((1.0 + (2.0 * a)) * b)) / 16.0))) <= 5e+230) {
tmp = x * (cos((t * ((b * fma(2.0, a, 1.0)) / 16.0))) * cos(((z * fma(y, 2.0, 1.0)) / (16.0 / t))));
} else {
tmp = x;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(x * cos(Float64(Float64(Float64(Float64(Float64(y * 2.0) + 1.0) * z) * t) / 16.0))) * cos(Float64(Float64(t * Float64(Float64(1.0 + Float64(2.0 * a)) * b)) / 16.0))) <= 5e+230) tmp = Float64(x * Float64(cos(Float64(t * Float64(Float64(b * fma(2.0, a, 1.0)) / 16.0))) * cos(Float64(Float64(z * fma(y, 2.0, 1.0)) / Float64(16.0 / t))))); else tmp = x; end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(x * N[Cos[N[(N[(N[(N[(N[(y * 2.0), $MachinePrecision] + 1.0), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(t * N[(N[(1.0 + N[(2.0 * a), $MachinePrecision]), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 5e+230], N[(x * N[(N[Cos[N[(t * N[(N[(b * N[(2.0 * a + 1.0), $MachinePrecision]), $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Cos[N[(N[(z * N[(y * 2.0 + 1.0), $MachinePrecision]), $MachinePrecision] / N[(16.0 / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(x \cdot \cos \left(\frac{\left(\left(y \cdot 2 + 1\right) \cdot z\right) \cdot t}{16}\right)\right) \cdot \cos \left(\frac{t \cdot \left(\left(1 + 2 \cdot a\right) \cdot b\right)}{16}\right) \leq 5 \cdot 10^{+230}:\\
\;\;\;\;x \cdot \left(\cos \left(t \cdot \frac{b \cdot \mathsf{fma}\left(2, a, 1\right)}{16}\right) \cdot \cos \left(\frac{z \cdot \mathsf{fma}\left(y, 2, 1\right)}{\frac{16}{t}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) < 5.0000000000000003e230Initial program 43.5%
associate-*l*43.5%
*-commutative43.5%
associate-*r*42.9%
*-commutative42.9%
associate-*l*43.8%
associate-/l*43.8%
fma-define43.8%
*-commutative43.8%
associate-*r/43.8%
*-commutative43.8%
associate-*r*42.8%
*-commutative42.8%
Simplified43.8%
associate-*r*43.5%
clear-num43.4%
un-div-inv43.8%
Applied egg-rr43.8%
if 5.0000000000000003e230 < (*.f64 (*.f64 x (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 y 2) 1) z) t) 16))) (cos.f64 (/.f64 (*.f64 (*.f64 (+.f64 (*.f64 a 2) 1) b) t) 16))) Initial program 2.3%
associate-*l*2.3%
*-commutative2.3%
associate-*r*3.2%
*-commutative3.2%
associate-*l*3.3%
associate-/l*3.5%
fma-define3.5%
*-commutative3.5%
associate-*r/3.5%
*-commutative3.5%
associate-*r*4.4%
*-commutative4.4%
Simplified3.5%
Taylor expanded in z around 0 6.7%
Taylor expanded in t around 0 12.3%
Final simplification29.3%
(FPCore (x y z t a b) :precision binary64 x)
double code(double x, double y, double z, double t, double a, double b) {
return x;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x;
}
def code(x, y, z, t, a, b): return x
function code(x, y, z, t, a, b) return x end
function tmp = code(x, y, z, t, a, b) tmp = x; end
code[x_, y_, z_, t_, a_, b_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 24.5%
associate-*l*24.5%
*-commutative24.5%
associate-*r*24.6%
*-commutative24.6%
associate-*l*25.1%
associate-/l*25.2%
fma-define25.2%
*-commutative25.2%
associate-*r/25.2%
*-commutative25.2%
associate-*r*25.1%
*-commutative25.1%
Simplified25.2%
Taylor expanded in z around 0 26.1%
Taylor expanded in t around 0 27.9%
Final simplification27.9%
(FPCore (x y z t a b) :precision binary64 (* x (cos (* (/ b 16.0) (/ t (+ (- 1.0 (* a 2.0)) (pow (* a 2.0) 2.0)))))))
double code(double x, double y, double z, double t, double a, double b) {
return x * cos(((b / 16.0) * (t / ((1.0 - (a * 2.0)) + pow((a * 2.0), 2.0)))));
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x * cos(((b / 16.0d0) * (t / ((1.0d0 - (a * 2.0d0)) + ((a * 2.0d0) ** 2.0d0)))))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x * Math.cos(((b / 16.0) * (t / ((1.0 - (a * 2.0)) + Math.pow((a * 2.0), 2.0)))));
}
def code(x, y, z, t, a, b): return x * math.cos(((b / 16.0) * (t / ((1.0 - (a * 2.0)) + math.pow((a * 2.0), 2.0)))))
function code(x, y, z, t, a, b) return Float64(x * cos(Float64(Float64(b / 16.0) * Float64(t / Float64(Float64(1.0 - Float64(a * 2.0)) + (Float64(a * 2.0) ^ 2.0)))))) end
function tmp = code(x, y, z, t, a, b) tmp = x * cos(((b / 16.0) * (t / ((1.0 - (a * 2.0)) + ((a * 2.0) ^ 2.0))))); end
code[x_, y_, z_, t_, a_, b_] := N[(x * N[Cos[N[(N[(b / 16.0), $MachinePrecision] * N[(t / N[(N[(1.0 - N[(a * 2.0), $MachinePrecision]), $MachinePrecision] + N[Power[N[(a * 2.0), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \cos \left(\frac{b}{16} \cdot \frac{t}{\left(1 - a \cdot 2\right) + {\left(a \cdot 2\right)}^{2}}\right)
\end{array}
herbie shell --seed 2024046
(FPCore (x y z t a b)
:name "Codec.Picture.Jpg.FastDct:referenceDct from JuicyPixels-3.2.6.1"
:precision binary64
:alt
(* x (cos (* (/ b 16.0) (/ t (+ (- 1.0 (* a 2.0)) (pow (* a 2.0) 2.0))))))
(* (* x (cos (/ (* (* (+ (* y 2.0) 1.0) z) t) 16.0))) (cos (/ (* (* (+ (* a 2.0) 1.0) b) t) 16.0))))