
(FPCore (x y) :precision binary64 (- x (/ y 200.0)))
double code(double x, double y) {
return x - (y / 200.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x - (y / 200.0d0)
end function
public static double code(double x, double y) {
return x - (y / 200.0);
}
def code(x, y): return x - (y / 200.0)
function code(x, y) return Float64(x - Float64(y / 200.0)) end
function tmp = code(x, y) tmp = x - (y / 200.0); end
code[x_, y_] := N[(x - N[(y / 200.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{y}{200}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (- x (/ y 200.0)))
double code(double x, double y) {
return x - (y / 200.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x - (y / 200.0d0)
end function
public static double code(double x, double y) {
return x - (y / 200.0);
}
def code(x, y): return x - (y / 200.0)
function code(x, y) return Float64(x - Float64(y / 200.0)) end
function tmp = code(x, y) tmp = x - (y / 200.0); end
code[x_, y_] := N[(x - N[(y / 200.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{y}{200}
\end{array}
(FPCore (x y) :precision binary64 (- x (/ y 200.0)))
double code(double x, double y) {
return x - (y / 200.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x - (y / 200.0d0)
end function
public static double code(double x, double y) {
return x - (y / 200.0);
}
def code(x, y): return x - (y / 200.0)
function code(x, y) return Float64(x - Float64(y / 200.0)) end
function tmp = code(x, y) tmp = x - (y / 200.0); end
code[x_, y_] := N[(x - N[(y / 200.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{y}{200}
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (x y) :precision binary64 (if (or (<= y -2.6e+35) (not (<= y 5.8e+25))) (/ (- y) 200.0) x))
double code(double x, double y) {
double tmp;
if ((y <= -2.6e+35) || !(y <= 5.8e+25)) {
tmp = -y / 200.0;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-2.6d+35)) .or. (.not. (y <= 5.8d+25))) then
tmp = -y / 200.0d0
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -2.6e+35) || !(y <= 5.8e+25)) {
tmp = -y / 200.0;
} else {
tmp = x;
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -2.6e+35) or not (y <= 5.8e+25): tmp = -y / 200.0 else: tmp = x return tmp
function code(x, y) tmp = 0.0 if ((y <= -2.6e+35) || !(y <= 5.8e+25)) tmp = Float64(Float64(-y) / 200.0); else tmp = x; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -2.6e+35) || ~((y <= 5.8e+25))) tmp = -y / 200.0; else tmp = x; end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -2.6e+35], N[Not[LessEqual[y, 5.8e+25]], $MachinePrecision]], N[((-y) / 200.0), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.6 \cdot 10^{+35} \lor \neg \left(y \leq 5.8 \cdot 10^{+25}\right):\\
\;\;\;\;\frac{-y}{200}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -2.60000000000000007e35 or 5.7999999999999998e25 < y Initial program 100.0%
sub-neg100.0%
distribute-frac-neg100.0%
neg-mul-1100.0%
associate-/l*99.7%
associate-/r/99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 83.9%
*-commutative83.9%
metadata-eval83.9%
metadata-eval83.9%
distribute-rgt-neg-in83.9%
div-inv84.1%
distribute-neg-frac84.1%
Applied egg-rr84.1%
if -2.60000000000000007e35 < y < 5.7999999999999998e25Initial program 100.0%
sub-neg100.0%
distribute-frac-neg100.0%
neg-mul-1100.0%
associate-/l*99.9%
associate-/r/100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 77.4%
Final simplification80.7%
(FPCore (x y) :precision binary64 (if (or (<= y -5e+38) (not (<= y 4.5e+23))) (* y -0.005) x))
double code(double x, double y) {
double tmp;
if ((y <= -5e+38) || !(y <= 4.5e+23)) {
tmp = y * -0.005;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-5d+38)) .or. (.not. (y <= 4.5d+23))) then
tmp = y * (-0.005d0)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -5e+38) || !(y <= 4.5e+23)) {
tmp = y * -0.005;
} else {
tmp = x;
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -5e+38) or not (y <= 4.5e+23): tmp = y * -0.005 else: tmp = x return tmp
function code(x, y) tmp = 0.0 if ((y <= -5e+38) || !(y <= 4.5e+23)) tmp = Float64(y * -0.005); else tmp = x; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -5e+38) || ~((y <= 4.5e+23))) tmp = y * -0.005; else tmp = x; end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -5e+38], N[Not[LessEqual[y, 4.5e+23]], $MachinePrecision]], N[(y * -0.005), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5 \cdot 10^{+38} \lor \neg \left(y \leq 4.5 \cdot 10^{+23}\right):\\
\;\;\;\;y \cdot -0.005\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -4.9999999999999997e38 or 4.49999999999999979e23 < y Initial program 100.0%
sub-neg100.0%
distribute-frac-neg100.0%
neg-mul-1100.0%
associate-/l*99.7%
associate-/r/99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 83.9%
if -4.9999999999999997e38 < y < 4.49999999999999979e23Initial program 100.0%
sub-neg100.0%
distribute-frac-neg100.0%
neg-mul-1100.0%
associate-/l*99.9%
associate-/r/100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 77.4%
Final simplification80.6%
(FPCore (x y) :precision binary64 (+ x (* y -0.005)))
double code(double x, double y) {
return x + (y * -0.005);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + (y * (-0.005d0))
end function
public static double code(double x, double y) {
return x + (y * -0.005);
}
def code(x, y): return x + (y * -0.005)
function code(x, y) return Float64(x + Float64(y * -0.005)) end
function tmp = code(x, y) tmp = x + (y * -0.005); end
code[x_, y_] := N[(x + N[(y * -0.005), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + y \cdot -0.005
\end{array}
Initial program 100.0%
sub-neg100.0%
distribute-frac-neg100.0%
neg-mul-1100.0%
associate-/l*99.8%
associate-/r/99.9%
metadata-eval99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x y) :precision binary64 x)
double code(double x, double y) {
return x;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x
end function
public static double code(double x, double y) {
return x;
}
def code(x, y): return x
function code(x, y) return x end
function tmp = code(x, y) tmp = x; end
code[x_, y_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 100.0%
sub-neg100.0%
distribute-frac-neg100.0%
neg-mul-1100.0%
associate-/l*99.8%
associate-/r/99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around inf 48.1%
Final simplification48.1%
herbie shell --seed 2023326
(FPCore (x y)
:name "Data.Colour.CIE:cieLAB from colour-2.3.3, D"
:precision binary64
(- x (/ y 200.0)))