Codec.Picture.Jpg.FastDct:referenceDct from JuicyPixels-3.2.6.1

Percentage Accurate: 27.2% → 30.6%

Time: 33.9s

Alternatives: 3

Speedup: 225.0×

Specification

Math

\[\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

(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))))

C

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));
}

Fortran

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

Java

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));
}

Python

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))

Julia

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

MATLAB

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

Wolfram

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]

Initial Program: 27.2% accurate, 1.0× speedup

Math

\[\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

(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))))

C

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));
}

Fortran

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

Java

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));
}

Python

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))

Julia

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

MATLAB

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

Wolfram

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]

Alternative 1: 30.6% accurate, 0.2× speedup

Math

\[\begin{array}{l} z_m = \left|z\right| \\ t_m = \left|t\right| \\ \begin{array}{l} t_1 := \log \left(\frac{\frac{16}{z_m}}{t_m}\right)\\ x \cdot \cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{{t_1}^{2}} \cdot \sqrt[3]{t_1}\right|}}\right) \end{array} \end{array} \]

FPCore

z_m = (fabs.f64 z)
t_m = (fabs.f64 t)
(FPCore (x y z_m t_m a b)
 :precision binary64
 (let* ((t_1 (log (/ (/ 16.0 z_m) t_m))))
   (*
    x
    (cos
     (/ (fma y 2.0 1.0) (exp (fabs (* (cbrt (pow t_1 2.0)) (cbrt t_1)))))))))

C

z_m = fabs(z);
t_m = fabs(t);
double code(double x, double y, double z_m, double t_m, double a, double b) {
	double t_1 = log(((16.0 / z_m) / t_m));
	return x * cos((fma(y, 2.0, 1.0) / exp(fabs((cbrt(pow(t_1, 2.0)) * cbrt(t_1))))));
}

Julia

z_m = abs(z)
t_m = abs(t)
function code(x, y, z_m, t_m, a, b)
	t_1 = log(Float64(Float64(16.0 / z_m) / t_m))
	return Float64(x * cos(Float64(fma(y, 2.0, 1.0) / exp(abs(Float64(cbrt((t_1 ^ 2.0)) * cbrt(t_1)))))))
end

Wolfram

z_m = N[Abs[z], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
code[x_, y_, z$95$m_, t$95$m_, a_, b_] := Block[{t$95$1 = N[Log[N[(N[(16.0 / z$95$m), $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]}, N[(x * N[Cos[N[(N[(y * 2.0 + 1.0), $MachinePrecision] / N[Exp[N[Abs[N[(N[Power[N[Power[t$95$1, 2.0], $MachinePrecision], 1/3], $MachinePrecision] * N[Power[t$95$1, 1/3], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]

Derivation

1. Initial program 27.8%

   \[\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) \]

3. Simplified 29.0%

   \[\leadsto \color{blue}{x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\frac{\frac{16}{t}}{z}}\right) \cdot \cos \left(\frac{t \cdot \mathsf{fma}\left(2, a, 1\right)}{16} \cdot b\right)\right)} \]
4. Add Preprocessing

5. Taylor expanded in t around 0 30.6%

   \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\frac{\frac{16}{t}}{z}}\right) \cdot \color{blue}{1}\right) \]
6. Step-by-step derivation

   1. add-exp-log 14.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\color{blue}{e^{\log \left(\frac{\frac{16}{t}}{z}\right)}}}\right) \cdot 1\right) \]

   2. associate-/l/ 14.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\log \color{blue}{\left(\frac{16}{z \cdot t}\right)}}}\right) \cdot 1\right) \]

   3. *-commutative 14.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\log \left(\frac{16}{\color{blue}{t \cdot z}}\right)}}\right) \cdot 1\right) \]

7. Applied egg-rr 14.2%

   \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\color{blue}{e^{\log \left(\frac{16}{t \cdot z}\right)}}}\right) \cdot 1\right) \]
8. Step-by-step derivation

   1. add-sqr-sqrt 12.8%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\sqrt{\log \left(\frac{16}{t \cdot z}\right)} \cdot \sqrt{\log \left(\frac{16}{t \cdot z}\right)}}}}\right) \cdot 1\right) \]

   2. sqrt-unprod 16.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\sqrt{\log \left(\frac{16}{t \cdot z}\right) \cdot \log \left(\frac{16}{t \cdot z}\right)}}}}\right) \cdot 1\right) \]

   3. pow2 16.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\sqrt{\color{blue}{{\log \left(\frac{16}{t \cdot z}\right)}^{2}}}}}\right) \cdot 1\right) \]

   4. associate-/r* 16.1%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\sqrt{{\log \color{blue}{\left(\frac{\frac{16}{t}}{z}\right)}}^{2}}}}\right) \cdot 1\right) \]

