
(FPCore (x y) :precision binary64 (/ (* x y) (* (* (+ x y) (+ x y)) (+ (+ x y) 1.0))))
double code(double x, double y) {
return (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0d0))
end function
public static double code(double x, double y) {
return (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0));
}
def code(x, y): return (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0))
function code(x, y) return Float64(Float64(x * y) / Float64(Float64(Float64(x + y) * Float64(x + y)) * Float64(Float64(x + y) + 1.0))) end
function tmp = code(x, y) tmp = (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0)); end
code[x_, y_] := N[(N[(x * y), $MachinePrecision] / N[(N[(N[(x + y), $MachinePrecision] * N[(x + y), $MachinePrecision]), $MachinePrecision] * N[(N[(x + y), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot y}{\left(\left(x + y\right) \cdot \left(x + y\right)\right) \cdot \left(\left(x + y\right) + 1\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 26 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (/ (* x y) (* (* (+ x y) (+ x y)) (+ (+ x y) 1.0))))
double code(double x, double y) {
return (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0d0))
end function
public static double code(double x, double y) {
return (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0));
}
def code(x, y): return (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0))
function code(x, y) return Float64(Float64(x * y) / Float64(Float64(Float64(x + y) * Float64(x + y)) * Float64(Float64(x + y) + 1.0))) end
function tmp = code(x, y) tmp = (x * y) / (((x + y) * (x + y)) * ((x + y) + 1.0)); end
code[x_, y_] := N[(N[(x * y), $MachinePrecision] / N[(N[(N[(x + y), $MachinePrecision] * N[(x + y), $MachinePrecision]), $MachinePrecision] * N[(N[(x + y), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot y}{\left(\left(x + y\right) \cdot \left(x + y\right)\right) \cdot \left(\left(x + y\right) + 1\right)}
\end{array}
(FPCore (x y) :precision binary64 (/ (* (/ y (+ (+ y x) 1.0)) (/ x (+ y x))) (+ y x)))
double code(double x, double y) {
return ((y / ((y + x) + 1.0)) * (x / (y + x))) / (y + x);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = ((y / ((y + x) + 1.0d0)) * (x / (y + x))) / (y + x)
end function
public static double code(double x, double y) {
return ((y / ((y + x) + 1.0)) * (x / (y + x))) / (y + x);
}
def code(x, y): return ((y / ((y + x) + 1.0)) * (x / (y + x))) / (y + x)
function code(x, y) return Float64(Float64(Float64(y / Float64(Float64(y + x) + 1.0)) * Float64(x / Float64(y + x))) / Float64(y + x)) end
function tmp = code(x, y) tmp = ((y / ((y + x) + 1.0)) * (x / (y + x))) / (y + x); end
code[x_, y_] := N[(N[(N[(y / N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] * N[(x / N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{y}{\left(y + x\right) + 1} \cdot \frac{x}{y + x}}{y + x}
\end{array}
Initial program 74.5%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (x y)
:precision binary64
(if (<= x -9.5e+157)
(* (/ (/ x (+ y x)) (+ y x)) (/ y x))
(if (<= x 1.8e-44)
(/ (* x (/ y (+ y x))) (* (+ y x) (+ (+ y x) 1.0)))
(/ (- -1.0) (* y (/ y x))))))
double code(double x, double y) {
double tmp;
if (x <= -9.5e+157) {
tmp = ((x / (y + x)) / (y + x)) * (y / x);
} else if (x <= 1.8e-44) {
tmp = (x * (y / (y + x))) / ((y + x) * ((y + x) + 1.0));
} else {
tmp = -(-1.0) / (y * (y / x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= (-9.5d+157)) then
tmp = ((x / (y + x)) / (y + x)) * (y / x)
else if (x <= 1.8d-44) then
tmp = (x * (y / (y + x))) / ((y + x) * ((y + x) + 1.0d0))
else
tmp = -(-1.0d0) / (y * (y / x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -9.5e+157) {
tmp = ((x / (y + x)) / (y + x)) * (y / x);
} else if (x <= 1.8e-44) {
tmp = (x * (y / (y + x))) / ((y + x) * ((y + x) + 1.0));
} else {
tmp = -(-1.0) / (y * (y / x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -9.5e+157: tmp = ((x / (y + x)) / (y + x)) * (y / x) elif x <= 1.8e-44: tmp = (x * (y / (y + x))) / ((y + x) * ((y + x) + 1.0)) else: tmp = -(-1.0) / (y * (y / x)) return tmp
function code(x, y) tmp = 0.0 if (x <= -9.5e+157) tmp = Float64(Float64(Float64(x / Float64(y + x)) / Float64(y + x)) * Float64(y / x)); elseif (x <= 1.8e-44) tmp = Float64(Float64(x * Float64(y / Float64(y + x))) / Float64(Float64(y + x) * Float64(Float64(y + x) + 1.0))); else tmp = Float64(Float64(-(-1.0)) / Float64(y * Float64(y / x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -9.5e+157) tmp = ((x / (y + x)) / (y + x)) * (y / x); elseif (x <= 1.8e-44) tmp = (x * (y / (y + x))) / ((y + x) * ((y + x) + 1.0)); else tmp = -(-1.0) / (y * (y / x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -9.5e+157], N[(N[(N[(x / N[(y + x), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision] * N[(y / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.8e-44], N[(N[(x * N[(y / N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y + x), $MachinePrecision] * N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[((--1.0) / N[(y * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -9.5 \cdot 10^{+157}:\\
\;\;\;\;\frac{\frac{x}{y + x}}{y + x} \cdot \frac{y}{x}\\
\mathbf{elif}\;x \leq 1.8 \cdot 10^{-44}:\\
\;\;\;\;\frac{x \cdot \frac{y}{y + x}}{\left(y + x\right) \cdot \left(\left(y + x\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{--1}{y \cdot \frac{y}{x}}\\
\end{array}
\end{array}
if x < -9.4999999999999996e157Initial program 70.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6491.3
Applied egg-rr91.3%
Taylor expanded in x around inf
Simplified91.3%
associate-*l/N/A
associate-*l*N/A
associate-/r*N/A
associate-*r/N/A
+-commutativeN/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64100.0
Applied egg-rr100.0%
if -9.4999999999999996e157 < x < 1.7999999999999999e-44Initial program 76.3%
*-commutativeN/A
associate-*l*N/A
times-fracN/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6497.4
Applied egg-rr97.4%
if 1.7999999999999999e-44 < x Initial program 72.5%
Taylor expanded in y around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6423.5
Simplified23.5%
*-lft-identityN/A
*-inversesN/A
times-fracN/A
associate-*r*N/A
times-fracN/A
clear-numN/A
un-div-invN/A
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6427.0
Applied egg-rr27.0%
associate-/l/N/A
frac-2negN/A
metadata-evalN/A
/-lowering-/.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
neg-lowering-neg.f6427.0
Applied egg-rr27.0%
Final simplification76.2%
(FPCore (x y)
:precision binary64
(if (<= x -1.8e+74)
(* (/ (/ x (+ y x)) (+ y x)) (/ y x))
(if (<= x -2.1e-298)
(* x (/ (/ y (+ y x)) (* (+ y x) (+ (+ y x) 1.0))))
(/ (/ x (+ y 1.0)) (+ y x)))))
double code(double x, double y) {
double tmp;
if (x <= -1.8e+74) {
tmp = ((x / (y + x)) / (y + x)) * (y / x);
} else if (x <= -2.1e-298) {
tmp = x * ((y / (y + x)) / ((y + x) * ((y + x) + 1.0)));
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= (-1.8d+74)) then
tmp = ((x / (y + x)) / (y + x)) * (y / x)
else if (x <= (-2.1d-298)) then
tmp = x * ((y / (y + x)) / ((y + x) * ((y + x) + 1.0d0)))
else
tmp = (x / (y + 1.0d0)) / (y + x)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -1.8e+74) {
tmp = ((x / (y + x)) / (y + x)) * (y / x);
} else if (x <= -2.1e-298) {
tmp = x * ((y / (y + x)) / ((y + x) * ((y + x) + 1.0)));
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -1.8e+74: tmp = ((x / (y + x)) / (y + x)) * (y / x) elif x <= -2.1e-298: tmp = x * ((y / (y + x)) / ((y + x) * ((y + x) + 1.0))) else: tmp = (x / (y + 1.0)) / (y + x) return tmp
function code(x, y) tmp = 0.0 if (x <= -1.8e+74) tmp = Float64(Float64(Float64(x / Float64(y + x)) / Float64(y + x)) * Float64(y / x)); elseif (x <= -2.1e-298) tmp = Float64(x * Float64(Float64(y / Float64(y + x)) / Float64(Float64(y + x) * Float64(Float64(y + x) + 1.0)))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / Float64(y + x)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -1.8e+74) tmp = ((x / (y + x)) / (y + x)) * (y / x); elseif (x <= -2.1e-298) tmp = x * ((y / (y + x)) / ((y + x) * ((y + x) + 1.0))); else tmp = (x / (y + 1.0)) / (y + x); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -1.8e+74], N[(N[(N[(x / N[(y + x), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision] * N[(y / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -2.1e-298], N[(x * N[(N[(y / N[(y + x), $MachinePrecision]), $MachinePrecision] / N[(N[(y + x), $MachinePrecision] * N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.8 \cdot 10^{+74}:\\
\;\;\;\;\frac{\frac{x}{y + x}}{y + x} \cdot \frac{y}{x}\\
\mathbf{elif}\;x \leq -2.1 \cdot 10^{-298}:\\
\;\;\;\;x \cdot \frac{\frac{y}{y + x}}{\left(y + x\right) \cdot \left(\left(y + x\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y + x}\\
\end{array}
\end{array}
if x < -1.79999999999999994e74Initial program 67.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6481.3
Applied egg-rr81.3%
Taylor expanded in x around inf
Simplified81.3%
associate-*l/N/A
associate-*l*N/A
associate-/r*N/A
associate-*r/N/A
+-commutativeN/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f6491.2
Applied egg-rr91.2%
if -1.79999999999999994e74 < x < -2.10000000000000005e-298Initial program 79.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6486.6
Applied egg-rr86.6%
associate-*l*N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6497.7
Applied egg-rr97.7%
if -2.10000000000000005e-298 < x Initial program 73.5%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6446.6
Simplified46.6%
Final simplification71.6%
(FPCore (x y)
:precision binary64
(if (<= x -9.5e+157)
(* (/ (/ x (+ y x)) (+ y x)) (/ y x))
(if (<= x -4e-296)
(* y (/ (/ x (* (+ y x) (+ (+ y x) 1.0))) (+ y x)))
(/ (/ x (+ y 1.0)) (+ y x)))))
double code(double x, double y) {
double tmp;
if (x <= -9.5e+157) {
tmp = ((x / (y + x)) / (y + x)) * (y / x);
} else if (x <= -4e-296) {
tmp = y * ((x / ((y + x) * ((y + x) + 1.0))) / (y + x));
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= (-9.5d+157)) then
tmp = ((x / (y + x)) / (y + x)) * (y / x)
else if (x <= (-4d-296)) then
tmp = y * ((x / ((y + x) * ((y + x) + 1.0d0))) / (y + x))
else
tmp = (x / (y + 1.0d0)) / (y + x)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -9.5e+157) {
tmp = ((x / (y + x)) / (y + x)) * (y / x);
} else if (x <= -4e-296) {
tmp = y * ((x / ((y + x) * ((y + x) + 1.0))) / (y + x));
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -9.5e+157: tmp = ((x / (y + x)) / (y + x)) * (y / x) elif x <= -4e-296: tmp = y * ((x / ((y + x) * ((y + x) + 1.0))) / (y + x)) else: tmp = (x / (y + 1.0)) / (y + x) return tmp
function code(x, y) tmp = 0.0 if (x <= -9.5e+157) tmp = Float64(Float64(Float64(x / Float64(y + x)) / Float64(y + x)) * Float64(y / x)); elseif (x <= -4e-296) tmp = Float64(y * Float64(Float64(x / Float64(Float64(y + x) * Float64(Float64(y + x) + 1.0))) / Float64(y + x))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / Float64(y + x)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -9.5e+157) tmp = ((x / (y + x)) / (y + x)) * (y / x); elseif (x <= -4e-296) tmp = y * ((x / ((y + x) * ((y + x) + 1.0))) / (y + x)); else tmp = (x / (y + 1.0)) / (y + x); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -9.5e+157], N[(N[(N[(x / N[(y + x), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision] * N[(y / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -4e-296], N[(y * N[(N[(x / N[(N[(y + x), $MachinePrecision] * N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -9.5 \cdot 10^{+157}:\\
\;\;\;\;\frac{\frac{x}{y + x}}{y + x} \cdot \frac{y}{x}\\
\mathbf{elif}\;x \leq -4 \cdot 10^{-296}:\\
\;\;\;\;y \cdot \frac{\frac{x}{\left(y + x\right) \cdot \left(\left(y + x\right) + 1\right)}}{y + x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y + x}\\
\end{array}
\end{array}
if x < -9.4999999999999996e157Initial program 70.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6491.3
Applied egg-rr91.3%
Taylor expanded in x around inf
Simplified91.3%
associate-*l/N/A
associate-*l*N/A
associate-/r*N/A
associate-*r/N/A
+-commutativeN/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64100.0
Applied egg-rr100.0%
if -9.4999999999999996e157 < x < -4e-296Initial program 76.7%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
*-commutativeN/A
clear-numN/A
frac-timesN/A
div-invN/A
*-lft-identityN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6499.4
Applied egg-rr99.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6492.4
Applied egg-rr92.4%
if -4e-296 < x Initial program 73.5%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6446.6
Simplified46.6%
Final simplification70.8%
(FPCore (x y)
:precision binary64
(let* ((t_0 (+ (+ y x) 1.0)))
(if (<= x -9.5e+157)
(/ (/ y t_0) (+ y x))
(if (<= x -1.1e-293)
(* y (/ (/ x (* (+ y x) t_0)) (+ y x)))
(/ (/ x (+ y 1.0)) (+ y x))))))
double code(double x, double y) {
double t_0 = (y + x) + 1.0;
double tmp;
if (x <= -9.5e+157) {
tmp = (y / t_0) / (y + x);
} else if (x <= -1.1e-293) {
tmp = y * ((x / ((y + x) * t_0)) / (y + x));
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
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 = (y + x) + 1.0d0
if (x <= (-9.5d+157)) then
tmp = (y / t_0) / (y + x)
else if (x <= (-1.1d-293)) then
tmp = y * ((x / ((y + x) * t_0)) / (y + x))
else
tmp = (x / (y + 1.0d0)) / (y + x)
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = (y + x) + 1.0;
double tmp;
if (x <= -9.5e+157) {
tmp = (y / t_0) / (y + x);
} else if (x <= -1.1e-293) {
tmp = y * ((x / ((y + x) * t_0)) / (y + x));
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
def code(x, y): t_0 = (y + x) + 1.0 tmp = 0 if x <= -9.5e+157: tmp = (y / t_0) / (y + x) elif x <= -1.1e-293: tmp = y * ((x / ((y + x) * t_0)) / (y + x)) else: tmp = (x / (y + 1.0)) / (y + x) return tmp
function code(x, y) t_0 = Float64(Float64(y + x) + 1.0) tmp = 0.0 if (x <= -9.5e+157) tmp = Float64(Float64(y / t_0) / Float64(y + x)); elseif (x <= -1.1e-293) tmp = Float64(y * Float64(Float64(x / Float64(Float64(y + x) * t_0)) / Float64(y + x))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / Float64(y + x)); end return tmp end
function tmp_2 = code(x, y) t_0 = (y + x) + 1.0; tmp = 0.0; if (x <= -9.5e+157) tmp = (y / t_0) / (y + x); elseif (x <= -1.1e-293) tmp = y * ((x / ((y + x) * t_0)) / (y + x)); else tmp = (x / (y + 1.0)) / (y + x); end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]}, If[LessEqual[x, -9.5e+157], N[(N[(y / t$95$0), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -1.1e-293], N[(y * N[(N[(x / N[(N[(y + x), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(y + x\right) + 1\\
\mathbf{if}\;x \leq -9.5 \cdot 10^{+157}:\\
\;\;\;\;\frac{\frac{y}{t\_0}}{y + x}\\
\mathbf{elif}\;x \leq -1.1 \cdot 10^{-293}:\\
\;\;\;\;y \cdot \frac{\frac{x}{\left(y + x\right) \cdot t\_0}}{y + x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y + x}\\
\end{array}
\end{array}
if x < -9.4999999999999996e157Initial program 70.4%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64100.0
Applied egg-rr100.0%
Taylor expanded in x around inf
Simplified100.0%
if -9.4999999999999996e157 < x < -1.1e-293Initial program 76.5%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
*-commutativeN/A
clear-numN/A
frac-timesN/A
div-invN/A
*-lft-identityN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6499.4
Applied egg-rr99.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
/-lowering-/.f64N/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6492.3
Applied egg-rr92.3%
if -1.1e-293 < x Initial program 73.8%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6447.0
Simplified47.0%
Final simplification70.8%
(FPCore (x y)
:precision binary64
(let* ((t_0 (+ (+ y x) 1.0)))
(if (<= x -1.85e+74)
(/ (* y (/ 1.0 t_0)) (+ y x))
(if (<= x -3.4e-159)
(* x (/ y (* t_0 (* (+ y x) (+ y x)))))
(/ (/ x (+ y 1.0)) y)))))
double code(double x, double y) {
double t_0 = (y + x) + 1.0;
double tmp;
if (x <= -1.85e+74) {
tmp = (y * (1.0 / t_0)) / (y + x);
} else if (x <= -3.4e-159) {
tmp = x * (y / (t_0 * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 1.0)) / y;
}
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 = (y + x) + 1.0d0
if (x <= (-1.85d+74)) then
tmp = (y * (1.0d0 / t_0)) / (y + x)
else if (x <= (-3.4d-159)) then
tmp = x * (y / (t_0 * ((y + x) * (y + x))))
else
tmp = (x / (y + 1.0d0)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = (y + x) + 1.0;
double tmp;
if (x <= -1.85e+74) {
tmp = (y * (1.0 / t_0)) / (y + x);
} else if (x <= -3.4e-159) {
tmp = x * (y / (t_0 * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 1.0)) / y;
}
return tmp;
}
def code(x, y): t_0 = (y + x) + 1.0 tmp = 0 if x <= -1.85e+74: tmp = (y * (1.0 / t_0)) / (y + x) elif x <= -3.4e-159: tmp = x * (y / (t_0 * ((y + x) * (y + x)))) else: tmp = (x / (y + 1.0)) / y return tmp
function code(x, y) t_0 = Float64(Float64(y + x) + 1.0) tmp = 0.0 if (x <= -1.85e+74) tmp = Float64(Float64(y * Float64(1.0 / t_0)) / Float64(y + x)); elseif (x <= -3.4e-159) tmp = Float64(x * Float64(y / Float64(t_0 * Float64(Float64(y + x) * Float64(y + x))))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / y); end return tmp end
function tmp_2 = code(x, y) t_0 = (y + x) + 1.0; tmp = 0.0; if (x <= -1.85e+74) tmp = (y * (1.0 / t_0)) / (y + x); elseif (x <= -3.4e-159) tmp = x * (y / (t_0 * ((y + x) * (y + x)))); else tmp = (x / (y + 1.0)) / y; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]}, If[LessEqual[x, -1.85e+74], N[(N[(y * N[(1.0 / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.4e-159], N[(x * N[(y / N[(t$95$0 * N[(N[(y + x), $MachinePrecision] * N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(y + x\right) + 1\\
\mathbf{if}\;x \leq -1.85 \cdot 10^{+74}:\\
\;\;\;\;\frac{y \cdot \frac{1}{t\_0}}{y + x}\\
\mathbf{elif}\;x \leq -3.4 \cdot 10^{-159}:\\
\;\;\;\;x \cdot \frac{y}{t\_0 \cdot \left(\left(y + x\right) \cdot \left(y + x\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y}\\
\end{array}
\end{array}
if x < -1.8500000000000001e74Initial program 67.8%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.9
Applied egg-rr99.9%
associate-*l/N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6499.9
Applied egg-rr99.9%
Taylor expanded in x around inf
Simplified91.1%
if -1.8500000000000001e74 < x < -3.39999999999999984e-159Initial program 84.5%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6493.5
Applied egg-rr93.5%
if -3.39999999999999984e-159 < x Initial program 72.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.0
Applied egg-rr77.0%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6452.5
Simplified52.5%
*-commutativeN/A
un-div-invN/A
distribute-rgt1-inN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6454.3
Applied egg-rr54.3%
Final simplification69.9%
(FPCore (x y) :precision binary64 (* (/ x (+ y x)) (/ (/ y (+ (+ y x) 1.0)) (+ y x))))
double code(double x, double y) {
return (x / (y + x)) * ((y / ((y + x) + 1.0)) / (y + x));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x / (y + x)) * ((y / ((y + x) + 1.0d0)) / (y + x))
end function
public static double code(double x, double y) {
return (x / (y + x)) * ((y / ((y + x) + 1.0)) / (y + x));
}
def code(x, y): return (x / (y + x)) * ((y / ((y + x) + 1.0)) / (y + x))
function code(x, y) return Float64(Float64(x / Float64(y + x)) * Float64(Float64(y / Float64(Float64(y + x) + 1.0)) / Float64(y + x))) end
function tmp = code(x, y) tmp = (x / (y + x)) * ((y / ((y + x) + 1.0)) / (y + x)); end
code[x_, y_] := N[(N[(x / N[(y + x), $MachinePrecision]), $MachinePrecision] * N[(N[(y / N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y + x} \cdot \frac{\frac{y}{\left(y + x\right) + 1}}{y + x}
\end{array}
Initial program 74.5%
times-fracN/A
associate-*l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (x y)
:precision binary64
(if (<= x -8500000000000.0)
(/ (/ y (+ (+ y x) 1.0)) (+ y x))
(if (<= x -3.4e-159)
(* x (/ y (* (+ y 1.0) (* (+ y x) (+ y x)))))
(/ (/ x (+ y 1.0)) y))))
double code(double x, double y) {
double tmp;
if (x <= -8500000000000.0) {
tmp = (y / ((y + x) + 1.0)) / (y + x);
} else if (x <= -3.4e-159) {
tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 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 <= (-8500000000000.0d0)) then
tmp = (y / ((y + x) + 1.0d0)) / (y + x)
else if (x <= (-3.4d-159)) then
tmp = x * (y / ((y + 1.0d0) * ((y + x) * (y + x))))
else
tmp = (x / (y + 1.0d0)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -8500000000000.0) {
tmp = (y / ((y + x) + 1.0)) / (y + x);
} else if (x <= -3.4e-159) {
tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 1.0)) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -8500000000000.0: tmp = (y / ((y + x) + 1.0)) / (y + x) elif x <= -3.4e-159: tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x)))) else: tmp = (x / (y + 1.0)) / y return tmp
function code(x, y) tmp = 0.0 if (x <= -8500000000000.0) tmp = Float64(Float64(y / Float64(Float64(y + x) + 1.0)) / Float64(y + x)); elseif (x <= -3.4e-159) tmp = Float64(x * Float64(y / Float64(Float64(y + 1.0) * Float64(Float64(y + x) * Float64(y + x))))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -8500000000000.0) tmp = (y / ((y + x) + 1.0)) / (y + x); elseif (x <= -3.4e-159) tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x)))); else tmp = (x / (y + 1.0)) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -8500000000000.0], N[(N[(y / N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.4e-159], N[(x * N[(y / N[(N[(y + 1.0), $MachinePrecision] * N[(N[(y + x), $MachinePrecision] * N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -8500000000000:\\
\;\;\;\;\frac{\frac{y}{\left(y + x\right) + 1}}{y + x}\\
\mathbf{elif}\;x \leq -3.4 \cdot 10^{-159}:\\
\;\;\;\;x \cdot \frac{y}{\left(y + 1\right) \cdot \left(\left(y + x\right) \cdot \left(y + x\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y}\\
\end{array}
\end{array}
if x < -8.5e12Initial program 72.3%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
Taylor expanded in x around inf
Simplified78.4%
if -8.5e12 < x < -3.39999999999999984e-159Initial program 84.0%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6495.3
Applied egg-rr95.3%
Taylor expanded in x around 0
+-commutativeN/A
+-lowering-+.f6490.7
Simplified90.7%
if -3.39999999999999984e-159 < x Initial program 72.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.0
Applied egg-rr77.0%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6452.5
Simplified52.5%
*-commutativeN/A
un-div-invN/A
distribute-rgt1-inN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6454.3
Applied egg-rr54.3%
Final simplification66.4%
(FPCore (x y)
:precision binary64
(if (<= x -1.05e+16)
(/ (* y (/ 1.0 (+ (+ y x) 1.0))) (+ y x))
(if (<= x -3.4e-159)
(* x (/ y (* (+ y 1.0) (* (+ y x) (+ y x)))))
(/ (/ x (+ y 1.0)) y))))
double code(double x, double y) {
double tmp;
if (x <= -1.05e+16) {
tmp = (y * (1.0 / ((y + x) + 1.0))) / (y + x);
} else if (x <= -3.4e-159) {
tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 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 <= (-1.05d+16)) then
tmp = (y * (1.0d0 / ((y + x) + 1.0d0))) / (y + x)
else if (x <= (-3.4d-159)) then
tmp = x * (y / ((y + 1.0d0) * ((y + x) * (y + x))))
else
tmp = (x / (y + 1.0d0)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -1.05e+16) {
tmp = (y * (1.0 / ((y + x) + 1.0))) / (y + x);
} else if (x <= -3.4e-159) {
tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 1.0)) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -1.05e+16: tmp = (y * (1.0 / ((y + x) + 1.0))) / (y + x) elif x <= -3.4e-159: tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x)))) else: tmp = (x / (y + 1.0)) / y return tmp
function code(x, y) tmp = 0.0 if (x <= -1.05e+16) tmp = Float64(Float64(y * Float64(1.0 / Float64(Float64(y + x) + 1.0))) / Float64(y + x)); elseif (x <= -3.4e-159) tmp = Float64(x * Float64(y / Float64(Float64(y + 1.0) * Float64(Float64(y + x) * Float64(y + x))))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -1.05e+16) tmp = (y * (1.0 / ((y + x) + 1.0))) / (y + x); elseif (x <= -3.4e-159) tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x)))); else tmp = (x / (y + 1.0)) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -1.05e+16], N[(N[(y * N[(1.0 / N[(N[(y + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -3.4e-159], N[(x * N[(y / N[(N[(y + 1.0), $MachinePrecision] * N[(N[(y + x), $MachinePrecision] * N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.05 \cdot 10^{+16}:\\
\;\;\;\;\frac{y \cdot \frac{1}{\left(y + x\right) + 1}}{y + x}\\
\mathbf{elif}\;x \leq -3.4 \cdot 10^{-159}:\\
\;\;\;\;x \cdot \frac{y}{\left(y + 1\right) \cdot \left(\left(y + x\right) \cdot \left(y + x\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y}\\
\end{array}
\end{array}
if x < -1.05e16Initial program 72.3%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
associate-*l/N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in x around inf
Simplified78.4%
if -1.05e16 < x < -3.39999999999999984e-159Initial program 84.0%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6495.3
Applied egg-rr95.3%
Taylor expanded in x around 0
+-commutativeN/A
+-lowering-+.f6490.7
Simplified90.7%
if -3.39999999999999984e-159 < x Initial program 72.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.0
Applied egg-rr77.0%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6452.5
Simplified52.5%
*-commutativeN/A
un-div-invN/A
distribute-rgt1-inN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6454.3
Applied egg-rr54.3%
Final simplification66.4%
(FPCore (x y)
:precision binary64
(if (<= x -6500000000000.0)
(/ (/ y (+ y x)) x)
(if (<= x -3.4e-159)
(* x (/ y (* (+ y 1.0) (* (+ y x) (+ y x)))))
(/ (/ x (+ y 1.0)) y))))
double code(double x, double y) {
double tmp;
if (x <= -6500000000000.0) {
tmp = (y / (y + x)) / x;
} else if (x <= -3.4e-159) {
tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 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 <= (-6500000000000.0d0)) then
tmp = (y / (y + x)) / x
else if (x <= (-3.4d-159)) then
tmp = x * (y / ((y + 1.0d0) * ((y + x) * (y + x))))
else
tmp = (x / (y + 1.0d0)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -6500000000000.0) {
tmp = (y / (y + x)) / x;
} else if (x <= -3.4e-159) {
tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x))));
} else {
tmp = (x / (y + 1.0)) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -6500000000000.0: tmp = (y / (y + x)) / x elif x <= -3.4e-159: tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x)))) else: tmp = (x / (y + 1.0)) / y return tmp
function code(x, y) tmp = 0.0 if (x <= -6500000000000.0) tmp = Float64(Float64(y / Float64(y + x)) / x); elseif (x <= -3.4e-159) tmp = Float64(x * Float64(y / Float64(Float64(y + 1.0) * Float64(Float64(y + x) * Float64(y + x))))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -6500000000000.0) tmp = (y / (y + x)) / x; elseif (x <= -3.4e-159) tmp = x * (y / ((y + 1.0) * ((y + x) * (y + x)))); else tmp = (x / (y + 1.0)) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -6500000000000.0], N[(N[(y / N[(y + x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], If[LessEqual[x, -3.4e-159], N[(x * N[(y / N[(N[(y + 1.0), $MachinePrecision] * N[(N[(y + x), $MachinePrecision] * N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -6500000000000:\\
\;\;\;\;\frac{\frac{y}{y + x}}{x}\\
\mathbf{elif}\;x \leq -3.4 \cdot 10^{-159}:\\
\;\;\;\;x \cdot \frac{y}{\left(y + 1\right) \cdot \left(\left(y + x\right) \cdot \left(y + x\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y}\\
\end{array}
\end{array}
if x < -6.5e12Initial program 72.3%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
Taylor expanded in x around inf
/-lowering-/.f6478.3
Simplified78.3%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6478.3
Applied egg-rr78.3%
if -6.5e12 < x < -3.39999999999999984e-159Initial program 84.0%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6495.3
Applied egg-rr95.3%
Taylor expanded in x around 0
+-commutativeN/A
+-lowering-+.f6490.7
Simplified90.7%
if -3.39999999999999984e-159 < x Initial program 72.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.0
Applied egg-rr77.0%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6452.5
Simplified52.5%
*-commutativeN/A
un-div-invN/A
distribute-rgt1-inN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6454.3
Applied egg-rr54.3%
Final simplification66.3%
(FPCore (x y) :precision binary64 (if (<= y 4.4e-71) (/ y (* (+ y x) (+ x 1.0))) (if (<= y 1.72e+37) (/ x (fma y y y)) (/ (/ x y) (+ y x)))))
double code(double x, double y) {
double tmp;
if (y <= 4.4e-71) {
tmp = y / ((y + x) * (x + 1.0));
} else if (y <= 1.72e+37) {
tmp = x / fma(y, y, y);
} else {
tmp = (x / y) / (y + x);
}
return tmp;
}
function code(x, y) tmp = 0.0 if (y <= 4.4e-71) tmp = Float64(y / Float64(Float64(y + x) * Float64(x + 1.0))); elseif (y <= 1.72e+37) tmp = Float64(x / fma(y, y, y)); else tmp = Float64(Float64(x / y) / Float64(y + x)); end return tmp end
code[x_, y_] := If[LessEqual[y, 4.4e-71], N[(y / N[(N[(y + x), $MachinePrecision] * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.72e+37], N[(x / N[(y * y + y), $MachinePrecision]), $MachinePrecision], N[(N[(x / y), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 4.4 \cdot 10^{-71}:\\
\;\;\;\;\frac{y}{\left(y + x\right) \cdot \left(x + 1\right)}\\
\mathbf{elif}\;y \leq 1.72 \cdot 10^{+37}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, y, y\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y}}{y + x}\\
\end{array}
\end{array}
if y < 4.39999999999999995e-71Initial program 74.1%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6460.2
Simplified60.2%
associate-/l/N/A
/-lowering-/.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6462.0
Applied egg-rr62.0%
if 4.39999999999999995e-71 < y < 1.72000000000000002e37Initial program 92.0%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f6446.9
Simplified46.9%
if 1.72000000000000002e37 < y Initial program 67.5%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around inf
/-lowering-/.f6474.8
Simplified74.8%
Final simplification63.2%
(FPCore (x y) :precision binary64 (if (<= y -4.6e-144) (/ y (* x x)) (if (<= y 2.2e-210) (/ y (+ y x)) (if (<= y 1.0) (/ x y) (/ x (* y y))))))
double code(double x, double y) {
double tmp;
if (y <= -4.6e-144) {
tmp = y / (x * x);
} else if (y <= 2.2e-210) {
tmp = y / (y + x);
} else if (y <= 1.0) {
tmp = x / y;
} else {
tmp = x / (y * y);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= (-4.6d-144)) then
tmp = y / (x * x)
else if (y <= 2.2d-210) then
tmp = y / (y + x)
else if (y <= 1.0d0) then
tmp = x / y
else
tmp = x / (y * y)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= -4.6e-144) {
tmp = y / (x * x);
} else if (y <= 2.2e-210) {
tmp = y / (y + x);
} else if (y <= 1.0) {
tmp = x / y;
} else {
tmp = x / (y * y);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= -4.6e-144: tmp = y / (x * x) elif y <= 2.2e-210: tmp = y / (y + x) elif y <= 1.0: tmp = x / y else: tmp = x / (y * y) return tmp
function code(x, y) tmp = 0.0 if (y <= -4.6e-144) tmp = Float64(y / Float64(x * x)); elseif (y <= 2.2e-210) tmp = Float64(y / Float64(y + x)); elseif (y <= 1.0) tmp = Float64(x / y); else tmp = Float64(x / Float64(y * y)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= -4.6e-144) tmp = y / (x * x); elseif (y <= 2.2e-210) tmp = y / (y + x); elseif (y <= 1.0) tmp = x / y; else tmp = x / (y * y); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, -4.6e-144], N[(y / N[(x * x), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.2e-210], N[(y / N[(y + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.0], N[(x / y), $MachinePrecision], N[(x / N[(y * y), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.6 \cdot 10^{-144}:\\
\;\;\;\;\frac{y}{x \cdot x}\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{-210}:\\
\;\;\;\;\frac{y}{y + x}\\
\mathbf{elif}\;y \leq 1:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot y}\\
\end{array}
\end{array}
if y < -4.6e-144Initial program 78.9%
Taylor expanded in x around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6432.8
Simplified32.8%
if -4.6e-144 < y < 2.19999999999999989e-210Initial program 69.1%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6488.5
Simplified88.5%
Taylor expanded in x around 0
Simplified73.8%
if 2.19999999999999989e-210 < y < 1Initial program 77.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6476.6
Applied egg-rr76.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6437.2
Simplified37.2%
Taylor expanded in y around 0
/-lowering-/.f6435.1
Simplified35.1%
if 1 < y Initial program 70.3%
Taylor expanded in y around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6467.4
Simplified67.4%
Final simplification50.7%
(FPCore (x y) :precision binary64 (if (<= y 8.2e-72) (/ y (* (+ y x) (+ x 1.0))) (if (<= y 2e+37) (/ x (fma y y y)) (/ (/ x y) y))))
double code(double x, double y) {
double tmp;
if (y <= 8.2e-72) {
tmp = y / ((y + x) * (x + 1.0));
} else if (y <= 2e+37) {
tmp = x / fma(y, y, y);
} else {
tmp = (x / y) / y;
}
return tmp;
}
function code(x, y) tmp = 0.0 if (y <= 8.2e-72) tmp = Float64(y / Float64(Float64(y + x) * Float64(x + 1.0))); elseif (y <= 2e+37) tmp = Float64(x / fma(y, y, y)); else tmp = Float64(Float64(x / y) / y); end return tmp end
code[x_, y_] := If[LessEqual[y, 8.2e-72], N[(y / N[(N[(y + x), $MachinePrecision] * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2e+37], N[(x / N[(y * y + y), $MachinePrecision]), $MachinePrecision], N[(N[(x / y), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 8.2 \cdot 10^{-72}:\\
\;\;\;\;\frac{y}{\left(y + x\right) \cdot \left(x + 1\right)}\\
\mathbf{elif}\;y \leq 2 \cdot 10^{+37}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, y, y\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y}}{y}\\
\end{array}
\end{array}
if y < 8.20000000000000007e-72Initial program 74.1%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6460.2
Simplified60.2%
associate-/l/N/A
/-lowering-/.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6462.0
Applied egg-rr62.0%
if 8.20000000000000007e-72 < y < 1.99999999999999991e37Initial program 92.0%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f6446.9
Simplified46.9%
if 1.99999999999999991e37 < y Initial program 67.5%
Taylor expanded in y around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6466.1
Simplified66.1%
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6474.3
Applied egg-rr74.3%
Final simplification63.1%
(FPCore (x y) :precision binary64 (if (<= y 4.4e-76) (/ (/ y (+ x 1.0)) (+ y x)) (/ (/ x (+ y 1.0)) (+ y x))))
double code(double x, double y) {
double tmp;
if (y <= 4.4e-76) {
tmp = (y / (x + 1.0)) / (y + x);
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 4.4d-76) then
tmp = (y / (x + 1.0d0)) / (y + x)
else
tmp = (x / (y + 1.0d0)) / (y + x)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 4.4e-76) {
tmp = (y / (x + 1.0)) / (y + x);
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 4.4e-76: tmp = (y / (x + 1.0)) / (y + x) else: tmp = (x / (y + 1.0)) / (y + x) return tmp
function code(x, y) tmp = 0.0 if (y <= 4.4e-76) tmp = Float64(Float64(y / Float64(x + 1.0)) / Float64(y + x)); else tmp = Float64(Float64(x / Float64(y + 1.0)) / Float64(y + x)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 4.4e-76) tmp = (y / (x + 1.0)) / (y + x); else tmp = (x / (y + 1.0)) / (y + x); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 4.4e-76], N[(N[(y / N[(x + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 4.4 \cdot 10^{-76}:\\
\;\;\;\;\frac{\frac{y}{x + 1}}{y + x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y + x}\\
\end{array}
\end{array}
if y < 4.39999999999999999e-76Initial program 73.8%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6459.7
Simplified59.7%
if 4.39999999999999999e-76 < y Initial program 76.1%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6464.5
Simplified64.5%
Final simplification61.2%
(FPCore (x y) :precision binary64 (if (<= y 3.8e-69) (/ (/ y (+ x 1.0)) x) (/ (/ x (+ y 1.0)) (+ y x))))
double code(double x, double y) {
double tmp;
if (y <= 3.8e-69) {
tmp = (y / (x + 1.0)) / x;
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 3.8d-69) then
tmp = (y / (x + 1.0d0)) / x
else
tmp = (x / (y + 1.0d0)) / (y + x)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 3.8e-69) {
tmp = (y / (x + 1.0)) / x;
} else {
tmp = (x / (y + 1.0)) / (y + x);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 3.8e-69: tmp = (y / (x + 1.0)) / x else: tmp = (x / (y + 1.0)) / (y + x) return tmp
function code(x, y) tmp = 0.0 if (y <= 3.8e-69) tmp = Float64(Float64(y / Float64(x + 1.0)) / x); else tmp = Float64(Float64(x / Float64(y + 1.0)) / Float64(y + x)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 3.8e-69) tmp = (y / (x + 1.0)) / x; else tmp = (x / (y + 1.0)) / (y + x); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 3.8e-69], N[(N[(y / N[(x + 1.0), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 3.8 \cdot 10^{-69}:\\
\;\;\;\;\frac{\frac{y}{x + 1}}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y + x}\\
\end{array}
\end{array}
if y < 3.7999999999999998e-69Initial program 74.4%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6460.6
Simplified60.6%
Taylor expanded in x around inf
Simplified60.0%
if 3.7999999999999998e-69 < y Initial program 74.8%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6467.7
Simplified67.7%
Final simplification62.2%
(FPCore (x y) :precision binary64 (if (<= y 3.8e-69) (/ (/ y (+ x 1.0)) x) (/ (/ x (+ y 1.0)) y)))
double code(double x, double y) {
double tmp;
if (y <= 3.8e-69) {
tmp = (y / (x + 1.0)) / x;
} else {
tmp = (x / (y + 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 <= 3.8d-69) then
tmp = (y / (x + 1.0d0)) / x
else
tmp = (x / (y + 1.0d0)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 3.8e-69) {
tmp = (y / (x + 1.0)) / x;
} else {
tmp = (x / (y + 1.0)) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 3.8e-69: tmp = (y / (x + 1.0)) / x else: tmp = (x / (y + 1.0)) / y return tmp
function code(x, y) tmp = 0.0 if (y <= 3.8e-69) tmp = Float64(Float64(y / Float64(x + 1.0)) / x); else tmp = Float64(Float64(x / Float64(y + 1.0)) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 3.8e-69) tmp = (y / (x + 1.0)) / x; else tmp = (x / (y + 1.0)) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 3.8e-69], N[(N[(y / N[(x + 1.0), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 3.8 \cdot 10^{-69}:\\
\;\;\;\;\frac{\frac{y}{x + 1}}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y}\\
\end{array}
\end{array}
if y < 3.7999999999999998e-69Initial program 74.4%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6460.6
Simplified60.6%
Taylor expanded in x around inf
Simplified60.0%
if 3.7999999999999998e-69 < y Initial program 74.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6479.9
Applied egg-rr79.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6461.4
Simplified61.4%
*-commutativeN/A
un-div-invN/A
distribute-rgt1-inN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6467.1
Applied egg-rr67.1%
Final simplification62.1%
(FPCore (x y) :precision binary64 (if (<= y 3.9e-69) (/ y (* (+ y x) (+ x 1.0))) (/ (/ x (+ y 1.0)) y)))
double code(double x, double y) {
double tmp;
if (y <= 3.9e-69) {
tmp = y / ((y + x) * (x + 1.0));
} else {
tmp = (x / (y + 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 <= 3.9d-69) then
tmp = y / ((y + x) * (x + 1.0d0))
else
tmp = (x / (y + 1.0d0)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 3.9e-69) {
tmp = y / ((y + x) * (x + 1.0));
} else {
tmp = (x / (y + 1.0)) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 3.9e-69: tmp = y / ((y + x) * (x + 1.0)) else: tmp = (x / (y + 1.0)) / y return tmp
function code(x, y) tmp = 0.0 if (y <= 3.9e-69) tmp = Float64(y / Float64(Float64(y + x) * Float64(x + 1.0))); else tmp = Float64(Float64(x / Float64(y + 1.0)) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 3.9e-69) tmp = y / ((y + x) * (x + 1.0)); else tmp = (x / (y + 1.0)) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 3.9e-69], N[(y / N[(N[(y + x), $MachinePrecision] * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(y + 1.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 3.9 \cdot 10^{-69}:\\
\;\;\;\;\frac{y}{\left(y + x\right) \cdot \left(x + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{y + 1}}{y}\\
\end{array}
\end{array}
if y < 3.89999999999999981e-69Initial program 74.4%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6460.6
Simplified60.6%
associate-/l/N/A
/-lowering-/.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6462.5
Applied egg-rr62.5%
if 3.89999999999999981e-69 < y Initial program 74.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6479.9
Applied egg-rr79.9%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6461.4
Simplified61.4%
*-commutativeN/A
un-div-invN/A
distribute-rgt1-inN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6467.1
Applied egg-rr67.1%
Final simplification63.8%
(FPCore (x y) :precision binary64 (if (<= y 2.2e-210) (/ y (+ y x)) (if (<= y 1.0) (/ x y) (/ x (* y y)))))
double code(double x, double y) {
double tmp;
if (y <= 2.2e-210) {
tmp = y / (y + x);
} else if (y <= 1.0) {
tmp = x / y;
} else {
tmp = x / (y * y);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 2.2d-210) then
tmp = y / (y + x)
else if (y <= 1.0d0) then
tmp = x / y
else
tmp = x / (y * y)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 2.2e-210) {
tmp = y / (y + x);
} else if (y <= 1.0) {
tmp = x / y;
} else {
tmp = x / (y * y);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 2.2e-210: tmp = y / (y + x) elif y <= 1.0: tmp = x / y else: tmp = x / (y * y) return tmp
function code(x, y) tmp = 0.0 if (y <= 2.2e-210) tmp = Float64(y / Float64(y + x)); elseif (y <= 1.0) tmp = Float64(x / y); else tmp = Float64(x / Float64(y * y)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 2.2e-210) tmp = y / (y + x); elseif (y <= 1.0) tmp = x / y; else tmp = x / (y * y); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 2.2e-210], N[(y / N[(y + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.0], N[(x / y), $MachinePrecision], N[(x / N[(y * y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2.2 \cdot 10^{-210}:\\
\;\;\;\;\frac{y}{y + x}\\
\mathbf{elif}\;y \leq 1:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot y}\\
\end{array}
\end{array}
if y < 2.19999999999999989e-210Initial program 75.1%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6459.3
Simplified59.3%
Taylor expanded in x around 0
Simplified35.5%
if 2.19999999999999989e-210 < y < 1Initial program 77.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6476.6
Applied egg-rr76.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6437.2
Simplified37.2%
Taylor expanded in y around 0
/-lowering-/.f6435.1
Simplified35.1%
if 1 < y Initial program 70.3%
Taylor expanded in y around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6467.4
Simplified67.4%
Final simplification42.5%
(FPCore (x y) :precision binary64 (if (<= y 2.3e-69) (/ y (* (+ y x) (+ x 1.0))) (/ x (fma y y y))))
double code(double x, double y) {
double tmp;
if (y <= 2.3e-69) {
tmp = y / ((y + x) * (x + 1.0));
} else {
tmp = x / fma(y, y, y);
}
return tmp;
}
function code(x, y) tmp = 0.0 if (y <= 2.3e-69) tmp = Float64(y / Float64(Float64(y + x) * Float64(x + 1.0))); else tmp = Float64(x / fma(y, y, y)); end return tmp end
code[x_, y_] := If[LessEqual[y, 2.3e-69], N[(y / N[(N[(y + x), $MachinePrecision] * N[(x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(y * y + y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2.3 \cdot 10^{-69}:\\
\;\;\;\;\frac{y}{\left(y + x\right) \cdot \left(x + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, y, y\right)}\\
\end{array}
\end{array}
if y < 2.3000000000000001e-69Initial program 74.4%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6460.6
Simplified60.6%
associate-/l/N/A
/-lowering-/.f64N/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f6462.5
Applied egg-rr62.5%
if 2.3000000000000001e-69 < y Initial program 74.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f6461.4
Simplified61.4%
Final simplification62.2%
(FPCore (x y) :precision binary64 (if (<= y 3.7e-69) (/ y (fma x x x)) (/ x (fma y y y))))
double code(double x, double y) {
double tmp;
if (y <= 3.7e-69) {
tmp = y / fma(x, x, x);
} else {
tmp = x / fma(y, y, y);
}
return tmp;
}
function code(x, y) tmp = 0.0 if (y <= 3.7e-69) tmp = Float64(y / fma(x, x, x)); else tmp = Float64(x / fma(y, y, y)); end return tmp end
code[x_, y_] := If[LessEqual[y, 3.7e-69], N[(y / N[(x * x + x), $MachinePrecision]), $MachinePrecision], N[(x / N[(y * y + y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 3.7 \cdot 10^{-69}:\\
\;\;\;\;\frac{y}{\mathsf{fma}\left(x, x, x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, y, y\right)}\\
\end{array}
\end{array}
if y < 3.7000000000000002e-69Initial program 74.4%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f6460.3
Simplified60.3%
if 3.7000000000000002e-69 < y Initial program 74.8%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f6461.4
Simplified61.4%
(FPCore (x y) :precision binary64 (if (<= x -8.2e+14) (/ y (* x x)) (/ x (fma y y y))))
double code(double x, double y) {
double tmp;
if (x <= -8.2e+14) {
tmp = y / (x * x);
} else {
tmp = x / fma(y, y, y);
}
return tmp;
}
function code(x, y) tmp = 0.0 if (x <= -8.2e+14) tmp = Float64(y / Float64(x * x)); else tmp = Float64(x / fma(y, y, y)); end return tmp end
code[x_, y_] := If[LessEqual[x, -8.2e+14], N[(y / N[(x * x), $MachinePrecision]), $MachinePrecision], N[(x / N[(y * y + y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -8.2 \cdot 10^{+14}:\\
\;\;\;\;\frac{y}{x \cdot x}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, y, y\right)}\\
\end{array}
\end{array}
if x < -8.2e14Initial program 72.3%
Taylor expanded in x around inf
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6473.7
Simplified73.7%
if -8.2e14 < x Initial program 75.1%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
accelerator-lowering-fma.f6454.3
Simplified54.3%
(FPCore (x y) :precision binary64 (if (<= x -6.1e-134) (/ y (+ y x)) (/ x y)))
double code(double x, double y) {
double tmp;
if (x <= -6.1e-134) {
tmp = y / (y + x);
} 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 <= (-6.1d-134)) then
tmp = y / (y + x)
else
tmp = x / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -6.1e-134) {
tmp = y / (y + x);
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -6.1e-134: tmp = y / (y + x) else: tmp = x / y return tmp
function code(x, y) tmp = 0.0 if (x <= -6.1e-134) tmp = Float64(y / Float64(y + x)); else tmp = Float64(x / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -6.1e-134) tmp = y / (y + x); else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -6.1e-134], N[(y / N[(y + x), $MachinePrecision]), $MachinePrecision], N[(x / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -6.1 \cdot 10^{-134}:\\
\;\;\;\;\frac{y}{y + x}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if x < -6.0999999999999996e-134Initial program 75.0%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.7
Simplified68.7%
Taylor expanded in x around 0
Simplified34.5%
if -6.0999999999999996e-134 < x Initial program 74.2%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6478.6
Applied egg-rr78.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6455.6
Simplified55.6%
Taylor expanded in y around 0
/-lowering-/.f6435.9
Simplified35.9%
Final simplification35.4%
(FPCore (x y) :precision binary64 (if (<= x -5.3e+14) (/ 1.0 x) (/ x y)))
double code(double x, double y) {
double tmp;
if (x <= -5.3e+14) {
tmp = 1.0 / x;
} 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 <= (-5.3d+14)) then
tmp = 1.0d0 / x
else
tmp = x / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= -5.3e+14) {
tmp = 1.0 / x;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y): tmp = 0 if x <= -5.3e+14: tmp = 1.0 / x else: tmp = x / y return tmp
function code(x, y) tmp = 0.0 if (x <= -5.3e+14) tmp = Float64(1.0 / x); else tmp = Float64(x / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= -5.3e+14) tmp = 1.0 / x; else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, -5.3e+14], N[(1.0 / x), $MachinePrecision], N[(x / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.3 \cdot 10^{+14}:\\
\;\;\;\;\frac{1}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if x < -5.3e14Initial program 72.3%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.7
Applied egg-rr99.7%
Taylor expanded in x around inf
/-lowering-/.f6478.3
Simplified78.3%
Taylor expanded in y around inf
/-lowering-/.f645.7
Simplified5.7%
if -5.3e14 < x Initial program 75.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6481.3
Applied egg-rr81.3%
Taylor expanded in x around 0
/-lowering-/.f64N/A
distribute-lft-inN/A
*-rgt-identityN/A
unpow2N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6454.2
Simplified54.2%
Taylor expanded in y around 0
/-lowering-/.f6431.9
Simplified31.9%
(FPCore (x y) :precision binary64 (/ 1.0 x))
double code(double x, double y) {
return 1.0 / x;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 / x
end function
public static double code(double x, double y) {
return 1.0 / x;
}
def code(x, y): return 1.0 / x
function code(x, y) return Float64(1.0 / x) end
function tmp = code(x, y) tmp = 1.0 / x; end
code[x_, y_] := N[(1.0 / x), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{x}
\end{array}
Initial program 74.5%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in x around inf
/-lowering-/.f6439.3
Simplified39.3%
Taylor expanded in y around inf
/-lowering-/.f644.2
Simplified4.2%
(FPCore (x y) :precision binary64 (/ 0.5 y))
double code(double x, double y) {
return 0.5 / y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 0.5d0 / y
end function
public static double code(double x, double y) {
return 0.5 / y;
}
def code(x, y): return 0.5 / y
function code(x, y) return Float64(0.5 / y) end
function tmp = code(x, y) tmp = 0.5 / y; end
code[x_, y_] := N[(0.5 / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5}{y}
\end{array}
Initial program 74.5%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6481.7
Applied egg-rr81.7%
Taylor expanded in y around 0
+-commutativeN/A
unpow2N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6456.4
Simplified56.4%
Taylor expanded in y around inf
/-lowering-/.f644.1
Simplified4.1%
(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 74.5%
times-fracN/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6499.8
Applied egg-rr99.8%
Taylor expanded in y around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6453.0
Simplified53.0%
Taylor expanded in x around 0
Simplified3.7%
(FPCore (x y) :precision binary64 (/ (/ (/ x (+ (+ y 1.0) x)) (+ y x)) (/ 1.0 (/ y (+ y x)))))
double code(double x, double y) {
return ((x / ((y + 1.0) + x)) / (y + x)) / (1.0 / (y / (y + x)));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = ((x / ((y + 1.0d0) + x)) / (y + x)) / (1.0d0 / (y / (y + x)))
end function
public static double code(double x, double y) {
return ((x / ((y + 1.0) + x)) / (y + x)) / (1.0 / (y / (y + x)));
}
def code(x, y): return ((x / ((y + 1.0) + x)) / (y + x)) / (1.0 / (y / (y + x)))
function code(x, y) return Float64(Float64(Float64(x / Float64(Float64(y + 1.0) + x)) / Float64(y + x)) / Float64(1.0 / Float64(y / Float64(y + x)))) end
function tmp = code(x, y) tmp = ((x / ((y + 1.0) + x)) / (y + x)) / (1.0 / (y / (y + x))); end
code[x_, y_] := N[(N[(N[(x / N[(N[(y + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision] / N[(y + x), $MachinePrecision]), $MachinePrecision] / N[(1.0 / N[(y / N[(y + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\frac{x}{\left(y + 1\right) + x}}{y + x}}{\frac{1}{\frac{y}{y + x}}}
\end{array}
herbie shell --seed 2024204
(FPCore (x y)
:name "Numeric.SpecFunctions:incompleteBetaApprox from math-functions-0.1.5.2, A"
:precision binary64
:alt
(! :herbie-platform default (/ (/ (/ x (+ (+ y 1) x)) (+ y x)) (/ 1 (/ y (+ y x)))))
(/ (* x y) (* (* (+ x y) (+ x y)) (+ (+ x y) 1.0))))