
(FPCore (a b c d e) :precision binary64 (+ (+ (+ (+ e d) c) b) a))
double code(double a, double b, double c, double d, double e) {
return (((e + d) + c) + b) + a;
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = (((e + d) + c) + b) + a
end function
public static double code(double a, double b, double c, double d, double e) {
return (((e + d) + c) + b) + a;
}
def code(a, b, c, d, e): return (((e + d) + c) + b) + a
function code(a, b, c, d, e) return Float64(Float64(Float64(Float64(e + d) + c) + b) + a) end
function tmp = code(a, b, c, d, e) tmp = (((e + d) + c) + b) + a; end
code[a_, b_, c_, d_, e_] := N[(N[(N[(N[(e + d), $MachinePrecision] + c), $MachinePrecision] + b), $MachinePrecision] + a), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(e + d\right) + c\right) + b\right) + a
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b c d e) :precision binary64 (+ (+ (+ (+ e d) c) b) a))
double code(double a, double b, double c, double d, double e) {
return (((e + d) + c) + b) + a;
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = (((e + d) + c) + b) + a
end function
public static double code(double a, double b, double c, double d, double e) {
return (((e + d) + c) + b) + a;
}
def code(a, b, c, d, e): return (((e + d) + c) + b) + a
function code(a, b, c, d, e) return Float64(Float64(Float64(Float64(e + d) + c) + b) + a) end
function tmp = code(a, b, c, d, e) tmp = (((e + d) + c) + b) + a; end
code[a_, b_, c_, d_, e_] := N[(N[(N[(N[(e + d), $MachinePrecision] + c), $MachinePrecision] + b), $MachinePrecision] + a), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(e + d\right) + c\right) + b\right) + a
\end{array}
(FPCore (a b c d e) :precision binary64 (+ (+ d c) (+ e (+ b a))))
double code(double a, double b, double c, double d, double e) {
return (d + c) + (e + (b + a));
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = (d + c) + (e + (b + a))
end function
public static double code(double a, double b, double c, double d, double e) {
return (d + c) + (e + (b + a));
}
def code(a, b, c, d, e): return (d + c) + (e + (b + a))
function code(a, b, c, d, e) return Float64(Float64(d + c) + Float64(e + Float64(b + a))) end
function tmp = code(a, b, c, d, e) tmp = (d + c) + (e + (b + a)); end
code[a_, b_, c_, d_, e_] := N[(N[(d + c), $MachinePrecision] + N[(e + N[(b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(d + c\right) + \left(e + \left(b + a\right)\right)
\end{array}
Initial program 99.4%
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6499.6
Applied rewrites99.6%
Final simplification99.6%
(FPCore (a b c d e) :precision binary64 (+ (+ c e) (+ d (+ b a))))
double code(double a, double b, double c, double d, double e) {
return (c + e) + (d + (b + a));
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = (c + e) + (d + (b + a))
end function
public static double code(double a, double b, double c, double d, double e) {
return (c + e) + (d + (b + a));
}
def code(a, b, c, d, e): return (c + e) + (d + (b + a))
function code(a, b, c, d, e) return Float64(Float64(c + e) + Float64(d + Float64(b + a))) end
function tmp = code(a, b, c, d, e) tmp = (c + e) + (d + (b + a)); end
code[a_, b_, c_, d_, e_] := N[(N[(c + e), $MachinePrecision] + N[(d + N[(b + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(c + e\right) + \left(d + \left(b + a\right)\right)
\end{array}
Initial program 99.4%
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
+-commutativeN/A
lift-+.f64N/A
lift-+.f64N/A
+-commutativeN/A
associate-+l+N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-+.f6499.6
Applied rewrites99.6%
Final simplification99.6%
(FPCore (a b c d e) :precision binary64 (+ (+ (+ c e) b) (+ d a)))
double code(double a, double b, double c, double d, double e) {
return ((c + e) + b) + (d + a);
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = ((c + e) + b) + (d + a)
end function
public static double code(double a, double b, double c, double d, double e) {
return ((c + e) + b) + (d + a);
}
def code(a, b, c, d, e): return ((c + e) + b) + (d + a)
function code(a, b, c, d, e) return Float64(Float64(Float64(c + e) + b) + Float64(d + a)) end
function tmp = code(a, b, c, d, e) tmp = ((c + e) + b) + (d + a); end
code[a_, b_, c_, d_, e_] := N[(N[(N[(c + e), $MachinePrecision] + b), $MachinePrecision] + N[(d + a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(c + e\right) + b\right) + \left(d + a\right)
\end{array}
Initial program 99.4%
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
lift-+.f64N/A
+-commutativeN/A
associate-+l+N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f6499.6
Applied rewrites99.6%
Final simplification99.6%
(FPCore (a b c d e) :precision binary64 (+ (+ (+ (+ c b) a) e) d))
double code(double a, double b, double c, double d, double e) {
return (((c + b) + a) + e) + d;
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = (((c + b) + a) + e) + d
end function
public static double code(double a, double b, double c, double d, double e) {
return (((c + b) + a) + e) + d;
}
def code(a, b, c, d, e): return (((c + b) + a) + e) + d
function code(a, b, c, d, e) return Float64(Float64(Float64(Float64(c + b) + a) + e) + d) end
function tmp = code(a, b, c, d, e) tmp = (((c + b) + a) + e) + d; end
code[a_, b_, c_, d_, e_] := N[(N[(N[(N[(c + b), $MachinePrecision] + a), $MachinePrecision] + e), $MachinePrecision] + d), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(c + b\right) + a\right) + e\right) + d
\end{array}
Initial program 99.4%
lift-+.f64N/A
lift-+.f64N/A
lift-+.f64N/A
associate-+l+N/A
associate-+l+N/A
lift-+.f64N/A
+-commutativeN/A
associate-+l+N/A
lower-+.f64N/A
lower-+.f64N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6499.6
Applied rewrites99.6%
Final simplification99.6%
(FPCore (a b c d e) :precision binary64 (+ (+ (+ d c) b) e))
double code(double a, double b, double c, double d, double e) {
return ((d + c) + b) + e;
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = ((d + c) + b) + e
end function
public static double code(double a, double b, double c, double d, double e) {
return ((d + c) + b) + e;
}
def code(a, b, c, d, e): return ((d + c) + b) + e
function code(a, b, c, d, e) return Float64(Float64(Float64(d + c) + b) + e) end
function tmp = code(a, b, c, d, e) tmp = ((d + c) + b) + e; end
code[a_, b_, c_, d_, e_] := N[(N[(N[(d + c), $MachinePrecision] + b), $MachinePrecision] + e), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(d + c\right) + b\right) + e
\end{array}
Initial program 99.4%
Taylor expanded in a around 0
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-+.f6425.7
Applied rewrites25.7%
Applied rewrites25.7%
Final simplification25.7%
(FPCore (a b c d e) :precision binary64 (+ (+ d c) e))
double code(double a, double b, double c, double d, double e) {
return (d + c) + e;
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = (d + c) + e
end function
public static double code(double a, double b, double c, double d, double e) {
return (d + c) + e;
}
def code(a, b, c, d, e): return (d + c) + e
function code(a, b, c, d, e) return Float64(Float64(d + c) + e) end
function tmp = code(a, b, c, d, e) tmp = (d + c) + e; end
code[a_, b_, c_, d_, e_] := N[(N[(d + c), $MachinePrecision] + e), $MachinePrecision]
\begin{array}{l}
\\
\left(d + c\right) + e
\end{array}
Initial program 99.4%
Taylor expanded in a around 0
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-+.f6425.7
Applied rewrites25.7%
Taylor expanded in b around 0
Applied rewrites23.2%
Applied rewrites23.2%
Final simplification23.2%
(FPCore (a b c d e) :precision binary64 (+ c e))
double code(double a, double b, double c, double d, double e) {
return c + e;
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = c + e
end function
public static double code(double a, double b, double c, double d, double e) {
return c + e;
}
def code(a, b, c, d, e): return c + e
function code(a, b, c, d, e) return Float64(c + e) end
function tmp = code(a, b, c, d, e) tmp = c + e; end
code[a_, b_, c_, d_, e_] := N[(c + e), $MachinePrecision]
\begin{array}{l}
\\
c + e
\end{array}
Initial program 99.4%
Taylor expanded in a around 0
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-+.f6425.7
Applied rewrites25.7%
Taylor expanded in b around 0
Applied rewrites23.2%
Taylor expanded in d around 0
Applied rewrites19.9%
Final simplification19.9%
(FPCore (a b c d e) :precision binary64 (+ (+ d (+ c (+ a b))) e))
double code(double a, double b, double c, double d, double e) {
return (d + (c + (a + b))) + e;
}
real(8) function code(a, b, c, d, e)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: d
real(8), intent (in) :: e
code = (d + (c + (a + b))) + e
end function
public static double code(double a, double b, double c, double d, double e) {
return (d + (c + (a + b))) + e;
}
def code(a, b, c, d, e): return (d + (c + (a + b))) + e
function code(a, b, c, d, e) return Float64(Float64(d + Float64(c + Float64(a + b))) + e) end
function tmp = code(a, b, c, d, e) tmp = (d + (c + (a + b))) + e; end
code[a_, b_, c_, d_, e_] := N[(N[(d + N[(c + N[(a + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + e), $MachinePrecision]
\begin{array}{l}
\\
\left(d + \left(c + \left(a + b\right)\right)\right) + e
\end{array}
herbie shell --seed 2024295
(FPCore (a b c d e)
:name "Expression 1, p15"
:precision binary64
:pre (and (and (and (and (and (and (and (and (and (<= 1.0 a) (<= a 2.0)) (<= 2.0 b)) (<= b 4.0)) (<= 4.0 c)) (<= c 8.0)) (<= 8.0 d)) (<= d 16.0)) (<= 16.0 e)) (<= e 32.0))
:alt
(! :herbie-platform default (+ (+ d (+ c (+ a b))) e))
(+ (+ (+ (+ e d) c) b) a))