9. Applied egg-rr 16.1%

   \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\sqrt{{\log \left(\frac{\frac{16}{t}}{z}\right)}^{2}}}}}\right) \cdot 1\right) \]
10. Step-by-step derivation

    1. unpow2 16.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\sqrt{\color{blue}{\log \left(\frac{\frac{16}{t}}{z}\right) \cdot \log \left(\frac{\frac{16}{t}}{z}\right)}}}}\right) \cdot 1\right) \]

    2. rem-sqrt-square 16.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\left|\log \left(\frac{\frac{16}{t}}{z}\right)\right|}}}\right) \cdot 1\right) \]

11. Simplified 16.1%

    \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\left|\log \left(\frac{\frac{16}{t}}{z}\right)\right|}}}\right) \cdot 1\right) \]
12. Step-by-step derivation

    1. diff-log 7.0%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\color{blue}{\log \left(\frac{16}{t}\right) - \log z}\right|}}\right) \cdot 1\right) \]

    2. add-sqr-sqrt 5.7%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\color{blue}{\sqrt{\log \left(\frac{16}{t}\right) - \log z} \cdot \sqrt{\log \left(\frac{16}{t}\right) - \log z}}\right|}}\right) \cdot 1\right) \]

    3. fabs-sqr 5.7%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\color{blue}{\left|\sqrt{\log \left(\frac{16}{t}\right) - \log z} \cdot \sqrt{\log \left(\frac{16}{t}\right) - \log z}\right|}\right|}}\right) \cdot 1\right) \]

    4. add-sqr-sqrt 7.0%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\left|\color{blue}{\log \left(\frac{16}{t}\right) - \log z}\right|\right|}}\right) \cdot 1\right) \]

    5. add-cbrt-cube 7.0%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\color{blue}{\sqrt[3]{\left(\left|\log \left(\frac{16}{t}\right) - \log z\right| \cdot \left|\log \left(\frac{16}{t}\right) - \log z\right|\right) \cdot \left|\log \left(\frac{16}{t}\right) - \log z\right|}}\right|}}\right) \cdot 1\right) \]

    6. cbrt-prod 7.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\color{blue}{\sqrt[3]{\left|\log \left(\frac{16}{t}\right) - \log z\right| \cdot \left|\log \left(\frac{16}{t}\right) - \log z\right|} \cdot \sqrt[3]{\left|\log \left(\frac{16}{t}\right) - \log z\right|}}\right|}}\right) \cdot 1\right) \]

    7. sqr-abs 7.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{\color{blue}{\left(\log \left(\frac{16}{t}\right) - \log z\right) \cdot \left(\log \left(\frac{16}{t}\right) - \log z\right)}} \cdot \sqrt[3]{\left|\log \left(\frac{16}{t}\right) - \log z\right|}\right|}}\right) \cdot 1\right) \]

    8. pow2 7.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{\color{blue}{{\left(\log \left(\frac{16}{t}\right) - \log z\right)}^{2}}} \cdot \sqrt[3]{\left|\log \left(\frac{16}{t}\right) - \log z\right|}\right|}}\right) \cdot 1\right) \]

    9. diff-log 7.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{{\color{blue}{\log \left(\frac{\frac{16}{t}}{z}\right)}}^{2}} \cdot \sqrt[3]{\left|\log \left(\frac{16}{t}\right) - \log z\right|}\right|}}\right) \cdot 1\right) \]

    10. associate-/r* 7.1%

        \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{{\log \color{blue}{\left(\frac{16}{t \cdot z}\right)}}^{2}} \cdot \sqrt[3]{\left|\log \left(\frac{16}{t}\right) - \log z\right|}\right|}}\right) \cdot 1\right) \]

    11. add-sqr-sqrt 5.6%

        \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{{\log \left(\frac{16}{t \cdot z}\right)}^{2}} \cdot \sqrt[3]{\left|\color{blue}{\sqrt{\log \left(\frac{16}{t}\right) - \log z} \cdot \sqrt{\log \left(\frac{16}{t}\right) - \log z}}\right|}\right|}}\right) \cdot 1\right) \]

13. Applied egg-rr 16.5%

    \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\color{blue}{\sqrt[3]{{\log \left(\frac{16}{t \cdot z}\right)}^{2}} \cdot \sqrt[3]{\log \left(\frac{16}{t \cdot z}\right)}}\right|}}\right) \cdot 1\right) \]
