
(FPCore (x y z) :precision binary64 (/ (* x (- y z)) y))
double code(double x, double y, double z) {
return (x * (y - z)) / y;
}
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 - z)) / y
end function
public static double code(double x, double y, double z) {
return (x * (y - z)) / y;
}
def code(x, y, z): return (x * (y - z)) / y
function code(x, y, z) return Float64(Float64(x * Float64(y - z)) / y) end
function tmp = code(x, y, z) tmp = (x * (y - z)) / y; end
code[x_, y_, z_] := N[(N[(x * N[(y - z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot \left(y - z\right)}{y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z) :precision binary64 (/ (* x (- y z)) y))
double code(double x, double y, double z) {
return (x * (y - z)) / y;
}
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 - z)) / y
end function
public static double code(double x, double y, double z) {
return (x * (y - z)) / y;
}
def code(x, y, z): return (x * (y - z)) / y
function code(x, y, z) return Float64(Float64(x * Float64(y - z)) / y) end
function tmp = code(x, y, z) tmp = (x * (y - z)) / y; end
code[x_, y_, z_] := N[(N[(x * N[(y - z), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot \left(y - z\right)}{y}
\end{array}
(FPCore (x y z) :precision binary64 (* x (- 1.0 (/ z y))))
double code(double x, double y, double z) {
return x * (1.0 - (z / y));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = x * (1.0d0 - (z / y))
end function
public static double code(double x, double y, double z) {
return x * (1.0 - (z / y));
}
def code(x, y, z): return x * (1.0 - (z / y))
function code(x, y, z) return Float64(x * Float64(1.0 - Float64(z / y))) end
function tmp = code(x, y, z) tmp = x * (1.0 - (z / y)); end
code[x_, y_, z_] := N[(x * N[(1.0 - N[(z / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 - \frac{z}{y}\right)
\end{array}
Initial program 84.3%
associate-*r/97.3%
div-sub97.3%
*-inverses97.3%
Simplified97.3%
Final simplification97.3%
(FPCore (x y z) :precision binary64 (if (<= y -6e+100) x (if (<= y 5.2e+31) (* x (/ (- z) y)) x)))
double code(double x, double y, double z) {
double tmp;
if (y <= -6e+100) {
tmp = x;
} else if (y <= 5.2e+31) {
tmp = x * (-z / y);
} else {
tmp = x;
}
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 <= (-6d+100)) then
tmp = x
else if (y <= 5.2d+31) then
tmp = x * (-z / y)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -6e+100) {
tmp = x;
} else if (y <= 5.2e+31) {
tmp = x * (-z / y);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -6e+100: tmp = x elif y <= 5.2e+31: tmp = x * (-z / y) else: tmp = x return tmp
function code(x, y, z) tmp = 0.0 if (y <= -6e+100) tmp = x; elseif (y <= 5.2e+31) tmp = Float64(x * Float64(Float64(-z) / y)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -6e+100) tmp = x; elseif (y <= 5.2e+31) tmp = x * (-z / y); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -6e+100], x, If[LessEqual[y, 5.2e+31], N[(x * N[((-z) / y), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6 \cdot 10^{+100}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 5.2 \cdot 10^{+31}:\\
\;\;\;\;x \cdot \frac{-z}{y}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -5.99999999999999971e100 or 5.2e31 < y Initial program 70.0%
associate-*r/99.9%
div-sub100.0%
*-inverses100.0%
Simplified100.0%
Taylor expanded in z around 0 82.6%
if -5.99999999999999971e100 < y < 5.2e31Initial program 93.2%
associate-*r/95.6%
div-sub95.6%
*-inverses95.6%
Simplified95.6%
Taylor expanded in z around inf 70.0%
mul-1-neg70.0%
distribute-frac-neg70.0%
Simplified70.0%
Final simplification74.8%
(FPCore (x y z) :precision binary64 (if (<= y -1.55e+100) x (if (<= y 3.8e+35) (* z (/ (- x) y)) x)))
double code(double x, double y, double z) {
double tmp;
if (y <= -1.55e+100) {
tmp = x;
} else if (y <= 3.8e+35) {
tmp = z * (-x / y);
} else {
tmp = x;
}
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.55d+100)) then
tmp = x
else if (y <= 3.8d+35) then
tmp = z * (-x / y)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -1.55e+100) {
tmp = x;
} else if (y <= 3.8e+35) {
tmp = z * (-x / y);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -1.55e+100: tmp = x elif y <= 3.8e+35: tmp = z * (-x / y) else: tmp = x return tmp
function code(x, y, z) tmp = 0.0 if (y <= -1.55e+100) tmp = x; elseif (y <= 3.8e+35) tmp = Float64(z * Float64(Float64(-x) / y)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -1.55e+100) tmp = x; elseif (y <= 3.8e+35) tmp = z * (-x / y); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -1.55e+100], x, If[LessEqual[y, 3.8e+35], N[(z * N[((-x) / y), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.55 \cdot 10^{+100}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 3.8 \cdot 10^{+35}:\\
\;\;\;\;z \cdot \frac{-x}{y}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -1.55000000000000003e100 or 3.8e35 < y Initial program 70.0%
associate-*r/99.9%
div-sub100.0%
*-inverses100.0%
Simplified100.0%
Taylor expanded in z around 0 82.6%
if -1.55000000000000003e100 < y < 3.8e35Initial program 93.2%
associate-*r/95.6%
div-sub95.6%
*-inverses95.6%
Simplified95.6%
Taylor expanded in z around inf 72.0%
mul-1-neg72.0%
*-commutative72.0%
associate-*r/70.6%
distribute-rgt-neg-in70.6%
Simplified70.6%
Final simplification75.2%
(FPCore (x y z) :precision binary64 (if (<= y -1.5e+100) x (if (<= y 6e+30) (/ z (/ (- y) x)) x)))
double code(double x, double y, double z) {
double tmp;
if (y <= -1.5e+100) {
tmp = x;
} else if (y <= 6e+30) {
tmp = z / (-y / x);
} else {
tmp = x;
}
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.5d+100)) then
tmp = x
else if (y <= 6d+30) then
tmp = z / (-y / x)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -1.5e+100) {
tmp = x;
} else if (y <= 6e+30) {
tmp = z / (-y / x);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -1.5e+100: tmp = x elif y <= 6e+30: tmp = z / (-y / x) else: tmp = x return tmp
function code(x, y, z) tmp = 0.0 if (y <= -1.5e+100) tmp = x; elseif (y <= 6e+30) tmp = Float64(z / Float64(Float64(-y) / x)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -1.5e+100) tmp = x; elseif (y <= 6e+30) tmp = z / (-y / x); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -1.5e+100], x, If[LessEqual[y, 6e+30], N[(z / N[((-y) / x), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.5 \cdot 10^{+100}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 6 \cdot 10^{+30}:\\
\;\;\;\;\frac{z}{\frac{-y}{x}}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -1.49999999999999993e100 or 5.99999999999999956e30 < y Initial program 70.0%
associate-*r/99.9%
div-sub100.0%
*-inverses100.0%
Simplified100.0%
Taylor expanded in z around 0 82.6%
if -1.49999999999999993e100 < y < 5.99999999999999956e30Initial program 93.2%
associate-*r/95.6%
div-sub95.6%
*-inverses95.6%
Simplified95.6%
Taylor expanded in z around inf 70.0%
mul-1-neg70.0%
distribute-frac-neg70.0%
Simplified70.0%
*-commutative70.0%
associate-*l/72.0%
distribute-lft-neg-in72.0%
distribute-rgt-neg-out72.0%
associate-/l*71.2%
frac-2neg71.2%
add-sqr-sqrt35.4%
sqrt-unprod29.0%
sqr-neg29.0%
sqrt-unprod1.0%
add-sqr-sqrt2.0%
add-sqr-sqrt1.0%
sqrt-unprod29.8%
sqr-neg29.8%
sqrt-unprod35.5%
add-sqr-sqrt71.2%
*-un-lft-identity71.2%
*-un-lft-identity71.2%
Applied egg-rr71.2%
Final simplification75.6%
(FPCore (x y z) :precision binary64 (if (<= y -2e+100) x (if (<= y 3.05e+31) (/ (* z (- x)) y) x)))
double code(double x, double y, double z) {
double tmp;
if (y <= -2e+100) {
tmp = x;
} else if (y <= 3.05e+31) {
tmp = (z * -x) / y;
} else {
tmp = x;
}
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 <= (-2d+100)) then
tmp = x
else if (y <= 3.05d+31) then
tmp = (z * -x) / y
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -2e+100) {
tmp = x;
} else if (y <= 3.05e+31) {
tmp = (z * -x) / y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -2e+100: tmp = x elif y <= 3.05e+31: tmp = (z * -x) / y else: tmp = x return tmp
function code(x, y, z) tmp = 0.0 if (y <= -2e+100) tmp = x; elseif (y <= 3.05e+31) tmp = Float64(Float64(z * Float64(-x)) / y); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -2e+100) tmp = x; elseif (y <= 3.05e+31) tmp = (z * -x) / y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -2e+100], x, If[LessEqual[y, 3.05e+31], N[(N[(z * (-x)), $MachinePrecision] / y), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2 \cdot 10^{+100}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 3.05 \cdot 10^{+31}:\\
\;\;\;\;\frac{z \cdot \left(-x\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -2.00000000000000003e100 or 3.05000000000000005e31 < y Initial program 70.0%
associate-*r/99.9%
div-sub100.0%
*-inverses100.0%
Simplified100.0%
Taylor expanded in z around 0 82.6%
if -2.00000000000000003e100 < y < 3.05000000000000005e31Initial program 93.2%
Taylor expanded in y around 0 72.0%
associate-*r*72.0%
neg-mul-172.0%
*-commutative72.0%
Simplified72.0%
Final simplification76.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 84.3%
associate-*r/97.3%
div-sub97.3%
*-inverses97.3%
Simplified97.3%
Taylor expanded in z around 0 47.9%
Final simplification47.9%
(FPCore (x y z) :precision binary64 (if (< z -2.060202331921739e+104) (- x (/ (* z x) y)) (if (< z 1.6939766013828526e+213) (/ x (/ y (- y z))) (* (- y z) (/ x y)))))
double code(double x, double y, double z) {
double tmp;
if (z < -2.060202331921739e+104) {
tmp = x - ((z * x) / y);
} else if (z < 1.6939766013828526e+213) {
tmp = x / (y / (y - z));
} else {
tmp = (y - z) * (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 (z < (-2.060202331921739d+104)) then
tmp = x - ((z * x) / y)
else if (z < 1.6939766013828526d+213) then
tmp = x / (y / (y - z))
else
tmp = (y - z) * (x / y)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z < -2.060202331921739e+104) {
tmp = x - ((z * x) / y);
} else if (z < 1.6939766013828526e+213) {
tmp = x / (y / (y - z));
} else {
tmp = (y - z) * (x / y);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z < -2.060202331921739e+104: tmp = x - ((z * x) / y) elif z < 1.6939766013828526e+213: tmp = x / (y / (y - z)) else: tmp = (y - z) * (x / y) return tmp
function code(x, y, z) tmp = 0.0 if (z < -2.060202331921739e+104) tmp = Float64(x - Float64(Float64(z * x) / y)); elseif (z < 1.6939766013828526e+213) tmp = Float64(x / Float64(y / Float64(y - z))); else tmp = Float64(Float64(y - z) * Float64(x / y)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z < -2.060202331921739e+104) tmp = x - ((z * x) / y); elseif (z < 1.6939766013828526e+213) tmp = x / (y / (y - z)); else tmp = (y - z) * (x / y); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Less[z, -2.060202331921739e+104], N[(x - N[(N[(z * x), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[Less[z, 1.6939766013828526e+213], N[(x / N[(y / N[(y - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y - z), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z < -2.060202331921739 \cdot 10^{+104}:\\
\;\;\;\;x - \frac{z \cdot x}{y}\\
\mathbf{elif}\;z < 1.6939766013828526 \cdot 10^{+213}:\\
\;\;\;\;\frac{x}{\frac{y}{y - z}}\\
\mathbf{else}:\\
\;\;\;\;\left(y - z\right) \cdot \frac{x}{y}\\
\end{array}
\end{array}
herbie shell --seed 2023309
(FPCore (x y z)
:name "Diagrams.Backend.Cairo.Internal:setTexture from diagrams-cairo-1.3.0.3"
:precision binary64
:herbie-target
(if (< z -2.060202331921739e+104) (- x (/ (* z x) y)) (if (< z 1.6939766013828526e+213) (/ x (/ y (- y z))) (* (- y z) (/ x y))))
(/ (* x (- y z)) y))