
(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 10 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%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
distribute-neg-out99.9%
fma-define100.0%
Simplified100.0%
(FPCore (x y z) :precision binary64 (if (<= z -4e+106) z (if (<= z -0.0068) (* x 3.0) (if (<= z 1.75e+92) (* 2.0 (+ x y)) z))))
double code(double x, double y, double z) {
double tmp;
if (z <= -4e+106) {
tmp = z;
} else if (z <= -0.0068) {
tmp = x * 3.0;
} else if (z <= 1.75e+92) {
tmp = 2.0 * (x + y);
} 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 (z <= (-4d+106)) then
tmp = z
else if (z <= (-0.0068d0)) then
tmp = x * 3.0d0
else if (z <= 1.75d+92) then
tmp = 2.0d0 * (x + y)
else
tmp = z
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -4e+106) {
tmp = z;
} else if (z <= -0.0068) {
tmp = x * 3.0;
} else if (z <= 1.75e+92) {
tmp = 2.0 * (x + y);
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -4e+106: tmp = z elif z <= -0.0068: tmp = x * 3.0 elif z <= 1.75e+92: tmp = 2.0 * (x + y) else: tmp = z return tmp
function code(x, y, z) tmp = 0.0 if (z <= -4e+106) tmp = z; elseif (z <= -0.0068) tmp = Float64(x * 3.0); elseif (z <= 1.75e+92) tmp = Float64(2.0 * Float64(x + y)); else tmp = z; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -4e+106) tmp = z; elseif (z <= -0.0068) tmp = x * 3.0; elseif (z <= 1.75e+92) tmp = 2.0 * (x + y); else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -4e+106], z, If[LessEqual[z, -0.0068], N[(x * 3.0), $MachinePrecision], If[LessEqual[z, 1.75e+92], N[(2.0 * N[(x + y), $MachinePrecision]), $MachinePrecision], z]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4 \cdot 10^{+106}:\\
\;\;\;\;z\\
\mathbf{elif}\;z \leq -0.0068:\\
\;\;\;\;x \cdot 3\\
\mathbf{elif}\;z \leq 1.75 \cdot 10^{+92}:\\
\;\;\;\;2 \cdot \left(x + y\right)\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if z < -4.00000000000000036e106 or 1.74999999999999993e92 < z Initial program 100.0%
+-commutative100.0%
associate-+l+100.0%
remove-double-neg100.0%
unsub-neg100.0%
+-commutative100.0%
+-commutative100.0%
associate-+l+100.0%
associate-+r+100.0%
associate-+r+100.0%
distribute-neg-in100.0%
distribute-neg-out100.0%
neg-mul-1100.0%
count-2100.0%
distribute-lft-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-rgt-out100.0%
distribute-neg-out100.0%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around 0 85.3%
Taylor expanded in z around inf 70.4%
if -4.00000000000000036e106 < z < -0.00679999999999999962Initial 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 57.4%
if -0.00679999999999999962 < z < 1.74999999999999993e92Initial 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 0 71.9%
Taylor expanded in z around 0 61.4%
+-commutative61.4%
Simplified61.4%
Final simplification63.7%
(FPCore (x y z) :precision binary64 (if (<= z -4.1e+105) z (if (<= z -0.0285) (* x 3.0) (if (<= z 8.6e+89) (* y 2.0) z))))
double code(double x, double y, double z) {
double tmp;
if (z <= -4.1e+105) {
tmp = z;
} else if (z <= -0.0285) {
tmp = x * 3.0;
} else if (z <= 8.6e+89) {
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 (z <= (-4.1d+105)) then
tmp = z
else if (z <= (-0.0285d0)) then
tmp = x * 3.0d0
else if (z <= 8.6d+89) 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 (z <= -4.1e+105) {
tmp = z;
} else if (z <= -0.0285) {
tmp = x * 3.0;
} else if (z <= 8.6e+89) {
tmp = y * 2.0;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -4.1e+105: tmp = z elif z <= -0.0285: tmp = x * 3.0 elif z <= 8.6e+89: tmp = y * 2.0 else: tmp = z return tmp
function code(x, y, z) tmp = 0.0 if (z <= -4.1e+105) tmp = z; elseif (z <= -0.0285) tmp = Float64(x * 3.0); elseif (z <= 8.6e+89) tmp = Float64(y * 2.0); else tmp = z; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -4.1e+105) tmp = z; elseif (z <= -0.0285) tmp = x * 3.0; elseif (z <= 8.6e+89) tmp = y * 2.0; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -4.1e+105], z, If[LessEqual[z, -0.0285], N[(x * 3.0), $MachinePrecision], If[LessEqual[z, 8.6e+89], N[(y * 2.0), $MachinePrecision], z]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.1 \cdot 10^{+105}:\\
\;\;\;\;z\\
\mathbf{elif}\;z \leq -0.0285:\\
\;\;\;\;x \cdot 3\\
\mathbf{elif}\;z \leq 8.6 \cdot 10^{+89}:\\
\;\;\;\;y \cdot 2\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if z < -4.1000000000000002e105 or 8.6000000000000003e89 < z Initial program 100.0%
+-commutative100.0%
associate-+l+100.0%
remove-double-neg100.0%
unsub-neg100.0%
+-commutative100.0%
+-commutative100.0%
associate-+l+100.0%
associate-+r+100.0%
associate-+r+100.0%
distribute-neg-in100.0%
distribute-neg-out100.0%
neg-mul-1100.0%
count-2100.0%
distribute-lft-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-rgt-out100.0%
distribute-neg-out100.0%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around 0 85.3%
Taylor expanded in z around inf 70.4%
if -4.1000000000000002e105 < z < -0.028500000000000001Initial 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 57.4%
if -0.028500000000000001 < z < 8.6000000000000003e89Initial 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%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
distribute-neg-out99.9%
fma-define100.0%
Simplified100.0%
Taylor expanded in x around 0 66.4%
Taylor expanded in z around 0 55.7%
*-commutative55.7%
Simplified55.7%
Final simplification60.5%
(FPCore (x y z) :precision binary64 (if (or (<= x -1.85e+21) (not (<= x 4.8e+47))) (- z (* x -3.0)) (+ z (* 2.0 (+ x y)))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -1.85e+21) || !(x <= 4.8e+47)) {
tmp = z - (x * -3.0);
} else {
tmp = z + (2.0 * (x + y));
}
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 ((x <= (-1.85d+21)) .or. (.not. (x <= 4.8d+47))) then
tmp = z - (x * (-3.0d0))
else
tmp = z + (2.0d0 * (x + y))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -1.85e+21) || !(x <= 4.8e+47)) {
tmp = z - (x * -3.0);
} else {
tmp = z + (2.0 * (x + y));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -1.85e+21) or not (x <= 4.8e+47): tmp = z - (x * -3.0) else: tmp = z + (2.0 * (x + y)) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -1.85e+21) || !(x <= 4.8e+47)) tmp = Float64(z - Float64(x * -3.0)); else tmp = Float64(z + Float64(2.0 * Float64(x + y))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -1.85e+21) || ~((x <= 4.8e+47))) tmp = z - (x * -3.0); else tmp = z + (2.0 * (x + y)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -1.85e+21], N[Not[LessEqual[x, 4.8e+47]], $MachinePrecision]], N[(z - N[(x * -3.0), $MachinePrecision]), $MachinePrecision], N[(z + N[(2.0 * N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.85 \cdot 10^{+21} \lor \neg \left(x \leq 4.8 \cdot 10^{+47}\right):\\
\;\;\;\;z - x \cdot -3\\
\mathbf{else}:\\
\;\;\;\;z + 2 \cdot \left(x + y\right)\\
\end{array}
\end{array}
if x < -1.85e21 or 4.80000000000000037e47 < x 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%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
distribute-neg-out99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around inf 83.5%
if -1.85e21 < x < 4.80000000000000037e47Initial 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 x around 0 94.7%
Final simplification89.9%
(FPCore (x y z) :precision binary64 (if (<= x -7.2e+20) (+ x (+ z (* x 2.0))) (if (<= x 3e+49) (+ z (* 2.0 (+ x y))) (- z (* x -3.0)))))
double code(double x, double y, double z) {
double tmp;
if (x <= -7.2e+20) {
tmp = x + (z + (x * 2.0));
} else if (x <= 3e+49) {
tmp = z + (2.0 * (x + y));
} 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 (x <= (-7.2d+20)) then
tmp = x + (z + (x * 2.0d0))
else if (x <= 3d+49) then
tmp = z + (2.0d0 * (x + y))
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 (x <= -7.2e+20) {
tmp = x + (z + (x * 2.0));
} else if (x <= 3e+49) {
tmp = z + (2.0 * (x + y));
} else {
tmp = z - (x * -3.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -7.2e+20: tmp = x + (z + (x * 2.0)) elif x <= 3e+49: tmp = z + (2.0 * (x + y)) else: tmp = z - (x * -3.0) return tmp
function code(x, y, z) tmp = 0.0 if (x <= -7.2e+20) tmp = Float64(x + Float64(z + Float64(x * 2.0))); elseif (x <= 3e+49) tmp = Float64(z + Float64(2.0 * Float64(x + y))); else tmp = Float64(z - Float64(x * -3.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -7.2e+20) tmp = x + (z + (x * 2.0)); elseif (x <= 3e+49) tmp = z + (2.0 * (x + y)); else tmp = z - (x * -3.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -7.2e+20], N[(x + N[(z + N[(x * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3e+49], N[(z + N[(2.0 * N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(z - N[(x * -3.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7.2 \cdot 10^{+20}:\\
\;\;\;\;x + \left(z + x \cdot 2\right)\\
\mathbf{elif}\;x \leq 3 \cdot 10^{+49}:\\
\;\;\;\;z + 2 \cdot \left(x + y\right)\\
\mathbf{else}:\\
\;\;\;\;z - x \cdot -3\\
\end{array}
\end{array}
if x < -7.2e20Initial program 99.8%
Taylor expanded in y around 0 83.8%
if -7.2e20 < x < 3.0000000000000002e49Initial 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 x around 0 94.7%
if 3.0000000000000002e49 < x 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%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
distribute-neg-out99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around inf 83.3%
Final simplification90.0%
(FPCore (x y z) :precision binary64 (if (or (<= x -7.2e+20) (not (<= x 9.5e+48))) (- z (* x -3.0)) (- z (* y -2.0))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -7.2e+20) || !(x <= 9.5e+48)) {
tmp = z - (x * -3.0);
} else {
tmp = z - (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 ((x <= (-7.2d+20)) .or. (.not. (x <= 9.5d+48))) then
tmp = z - (x * (-3.0d0))
else
tmp = z - (y * (-2.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -7.2e+20) || !(x <= 9.5e+48)) {
tmp = z - (x * -3.0);
} else {
tmp = z - (y * -2.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -7.2e+20) or not (x <= 9.5e+48): tmp = z - (x * -3.0) else: tmp = z - (y * -2.0) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -7.2e+20) || !(x <= 9.5e+48)) tmp = Float64(z - Float64(x * -3.0)); else tmp = Float64(z - Float64(y * -2.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -7.2e+20) || ~((x <= 9.5e+48))) tmp = z - (x * -3.0); else tmp = z - (y * -2.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -7.2e+20], N[Not[LessEqual[x, 9.5e+48]], $MachinePrecision]], N[(z - N[(x * -3.0), $MachinePrecision]), $MachinePrecision], N[(z - N[(y * -2.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7.2 \cdot 10^{+20} \lor \neg \left(x \leq 9.5 \cdot 10^{+48}\right):\\
\;\;\;\;z - x \cdot -3\\
\mathbf{else}:\\
\;\;\;\;z - y \cdot -2\\
\end{array}
\end{array}
if x < -7.2e20 or 9.4999999999999997e48 < x 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%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
distribute-neg-out99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around inf 83.5%
if -7.2e20 < x < 9.4999999999999997e48Initial program 100.0%
+-commutative100.0%
associate-+l+100.0%
remove-double-neg100.0%
unsub-neg100.0%
+-commutative100.0%
+-commutative100.0%
associate-+l+100.0%
associate-+r+100.0%
associate-+r+100.0%
distribute-neg-in100.0%
distribute-neg-out100.0%
neg-mul-1100.0%
count-2100.0%
distribute-lft-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-rgt-out100.0%
distribute-neg-out100.0%
fma-define100.0%
Simplified100.0%
Taylor expanded in x around 0 93.9%
Final simplification89.5%
(FPCore (x y z) :precision binary64 (if (or (<= y -2.7e+99) (not (<= y 1.2e+131))) (* 2.0 (+ x y)) (- z (* x -3.0))))
double code(double x, double y, double z) {
double tmp;
if ((y <= -2.7e+99) || !(y <= 1.2e+131)) {
tmp = 2.0 * (x + y);
} 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 <= (-2.7d+99)) .or. (.not. (y <= 1.2d+131))) then
tmp = 2.0d0 * (x + y)
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 <= -2.7e+99) || !(y <= 1.2e+131)) {
tmp = 2.0 * (x + y);
} else {
tmp = z - (x * -3.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (y <= -2.7e+99) or not (y <= 1.2e+131): tmp = 2.0 * (x + y) else: tmp = z - (x * -3.0) return tmp
function code(x, y, z) tmp = 0.0 if ((y <= -2.7e+99) || !(y <= 1.2e+131)) tmp = Float64(2.0 * Float64(x + y)); else tmp = Float64(z - Float64(x * -3.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((y <= -2.7e+99) || ~((y <= 1.2e+131))) tmp = 2.0 * (x + y); else tmp = z - (x * -3.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[y, -2.7e+99], N[Not[LessEqual[y, 1.2e+131]], $MachinePrecision]], N[(2.0 * N[(x + y), $MachinePrecision]), $MachinePrecision], N[(z - N[(x * -3.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.7 \cdot 10^{+99} \lor \neg \left(y \leq 1.2 \cdot 10^{+131}\right):\\
\;\;\;\;2 \cdot \left(x + y\right)\\
\mathbf{else}:\\
\;\;\;\;z - x \cdot -3\\
\end{array}
\end{array}
if y < -2.69999999999999989e99 or 1.2e131 < 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 x around 0 88.5%
Taylor expanded in z around 0 73.4%
+-commutative73.4%
Simplified73.4%
if -2.69999999999999989e99 < y < 1.2e131Initial 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%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
distribute-neg-out99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around inf 82.7%
Final simplification79.3%
(FPCore (x y z) :precision binary64 (if (or (<= x -0.00082) (not (<= x 85000.0))) (* x 3.0) z))
double code(double x, double y, double z) {
double tmp;
if ((x <= -0.00082) || !(x <= 85000.0)) {
tmp = x * 3.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 ((x <= (-0.00082d0)) .or. (.not. (x <= 85000.0d0))) then
tmp = x * 3.0d0
else
tmp = z
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -0.00082) || !(x <= 85000.0)) {
tmp = x * 3.0;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -0.00082) or not (x <= 85000.0): tmp = x * 3.0 else: tmp = z return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -0.00082) || !(x <= 85000.0)) tmp = Float64(x * 3.0); else tmp = z; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -0.00082) || ~((x <= 85000.0))) tmp = x * 3.0; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -0.00082], N[Not[LessEqual[x, 85000.0]], $MachinePrecision]], N[(x * 3.0), $MachinePrecision], z]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.00082 \lor \neg \left(x \leq 85000\right):\\
\;\;\;\;x \cdot 3\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if x < -8.1999999999999998e-4 or 85000 < x Initial program 99.9%
associate-+l+99.8%
associate-+l+99.9%
+-commutative99.9%
count-299.9%
+-commutative99.9%
+-commutative99.9%
Simplified99.9%
Taylor expanded in x around inf 59.0%
if -8.1999999999999998e-4 < x < 85000Initial program 100.0%
+-commutative100.0%
associate-+l+100.0%
remove-double-neg100.0%
unsub-neg100.0%
+-commutative100.0%
+-commutative100.0%
associate-+l+100.0%
associate-+r+100.0%
associate-+r+100.0%
distribute-neg-in100.0%
distribute-neg-out100.0%
neg-mul-1100.0%
count-2100.0%
distribute-lft-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-rgt-out100.0%
distribute-neg-out100.0%
fma-define100.0%
Simplified100.0%
Taylor expanded in x around 0 96.7%
Taylor expanded in z around inf 43.0%
Final simplification50.8%
(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%
+-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%
neg-mul-199.9%
count-299.9%
distribute-lft-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-rgt-out99.9%
distribute-neg-out99.9%
fma-define100.0%
Simplified100.0%
Taylor expanded in x around 0 69.1%
Taylor expanded in z around inf 32.3%
herbie shell --seed 2024146
(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))