
(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 12 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 (+ e (+ d (fabs (+ a (+ b c))))))
double code(double a, double b, double c, double d, double e) {
return e + (d + fabs((a + (b + c))));
}
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 + abs((a + (b + c))))
end function
public static double code(double a, double b, double c, double d, double e) {
return e + (d + Math.abs((a + (b + c))));
}
def code(a, b, c, d, e): return e + (d + math.fabs((a + (b + c))))
function code(a, b, c, d, e) return Float64(e + Float64(d + abs(Float64(a + Float64(b + c))))) end
function tmp = code(a, b, c, d, e) tmp = e + (d + abs((a + (b + c)))); end
code[a_, b_, c_, d_, e_] := N[(e + N[(d + N[Abs[N[(a + N[(b + c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e + \left(d + \left|a + \left(b + c\right)\right|\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
add-sqr-sqrt99.5%
sqrt-unprod99.6%
pow299.6%
Applied egg-rr99.6%
unpow299.6%
rem-sqrt-square99.6%
associate-+r+99.6%
+-commutative99.6%
+-commutative99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (a b c d e) :precision binary64 (+ a (+ b (+ d (+ e c)))))
double code(double a, double b, double c, double d, double e) {
return a + (b + (d + (e + c)));
}
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 = a + (b + (d + (e + c)))
end function
public static double code(double a, double b, double c, double d, double e) {
return a + (b + (d + (e + c)));
}
def code(a, b, c, d, e): return a + (b + (d + (e + c)))
function code(a, b, c, d, e) return Float64(a + Float64(b + Float64(d + Float64(e + c)))) end
function tmp = code(a, b, c, d, e) tmp = a + (b + (d + (e + c))); end
code[a_, b_, c_, d_, e_] := N[(a + N[(b + N[(d + N[(e + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
a + \left(b + \left(d + \left(e + c\right)\right)\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in e around 0 99.3%
associate-+r+99.4%
+-commutative99.4%
associate-+r+99.4%
Simplified99.4%
Final simplification99.4%
(FPCore (a b c d e) :precision binary64 (+ a (+ c (+ e (+ d b)))))
double code(double a, double b, double c, double d, double e) {
return a + (c + (e + (d + b)));
}
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 = a + (c + (e + (d + b)))
end function
public static double code(double a, double b, double c, double d, double e) {
return a + (c + (e + (d + b)));
}
def code(a, b, c, d, e): return a + (c + (e + (d + b)))
function code(a, b, c, d, e) return Float64(a + Float64(c + Float64(e + Float64(d + b)))) end
function tmp = code(a, b, c, d, e) tmp = a + (c + (e + (d + b))); end
code[a_, b_, c_, d_, e_] := N[(a + N[(c + N[(e + N[(d + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
a + \left(c + \left(e + \left(d + b\right)\right)\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
add-sqr-sqrt99.5%
sqrt-unprod99.6%
pow299.6%
Applied egg-rr99.6%
unpow299.6%
rem-sqrt-square99.6%
associate-+r+99.6%
+-commutative99.6%
+-commutative99.6%
Simplified99.6%
Taylor expanded in e around 0 99.6%
associate-+r+99.6%
rem-square-sqrt99.6%
fabs-sqr99.6%
rem-square-sqrt99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in d around 0 99.3%
associate-+r+99.5%
+-commutative99.5%
associate-+l+99.3%
+-commutative99.3%
+-commutative99.3%
associate-+l+99.4%
+-commutative99.4%
Simplified99.4%
Final simplification99.4%
(FPCore (a b c d e) :precision binary64 (+ a (+ e (+ b (+ d c)))))
double code(double a, double b, double c, double d, double e) {
return a + (e + (b + (d + c)));
}
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 = a + (e + (b + (d + c)))
end function
public static double code(double a, double b, double c, double d, double e) {
return a + (e + (b + (d + c)));
}
def code(a, b, c, d, e): return a + (e + (b + (d + c)))
function code(a, b, c, d, e) return Float64(a + Float64(e + Float64(b + Float64(d + c)))) end
function tmp = code(a, b, c, d, e) tmp = a + (e + (b + (d + c))); end
code[a_, b_, c_, d_, e_] := N[(a + N[(e + N[(b + N[(d + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
a + \left(e + \left(b + \left(d + c\right)\right)\right)
\end{array}
Initial program 99.3%
*-un-lft-identity99.3%
*-commutative99.3%
associate-+l+99.4%
associate-+l+99.5%
Applied egg-rr99.5%
expm1-log1p-u98.7%
expm1-udef98.7%
+-commutative98.7%
*-rgt-identity98.7%
+-commutative98.7%
Applied egg-rr98.7%
expm1-def98.7%
expm1-log1p99.5%
Simplified99.5%
Final simplification99.5%
(FPCore (a b c d e) :precision binary64 (+ b (+ e (+ a (+ d c)))))
double code(double a, double b, double c, double d, double e) {
return b + (e + (a + (d + c)));
}
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 = b + (e + (a + (d + c)))
end function
public static double code(double a, double b, double c, double d, double e) {
return b + (e + (a + (d + c)));
}
def code(a, b, c, d, e): return b + (e + (a + (d + c)))
function code(a, b, c, d, e) return Float64(b + Float64(e + Float64(a + Float64(d + c)))) end
function tmp = code(a, b, c, d, e) tmp = b + (e + (a + (d + c))); end
code[a_, b_, c_, d_, e_] := N[(b + N[(e + N[(a + N[(d + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
b + \left(e + \left(a + \left(d + c\right)\right)\right)
\end{array}
Initial program 99.3%
*-un-lft-identity99.3%
*-commutative99.3%
associate-+l+99.4%
associate-+l+99.5%
Applied egg-rr99.5%
expm1-log1p-u98.7%
expm1-udef98.7%
+-commutative98.7%
*-rgt-identity98.7%
+-commutative98.7%
Applied egg-rr98.7%
expm1-def98.7%
expm1-log1p99.5%
Simplified99.5%
Taylor expanded in a around 0 99.3%
+-commutative99.3%
+-commutative99.3%
+-commutative99.3%
associate-+r+99.4%
associate-+l+99.4%
associate-+l+99.5%
+-commutative99.5%
Simplified99.5%
Final simplification99.5%
(FPCore (a b c d e) :precision binary64 (+ d (+ e (+ a (+ b c)))))
double code(double a, double b, double c, double d, double e) {
return d + (e + (a + (b + c)));
}
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 + (e + (a + (b + c)))
end function
public static double code(double a, double b, double c, double d, double e) {
return d + (e + (a + (b + c)));
}
def code(a, b, c, d, e): return d + (e + (a + (b + c)))
function code(a, b, c, d, e) return Float64(d + Float64(e + Float64(a + Float64(b + c)))) end
function tmp = code(a, b, c, d, e) tmp = d + (e + (a + (b + c))); end
code[a_, b_, c_, d_, e_] := N[(d + N[(e + N[(a + N[(b + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
d + \left(e + \left(a + \left(b + c\right)\right)\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
add-sqr-sqrt99.5%
sqrt-unprod99.6%
pow299.6%
Applied egg-rr99.6%
unpow299.6%
rem-sqrt-square99.6%
associate-+r+99.6%
+-commutative99.6%
+-commutative99.6%
Simplified99.6%
Taylor expanded in e around 0 99.6%
associate-+r+99.6%
rem-square-sqrt99.6%
fabs-sqr99.6%
rem-square-sqrt99.6%
associate-+r+99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (a b c d e) :precision binary64 (+ b (+ c (+ e d))))
double code(double a, double b, double c, double d, double e) {
return b + (c + (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 = b + (c + (e + d))
end function
public static double code(double a, double b, double c, double d, double e) {
return b + (c + (e + d));
}
def code(a, b, c, d, e): return b + (c + (e + d))
function code(a, b, c, d, e) return Float64(b + Float64(c + Float64(e + d))) end
function tmp = code(a, b, c, d, e) tmp = b + (c + (e + d)); end
code[a_, b_, c_, d_, e_] := N[(b + N[(c + N[(e + d), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
b + \left(c + \left(e + d\right)\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in a around 0 25.7%
Final simplification25.7%
(FPCore (a b c d e) :precision binary64 (+ b (+ d (+ e c))))
double code(double a, double b, double c, double d, double e) {
return b + (d + (e + c));
}
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 = b + (d + (e + c))
end function
public static double code(double a, double b, double c, double d, double e) {
return b + (d + (e + c));
}
def code(a, b, c, d, e): return b + (d + (e + c))
function code(a, b, c, d, e) return Float64(b + Float64(d + Float64(e + c))) end
function tmp = code(a, b, c, d, e) tmp = b + (d + (e + c)); end
code[a_, b_, c_, d_, e_] := N[(b + N[(d + N[(e + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
b + \left(d + \left(e + c\right)\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in a around 0 25.7%
associate-+r+25.7%
+-commutative25.7%
associate-+r+25.7%
Simplified25.7%
Final simplification25.7%
(FPCore (a b c d e) :precision binary64 (+ d (+ e (+ b c))))
double code(double a, double b, double c, double d, double e) {
return d + (e + (b + c));
}
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 + (e + (b + c))
end function
public static double code(double a, double b, double c, double d, double e) {
return d + (e + (b + c));
}
def code(a, b, c, d, e): return d + (e + (b + c))
function code(a, b, c, d, e) return Float64(d + Float64(e + Float64(b + c))) end
function tmp = code(a, b, c, d, e) tmp = d + (e + (b + c)); end
code[a_, b_, c_, d_, e_] := N[(d + N[(e + N[(b + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
d + \left(e + \left(b + c\right)\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
add-sqr-sqrt99.5%
sqrt-unprod99.6%
pow299.6%
Applied egg-rr99.6%
unpow299.6%
rem-sqrt-square99.6%
associate-+r+99.6%
+-commutative99.6%
+-commutative99.6%
Simplified99.6%
Taylor expanded in e around 0 99.6%
associate-+r+99.6%
rem-square-sqrt99.6%
fabs-sqr99.6%
rem-square-sqrt99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in a around 0 25.7%
Final simplification25.7%
(FPCore (a b c d e) :precision binary64 (+ c (+ e d)))
double code(double a, double b, double c, double d, double e) {
return c + (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 + (e + d)
end function
public static double code(double a, double b, double c, double d, double e) {
return c + (e + d);
}
def code(a, b, c, d, e): return c + (e + d)
function code(a, b, c, d, e) return Float64(c + Float64(e + d)) end
function tmp = code(a, b, c, d, e) tmp = c + (e + d); end
code[a_, b_, c_, d_, e_] := N[(c + N[(e + d), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
c + \left(e + d\right)
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in a around 0 25.7%
Taylor expanded in b around 0 23.2%
Final simplification23.2%
(FPCore (a b c d e) :precision binary64 (+ e d))
double code(double a, double b, double c, double d, double e) {
return 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 = e + d
end function
public static double code(double a, double b, double c, double d, double e) {
return e + d;
}
def code(a, b, c, d, e): return e + d
function code(a, b, c, d, e) return Float64(e + d) end
function tmp = code(a, b, c, d, e) tmp = e + d; end
code[a_, b_, c_, d_, e_] := N[(e + d), $MachinePrecision]
\begin{array}{l}
\\
e + d
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in a around 0 25.7%
Taylor expanded in b around 0 23.2%
Taylor expanded in c around 0 21.1%
Final simplification21.1%
(FPCore (a b c d e) :precision binary64 e)
double code(double a, double b, double c, double d, double e) {
return 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 = e
end function
public static double code(double a, double b, double c, double d, double e) {
return e;
}
def code(a, b, c, d, e): return e
function code(a, b, c, d, e) return e end
function tmp = code(a, b, c, d, e) tmp = e; end
code[a_, b_, c_, d_, e_] := e
\begin{array}{l}
\\
e
\end{array}
Initial program 99.3%
associate-+l+99.5%
associate-+l+99.6%
associate-+r+99.6%
associate-+r+99.6%
Simplified99.6%
Taylor expanded in e around inf 18.9%
Final simplification18.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 2023306
(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))
:herbie-target
(+ (+ d (+ c (+ a b))) e)
(+ (+ (+ (+ e d) c) b) a))