
(FPCore (x y) :precision binary64 (+ x (* (- 1.0 x) (- 1.0 y))))
double code(double x, double y) {
return x + ((1.0 - x) * (1.0 - y));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + ((1.0d0 - x) * (1.0d0 - y))
end function
public static double code(double x, double y) {
return x + ((1.0 - x) * (1.0 - y));
}
def code(x, y): return x + ((1.0 - x) * (1.0 - y))
function code(x, y) return Float64(x + Float64(Float64(1.0 - x) * Float64(1.0 - y))) end
function tmp = code(x, y) tmp = x + ((1.0 - x) * (1.0 - y)); end
code[x_, y_] := N[(x + N[(N[(1.0 - x), $MachinePrecision] * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(1 - x\right) \cdot \left(1 - y\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (+ x (* (- 1.0 x) (- 1.0 y))))
double code(double x, double y) {
return x + ((1.0 - x) * (1.0 - y));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + ((1.0d0 - x) * (1.0d0 - y))
end function
public static double code(double x, double y) {
return x + ((1.0 - x) * (1.0 - y));
}
def code(x, y): return x + ((1.0 - x) * (1.0 - y))
function code(x, y) return Float64(x + Float64(Float64(1.0 - x) * Float64(1.0 - y))) end
function tmp = code(x, y) tmp = x + ((1.0 - x) * (1.0 - y)); end
code[x_, y_] := N[(x + N[(N[(1.0 - x), $MachinePrecision] * N[(1.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(1 - x\right) \cdot \left(1 - y\right)
\end{array}
(FPCore (x y) :precision binary64 (- (+ 1.0 (* x y)) y))
double code(double x, double y) {
return (1.0 + (x * y)) - y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.0d0 + (x * y)) - y
end function
public static double code(double x, double y) {
return (1.0 + (x * y)) - y;
}
def code(x, y): return (1.0 + (x * y)) - y
function code(x, y) return Float64(Float64(1.0 + Float64(x * y)) - y) end
function tmp = code(x, y) tmp = (1.0 + (x * y)) - y; end
code[x_, y_] := N[(N[(1.0 + N[(x * y), $MachinePrecision]), $MachinePrecision] - y), $MachinePrecision]
\begin{array}{l}
\\
\left(1 + x \cdot y\right) - y
\end{array}
Initial program 77.4%
+-commutative77.4%
remove-double-neg77.4%
unsub-neg77.4%
sub-neg77.4%
+-commutative77.4%
distribute-rgt-in77.4%
*-lft-identity77.4%
associate--l+87.3%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0 100.0%
Final simplification100.0%
(FPCore (x y)
:precision binary64
(if (<= x -2.85e+118)
(* x y)
(if (<= x -1.05e-292)
1.0
(if (<= x 9e-267)
(- y)
(if (<= x 2.3e-228) 1.0 (if (<= x 1.0) (- y) (* x y)))))))
double code(double x, double y) {
double tmp;
if (x <= -2.85e+118) {
tmp = x * y;
} else if (x <= -1.05e-292) {
tmp = 1.0;
} else if (x <= 9e-267) {
tmp = -y;
} else if (x <= 2.3e-228) {
tmp = 1.0;
} else if (x <= 1.0) {
tmp = -y;
} else {
tmp = x * y;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= (-2.85d+118)) then
tmp = x * y
else if (x <= (-1.05d-292)) then
tmp = 1.0d0
else if (x <= 9d-267) then
tmp = -y
else if (x <= 2.3d-228) then
tmp = 1.0d0
else if (x <= 1.0d0) then
tmp = -y
else
tmp = x * y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -2.85e+118) {
tmp = x * y;
} else if (x <= -1.05e-292) {
tmp = 1.0;
} else if (x <= 9e-267) {
tmp = -y;
} else if (x <= 2.3e-228) {
tmp = 1.0;
} else if (x <= 1.0) {
tmp = -y;
} else {
tmp = x * y;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -2.85e+118: tmp = x * y elif x <= -1.05e-292: tmp = 1.0 elif x <= 9e-267: tmp = -y elif x <= 2.3e-228: tmp = 1.0 elif x <= 1.0: tmp = -y else: tmp = x * y return tmp
function code(x, y) tmp = 0.0 if (x <= -2.85e+118) tmp = Float64(x * y); elseif (x <= -1.05e-292) tmp = 1.0; elseif (x <= 9e-267) tmp = Float64(-y); elseif (x <= 2.3e-228) tmp = 1.0; elseif (x <= 1.0) tmp = Float64(-y); else tmp = Float64(x * y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -2.85e+118) tmp = x * y; elseif (x <= -1.05e-292) tmp = 1.0; elseif (x <= 9e-267) tmp = -y; elseif (x <= 2.3e-228) tmp = 1.0; elseif (x <= 1.0) tmp = -y; else tmp = x * y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -2.85e+118], N[(x * y), $MachinePrecision], If[LessEqual[x, -1.05e-292], 1.0, If[LessEqual[x, 9e-267], (-y), If[LessEqual[x, 2.3e-228], 1.0, If[LessEqual[x, 1.0], (-y), N[(x * y), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.85 \cdot 10^{+118}:\\
\;\;\;\;x \cdot y\\
\mathbf{elif}\;x \leq -1.05 \cdot 10^{-292}:\\
\;\;\;\;1\\
\mathbf{elif}\;x \leq 9 \cdot 10^{-267}:\\
\;\;\;\;-y\\
\mathbf{elif}\;x \leq 2.3 \cdot 10^{-228}:\\
\;\;\;\;1\\
\mathbf{elif}\;x \leq 1:\\
\;\;\;\;-y\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}
\end{array}
if x < -2.85000000000000001e118 or 1 < x Initial program 58.4%
+-commutative58.4%
remove-double-neg58.4%
unsub-neg58.4%
sub-neg58.4%
+-commutative58.4%
distribute-rgt-in58.4%
*-lft-identity58.4%
associate--l+79.8%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf 78.0%
*-commutative78.0%
Simplified78.0%
if -2.85000000000000001e118 < x < -1.04999999999999994e-292 or 8.9999999999999999e-267 < x < 2.2999999999999999e-228Initial program 88.6%
+-commutative88.6%
remove-double-neg88.6%
unsub-neg88.6%
sub-neg88.6%
+-commutative88.6%
distribute-rgt-in88.6%
*-lft-identity88.6%
associate--l+88.7%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in y around 0 58.4%
if -1.04999999999999994e-292 < x < 8.9999999999999999e-267 or 2.2999999999999999e-228 < x < 1Initial program 100.0%
+-commutative100.0%
remove-double-neg100.0%
unsub-neg100.0%
sub-neg100.0%
+-commutative100.0%
distribute-rgt-in100.0%
*-lft-identity100.0%
associate--l+100.0%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in y around inf 69.2%
Taylor expanded in x around 0 66.6%
neg-mul-166.6%
Simplified66.6%
Final simplification69.4%
(FPCore (x y) :precision binary64 (if (or (<= (- 1.0 y) -1000000000000.0) (not (<= (- 1.0 y) 2.0))) (* y (+ x -1.0)) (+ 1.0 (* x y))))
double code(double x, double y) {
double tmp;
if (((1.0 - y) <= -1000000000000.0) || !((1.0 - y) <= 2.0)) {
tmp = y * (x + -1.0);
} else {
tmp = 1.0 + (x * y);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (((1.0d0 - y) <= (-1000000000000.0d0)) .or. (.not. ((1.0d0 - y) <= 2.0d0))) then
tmp = y * (x + (-1.0d0))
else
tmp = 1.0d0 + (x * y)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (((1.0 - y) <= -1000000000000.0) || !((1.0 - y) <= 2.0)) {
tmp = y * (x + -1.0);
} else {
tmp = 1.0 + (x * y);
}
return tmp;
}
def code(x, y): tmp = 0 if ((1.0 - y) <= -1000000000000.0) or not ((1.0 - y) <= 2.0): tmp = y * (x + -1.0) else: tmp = 1.0 + (x * y) return tmp
function code(x, y) tmp = 0.0 if ((Float64(1.0 - y) <= -1000000000000.0) || !(Float64(1.0 - y) <= 2.0)) tmp = Float64(y * Float64(x + -1.0)); else tmp = Float64(1.0 + Float64(x * y)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (((1.0 - y) <= -1000000000000.0) || ~(((1.0 - y) <= 2.0))) tmp = y * (x + -1.0); else tmp = 1.0 + (x * y); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[N[(1.0 - y), $MachinePrecision], -1000000000000.0], N[Not[LessEqual[N[(1.0 - y), $MachinePrecision], 2.0]], $MachinePrecision]], N[(y * N[(x + -1.0), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(x * y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;1 - y \leq -1000000000000 \lor \neg \left(1 - y \leq 2\right):\\
\;\;\;\;y \cdot \left(x + -1\right)\\
\mathbf{else}:\\
\;\;\;\;1 + x \cdot y\\
\end{array}
\end{array}
if (-.f64 #s(literal 1 binary64) y) < -1e12 or 2 < (-.f64 #s(literal 1 binary64) y) Initial program 100.0%
+-commutative100.0%
remove-double-neg100.0%
unsub-neg100.0%
sub-neg100.0%
+-commutative100.0%
distribute-rgt-in100.0%
*-lft-identity100.0%
associate--l+100.0%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in y around inf 99.6%
if -1e12 < (-.f64 #s(literal 1 binary64) y) < 2Initial program 49.7%
+-commutative49.7%
remove-double-neg49.7%
unsub-neg49.7%
sub-neg49.7%
+-commutative49.7%
distribute-rgt-in49.7%
*-lft-identity49.7%
associate--l+71.8%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in x around inf 99.4%
*-commutative99.4%
Simplified99.4%
Final simplification99.5%
(FPCore (x y) :precision binary64 (if (or (<= (- 1.0 y) -1000000000000.0) (not (<= (- 1.0 y) 2.0))) (- (* x y) y) (+ 1.0 (* x y))))
double code(double x, double y) {
double tmp;
if (((1.0 - y) <= -1000000000000.0) || !((1.0 - y) <= 2.0)) {
tmp = (x * y) - y;
} else {
tmp = 1.0 + (x * y);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (((1.0d0 - y) <= (-1000000000000.0d0)) .or. (.not. ((1.0d0 - y) <= 2.0d0))) then
tmp = (x * y) - y
else
tmp = 1.0d0 + (x * y)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (((1.0 - y) <= -1000000000000.0) || !((1.0 - y) <= 2.0)) {
tmp = (x * y) - y;
} else {
tmp = 1.0 + (x * y);
}
return tmp;
}
def code(x, y): tmp = 0 if ((1.0 - y) <= -1000000000000.0) or not ((1.0 - y) <= 2.0): tmp = (x * y) - y else: tmp = 1.0 + (x * y) return tmp
function code(x, y) tmp = 0.0 if ((Float64(1.0 - y) <= -1000000000000.0) || !(Float64(1.0 - y) <= 2.0)) tmp = Float64(Float64(x * y) - y); else tmp = Float64(1.0 + Float64(x * y)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (((1.0 - y) <= -1000000000000.0) || ~(((1.0 - y) <= 2.0))) tmp = (x * y) - y; else tmp = 1.0 + (x * y); end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[N[(1.0 - y), $MachinePrecision], -1000000000000.0], N[Not[LessEqual[N[(1.0 - y), $MachinePrecision], 2.0]], $MachinePrecision]], N[(N[(x * y), $MachinePrecision] - y), $MachinePrecision], N[(1.0 + N[(x * y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;1 - y \leq -1000000000000 \lor \neg \left(1 - y \leq 2\right):\\
\;\;\;\;x \cdot y - y\\
\mathbf{else}:\\
\;\;\;\;1 + x \cdot y\\
\end{array}
\end{array}
if (-.f64 #s(literal 1 binary64) y) < -1e12 or 2 < (-.f64 #s(literal 1 binary64) y) Initial program 100.0%
+-commutative100.0%
remove-double-neg100.0%
unsub-neg100.0%
sub-neg100.0%
+-commutative100.0%
distribute-rgt-in100.0%
*-lft-identity100.0%
associate--l+100.0%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in y around inf 99.6%
sub-neg99.6%
metadata-eval99.6%
distribute-rgt-in99.7%
*-commutative99.7%
neg-mul-199.7%
Applied egg-rr99.7%
if -1e12 < (-.f64 #s(literal 1 binary64) y) < 2Initial program 49.7%
+-commutative49.7%
remove-double-neg49.7%
unsub-neg49.7%
sub-neg49.7%
+-commutative49.7%
distribute-rgt-in49.7%
*-lft-identity49.7%
associate--l+71.8%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in x around inf 99.4%
*-commutative99.4%
Simplified99.4%
Final simplification99.5%
(FPCore (x y) :precision binary64 (if (or (<= y -6.5e-69) (not (<= y 6.6e-15))) (* y (+ x -1.0)) (- 1.0 y)))
double code(double x, double y) {
double tmp;
if ((y <= -6.5e-69) || !(y <= 6.6e-15)) {
tmp = y * (x + -1.0);
} else {
tmp = 1.0 - y;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-6.5d-69)) .or. (.not. (y <= 6.6d-15))) then
tmp = y * (x + (-1.0d0))
else
tmp = 1.0d0 - y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -6.5e-69) || !(y <= 6.6e-15)) {
tmp = y * (x + -1.0);
} else {
tmp = 1.0 - y;
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -6.5e-69) or not (y <= 6.6e-15): tmp = y * (x + -1.0) else: tmp = 1.0 - y return tmp
function code(x, y) tmp = 0.0 if ((y <= -6.5e-69) || !(y <= 6.6e-15)) tmp = Float64(y * Float64(x + -1.0)); else tmp = Float64(1.0 - y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -6.5e-69) || ~((y <= 6.6e-15))) tmp = y * (x + -1.0); else tmp = 1.0 - y; end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -6.5e-69], N[Not[LessEqual[y, 6.6e-15]], $MachinePrecision]], N[(y * N[(x + -1.0), $MachinePrecision]), $MachinePrecision], N[(1.0 - y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.5 \cdot 10^{-69} \lor \neg \left(y \leq 6.6 \cdot 10^{-15}\right):\\
\;\;\;\;y \cdot \left(x + -1\right)\\
\mathbf{else}:\\
\;\;\;\;1 - y\\
\end{array}
\end{array}
if y < -6.49999999999999951e-69 or 6.6e-15 < y Initial program 94.6%
+-commutative94.6%
remove-double-neg94.6%
unsub-neg94.6%
sub-neg94.6%
+-commutative94.6%
distribute-rgt-in94.6%
*-lft-identity94.6%
associate--l+100.0%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in y around inf 96.1%
if -6.49999999999999951e-69 < y < 6.6e-15Initial program 50.5%
+-commutative50.5%
remove-double-neg50.5%
unsub-neg50.5%
sub-neg50.5%
+-commutative50.5%
distribute-rgt-in50.5%
*-lft-identity50.5%
associate--l+67.5%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in x around 0 79.9%
neg-mul-179.9%
unsub-neg79.9%
Simplified79.9%
Final simplification89.8%
(FPCore (x y) :precision binary64 (if (or (<= x -2.85e+118) (not (<= x 1.22e+64))) (* x y) (- 1.0 y)))
double code(double x, double y) {
double tmp;
if ((x <= -2.85e+118) || !(x <= 1.22e+64)) {
tmp = x * y;
} else {
tmp = 1.0 - y;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((x <= (-2.85d+118)) .or. (.not. (x <= 1.22d+64))) then
tmp = x * y
else
tmp = 1.0d0 - y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((x <= -2.85e+118) || !(x <= 1.22e+64)) {
tmp = x * y;
} else {
tmp = 1.0 - y;
}
return tmp;
}
def code(x, y): tmp = 0 if (x <= -2.85e+118) or not (x <= 1.22e+64): tmp = x * y else: tmp = 1.0 - y return tmp
function code(x, y) tmp = 0.0 if ((x <= -2.85e+118) || !(x <= 1.22e+64)) tmp = Float64(x * y); else tmp = Float64(1.0 - y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((x <= -2.85e+118) || ~((x <= 1.22e+64))) tmp = x * y; else tmp = 1.0 - y; end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[x, -2.85e+118], N[Not[LessEqual[x, 1.22e+64]], $MachinePrecision]], N[(x * y), $MachinePrecision], N[(1.0 - y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.85 \cdot 10^{+118} \lor \neg \left(x \leq 1.22 \cdot 10^{+64}\right):\\
\;\;\;\;x \cdot y\\
\mathbf{else}:\\
\;\;\;\;1 - y\\
\end{array}
\end{array}
if x < -2.85000000000000001e118 or 1.21999999999999994e64 < x Initial program 59.0%
+-commutative59.0%
remove-double-neg59.0%
unsub-neg59.0%
sub-neg59.0%
+-commutative59.0%
distribute-rgt-in59.0%
*-lft-identity59.0%
associate--l+84.1%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf 84.1%
*-commutative84.1%
Simplified84.1%
if -2.85000000000000001e118 < x < 1.21999999999999994e64Initial program 89.4%
+-commutative89.4%
remove-double-neg89.4%
unsub-neg89.4%
sub-neg89.4%
+-commutative89.4%
distribute-rgt-in89.4%
*-lft-identity89.4%
associate--l+89.4%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in x around 0 87.2%
neg-mul-187.2%
unsub-neg87.2%
Simplified87.2%
Final simplification86.0%
(FPCore (x y) :precision binary64 (if (or (<= y -1.0) (not (<= y 2.05e-10))) (- y) 1.0))
double code(double x, double y) {
double tmp;
if ((y <= -1.0) || !(y <= 2.05e-10)) {
tmp = -y;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if ((y <= (-1.0d0)) .or. (.not. (y <= 2.05d-10))) then
tmp = -y
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if ((y <= -1.0) || !(y <= 2.05e-10)) {
tmp = -y;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if (y <= -1.0) or not (y <= 2.05e-10): tmp = -y else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if ((y <= -1.0) || !(y <= 2.05e-10)) tmp = Float64(-y); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if ((y <= -1.0) || ~((y <= 2.05e-10))) tmp = -y; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[Or[LessEqual[y, -1.0], N[Not[LessEqual[y, 2.05e-10]], $MachinePrecision]], (-y), 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1 \lor \neg \left(y \leq 2.05 \cdot 10^{-10}\right):\\
\;\;\;\;-y\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -1 or 2.0499999999999999e-10 < y Initial program 99.6%
+-commutative99.6%
remove-double-neg99.6%
unsub-neg99.6%
sub-neg99.6%
+-commutative99.6%
distribute-rgt-in99.6%
*-lft-identity99.6%
associate--l+100.0%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
distribute-lft-in100.0%
*-commutative100.0%
mul-1-neg100.0%
Applied egg-rr100.0%
Taylor expanded in y around inf 99.6%
Taylor expanded in x around 0 45.4%
neg-mul-145.4%
Simplified45.4%
if -1 < y < 2.0499999999999999e-10Initial program 49.7%
+-commutative49.7%
remove-double-neg49.7%
unsub-neg49.7%
sub-neg49.7%
+-commutative49.7%
distribute-rgt-in49.7%
*-lft-identity49.7%
associate--l+71.5%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in y around 0 74.9%
Final simplification58.5%
(FPCore (x y) :precision binary64 (+ 1.0 (* y (+ x -1.0))))
double code(double x, double y) {
return 1.0 + (y * (x + -1.0));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 + (y * (x + (-1.0d0)))
end function
public static double code(double x, double y) {
return 1.0 + (y * (x + -1.0));
}
def code(x, y): return 1.0 + (y * (x + -1.0))
function code(x, y) return Float64(1.0 + Float64(y * Float64(x + -1.0))) end
function tmp = code(x, y) tmp = 1.0 + (y * (x + -1.0)); end
code[x_, y_] := N[(1.0 + N[(y * N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + y \cdot \left(x + -1\right)
\end{array}
Initial program 77.4%
+-commutative77.4%
remove-double-neg77.4%
unsub-neg77.4%
sub-neg77.4%
+-commutative77.4%
distribute-rgt-in77.4%
*-lft-identity77.4%
associate--l+87.3%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (x y) :precision binary64 1.0)
double code(double x, double y) {
return 1.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0
end function
public static double code(double x, double y) {
return 1.0;
}
def code(x, y): return 1.0
function code(x, y) return 1.0 end
function tmp = code(x, y) tmp = 1.0; end
code[x_, y_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 77.4%
+-commutative77.4%
remove-double-neg77.4%
unsub-neg77.4%
sub-neg77.4%
+-commutative77.4%
distribute-rgt-in77.4%
*-lft-identity77.4%
associate--l+87.3%
associate--l-100.0%
sub-neg100.0%
+-inverses100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
neg-sub0100.0%
associate--r-100.0%
metadata-eval100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in y around 0 34.8%
Final simplification34.8%
(FPCore (x y) :precision binary64 (- (* y x) (- y 1.0)))
double code(double x, double y) {
return (y * x) - (y - 1.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (y * x) - (y - 1.0d0)
end function
public static double code(double x, double y) {
return (y * x) - (y - 1.0);
}
def code(x, y): return (y * x) - (y - 1.0)
function code(x, y) return Float64(Float64(y * x) - Float64(y - 1.0)) end
function tmp = code(x, y) tmp = (y * x) - (y - 1.0); end
code[x_, y_] := N[(N[(y * x), $MachinePrecision] - N[(y - 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
y \cdot x - \left(y - 1\right)
\end{array}
herbie shell --seed 2024073
(FPCore (x y)
:name "Graphics.Rendering.Chart.Plot.Vectors:renderPlotVectors from Chart-1.5.3"
:precision binary64
:alt
(- (* y x) (- y 1.0))
(+ x (* (- 1.0 x) (- 1.0 y))))