14. Step-by-step derivation

    1. associate-/l/ 16.4%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{{\log \color{blue}{\left(\frac{\frac{16}{z}}{t}\right)}}^{2}} \cdot \sqrt[3]{\log \left(\frac{16}{t \cdot z}\right)}\right|}}\right) \cdot 1\right) \]

    2. associate-/l/ 16.4%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{{\log \left(\frac{\frac{16}{z}}{t}\right)}^{2}} \cdot \sqrt[3]{\log \color{blue}{\left(\frac{\frac{16}{z}}{t}\right)}}\right|}}\right) \cdot 1\right) \]

15. Simplified 16.4%

    \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\color{blue}{\sqrt[3]{{\log \left(\frac{\frac{16}{z}}{t}\right)}^{2}} \cdot \sqrt[3]{\log \left(\frac{\frac{16}{z}}{t}\right)}}\right|}}\right) \cdot 1\right) \]

16. Final simplification 16.4%

    \[\leadsto x \cdot \cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\left|\sqrt[3]{{\log \left(\frac{\frac{16}{z}}{t}\right)}^{2}} \cdot \sqrt[3]{\log \left(\frac{\frac{16}{z}}{t}\right)}\right|}}\right) \]
17. Add Preprocessing

Alternative 2: 30.6% accurate, 0.5× speedup

Math

\[\begin{array}{l} z_m = \left|z\right| \\ t_m = \left|t\right| \\ x \cdot \cos \left(\frac{y \cdot 2}{e^{\left|\log \left(\frac{16}{z_m \cdot t_m}\right)\right|}}\right) \end{array} \]

FPCore

z_m = (fabs.f64 z)
t_m = (fabs.f64 t)
(FPCore (x y z_m t_m a b)
 :precision binary64
 (* x (cos (/ (* y 2.0) (exp (fabs (log (/ 16.0 (* z_m t_m)))))))))

C

z_m = fabs(z);
t_m = fabs(t);
double code(double x, double y, double z_m, double t_m, double a, double b) {
	return x * cos(((y * 2.0) / exp(fabs(log((16.0 / (z_m * t_m)))))));
}

Fortran

z_m = abs(z)
t_m = abs(t)
real(8) function code(x, y, z_m, t_m, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z_m
    real(8), intent (in) :: t_m
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = x * cos(((y * 2.0d0) / exp(abs(log((16.0d0 / (z_m * t_m)))))))
end function

Java

z_m = Math.abs(z);
t_m = Math.abs(t);
public static double code(double x, double y, double z_m, double t_m, double a, double b) {
	return x * Math.cos(((y * 2.0) / Math.exp(Math.abs(Math.log((16.0 / (z_m * t_m)))))));
}

Python

z_m = math.fabs(z)
t_m = math.fabs(t)
def code(x, y, z_m, t_m, a, b):
	return x * math.cos(((y * 2.0) / math.exp(math.fabs(math.log((16.0 / (z_m * t_m)))))))

Julia

z_m = abs(z)
t_m = abs(t)
function code(x, y, z_m, t_m, a, b)
	return Float64(x * cos(Float64(Float64(y * 2.0) / exp(abs(log(Float64(16.0 / Float64(z_m * t_m))))))))
end

MATLAB

z_m = abs(z);
t_m = abs(t);
function tmp = code(x, y, z_m, t_m, a, b)
	tmp = x * cos(((y * 2.0) / exp(abs(log((16.0 / (z_m * t_m)))))));
end

Wolfram

z_m = N[Abs[z], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
code[x_, y_, z$95$m_, t$95$m_, a_, b_] := N[(x * N[Cos[N[(N[(y * 2.0), $MachinePrecision] / N[Exp[N[Abs[N[Log[N[(16.0 / N[(z$95$m * t$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 27.8%

   \[\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) \]

3. Simplified 29.0%

   \[\leadsto \color{blue}{x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\frac{\frac{16}{t}}{z}}\right) \cdot \cos \left(\frac{t \cdot \mathsf{fma}\left(2, a, 1\right)}{16} \cdot b\right)\right)} \]
4. Add Preprocessing

5. Taylor expanded in t around 0 30.6%

   \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\frac{\frac{16}{t}}{z}}\right) \cdot \color{blue}{1}\right) \]
6. Step-by-step derivation

   1. add-exp-log 14.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\color{blue}{e^{\log \left(\frac{\frac{16}{t}}{z}\right)}}}\right) \cdot 1\right) \]

   2. associate-/l/ 14.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\log \color{blue}{\left(\frac{16}{z \cdot t}\right)}}}\right) \cdot 1\right) \]

   3. *-commutative 14.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\log \left(\frac{16}{\color{blue}{t \cdot z}}\right)}}\right) \cdot 1\right) \]

