
(FPCore (x y) :precision binary64 (/ (- x y) (- 1.0 y)))
double code(double x, double y) {
return (x - y) / (1.0 - y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x - y) / (1.0d0 - y)
end function
public static double code(double x, double y) {
return (x - y) / (1.0 - y);
}
def code(x, y): return (x - y) / (1.0 - y)
function code(x, y) return Float64(Float64(x - y) / Float64(1.0 - y)) end
function tmp = code(x, y) tmp = (x - y) / (1.0 - y); end
code[x_, y_] := N[(N[(x - y), $MachinePrecision] / N[(1.0 - y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - y}{1 - y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (/ (- x y) (- 1.0 y)))
double code(double x, double y) {
return (x - y) / (1.0 - y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x - y) / (1.0d0 - y)
end function
public static double code(double x, double y) {
return (x - y) / (1.0 - y);
}
def code(x, y): return (x - y) / (1.0 - y)
function code(x, y) return Float64(Float64(x - y) / Float64(1.0 - y)) end
function tmp = code(x, y) tmp = (x - y) / (1.0 - y); end
code[x_, y_] := N[(N[(x - y), $MachinePrecision] / N[(1.0 - y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - y}{1 - y}
\end{array}
(FPCore (x y) :precision binary64 (/ (- x y) (- 1.0 y)))
double code(double x, double y) {
return (x - y) / (1.0 - y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x - y) / (1.0d0 - y)
end function
public static double code(double x, double y) {
return (x - y) / (1.0 - y);
}
def code(x, y): return (x - y) / (1.0 - y)
function code(x, y) return Float64(Float64(x - y) / Float64(1.0 - y)) end
function tmp = code(x, y) tmp = (x - y) / (1.0 - y); end
code[x_, y_] := N[(N[(x - y), $MachinePrecision] / N[(1.0 - y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - y}{1 - y}
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (x y)
:precision binary64
(let* ((t_0 (/ (- x) y)))
(if (<= y -2.4e+71)
1.0
(if (<= y -5.2e+52)
t_0
(if (<= y -1.0)
1.0
(if (<= y -1.02e-63)
(- y)
(if (<= y 1.0) x (if (<= y 1.65e+50) t_0 1.0))))))))
double code(double x, double y) {
double t_0 = -x / y;
double tmp;
if (y <= -2.4e+71) {
tmp = 1.0;
} else if (y <= -5.2e+52) {
tmp = t_0;
} else if (y <= -1.0) {
tmp = 1.0;
} else if (y <= -1.02e-63) {
tmp = -y;
} else if (y <= 1.0) {
tmp = x;
} else if (y <= 1.65e+50) {
tmp = t_0;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = -x / y
if (y <= (-2.4d+71)) then
tmp = 1.0d0
else if (y <= (-5.2d+52)) then
tmp = t_0
else if (y <= (-1.0d0)) then
tmp = 1.0d0
else if (y <= (-1.02d-63)) then
tmp = -y
else if (y <= 1.0d0) then
tmp = x
else if (y <= 1.65d+50) then
tmp = t_0
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = -x / y;
double tmp;
if (y <= -2.4e+71) {
tmp = 1.0;
} else if (y <= -5.2e+52) {
tmp = t_0;
} else if (y <= -1.0) {
tmp = 1.0;
} else if (y <= -1.02e-63) {
tmp = -y;
} else if (y <= 1.0) {
tmp = x;
} else if (y <= 1.65e+50) {
tmp = t_0;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): t_0 = -x / y tmp = 0 if y <= -2.4e+71: tmp = 1.0 elif y <= -5.2e+52: tmp = t_0 elif y <= -1.0: tmp = 1.0 elif y <= -1.02e-63: tmp = -y elif y <= 1.0: tmp = x elif y <= 1.65e+50: tmp = t_0 else: tmp = 1.0 return tmp
function code(x, y) t_0 = Float64(Float64(-x) / y) tmp = 0.0 if (y <= -2.4e+71) tmp = 1.0; elseif (y <= -5.2e+52) tmp = t_0; elseif (y <= -1.0) tmp = 1.0; elseif (y <= -1.02e-63) tmp = Float64(-y); elseif (y <= 1.0) tmp = x; elseif (y <= 1.65e+50) tmp = t_0; else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) t_0 = -x / y; tmp = 0.0; if (y <= -2.4e+71) tmp = 1.0; elseif (y <= -5.2e+52) tmp = t_0; elseif (y <= -1.0) tmp = 1.0; elseif (y <= -1.02e-63) tmp = -y; elseif (y <= 1.0) tmp = x; elseif (y <= 1.65e+50) tmp = t_0; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[((-x) / y), $MachinePrecision]}, If[LessEqual[y, -2.4e+71], 1.0, If[LessEqual[y, -5.2e+52], t$95$0, If[LessEqual[y, -1.0], 1.0, If[LessEqual[y, -1.02e-63], (-y), If[LessEqual[y, 1.0], x, If[LessEqual[y, 1.65e+50], t$95$0, 1.0]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-x}{y}\\
\mathbf{if}\;y \leq -2.4 \cdot 10^{+71}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -5.2 \cdot 10^{+52}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y \leq -1:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -1.02 \cdot 10^{-63}:\\
\;\;\;\;-y\\
\mathbf{elif}\;y \leq 1:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 1.65 \cdot 10^{+50}:\\
\;\;\;\;t_0\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -2.39999999999999981e71 or -5.2e52 < y < -1 or 1.65e50 < y Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 83.6%
if -2.39999999999999981e71 < y < -5.2e52 or 1 < y < 1.65e50Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 84.6%
sub-neg84.6%
metadata-eval84.6%
neg-mul-184.6%
distribute-neg-frac84.6%
+-commutative84.6%
Simplified84.6%
frac-2neg84.6%
+-commutative84.6%
metadata-eval84.6%
sub-neg84.6%
div-inv84.4%
remove-double-neg84.4%
sub-neg84.4%
metadata-eval84.4%
distribute-neg-in84.4%
neg-mul-184.4%
*-commutative84.4%
metadata-eval84.4%
fma-def84.4%
Applied egg-rr84.4%
associate-*r/84.6%
*-rgt-identity84.6%
fma-udef84.6%
*-commutative84.6%
neg-mul-184.6%
+-commutative84.6%
unsub-neg84.6%
Simplified84.6%
Taylor expanded in y around inf 79.3%
associate-*r/79.3%
neg-mul-179.3%
Simplified79.3%
if -1 < y < -1.01999999999999997e-63Initial program 99.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
times-frac99.8%
metadata-eval99.8%
*-lft-identity99.8%
sub-neg99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 62.4%
Taylor expanded in y around 0 54.3%
neg-mul-154.3%
Simplified54.3%
if -1.01999999999999997e-63 < y < 1Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around 0 82.1%
Final simplification80.6%
(FPCore (x y)
:precision binary64
(if (<= y -4e+71)
1.0
(if (<= y -2.85e+51)
(/ (- x) y)
(if (or (<= y -1e-63) (not (<= y 1.8e+50)))
(/ y (+ y -1.0))
(/ x (- 1.0 y))))))
double code(double x, double y) {
double tmp;
if (y <= -4e+71) {
tmp = 1.0;
} else if (y <= -2.85e+51) {
tmp = -x / y;
} else if ((y <= -1e-63) || !(y <= 1.8e+50)) {
tmp = y / (y + -1.0);
} else {
tmp = x / (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 <= (-4d+71)) then
tmp = 1.0d0
else if (y <= (-2.85d+51)) then
tmp = -x / y
else if ((y <= (-1d-63)) .or. (.not. (y <= 1.8d+50))) then
tmp = y / (y + (-1.0d0))
else
tmp = x / (1.0d0 - y)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -4e+71) {
tmp = 1.0;
} else if (y <= -2.85e+51) {
tmp = -x / y;
} else if ((y <= -1e-63) || !(y <= 1.8e+50)) {
tmp = y / (y + -1.0);
} else {
tmp = x / (1.0 - y);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -4e+71: tmp = 1.0 elif y <= -2.85e+51: tmp = -x / y elif (y <= -1e-63) or not (y <= 1.8e+50): tmp = y / (y + -1.0) else: tmp = x / (1.0 - y) return tmp
function code(x, y) tmp = 0.0 if (y <= -4e+71) tmp = 1.0; elseif (y <= -2.85e+51) tmp = Float64(Float64(-x) / y); elseif ((y <= -1e-63) || !(y <= 1.8e+50)) tmp = Float64(y / Float64(y + -1.0)); else tmp = Float64(x / Float64(1.0 - y)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -4e+71) tmp = 1.0; elseif (y <= -2.85e+51) tmp = -x / y; elseif ((y <= -1e-63) || ~((y <= 1.8e+50))) tmp = y / (y + -1.0); else tmp = x / (1.0 - y); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -4e+71], 1.0, If[LessEqual[y, -2.85e+51], N[((-x) / y), $MachinePrecision], If[Or[LessEqual[y, -1e-63], N[Not[LessEqual[y, 1.8e+50]], $MachinePrecision]], N[(y / N[(y + -1.0), $MachinePrecision]), $MachinePrecision], N[(x / N[(1.0 - y), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4 \cdot 10^{+71}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -2.85 \cdot 10^{+51}:\\
\;\;\;\;\frac{-x}{y}\\
\mathbf{elif}\;y \leq -1 \cdot 10^{-63} \lor \neg \left(y \leq 1.8 \cdot 10^{+50}\right):\\
\;\;\;\;\frac{y}{y + -1}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{1 - y}\\
\end{array}
\end{array}
if y < -4.0000000000000002e71Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 82.2%
if -4.0000000000000002e71 < y < -2.8500000000000001e51Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 90.4%
sub-neg90.4%
metadata-eval90.4%
neg-mul-190.4%
distribute-neg-frac90.4%
+-commutative90.4%
Simplified90.4%
frac-2neg90.4%
+-commutative90.4%
metadata-eval90.4%
sub-neg90.4%
div-inv90.4%
remove-double-neg90.4%
sub-neg90.4%
metadata-eval90.4%
distribute-neg-in90.4%
neg-mul-190.4%
*-commutative90.4%
metadata-eval90.4%
fma-def90.4%
Applied egg-rr90.4%
associate-*r/90.4%
*-rgt-identity90.4%
fma-udef90.4%
*-commutative90.4%
neg-mul-190.4%
+-commutative90.4%
unsub-neg90.4%
Simplified90.4%
Taylor expanded in y around inf 90.4%
associate-*r/90.4%
neg-mul-190.4%
Simplified90.4%
if -2.8500000000000001e51 < y < -1.00000000000000007e-63 or 1.79999999999999993e50 < y Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 80.2%
if -1.00000000000000007e-63 < y < 1.79999999999999993e50Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 83.1%
sub-neg83.1%
metadata-eval83.1%
neg-mul-183.1%
distribute-neg-frac83.1%
+-commutative83.1%
Simplified83.1%
frac-2neg83.1%
+-commutative83.1%
metadata-eval83.1%
sub-neg83.1%
div-inv83.0%
remove-double-neg83.0%
sub-neg83.0%
metadata-eval83.0%
distribute-neg-in83.0%
neg-mul-183.0%
*-commutative83.0%
metadata-eval83.0%
fma-def83.0%
Applied egg-rr83.0%
associate-*r/83.1%
*-rgt-identity83.1%
fma-udef83.1%
*-commutative83.1%
neg-mul-183.1%
+-commutative83.1%
unsub-neg83.1%
Simplified83.1%
Final simplification82.1%
(FPCore (x y)
:precision binary64
(if (<= y -2.7e+72)
1.0
(if (<= y -3e+53)
(/ (- x) y)
(if (<= y -0.0023) 1.0 (if (<= y 2.2e+50) (/ x (- 1.0 y)) 1.0)))))
double code(double x, double y) {
double tmp;
if (y <= -2.7e+72) {
tmp = 1.0;
} else if (y <= -3e+53) {
tmp = -x / y;
} else if (y <= -0.0023) {
tmp = 1.0;
} else if (y <= 2.2e+50) {
tmp = x / (1.0 - 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 <= (-2.7d+72)) then
tmp = 1.0d0
else if (y <= (-3d+53)) then
tmp = -x / y
else if (y <= (-0.0023d0)) then
tmp = 1.0d0
else if (y <= 2.2d+50) then
tmp = x / (1.0d0 - y)
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -2.7e+72) {
tmp = 1.0;
} else if (y <= -3e+53) {
tmp = -x / y;
} else if (y <= -0.0023) {
tmp = 1.0;
} else if (y <= 2.2e+50) {
tmp = x / (1.0 - y);
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -2.7e+72: tmp = 1.0 elif y <= -3e+53: tmp = -x / y elif y <= -0.0023: tmp = 1.0 elif y <= 2.2e+50: tmp = x / (1.0 - y) else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if (y <= -2.7e+72) tmp = 1.0; elseif (y <= -3e+53) tmp = Float64(Float64(-x) / y); elseif (y <= -0.0023) tmp = 1.0; elseif (y <= 2.2e+50) tmp = Float64(x / Float64(1.0 - y)); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -2.7e+72) tmp = 1.0; elseif (y <= -3e+53) tmp = -x / y; elseif (y <= -0.0023) tmp = 1.0; elseif (y <= 2.2e+50) tmp = x / (1.0 - y); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -2.7e+72], 1.0, If[LessEqual[y, -3e+53], N[((-x) / y), $MachinePrecision], If[LessEqual[y, -0.0023], 1.0, If[LessEqual[y, 2.2e+50], N[(x / N[(1.0 - y), $MachinePrecision]), $MachinePrecision], 1.0]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.7 \cdot 10^{+72}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -3 \cdot 10^{+53}:\\
\;\;\;\;\frac{-x}{y}\\
\mathbf{elif}\;y \leq -0.0023:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{+50}:\\
\;\;\;\;\frac{x}{1 - y}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -2.7000000000000001e72 or -2.99999999999999998e53 < y < -0.0023 or 2.20000000000000017e50 < y Initial program 99.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
sub-neg99.9%
+-commutative99.9%
neg-sub099.9%
associate-+l-99.9%
sub0-neg99.9%
neg-mul-199.9%
times-frac99.9%
metadata-eval99.9%
*-lft-identity99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 82.4%
if -2.7000000000000001e72 < y < -2.99999999999999998e53Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 90.4%
sub-neg90.4%
metadata-eval90.4%
neg-mul-190.4%
distribute-neg-frac90.4%
+-commutative90.4%
Simplified90.4%
frac-2neg90.4%
+-commutative90.4%
metadata-eval90.4%
sub-neg90.4%
div-inv90.4%
remove-double-neg90.4%
sub-neg90.4%
metadata-eval90.4%
distribute-neg-in90.4%
neg-mul-190.4%
*-commutative90.4%
metadata-eval90.4%
fma-def90.4%
Applied egg-rr90.4%
associate-*r/90.4%
*-rgt-identity90.4%
fma-udef90.4%
*-commutative90.4%
neg-mul-190.4%
+-commutative90.4%
unsub-neg90.4%
Simplified90.4%
Taylor expanded in y around inf 90.4%
associate-*r/90.4%
neg-mul-190.4%
Simplified90.4%
if -0.0023 < y < 2.20000000000000017e50Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 78.1%
sub-neg78.1%
metadata-eval78.1%
neg-mul-178.1%
distribute-neg-frac78.1%
+-commutative78.1%
Simplified78.1%
frac-2neg78.1%
+-commutative78.1%
metadata-eval78.1%
sub-neg78.1%
div-inv78.0%
remove-double-neg78.0%
sub-neg78.0%
metadata-eval78.0%
distribute-neg-in78.0%
neg-mul-178.0%
*-commutative78.0%
metadata-eval78.0%
fma-def78.0%
Applied egg-rr78.0%
associate-*r/78.1%
*-rgt-identity78.1%
fma-udef78.1%
*-commutative78.1%
neg-mul-178.1%
+-commutative78.1%
unsub-neg78.1%
Simplified78.1%
Final simplification80.3%
(FPCore (x y)
:precision binary64
(let* ((t_0 (+ 1.0 (/ (- 1.0 x) y))))
(if (<= y -11000.0)
t_0
(if (<= y -1.05e-63)
(* y (/ 1.0 (+ y -1.0)))
(if (<= y 10000000000000.0) (/ x (- 1.0 y)) t_0)))))
double code(double x, double y) {
double t_0 = 1.0 + ((1.0 - x) / y);
double tmp;
if (y <= -11000.0) {
tmp = t_0;
} else if (y <= -1.05e-63) {
tmp = y * (1.0 / (y + -1.0));
} else if (y <= 10000000000000.0) {
tmp = x / (1.0 - y);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 + ((1.0d0 - x) / y)
if (y <= (-11000.0d0)) then
tmp = t_0
else if (y <= (-1.05d-63)) then
tmp = y * (1.0d0 / (y + (-1.0d0)))
else if (y <= 10000000000000.0d0) then
tmp = x / (1.0d0 - y)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = 1.0 + ((1.0 - x) / y);
double tmp;
if (y <= -11000.0) {
tmp = t_0;
} else if (y <= -1.05e-63) {
tmp = y * (1.0 / (y + -1.0));
} else if (y <= 10000000000000.0) {
tmp = x / (1.0 - y);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y): t_0 = 1.0 + ((1.0 - x) / y) tmp = 0 if y <= -11000.0: tmp = t_0 elif y <= -1.05e-63: tmp = y * (1.0 / (y + -1.0)) elif y <= 10000000000000.0: tmp = x / (1.0 - y) else: tmp = t_0 return tmp
function code(x, y) t_0 = Float64(1.0 + Float64(Float64(1.0 - x) / y)) tmp = 0.0 if (y <= -11000.0) tmp = t_0; elseif (y <= -1.05e-63) tmp = Float64(y * Float64(1.0 / Float64(y + -1.0))); elseif (y <= 10000000000000.0) tmp = Float64(x / Float64(1.0 - y)); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y) t_0 = 1.0 + ((1.0 - x) / y); tmp = 0.0; if (y <= -11000.0) tmp = t_0; elseif (y <= -1.05e-63) tmp = y * (1.0 / (y + -1.0)); elseif (y <= 10000000000000.0) tmp = x / (1.0 - y); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(1.0 + N[(N[(1.0 - x), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -11000.0], t$95$0, If[LessEqual[y, -1.05e-63], N[(y * N[(1.0 / N[(y + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 10000000000000.0], N[(x / N[(1.0 - y), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 + \frac{1 - x}{y}\\
\mathbf{if}\;y \leq -11000:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y \leq -1.05 \cdot 10^{-63}:\\
\;\;\;\;y \cdot \frac{1}{y + -1}\\
\mathbf{elif}\;y \leq 10000000000000:\\
\;\;\;\;\frac{x}{1 - y}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if y < -11000 or 1e13 < y Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 99.6%
+-commutative99.6%
associate-+r+99.6%
mul-1-neg99.6%
unsub-neg99.6%
div-sub99.6%
unsub-neg99.6%
mul-1-neg99.6%
+-commutative99.6%
metadata-eval99.6%
distribute-lft-in99.6%
metadata-eval99.6%
sub-neg99.6%
associate-*r/99.6%
+-commutative99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
distribute-lft-in99.6%
metadata-eval99.6%
+-commutative99.6%
mul-1-neg99.6%
unsub-neg99.6%
Simplified99.6%
if -11000 < y < -1.05e-63Initial program 99.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
times-frac99.8%
metadata-eval99.8%
*-lft-identity99.8%
sub-neg99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 62.4%
clear-num62.3%
associate-/r/62.5%
sub-neg62.5%
metadata-eval62.5%
Applied egg-rr62.5%
if -1.05e-63 < y < 1e13Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 83.5%
sub-neg83.5%
metadata-eval83.5%
neg-mul-183.5%
distribute-neg-frac83.5%
+-commutative83.5%
Simplified83.5%
frac-2neg83.5%
+-commutative83.5%
metadata-eval83.5%
sub-neg83.5%
div-inv83.5%
remove-double-neg83.5%
sub-neg83.5%
metadata-eval83.5%
distribute-neg-in83.5%
neg-mul-183.5%
*-commutative83.5%
metadata-eval83.5%
fma-def83.5%
Applied egg-rr83.5%
associate-*r/83.5%
*-rgt-identity83.5%
fma-udef83.5%
*-commutative83.5%
neg-mul-183.5%
+-commutative83.5%
unsub-neg83.5%
Simplified83.5%
Final simplification90.2%
(FPCore (x y)
:precision binary64
(let* ((t_0 (- 1.0 (/ x y))))
(if (<= y -15600000000.0)
t_0
(if (<= y -6e-64)
(/ y (+ y -1.0))
(if (<= y 10000000000000.0) (/ x (- 1.0 y)) t_0)))))
double code(double x, double y) {
double t_0 = 1.0 - (x / y);
double tmp;
if (y <= -15600000000.0) {
tmp = t_0;
} else if (y <= -6e-64) {
tmp = y / (y + -1.0);
} else if (y <= 10000000000000.0) {
tmp = x / (1.0 - y);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 - (x / y)
if (y <= (-15600000000.0d0)) then
tmp = t_0
else if (y <= (-6d-64)) then
tmp = y / (y + (-1.0d0))
else if (y <= 10000000000000.0d0) then
tmp = x / (1.0d0 - y)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = 1.0 - (x / y);
double tmp;
if (y <= -15600000000.0) {
tmp = t_0;
} else if (y <= -6e-64) {
tmp = y / (y + -1.0);
} else if (y <= 10000000000000.0) {
tmp = x / (1.0 - y);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y): t_0 = 1.0 - (x / y) tmp = 0 if y <= -15600000000.0: tmp = t_0 elif y <= -6e-64: tmp = y / (y + -1.0) elif y <= 10000000000000.0: tmp = x / (1.0 - y) else: tmp = t_0 return tmp
function code(x, y) t_0 = Float64(1.0 - Float64(x / y)) tmp = 0.0 if (y <= -15600000000.0) tmp = t_0; elseif (y <= -6e-64) tmp = Float64(y / Float64(y + -1.0)); elseif (y <= 10000000000000.0) tmp = Float64(x / Float64(1.0 - y)); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y) t_0 = 1.0 - (x / y); tmp = 0.0; if (y <= -15600000000.0) tmp = t_0; elseif (y <= -6e-64) tmp = y / (y + -1.0); elseif (y <= 10000000000000.0) tmp = x / (1.0 - y); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(1.0 - N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -15600000000.0], t$95$0, If[LessEqual[y, -6e-64], N[(y / N[(y + -1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 10000000000000.0], N[(x / N[(1.0 - y), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 - \frac{x}{y}\\
\mathbf{if}\;y \leq -15600000000:\\
\;\;\;\;t_0\\
\mathbf{elif}\;y \leq -6 \cdot 10^{-64}:\\
\;\;\;\;\frac{y}{y + -1}\\
\mathbf{elif}\;y \leq 10000000000000:\\
\;\;\;\;\frac{x}{1 - y}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if y < -1.56e10 or 1e13 < y Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 100.0%
+-commutative100.0%
associate-+r+100.0%
mul-1-neg100.0%
unsub-neg100.0%
div-sub100.0%
unsub-neg100.0%
mul-1-neg100.0%
+-commutative100.0%
metadata-eval100.0%
distribute-lft-in100.0%
metadata-eval100.0%
sub-neg100.0%
associate-*r/100.0%
+-commutative100.0%
associate-*r/100.0%
sub-neg100.0%
metadata-eval100.0%
distribute-lft-in100.0%
metadata-eval100.0%
+-commutative100.0%
mul-1-neg100.0%
unsub-neg100.0%
Simplified100.0%
Taylor expanded in x around inf 99.9%
neg-mul-199.9%
distribute-neg-frac99.9%
Simplified99.9%
if -1.56e10 < y < -6.0000000000000001e-64Initial program 99.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
times-frac99.8%
metadata-eval99.8%
*-lft-identity99.8%
sub-neg99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 62.8%
if -6.0000000000000001e-64 < y < 1e13Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 83.5%
sub-neg83.5%
metadata-eval83.5%
neg-mul-183.5%
distribute-neg-frac83.5%
+-commutative83.5%
Simplified83.5%
frac-2neg83.5%
+-commutative83.5%
metadata-eval83.5%
sub-neg83.5%
div-inv83.5%
remove-double-neg83.5%
sub-neg83.5%
metadata-eval83.5%
distribute-neg-in83.5%
neg-mul-183.5%
*-commutative83.5%
metadata-eval83.5%
fma-def83.5%
Applied egg-rr83.5%
associate-*r/83.5%
*-rgt-identity83.5%
fma-udef83.5%
*-commutative83.5%
neg-mul-183.5%
+-commutative83.5%
unsub-neg83.5%
Simplified83.5%
Final simplification90.1%
(FPCore (x y) :precision binary64 (if (<= y -1.0) 1.0 (if (<= y -5e-64) (- y) (if (<= y 1.0) x 1.0))))
double code(double x, double y) {
double tmp;
if (y <= -1.0) {
tmp = 1.0;
} else if (y <= -5e-64) {
tmp = -y;
} else if (y <= 1.0) {
tmp = x;
} 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)) then
tmp = 1.0d0
else if (y <= (-5d-64)) then
tmp = -y
else if (y <= 1.0d0) then
tmp = x
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -1.0) {
tmp = 1.0;
} else if (y <= -5e-64) {
tmp = -y;
} else if (y <= 1.0) {
tmp = x;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -1.0: tmp = 1.0 elif y <= -5e-64: tmp = -y elif y <= 1.0: tmp = x else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if (y <= -1.0) tmp = 1.0; elseif (y <= -5e-64) tmp = Float64(-y); elseif (y <= 1.0) tmp = x; else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -1.0) tmp = 1.0; elseif (y <= -5e-64) tmp = -y; elseif (y <= 1.0) tmp = x; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -1.0], 1.0, If[LessEqual[y, -5e-64], (-y), If[LessEqual[y, 1.0], x, 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -5 \cdot 10^{-64}:\\
\;\;\;\;-y\\
\mathbf{elif}\;y \leq 1:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -1 or 1 < y Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 74.1%
if -1 < y < -5.00000000000000033e-64Initial program 99.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
neg-mul-199.8%
times-frac99.8%
metadata-eval99.8%
*-lft-identity99.8%
sub-neg99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 62.4%
Taylor expanded in y around 0 54.3%
neg-mul-154.3%
Simplified54.3%
if -5.00000000000000033e-64 < y < 1Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around 0 82.1%
Final simplification76.0%
(FPCore (x y) :precision binary64 (if (<= y -0.0023) 1.0 (if (<= y 1.0) x 1.0)))
double code(double x, double y) {
double tmp;
if (y <= -0.0023) {
tmp = 1.0;
} else if (y <= 1.0) {
tmp = x;
} 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 <= (-0.0023d0)) then
tmp = 1.0d0
else if (y <= 1.0d0) then
tmp = x
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -0.0023) {
tmp = 1.0;
} else if (y <= 1.0) {
tmp = x;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -0.0023: tmp = 1.0 elif y <= 1.0: tmp = x else: tmp = 1.0 return tmp
function code(x, y) tmp = 0.0 if (y <= -0.0023) tmp = 1.0; elseif (y <= 1.0) tmp = x; else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -0.0023) tmp = 1.0; elseif (y <= 1.0) tmp = x; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -0.0023], 1.0, If[LessEqual[y, 1.0], x, 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -0.0023:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq 1:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -0.0023 or 1 < y Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 73.3%
if -0.0023 < y < 1Initial program 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around 0 75.6%
Final simplification74.4%
(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 100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
*-lft-identity100.0%
sub-neg100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 40.3%
Final simplification40.3%
herbie shell --seed 2023252
(FPCore (x y)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, C"
:precision binary64
(/ (- x y) (- 1.0 y)))