
(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 12 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
(let* ((t_0 (+ x (* y 2.0))))
(if (<= z -5.2e+79)
z
(if (<= z -1.65e-39)
t_0
(if (<= z -1.16e-125) (* x 3.0) (if (<= z 2.1e+56) t_0 z))))))
double code(double x, double y, double z) {
double t_0 = x + (y * 2.0);
double tmp;
if (z <= -5.2e+79) {
tmp = z;
} else if (z <= -1.65e-39) {
tmp = t_0;
} else if (z <= -1.16e-125) {
tmp = x * 3.0;
} else if (z <= 2.1e+56) {
tmp = t_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) :: t_0
real(8) :: tmp
t_0 = x + (y * 2.0d0)
if (z <= (-5.2d+79)) then
tmp = z
else if (z <= (-1.65d-39)) then
tmp = t_0
else if (z <= (-1.16d-125)) then
tmp = x * 3.0d0
else if (z <= 2.1d+56) then
tmp = t_0
else
tmp = z
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = x + (y * 2.0);
double tmp;
if (z <= -5.2e+79) {
tmp = z;
} else if (z <= -1.65e-39) {
tmp = t_0;
} else if (z <= -1.16e-125) {
tmp = x * 3.0;
} else if (z <= 2.1e+56) {
tmp = t_0;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z): t_0 = x + (y * 2.0) tmp = 0 if z <= -5.2e+79: tmp = z elif z <= -1.65e-39: tmp = t_0 elif z <= -1.16e-125: tmp = x * 3.0 elif z <= 2.1e+56: tmp = t_0 else: tmp = z return tmp
function code(x, y, z) t_0 = Float64(x + Float64(y * 2.0)) tmp = 0.0 if (z <= -5.2e+79) tmp = z; elseif (z <= -1.65e-39) tmp = t_0; elseif (z <= -1.16e-125) tmp = Float64(x * 3.0); elseif (z <= 2.1e+56) tmp = t_0; else tmp = z; end return tmp end
function tmp_2 = code(x, y, z) t_0 = x + (y * 2.0); tmp = 0.0; if (z <= -5.2e+79) tmp = z; elseif (z <= -1.65e-39) tmp = t_0; elseif (z <= -1.16e-125) tmp = x * 3.0; elseif (z <= 2.1e+56) tmp = t_0; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x + N[(y * 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -5.2e+79], z, If[LessEqual[z, -1.65e-39], t$95$0, If[LessEqual[z, -1.16e-125], N[(x * 3.0), $MachinePrecision], If[LessEqual[z, 2.1e+56], t$95$0, z]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + y \cdot 2\\
\mathbf{if}\;z \leq -5.2 \cdot 10^{+79}:\\
\;\;\;\;z\\
\mathbf{elif}\;z \leq -1.65 \cdot 10^{-39}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -1.16 \cdot 10^{-125}:\\
\;\;\;\;x \cdot 3\\
\mathbf{elif}\;z \leq 2.1 \cdot 10^{+56}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if z < -5.20000000000000029e79 or 2.10000000000000017e56 < z 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 z around inf 66.9%
if -5.20000000000000029e79 < z < -1.64999999999999992e-39 or -1.15999999999999995e-125 < z < 2.10000000000000017e56Initial 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 0 90.6%
Taylor expanded in x around 0 59.5%
if -1.64999999999999992e-39 < z < -1.15999999999999995e-125Initial program 99.6%
associate-+l+99.7%
associate-+l+99.7%
+-commutative99.7%
count-299.7%
+-commutative99.7%
+-commutative99.7%
Simplified99.7%
Taylor expanded in x around inf 70.6%
Final simplification63.2%
(FPCore (x y z)
:precision binary64
(if (<= z -2.25e+78)
z
(if (<= z -8.8e-41)
(* y 2.0)
(if (<= z -8e-126) (* x 3.0) (if (<= z 3.6e+24) (* y 2.0) z)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -2.25e+78) {
tmp = z;
} else if (z <= -8.8e-41) {
tmp = y * 2.0;
} else if (z <= -8e-126) {
tmp = x * 3.0;
} else if (z <= 3.6e+24) {
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 <= (-2.25d+78)) then
tmp = z
else if (z <= (-8.8d-41)) then
tmp = y * 2.0d0
else if (z <= (-8d-126)) then
tmp = x * 3.0d0
else if (z <= 3.6d+24) 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 <= -2.25e+78) {
tmp = z;
} else if (z <= -8.8e-41) {
tmp = y * 2.0;
} else if (z <= -8e-126) {
tmp = x * 3.0;
} else if (z <= 3.6e+24) {
tmp = y * 2.0;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -2.25e+78: tmp = z elif z <= -8.8e-41: tmp = y * 2.0 elif z <= -8e-126: tmp = x * 3.0 elif z <= 3.6e+24: tmp = y * 2.0 else: tmp = z return tmp
function code(x, y, z) tmp = 0.0 if (z <= -2.25e+78) tmp = z; elseif (z <= -8.8e-41) tmp = Float64(y * 2.0); elseif (z <= -8e-126) tmp = Float64(x * 3.0); elseif (z <= 3.6e+24) tmp = Float64(y * 2.0); else tmp = z; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -2.25e+78) tmp = z; elseif (z <= -8.8e-41) tmp = y * 2.0; elseif (z <= -8e-126) tmp = x * 3.0; elseif (z <= 3.6e+24) tmp = y * 2.0; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -2.25e+78], z, If[LessEqual[z, -8.8e-41], N[(y * 2.0), $MachinePrecision], If[LessEqual[z, -8e-126], N[(x * 3.0), $MachinePrecision], If[LessEqual[z, 3.6e+24], N[(y * 2.0), $MachinePrecision], z]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.25 \cdot 10^{+78}:\\
\;\;\;\;z\\
\mathbf{elif}\;z \leq -8.8 \cdot 10^{-41}:\\
\;\;\;\;y \cdot 2\\
\mathbf{elif}\;z \leq -8 \cdot 10^{-126}:\\
\;\;\;\;x \cdot 3\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{+24}:\\
\;\;\;\;y \cdot 2\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if z < -2.25e78 or 3.59999999999999983e24 < z 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 z around inf 64.5%
if -2.25e78 < z < -8.7999999999999999e-41 or -7.9999999999999996e-126 < z < 3.59999999999999983e24Initial 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 y around inf 55.5%
if -8.7999999999999999e-41 < z < -7.9999999999999996e-126Initial program 99.6%
associate-+l+99.7%
associate-+l+99.7%
+-commutative99.7%
count-299.7%
+-commutative99.7%
+-commutative99.7%
Simplified99.7%
Taylor expanded in x around inf 70.6%
Final simplification60.4%
(FPCore (x y z) :precision binary64 (if (<= x -5.5e+124) (+ (+ z x) (* x 2.0)) (if (<= x 4.3e+57) (- z (* y -2.0)) (- (* x 3.0) (* y -2.0)))))
double code(double x, double y, double z) {
double tmp;
if (x <= -5.5e+124) {
tmp = (z + x) + (x * 2.0);
} else if (x <= 4.3e+57) {
tmp = z - (y * -2.0);
} else {
tmp = (x * 3.0) - (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 <= (-5.5d+124)) then
tmp = (z + x) + (x * 2.0d0)
else if (x <= 4.3d+57) then
tmp = z - (y * (-2.0d0))
else
tmp = (x * 3.0d0) - (y * (-2.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -5.5e+124) {
tmp = (z + x) + (x * 2.0);
} else if (x <= 4.3e+57) {
tmp = z - (y * -2.0);
} else {
tmp = (x * 3.0) - (y * -2.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -5.5e+124: tmp = (z + x) + (x * 2.0) elif x <= 4.3e+57: tmp = z - (y * -2.0) else: tmp = (x * 3.0) - (y * -2.0) return tmp
function code(x, y, z) tmp = 0.0 if (x <= -5.5e+124) tmp = Float64(Float64(z + x) + Float64(x * 2.0)); elseif (x <= 4.3e+57) tmp = Float64(z - Float64(y * -2.0)); else tmp = Float64(Float64(x * 3.0) - Float64(y * -2.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -5.5e+124) tmp = (z + x) + (x * 2.0); elseif (x <= 4.3e+57) tmp = z - (y * -2.0); else tmp = (x * 3.0) - (y * -2.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -5.5e+124], N[(N[(z + x), $MachinePrecision] + N[(x * 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 4.3e+57], N[(z - N[(y * -2.0), $MachinePrecision]), $MachinePrecision], N[(N[(x * 3.0), $MachinePrecision] - N[(y * -2.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.5 \cdot 10^{+124}:\\
\;\;\;\;\left(z + x\right) + x \cdot 2\\
\mathbf{elif}\;x \leq 4.3 \cdot 10^{+57}:\\
\;\;\;\;z - y \cdot -2\\
\mathbf{else}:\\
\;\;\;\;x \cdot 3 - y \cdot -2\\
\end{array}
\end{array}
if x < -5.49999999999999977e124Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
+-commutative99.8%
count-299.8%
+-commutative99.8%
+-commutative99.8%
Simplified99.8%
Taylor expanded in x around inf 85.9%
if -5.49999999999999977e124 < x < 4.30000000000000033e57Initial 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 90.2%
if 4.30000000000000033e57 < x Initial program 99.8%
+-commutative99.8%
associate-+l+99.8%
remove-double-neg99.8%
unsub-neg99.8%
+-commutative99.8%
+-commutative99.8%
associate-+l+99.8%
associate-+r+99.8%
associate-+r+99.8%
distribute-neg-in99.8%
distribute-neg-out99.8%
neg-mul-199.8%
count-299.8%
distribute-lft-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
distribute-rgt-out99.8%
distribute-neg-out99.8%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around 0 99.8%
Taylor expanded in z around 0 91.6%
Final simplification89.9%
(FPCore (x y z) :precision binary64 (if (<= x -4e+123) (+ (+ z x) (* x 2.0)) (if (<= x 2.5e+59) (- z (* y -2.0)) (+ x (* 2.0 (+ x y))))))
double code(double x, double y, double z) {
double tmp;
if (x <= -4e+123) {
tmp = (z + x) + (x * 2.0);
} else if (x <= 2.5e+59) {
tmp = z - (y * -2.0);
} else {
tmp = x + (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 <= (-4d+123)) then
tmp = (z + x) + (x * 2.0d0)
else if (x <= 2.5d+59) then
tmp = z - (y * (-2.0d0))
else
tmp = x + (2.0d0 * (x + y))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -4e+123) {
tmp = (z + x) + (x * 2.0);
} else if (x <= 2.5e+59) {
tmp = z - (y * -2.0);
} else {
tmp = x + (2.0 * (x + y));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -4e+123: tmp = (z + x) + (x * 2.0) elif x <= 2.5e+59: tmp = z - (y * -2.0) else: tmp = x + (2.0 * (x + y)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= -4e+123) tmp = Float64(Float64(z + x) + Float64(x * 2.0)); elseif (x <= 2.5e+59) tmp = Float64(z - Float64(y * -2.0)); else tmp = Float64(x + Float64(2.0 * Float64(x + y))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -4e+123) tmp = (z + x) + (x * 2.0); elseif (x <= 2.5e+59) tmp = z - (y * -2.0); else tmp = x + (2.0 * (x + y)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -4e+123], N[(N[(z + x), $MachinePrecision] + N[(x * 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 2.5e+59], N[(z - N[(y * -2.0), $MachinePrecision]), $MachinePrecision], N[(x + N[(2.0 * N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -4 \cdot 10^{+123}:\\
\;\;\;\;\left(z + x\right) + x \cdot 2\\
\mathbf{elif}\;x \leq 2.5 \cdot 10^{+59}:\\
\;\;\;\;z - y \cdot -2\\
\mathbf{else}:\\
\;\;\;\;x + 2 \cdot \left(x + y\right)\\
\end{array}
\end{array}
if x < -3.99999999999999991e123Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
+-commutative99.8%
count-299.8%
+-commutative99.8%
+-commutative99.8%
Simplified99.8%
Taylor expanded in x around inf 85.9%
if -3.99999999999999991e123 < x < 2.4999999999999999e59Initial 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 90.2%
if 2.4999999999999999e59 < x Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
+-commutative99.8%
count-299.8%
+-commutative99.8%
+-commutative99.8%
Simplified99.8%
Taylor expanded in z around 0 91.6%
Final simplification89.9%
(FPCore (x y z) :precision binary64 (if (<= x -1.55e+127) (- z (* x -3.0)) (if (<= x 3.8e+59) (- z (* y -2.0)) (+ x (* 2.0 (+ x y))))))
double code(double x, double y, double z) {
double tmp;
if (x <= -1.55e+127) {
tmp = z - (x * -3.0);
} else if (x <= 3.8e+59) {
tmp = z - (y * -2.0);
} else {
tmp = x + (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.55d+127)) then
tmp = z - (x * (-3.0d0))
else if (x <= 3.8d+59) then
tmp = z - (y * (-2.0d0))
else
tmp = x + (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.55e+127) {
tmp = z - (x * -3.0);
} else if (x <= 3.8e+59) {
tmp = z - (y * -2.0);
} else {
tmp = x + (2.0 * (x + y));
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -1.55e+127: tmp = z - (x * -3.0) elif x <= 3.8e+59: tmp = z - (y * -2.0) else: tmp = x + (2.0 * (x + y)) return tmp
function code(x, y, z) tmp = 0.0 if (x <= -1.55e+127) tmp = Float64(z - Float64(x * -3.0)); elseif (x <= 3.8e+59) tmp = Float64(z - Float64(y * -2.0)); else tmp = Float64(x + Float64(2.0 * Float64(x + y))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -1.55e+127) tmp = z - (x * -3.0); elseif (x <= 3.8e+59) tmp = z - (y * -2.0); else tmp = x + (2.0 * (x + y)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -1.55e+127], N[(z - N[(x * -3.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.8e+59], N[(z - N[(y * -2.0), $MachinePrecision]), $MachinePrecision], N[(x + N[(2.0 * N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.55 \cdot 10^{+127}:\\
\;\;\;\;z - x \cdot -3\\
\mathbf{elif}\;x \leq 3.8 \cdot 10^{+59}:\\
\;\;\;\;z - y \cdot -2\\
\mathbf{else}:\\
\;\;\;\;x + 2 \cdot \left(x + y\right)\\
\end{array}
\end{array}
if x < -1.5500000000000001e127Initial program 99.8%
+-commutative99.8%
associate-+l+99.8%
remove-double-neg99.8%
unsub-neg99.8%
+-commutative99.8%
+-commutative99.8%
associate-+l+99.8%
associate-+r+99.8%
associate-+r+99.7%
distribute-neg-in99.7%
distribute-neg-out99.7%
neg-mul-199.7%
count-299.7%
distribute-lft-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
distribute-rgt-out99.7%
distribute-neg-out99.7%
fma-define99.8%
Simplified99.8%
Taylor expanded in x around inf 85.8%
*-commutative85.8%
Simplified85.8%
if -1.5500000000000001e127 < x < 3.8000000000000001e59Initial 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 90.2%
if 3.8000000000000001e59 < x Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
+-commutative99.8%
count-299.8%
+-commutative99.8%
+-commutative99.8%
Simplified99.8%
Taylor expanded in z around 0 91.6%
Final simplification89.9%
(FPCore (x y z) :precision binary64 (if (or (<= x -4e+123) (not (<= x 2.9e+60))) (- z (* x -3.0)) (- z (* y -2.0))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -4e+123) || !(x <= 2.9e+60)) {
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 <= (-4d+123)) .or. (.not. (x <= 2.9d+60))) 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 <= -4e+123) || !(x <= 2.9e+60)) {
tmp = z - (x * -3.0);
} else {
tmp = z - (y * -2.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -4e+123) or not (x <= 2.9e+60): tmp = z - (x * -3.0) else: tmp = z - (y * -2.0) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -4e+123) || !(x <= 2.9e+60)) 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 <= -4e+123) || ~((x <= 2.9e+60))) tmp = z - (x * -3.0); else tmp = z - (y * -2.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -4e+123], N[Not[LessEqual[x, 2.9e+60]], $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 -4 \cdot 10^{+123} \lor \neg \left(x \leq 2.9 \cdot 10^{+60}\right):\\
\;\;\;\;z - x \cdot -3\\
\mathbf{else}:\\
\;\;\;\;z - y \cdot -2\\
\end{array}
\end{array}
if x < -3.99999999999999991e123 or 2.9e60 < x Initial program 99.8%
+-commutative99.8%
associate-+l+99.8%
remove-double-neg99.8%
unsub-neg99.8%
+-commutative99.8%
+-commutative99.8%
associate-+l+99.8%
associate-+r+99.8%
associate-+r+99.8%
distribute-neg-in99.8%
distribute-neg-out99.8%
neg-mul-199.8%
count-299.8%
distribute-lft-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
distribute-rgt-out99.8%
distribute-neg-out99.8%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around inf 80.0%
*-commutative80.0%
Simplified80.0%
if -3.99999999999999991e123 < x < 2.9e60Initial 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 90.3%
Final simplification86.7%
(FPCore (x y z) :precision binary64 (if (or (<= x -4e+192) (not (<= x 2.5e+167))) (* x 3.0) (- z (* y -2.0))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -4e+192) || !(x <= 2.5e+167)) {
tmp = 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 <= (-4d+192)) .or. (.not. (x <= 2.5d+167))) then
tmp = 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 <= -4e+192) || !(x <= 2.5e+167)) {
tmp = x * 3.0;
} else {
tmp = z - (y * -2.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -4e+192) or not (x <= 2.5e+167): tmp = x * 3.0 else: tmp = z - (y * -2.0) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -4e+192) || !(x <= 2.5e+167)) tmp = 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 <= -4e+192) || ~((x <= 2.5e+167))) tmp = x * 3.0; else tmp = z - (y * -2.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -4e+192], N[Not[LessEqual[x, 2.5e+167]], $MachinePrecision]], N[(x * 3.0), $MachinePrecision], N[(z - N[(y * -2.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -4 \cdot 10^{+192} \lor \neg \left(x \leq 2.5 \cdot 10^{+167}\right):\\
\;\;\;\;x \cdot 3\\
\mathbf{else}:\\
\;\;\;\;z - y \cdot -2\\
\end{array}
\end{array}
if x < -4.00000000000000016e192 or 2.4999999999999998e167 < x Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
+-commutative99.8%
count-299.8%
+-commutative99.8%
+-commutative99.8%
Simplified99.8%
Taylor expanded in x around inf 78.4%
if -4.00000000000000016e192 < x < 2.4999999999999998e167Initial program 100.0%
+-commutative100.0%
associate-+l+100.0%
remove-double-neg100.0%
unsub-neg100.0%
+-commutative100.0%
+-commutative100.0%
associate-+l+99.9%
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 84.5%
Final simplification83.1%
(FPCore (x y z) :precision binary64 (if (<= z -4.8e+76) z (if (<= z 3.8e+24) (* y 2.0) z)))
double code(double x, double y, double z) {
double tmp;
if (z <= -4.8e+76) {
tmp = z;
} else if (z <= 3.8e+24) {
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.8d+76)) then
tmp = z
else if (z <= 3.8d+24) 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.8e+76) {
tmp = z;
} else if (z <= 3.8e+24) {
tmp = y * 2.0;
} else {
tmp = z;
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -4.8e+76: tmp = z elif z <= 3.8e+24: tmp = y * 2.0 else: tmp = z return tmp
function code(x, y, z) tmp = 0.0 if (z <= -4.8e+76) tmp = z; elseif (z <= 3.8e+24) 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.8e+76) tmp = z; elseif (z <= 3.8e+24) tmp = y * 2.0; else tmp = z; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -4.8e+76], z, If[LessEqual[z, 3.8e+24], N[(y * 2.0), $MachinePrecision], z]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.8 \cdot 10^{+76}:\\
\;\;\;\;z\\
\mathbf{elif}\;z \leq 3.8 \cdot 10^{+24}:\\
\;\;\;\;y \cdot 2\\
\mathbf{else}:\\
\;\;\;\;z\\
\end{array}
\end{array}
if z < -4.8e76 or 3.80000000000000015e24 < z 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 z around inf 64.5%
if -4.8e76 < z < 3.80000000000000015e24Initial 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 y around inf 50.3%
Final simplification56.3%
(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 33.1%
(FPCore (x y z) :precision binary64 x)
double code(double x, double y, double z) {
return 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
end function
public static double code(double x, double y, double z) {
return x;
}
def code(x, y, z): return x
function code(x, y, z) return x end
function tmp = code(x, y, z) tmp = x; end
code[x_, y_, z_] := x
\begin{array}{l}
\\
x
\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 0 67.9%
Taylor expanded in x around 0 42.9%
Taylor expanded in x around inf 7.6%
herbie shell --seed 2024135
(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))