7. Applied egg-rr 14.2%

   \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{\color{blue}{e^{\log \left(\frac{16}{t \cdot z}\right)}}}\right) \cdot 1\right) \]
8. Step-by-step derivation

   1. add-sqr-sqrt 12.8%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\sqrt{\log \left(\frac{16}{t \cdot z}\right)} \cdot \sqrt{\log \left(\frac{16}{t \cdot z}\right)}}}}\right) \cdot 1\right) \]

   2. sqrt-unprod 16.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\sqrt{\log \left(\frac{16}{t \cdot z}\right) \cdot \log \left(\frac{16}{t \cdot z}\right)}}}}\right) \cdot 1\right) \]

   3. pow2 16.2%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\sqrt{\color{blue}{{\log \left(\frac{16}{t \cdot z}\right)}^{2}}}}}\right) \cdot 1\right) \]

   4. associate-/r* 16.1%

      \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\sqrt{{\log \color{blue}{\left(\frac{\frac{16}{t}}{z}\right)}}^{2}}}}\right) \cdot 1\right) \]

9. Applied egg-rr 16.1%

   \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\sqrt{{\log \left(\frac{\frac{16}{t}}{z}\right)}^{2}}}}}\right) \cdot 1\right) \]
10. Step-by-step derivation

    1. unpow2 16.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\sqrt{\color{blue}{\log \left(\frac{\frac{16}{t}}{z}\right) \cdot \log \left(\frac{\frac{16}{t}}{z}\right)}}}}\right) \cdot 1\right) \]

    2. rem-sqrt-square 16.1%

       \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\left|\log \left(\frac{\frac{16}{t}}{z}\right)\right|}}}\right) \cdot 1\right) \]

11. Simplified 16.1%

    \[\leadsto x \cdot \left(\cos \left(\frac{\mathsf{fma}\left(y, 2, 1\right)}{e^{\color{blue}{\left|\log \left(\frac{\frac{16}{t}}{z}\right)\right|}}}\right) \cdot 1\right) \]

12. Taylor expanded in y around inf 16.2%

    \[\leadsto x \cdot \left(\cos \color{blue}{\left(2 \cdot \frac{y}{e^{\left|\log \left(\frac{16}{t \cdot z}\right)\right|}}\right)} \cdot 1\right) \]
13. Step-by-step derivation

    1. associate-*r/ 16.2%

       \[\leadsto x \cdot \left(\cos \color{blue}{\left(\frac{2 \cdot y}{e^{\left|\log \left(\frac{16}{t \cdot z}\right)\right|}}\right)} \cdot 1\right) \]

    2. rem-square-sqrt 12.8%

       \[\leadsto x \cdot \left(\cos \left(\frac{2 \cdot y}{e^{\left|\color{blue}{\sqrt{\log \left(\frac{16}{t \cdot z}\right)} \cdot \sqrt{\log \left(\frac{16}{t \cdot z}\right)}}\right|}}\right) \cdot 1\right) \]

    3. fabs-sqr 12.8%

       \[\leadsto x \cdot \left(\cos \left(\frac{2 \cdot y}{e^{\color{blue}{\sqrt{\log \left(\frac{16}{t \cdot z}\right)} \cdot \sqrt{\log \left(\frac{16}{t \cdot z}\right)}}}}\right) \cdot 1\right) \]

    4. rem-square-sqrt 14.3%

       \[\leadsto x \cdot \left(\cos \left(\frac{2 \cdot y}{e^{\color{blue}{\log \left(\frac{16}{t \cdot z}\right)}}}\right) \cdot 1\right) \]

    5. associate-/r* 14.4%

       \[\leadsto x \cdot \left(\cos \left(\frac{2 \cdot y}{e^{\log \color{blue}{\left(\frac{\frac{16}{t}}{z}\right)}}}\right) \cdot 1\right) \]

    6. log-div 6.4%

       \[\leadsto x \cdot \left(\cos \left(\frac{2 \cdot y}{e^{\color{blue}{\log \left(\frac{16}{t}\right) - \log z}}}\right) \cdot 1\right) \]

    7. *-commutative 6.4%

       \[\leadsto x \cdot \left(\cos \left(\frac{\color{blue}{y \cdot 2}}{e^{\log \left(\frac{16}{t}\right) - \log z}}\right) \cdot 1\right) \]

    8. log-div 14.4%

       \[\leadsto x \cdot \left(\cos \left(\frac{y \cdot 2}{e^{\color{blue}{\log \left(\frac{\frac{16}{t}}{z}\right)}}}\right) \cdot 1\right) \]

    9. associate-/r* 14.3%

       \[\leadsto x \cdot \left(\cos \left(\frac{y \cdot 2}{e^{\log \color{blue}{\left(\frac{16}{t \cdot z}\right)}}}\right) \cdot 1\right) \]

    10. rem-square-sqrt 12.8%

        \[\leadsto x \cdot \left(\cos \left(\frac{y \cdot 2}{e^{\color{blue}{\sqrt{\log \left(\frac{16}{t \cdot z}\right)} \cdot \sqrt{\log \left(\frac{16}{t \cdot z}\right)}}}}\right) \cdot 1\right) \]

    11. fabs-sqr 12.8%

        \[\leadsto x \cdot \left(\cos \left(\frac{y \cdot 2}{e^{\color{blue}{\left|\sqrt{\log \left(\frac{16}{t \cdot z}\right)} \cdot \sqrt{\log \left(\frac{16}{t \cdot z}\right)}\right|}}}\right) \cdot 1\right) \]

    12. rem-square-sqrt 16.2%

        \[\leadsto x \cdot \left(\cos \left(\frac{y \cdot 2}{e^{\left|\color{blue}{\log \left(\frac{16}{t \cdot z}\right)}\right|}}\right) \cdot 1\right) \]

    13. *-commutative 16.2%

        \[\leadsto x \cdot \left(\cos \left(\frac{y \cdot 2}{e^{\left|\log \left(\frac{16}{\color{blue}{z \cdot t}}\right)\right|}}\right) \cdot 1\right) \]

