
(FPCore (x y) :precision binary64 (* (* (* x 3.0) y) y))
double code(double x, double y) {
return ((x * 3.0) * y) * y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = ((x * 3.0d0) * y) * y
end function
public static double code(double x, double y) {
return ((x * 3.0) * y) * y;
}
def code(x, y): return ((x * 3.0) * y) * y
function code(x, y) return Float64(Float64(Float64(x * 3.0) * y) * y) end
function tmp = code(x, y) tmp = ((x * 3.0) * y) * y; end
code[x_, y_] := N[(N[(N[(x * 3.0), $MachinePrecision] * y), $MachinePrecision] * y), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot 3\right) \cdot y\right) \cdot y
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (* (* (* x 3.0) y) y))
double code(double x, double y) {
return ((x * 3.0) * y) * y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = ((x * 3.0d0) * y) * y
end function
public static double code(double x, double y) {
return ((x * 3.0) * y) * y;
}
def code(x, y): return ((x * 3.0) * y) * y
function code(x, y) return Float64(Float64(Float64(x * 3.0) * y) * y) end
function tmp = code(x, y) tmp = ((x * 3.0) * y) * y; end
code[x_, y_] := N[(N[(N[(x * 3.0), $MachinePrecision] * y), $MachinePrecision] * y), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot 3\right) \cdot y\right) \cdot y
\end{array}
(FPCore (x y) :precision binary64 (* (* y (* 3.0 x)) y))
double code(double x, double y) {
return (y * (3.0 * x)) * y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (y * (3.0d0 * x)) * y
end function
public static double code(double x, double y) {
return (y * (3.0 * x)) * y;
}
def code(x, y): return (y * (3.0 * x)) * y
function code(x, y) return Float64(Float64(y * Float64(3.0 * x)) * y) end
function tmp = code(x, y) tmp = (y * (3.0 * x)) * y; end
code[x_, y_] := N[(N[(y * N[(3.0 * x), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision]
\begin{array}{l}
\\
\left(y \cdot \left(3 \cdot x\right)\right) \cdot y
\end{array}
Initial program 99.7%
Final simplification99.7%
(FPCore (x y) :precision binary64 (* (* (* y x) 3.0) y))
double code(double x, double y) {
return ((y * x) * 3.0) * y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = ((y * x) * 3.0d0) * y
end function
public static double code(double x, double y) {
return ((y * x) * 3.0) * y;
}
def code(x, y): return ((y * x) * 3.0) * y
function code(x, y) return Float64(Float64(Float64(y * x) * 3.0) * y) end
function tmp = code(x, y) tmp = ((y * x) * 3.0) * y; end
code[x_, y_] := N[(N[(N[(y * x), $MachinePrecision] * 3.0), $MachinePrecision] * y), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(y \cdot x\right) \cdot 3\right) \cdot y
\end{array}
Initial program 99.7%
Taylor expanded in y around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.6
Applied rewrites99.6%
(FPCore (x y) :precision binary64 (* (* (* y y) x) 3.0))
double code(double x, double y) {
return ((y * y) * x) * 3.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = ((y * y) * x) * 3.0d0
end function
public static double code(double x, double y) {
return ((y * y) * x) * 3.0;
}
def code(x, y): return ((y * y) * x) * 3.0
function code(x, y) return Float64(Float64(Float64(y * y) * x) * 3.0) end
function tmp = code(x, y) tmp = ((y * y) * x) * 3.0; end
code[x_, y_] := N[(N[(N[(y * y), $MachinePrecision] * x), $MachinePrecision] * 3.0), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(y \cdot y\right) \cdot x\right) \cdot 3
\end{array}
Initial program 99.7%
Taylor expanded in y around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6488.2
Applied rewrites88.2%
(FPCore (x y) :precision binary64 (* y (* 3.0 x)))
double code(double x, double y) {
return y * (3.0 * x);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = y * (3.0d0 * x)
end function
public static double code(double x, double y) {
return y * (3.0 * x);
}
def code(x, y): return y * (3.0 * x)
function code(x, y) return Float64(y * Float64(3.0 * x)) end
function tmp = code(x, y) tmp = y * (3.0 * x); end
code[x_, y_] := N[(y * N[(3.0 * x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
y \cdot \left(3 \cdot x\right)
\end{array}
Initial program 99.7%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
pow2N/A
pow-to-expN/A
*-commutativeN/A
count-2N/A
flip-+N/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
flip-+N/A
log-prodN/A
sqr-powN/A
metadata-evalN/A
unpow1N/A
rem-exp-logN/A
lift-*.f6421.7
Applied rewrites21.7%
Final simplification21.7%
(FPCore (x y) :precision binary64 (* (* y x) 3.0))
double code(double x, double y) {
return (y * x) * 3.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (y * x) * 3.0d0
end function
public static double code(double x, double y) {
return (y * x) * 3.0;
}
def code(x, y): return (y * x) * 3.0
function code(x, y) return Float64(Float64(y * x) * 3.0) end
function tmp = code(x, y) tmp = (y * x) * 3.0; end
code[x_, y_] := N[(N[(y * x), $MachinePrecision] * 3.0), $MachinePrecision]
\begin{array}{l}
\\
\left(y \cdot x\right) \cdot 3
\end{array}
Initial program 99.7%
Taylor expanded in y around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6488.2
Applied rewrites88.2%
Applied rewrites21.7%
(FPCore (x y) :precision binary64 (* 3.0 x))
double code(double x, double y) {
return 3.0 * x;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 3.0d0 * x
end function
public static double code(double x, double y) {
return 3.0 * x;
}
def code(x, y): return 3.0 * x
function code(x, y) return Float64(3.0 * x) end
function tmp = code(x, y) tmp = 3.0 * x; end
code[x_, y_] := N[(3.0 * x), $MachinePrecision]
\begin{array}{l}
\\
3 \cdot x
\end{array}
Initial program 99.7%
Applied rewrites4.9%
Final simplification4.9%
(FPCore (x y) :precision binary64 (* (* x (* 3.0 y)) y))
double code(double x, double y) {
return (x * (3.0 * y)) * y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x * (3.0d0 * y)) * y
end function
public static double code(double x, double y) {
return (x * (3.0 * y)) * y;
}
def code(x, y): return (x * (3.0 * y)) * y
function code(x, y) return Float64(Float64(x * Float64(3.0 * y)) * y) end
function tmp = code(x, y) tmp = (x * (3.0 * y)) * y; end
code[x_, y_] := N[(N[(x * N[(3.0 * y), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot \left(3 \cdot y\right)\right) \cdot y
\end{array}
herbie shell --seed 2024277
(FPCore (x y)
:name "Diagrams.Segment:$catParam from diagrams-lib-1.3.0.3, B"
:precision binary64
:alt
(! :herbie-platform default (* (* x (* 3 y)) y))
(* (* (* x 3.0) y) y))