
(FPCore (x y z) :precision binary64 (+ (+ (+ (+ (+ x y) y) x) z) x))
double code(double x, double y, double z) {
return ((((x + y) + y) + x) + z) + x;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((((x + y) + y) + x) + z) + x
end function
public static double code(double x, double y, double z) {
return ((((x + y) + y) + x) + z) + x;
}
def code(x, y, z): return ((((x + y) + y) + x) + z) + x
function code(x, y, z) return Float64(Float64(Float64(Float64(Float64(x + y) + y) + x) + z) + x) end
function tmp = code(x, y, z) tmp = ((((x + y) + y) + x) + z) + x; end
code[x_, y_, z_] := N[(N[(N[(N[(N[(x + y), $MachinePrecision] + y), $MachinePrecision] + x), $MachinePrecision] + z), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(x + y\right) + y\right) + x\right) + z\right) + x
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z) :precision binary64 (+ (+ (+ (+ (+ x y) y) x) z) x))
double code(double x, double y, double z) {
return ((((x + y) + y) + x) + z) + x;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((((x + y) + y) + x) + z) + x
end function
public static double code(double x, double y, double z) {
return ((((x + y) + y) + x) + z) + x;
}
def code(x, y, z): return ((((x + y) + y) + x) + z) + x
function code(x, y, z) return Float64(Float64(Float64(Float64(Float64(x + y) + y) + x) + z) + x) end
function tmp = code(x, y, z) tmp = ((((x + y) + y) + x) + z) + x; end
code[x_, y_, z_] := N[(N[(N[(N[(N[(x + y), $MachinePrecision] + y), $MachinePrecision] + x), $MachinePrecision] + z), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(x + y\right) + y\right) + x\right) + z\right) + x
\end{array}
(FPCore (x y z) :precision binary64 (- z (fma x -3.0 (* y -2.0))))
double code(double x, double y, double z) {
return z - fma(x, -3.0, (y * -2.0));
}
function code(x, y, z) return Float64(z - fma(x, -3.0, Float64(y * -2.0))) end
code[x_, y_, z_] := N[(z - N[(x * -3.0 + N[(y * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
z - \mathsf{fma}\left(x, -3, y \cdot -2\right)
\end{array}
Initial program 99.9%
+-commutative99.9%
associate-+l+99.9%
remove-double-neg99.9%
unsub-neg99.9%
+-commutative99.9%
+-commutative99.9%
associate-+l+99.9%
associate-+r+99.9%
associate-+r+99.9%
distribute-neg-in99.9%
distribute-neg-out99.9%
distribute-neg-out99.9%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
fma-def99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x y z)
:precision binary64
(if (<= y -9e+149)
(* y 2.0)
(if (<= y -6.2e+47)
(* x 3.0)
(if (<= y -1.3e-79)
z
(if (<= y 1.04e-243)
(* x 3.0)
(if (<= y 2.8e-183) z (if (<= y 1.2e+87) (* x 3.0) (* y 2.0))))))))
double code(double x, double y, double z) {
double tmp;
if (y <= -9e+149) {
tmp = y * 2.0;
} else if (y <= -6.2e+47) {
tmp = x * 3.0;
} else if (y <= -1.3e-79) {
tmp = z;
} else if (y <= 1.04e-243) {
tmp = x * 3.0;
} else if (y <= 2.8e-183) {
tmp = z;
} else if (y <= 1.2e+87) {
tmp = x * 3.0;
} else {
tmp = y * 2.0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (y <= (-9d+149)) then
tmp = y * 2.0d0
else if (y <= (-6.2d+47)) then
tmp = x * 3.0d0
else if (y <= (-1.3d-79)) then
tmp = z
else if (y <= 1.04d-243) then
tmp = x * 3.0d0
else if (y <= 2.8d-183) then
tmp = z
else if (y <= 1.2d+87) then
tmp = x * 3.0d0
else
tmp = y * 2.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -9e+149) {
tmp = y * 2.0;
} else if (y <= -6.2e+47) {
tmp = x * 3.0;
} else if (y <= -1.3e-79) {
tmp = z;
} else if (y <= 1.04e-243) {
tmp = x * 3.0;
} else if (y <= 2.8e-183) {
tmp = z;
} else if (y <= 1.2e+87) {
tmp = x * 3.0;
} else {
tmp = y * 2.0;
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -9e+149: tmp = y * 2.0 elif y <= -6.2e+47: tmp = x * 3.0 elif y <= -1.3e-79: tmp = z elif y <= 1.04e-243: tmp = x * 3.0 elif y <= 2.8e-183: tmp = z elif y <= 1.2e+87: tmp = x * 3.0 else: tmp = y * 2.0 return tmp
function code(x, y, z) tmp = 0.0 if (y <= -9e+149) tmp = Float64(y * 2.0); elseif (y <= -6.2e+47) tmp = Float64(x * 3.0); elseif (y <= -1.3e-79) tmp = z; elseif (y <= 1.04e-243) tmp = Float64(x * 3.0); elseif (y <= 2.8e-183) tmp = z; elseif (y <= 1.2e+87) tmp = Float64(x * 3.0); else tmp = Float64(y * 2.0); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -9e+149) tmp = y * 2.0; elseif (y <= -6.2e+47) tmp = x * 3.0; elseif (y <= -1.3e-79) tmp = z; elseif (y <= 1.04e-243) tmp = x * 3.0; elseif (y <= 2.8e-183) tmp = z; elseif (y <= 1.2e+87) tmp = x * 3.0; else tmp = y * 2.0; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -9e+149], N[(y * 2.0), $MachinePrecision], If[LessEqual[y, -6.2e+47], N[(x * 3.0), $MachinePrecision], If[LessEqual[y, -1.3e-79], z, If[LessEqual[y, 1.04e-243], N[(x * 3.0), $MachinePrecision], If[LessEqual[y, 2.8e-183], z, If[LessEqual[y, 1.2e+87], N[(x * 3.0), $MachinePrecision], N[(y * 2.0), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9 \cdot 10^{+149}:\\
\;\;\;\;y \cdot 2\\
\mathbf{elif}\;y \leq -6.2 \cdot 10^{+47}:\\
\;\;\;\;x \cdot 3\\
\mathbf{elif}\;y \leq -1.3 \cdot 10^{-79}:\\
\;\;\;\;z\\
\mathbf{elif}\;y \leq 1.04 \cdot 10^{-243}:\\
\;\;\;\;x \cdot 3\\
\mathbf{elif}\;y \leq 2.8 \cdot 10^{-183}:\\
\;\;\;\;z\\
\mathbf{elif}\;y \leq 1.2 \cdot 10^{+87}:\\
\;\;\;\;x \cdot 3\\
\mathbf{else}:\\
\;\;\;\;y \cdot 2\\
\end{array}
\end{array}
if y < -8.99999999999999965e149 or 1.19999999999999991e87 < y Initial program 100.0%
associate-+l+100.0%
associate-+l+100.0%
+-commutative100.0%
count-2100.0%
+-commutative100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in y around inf 79.8%
if -8.99999999999999965e149 < y < -6.2000000000000001e47 or -1.29999999999999997e-79 < y < 1.0400000000000001e-243 or 2.79999999999999985e-183 < y < 1.19999999999999991e87Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
count-299.9%
+-commutative99.9%
+-commutative99.9%
Simplified99.9%
Taylor expanded in x around inf 54.6%
if -6.2000000000000001e47 < y < -1.29999999999999997e-79 or 1.0400000000000001e-243 < y < 2.79999999999999985e-183Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
count-299.9%
+-commutative99.9%
+-commutative99.9%
Simplified99.9%
Taylor expanded in z around inf 66.1%
Final simplification64.8%
(FPCore (x y z) :precision binary64 (if (or (<= y -1.15e+150) (not (<= y 1.8e+87))) (* y 2.0) (- z (* x -3.0))))
double code(double x, double y, double z) {
double tmp;
if ((y <= -1.15e+150) || !(y <= 1.8e+87)) {
tmp = y * 2.0;
} else {
tmp = z - (x * -3.0);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((y <= (-1.15d+150)) .or. (.not. (y <= 1.8d+87))) then
tmp = y * 2.0d0
else
tmp = z - (x * (-3.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((y <= -1.15e+150) || !(y <= 1.8e+87)) {
tmp = y * 2.0;
} else {
tmp = z - (x * -3.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (y <= -1.15e+150) or not (y <= 1.8e+87): tmp = y * 2.0 else: tmp = z - (x * -3.0) return tmp
function code(x, y, z) tmp = 0.0 if ((y <= -1.15e+150) || !(y <= 1.8e+87)) tmp = Float64(y * 2.0); else tmp = Float64(z - Float64(x * -3.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((y <= -1.15e+150) || ~((y <= 1.8e+87))) tmp = y * 2.0; else tmp = z - (x * -3.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[y, -1.15e+150], N[Not[LessEqual[y, 1.8e+87]], $MachinePrecision]], N[(y * 2.0), $MachinePrecision], N[(z - N[(x * -3.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.15 \cdot 10^{+150} \lor \neg \left(y \leq 1.8 \cdot 10^{+87}\right):\\
\;\;\;\;y \cdot 2\\
\mathbf{else}:\\
\;\;\;\;z - x \cdot -3\\
\end{array}
\end{array}
if y < -1.15000000000000001e150 or 1.79999999999999997e87 < y Initial program 100.0%
associate-+l+100.0%
associate-+l+100.0%
+-commutative100.0%
count-2100.0%
+-commutative100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in y around inf 79.8%
if -1.15000000000000001e150 < y < 1.79999999999999997e87Initial program 99.9%
+-commutative99.9%
associate-+l+99.9%
remove-double-neg99.9%
unsub-neg99.9%
+-commutative99.9%
+-commutative99.9%
associate-+l+99.9%
associate-+r+99.9%
associate-+r+99.9%
distribute-neg-in99.9%
distribute-neg-out99.9%
distribute-neg-out99.9%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
fma-def99.9%
Simplified99.9%
Taylor expanded in y around 0 88.2%
Final simplification85.5%
(FPCore (x y z) :precision binary64 (if (or (<= y -9e+149) (not (<= y 4.5e+68))) (- z (* y -2.0)) (- z (* x -3.0))))
double code(double x, double y, double z) {
double tmp;
if ((y <= -9e+149) || !(y <= 4.5e+68)) {
tmp = z - (y * -2.0);
} else {
tmp = z - (x * -3.0);
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((y <= (-9d+149)) .or. (.not. (y <= 4.5d+68))) then
tmp = z - (y * (-2.0d0))
else
tmp = z - (x * (-3.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((y <= -9e+149) || !(y <= 4.5e+68)) {
tmp = z - (y * -2.0);
} else {
tmp = z - (x * -3.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (y <= -9e+149) or not (y <= 4.5e+68): tmp = z - (y * -2.0) else: tmp = z - (x * -3.0) return tmp
function code(x, y, z) tmp = 0.0 if ((y <= -9e+149) || !(y <= 4.5e+68)) tmp = Float64(z - Float64(y * -2.0)); else tmp = Float64(z - Float64(x * -3.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((y <= -9e+149) || ~((y <= 4.5e+68))) tmp = z - (y * -2.0); else tmp = z - (x * -3.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[y, -9e+149], N[Not[LessEqual[y, 4.5e+68]], $MachinePrecision]], N[(z - N[(y * -2.0), $MachinePrecision]), $MachinePrecision], N[(z - N[(x * -3.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9 \cdot 10^{+149} \lor \neg \left(y \leq 4.5 \cdot 10^{+68}\right):\\
\;\;\;\;z - y \cdot -2\\
\mathbf{else}:\\
\;\;\;\;z - x \cdot -3\\
\end{array}
\end{array}
if y < -8.99999999999999965e149 or 4.5000000000000003e68 < y Initial program 100.0%
+-commutative100.0%
associate-+l+99.9%
remove-double-neg99.9%
unsub-neg99.9%
+-commutative99.9%
+-commutative99.9%
associate-+l+99.9%
associate-+r+99.9%
associate-+r+99.9%
distribute-neg-in99.9%
distribute-neg-out99.9%
distribute-neg-out99.9%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
fma-def100.0%
Simplified100.0%
Taylor expanded in x around 0 94.0%
if -8.99999999999999965e149 < y < 4.5000000000000003e68Initial program 99.9%
+-commutative99.9%
associate-+l+99.9%
remove-double-neg99.9%
unsub-neg99.9%
+-commutative99.9%
+-commutative99.9%
associate-+l+99.9%
associate-+r+99.9%
associate-+r+99.9%
distribute-neg-in99.9%
distribute-neg-out99.9%
distribute-neg-out99.9%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
fma-def99.9%
Simplified99.9%
Taylor expanded in y around 0 88.1%
Final simplification90.0%
(FPCore (x y z) :precision binary64 (if (or (<= y -9.6e+42) (not (<= y 8.6e+86))) (* y 2.0) z))
double code(double x, double y, double z) {
double tmp;
if ((y <= -9.6e+42) || !(y <= 8.6e+86)) {
tmp = y * 2.0;
} else {
tmp = z;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((y <= (-9.6d+42)) .or. (.not. (y <= 8.6d+86))) then
tmp = y * 2.0d0
else
tmp = z
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((y <= -9.6e+42) || !(y <= 8.6e+86)) {
tmp = y * 2.0;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (y <= -9.6e+42) or not (y <= 8.6e+86): tmp = y * 2.0 else: tmp = z return tmp
function code(x, y, z) tmp = 0.0 if ((y <= -9.6e+42) || !(y <= 8.6e+86)) tmp = Float64(y * 2.0); else tmp = z; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((y <= -9.6e+42) || ~((y <= 8.6e+86))) tmp = y * 2.0; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[y, -9.6e+42], N[Not[LessEqual[y, 8.6e+86]], $MachinePrecision]], N[(y * 2.0), $MachinePrecision], z]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.6 \cdot 10^{+42} \lor \neg \left(y \leq 8.6 \cdot 10^{+86}\right):\\
\;\;\;\;y \cdot 2\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if y < -9.5999999999999994e42 or 8.6000000000000004e86 < y Initial program 99.9%
associate-+l+99.9%
associate-+l+100.0%
+-commutative100.0%
count-2100.0%
+-commutative100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in y around inf 71.2%
if -9.5999999999999994e42 < y < 8.6000000000000004e86Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
count-299.9%
+-commutative99.9%
+-commutative99.9%
Simplified99.9%
Taylor expanded in z around inf 46.1%
Final simplification55.5%
(FPCore (x y z) :precision binary64 (+ (* 2.0 (+ x y)) (+ z x)))
double code(double x, double y, double z) {
return (2.0 * (x + y)) + (z + x);
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = (2.0d0 * (x + y)) + (z + x)
end function
public static double code(double x, double y, double z) {
return (2.0 * (x + y)) + (z + x);
}
def code(x, y, z): return (2.0 * (x + y)) + (z + x)
function code(x, y, z) return Float64(Float64(2.0 * Float64(x + y)) + Float64(z + x)) end
function tmp = code(x, y, z) tmp = (2.0 * (x + y)) + (z + x); end
code[x_, y_, z_] := N[(N[(2.0 * N[(x + y), $MachinePrecision]), $MachinePrecision] + N[(z + x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
2 \cdot \left(x + y\right) + \left(z + x\right)
\end{array}
Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
count-299.9%
+-commutative99.9%
+-commutative99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x y z) :precision binary64 z)
double code(double x, double y, double z) {
return z;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = z
end function
public static double code(double x, double y, double z) {
return z;
}
def code(x, y, z): return z
function code(x, y, z) return z end
function tmp = code(x, y, z) tmp = z; end
code[x_, y_, z_] := z
\begin{array}{l}
\\
z
\end{array}
Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
count-299.9%
+-commutative99.9%
+-commutative99.9%
Simplified99.9%
Taylor expanded in z around inf 35.1%
Final simplification35.1%
herbie shell --seed 2024024
(FPCore (x y z)
:name "Graphics.Rendering.Plot.Render.Plot.Legend:renderLegendInside from plot-0.2.3.4"
:precision binary64
(+ (+ (+ (+ (+ x y) y) x) z) x))