14. Simplified 16.2%

    \[\leadsto x \cdot \left(\cos \color{blue}{\left(\frac{y \cdot 2}{e^{\left|\log \left(\frac{16}{z \cdot t}\right)\right|}}\right)} \cdot 1\right) \]

15. Final simplification 16.2%

    \[\leadsto x \cdot \cos \left(\frac{y \cdot 2}{e^{\left|\log \left(\frac{16}{z \cdot t}\right)\right|}}\right) \]
16. Add Preprocessing

Alternative 3: 30.3% accurate, 225.0× speedup

Math

\[\begin{array}{l} z_m = \left|z\right| \\ t_m = \left|t\right| \\ x \end{array} \]

FPCore

z_m = (fabs.f64 z)
t_m = (fabs.f64 t)
(FPCore (x y z_m t_m a b) :precision binary64 x)

C

z_m = fabs(z);
t_m = fabs(t);
double code(double x, double y, double z_m, double t_m, double a, double b) {
	return x;
}

Fortran

z_m = abs(z)
t_m = abs(t)
real(8) function code(x, y, z_m, t_m, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z_m
    real(8), intent (in) :: t_m
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = x
end function

Java

z_m = Math.abs(z);
t_m = Math.abs(t);
public static double code(double x, double y, double z_m, double t_m, double a, double b) {
	return x;
}

Python

z_m = math.fabs(z)
t_m = math.fabs(t)
def code(x, y, z_m, t_m, a, b):
	return x

Julia

z_m = abs(z)
t_m = abs(t)
function code(x, y, z_m, t_m, a, b)
	return x
end

MATLAB

z_m = abs(z);
t_m = abs(t);
function tmp = code(x, y, z_m, t_m, a, b)
	tmp = x;
end

Wolfram

z_m = N[Abs[z], $MachinePrecision]
t_m = N[Abs[t], $MachinePrecision]
code[x_, y_, z$95$m_, t$95$m_, a_, b_] := x

Derivation

1. Initial program 27.8%

   \[\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) \]
2. Add Preprocessing

3. Taylor expanded in y around 0 27.6%

   \[\leadsto \left(x \cdot \cos \left(\frac{\color{blue}{t \cdot z}}{16}\right)\right) \cdot \cos \left(\frac{\left(\left(a \cdot 2 + 1\right) \cdot b\right) \cdot t}{16}\right) \]

4. Taylor expanded in t around 0 31.4%

   \[\leadsto \color{blue}{x} \]

5. Final simplification 31.4%

   \[\leadsto x \]
6. Add Preprocessing

Developer target: 29.9% accurate, 1.0× speedup

Math

\[\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} \]

FPCore

(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)))))))

C

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)))));
}

Fortran

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

Java

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)))));
}

Python

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)))))

Julia

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

MATLAB

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

Wolfram

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]

Reproduce

herbie shell --seed 2024012 
(FPCore (x y z t a b)
  :name "Codec.Picture.Jpg.FastDct:referenceDct from JuicyPixels-3.2.6.1"
  :precision binary64

  :herbie-target
  (* 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))))