
(FPCore (x y z t) :precision binary64 (+ (/ x y) (/ (+ 2.0 (* (* z 2.0) (- 1.0 t))) (* t z))))
double code(double x, double y, double z, double t) {
return (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x / y) + ((2.0d0 + ((z * 2.0d0) * (1.0d0 - t))) / (t * z))
end function
public static double code(double x, double y, double z, double t) {
return (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z));
}
def code(x, y, z, t): return (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z))
function code(x, y, z, t) return Float64(Float64(x / y) + Float64(Float64(2.0 + Float64(Float64(z * 2.0) * Float64(1.0 - t))) / Float64(t * z))) end
function tmp = code(x, y, z, t) tmp = (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z)); end
code[x_, y_, z_, t_] := N[(N[(x / y), $MachinePrecision] + N[(N[(2.0 + N[(N[(z * 2.0), $MachinePrecision] * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y} + \frac{2 + \left(z \cdot 2\right) \cdot \left(1 - t\right)}{t \cdot z}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (+ (/ x y) (/ (+ 2.0 (* (* z 2.0) (- 1.0 t))) (* t z))))
double code(double x, double y, double z, double t) {
return (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x / y) + ((2.0d0 + ((z * 2.0d0) * (1.0d0 - t))) / (t * z))
end function
public static double code(double x, double y, double z, double t) {
return (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z));
}
def code(x, y, z, t): return (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z))
function code(x, y, z, t) return Float64(Float64(x / y) + Float64(Float64(2.0 + Float64(Float64(z * 2.0) * Float64(1.0 - t))) / Float64(t * z))) end
function tmp = code(x, y, z, t) tmp = (x / y) + ((2.0 + ((z * 2.0) * (1.0 - t))) / (t * z)); end
code[x_, y_, z_, t_] := N[(N[(x / y), $MachinePrecision] + N[(N[(2.0 + N[(N[(z * 2.0), $MachinePrecision] * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y} + \frac{2 + \left(z \cdot 2\right) \cdot \left(1 - t\right)}{t \cdot z}
\end{array}
(FPCore (x y z t) :precision binary64 (let* ((t_1 (+ (/ (- 2.0 (* (- t 1.0) (* z 2.0))) (* t z)) (/ x y)))) (if (<= t_1 INFINITY) t_1 (+ (- -2.0 (/ -2.0 t)) (/ x y)))))
double code(double x, double y, double z, double t) {
double t_1 = ((2.0 - ((t - 1.0) * (z * 2.0))) / (t * z)) + (x / y);
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = (-2.0 - (-2.0 / t)) + (x / y);
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double t_1 = ((2.0 - ((t - 1.0) * (z * 2.0))) / (t * z)) + (x / y);
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = (-2.0 - (-2.0 / t)) + (x / y);
}
return tmp;
}
def code(x, y, z, t): t_1 = ((2.0 - ((t - 1.0) * (z * 2.0))) / (t * z)) + (x / y) tmp = 0 if t_1 <= math.inf: tmp = t_1 else: tmp = (-2.0 - (-2.0 / t)) + (x / y) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(Float64(2.0 - Float64(Float64(t - 1.0) * Float64(z * 2.0))) / Float64(t * z)) + Float64(x / y)) tmp = 0.0 if (t_1 <= Inf) tmp = t_1; else tmp = Float64(Float64(-2.0 - Float64(-2.0 / t)) + Float64(x / y)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = ((2.0 - ((t - 1.0) * (z * 2.0))) / (t * z)) + (x / y); tmp = 0.0; if (t_1 <= Inf) tmp = t_1; else tmp = (-2.0 - (-2.0 / t)) + (x / y); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(N[(2.0 - N[(N[(t - 1.0), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(N[(-2.0 - N[(-2.0 / t), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2 - \left(t - 1\right) \cdot \left(z \cdot 2\right)}{t \cdot z} + \frac{x}{y}\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(-2 - \frac{-2}{t}\right) + \frac{x}{y}\\
\end{array}
\end{array}
if (+.f64 (/.f64 x y) (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z))) < +inf.0Initial program 99.9%
if +inf.0 < (+.f64 (/.f64 x y) (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z))) Initial program 0.0%
Taylor expanded in z around inf
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
metadata-evalN/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6494.6
Applied rewrites94.6%
Final simplification99.1%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ -2.0 (/ x y)))
(t_2 (/ (- 2.0 (* (- t 1.0) (* z 2.0))) (* t z)))
(t_3 (- (/ 2.0 t) 2.0)))
(if (<= t_2 -1e+287)
(/ (/ 2.0 z) t)
(if (<= t_2 -2e+88)
t_3
(if (<= t_2 500000.0)
t_1
(if (<= t_2 1e+291)
t_3
(if (<= t_2 INFINITY) (/ 2.0 (* t z)) t_1)))))))
double code(double x, double y, double z, double t) {
double t_1 = -2.0 + (x / y);
double t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_3 = (2.0 / t) - 2.0;
double tmp;
if (t_2 <= -1e+287) {
tmp = (2.0 / z) / t;
} else if (t_2 <= -2e+88) {
tmp = t_3;
} else if (t_2 <= 500000.0) {
tmp = t_1;
} else if (t_2 <= 1e+291) {
tmp = t_3;
} else if (t_2 <= ((double) INFINITY)) {
tmp = 2.0 / (t * z);
} else {
tmp = t_1;
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double t_1 = -2.0 + (x / y);
double t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_3 = (2.0 / t) - 2.0;
double tmp;
if (t_2 <= -1e+287) {
tmp = (2.0 / z) / t;
} else if (t_2 <= -2e+88) {
tmp = t_3;
} else if (t_2 <= 500000.0) {
tmp = t_1;
} else if (t_2 <= 1e+291) {
tmp = t_3;
} else if (t_2 <= Double.POSITIVE_INFINITY) {
tmp = 2.0 / (t * z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = -2.0 + (x / y) t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z) t_3 = (2.0 / t) - 2.0 tmp = 0 if t_2 <= -1e+287: tmp = (2.0 / z) / t elif t_2 <= -2e+88: tmp = t_3 elif t_2 <= 500000.0: tmp = t_1 elif t_2 <= 1e+291: tmp = t_3 elif t_2 <= math.inf: tmp = 2.0 / (t * z) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(-2.0 + Float64(x / y)) t_2 = Float64(Float64(2.0 - Float64(Float64(t - 1.0) * Float64(z * 2.0))) / Float64(t * z)) t_3 = Float64(Float64(2.0 / t) - 2.0) tmp = 0.0 if (t_2 <= -1e+287) tmp = Float64(Float64(2.0 / z) / t); elseif (t_2 <= -2e+88) tmp = t_3; elseif (t_2 <= 500000.0) tmp = t_1; elseif (t_2 <= 1e+291) tmp = t_3; elseif (t_2 <= Inf) tmp = Float64(2.0 / Float64(t * z)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = -2.0 + (x / y); t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z); t_3 = (2.0 / t) - 2.0; tmp = 0.0; if (t_2 <= -1e+287) tmp = (2.0 / z) / t; elseif (t_2 <= -2e+88) tmp = t_3; elseif (t_2 <= 500000.0) tmp = t_1; elseif (t_2 <= 1e+291) tmp = t_3; elseif (t_2 <= Inf) tmp = 2.0 / (t * z); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(2.0 - N[(N[(t - 1.0), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision]}, If[LessEqual[t$95$2, -1e+287], N[(N[(2.0 / z), $MachinePrecision] / t), $MachinePrecision], If[LessEqual[t$95$2, -2e+88], t$95$3, If[LessEqual[t$95$2, 500000.0], t$95$1, If[LessEqual[t$95$2, 1e+291], t$95$3, If[LessEqual[t$95$2, Infinity], N[(2.0 / N[(t * z), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := -2 + \frac{x}{y}\\
t_2 := \frac{2 - \left(t - 1\right) \cdot \left(z \cdot 2\right)}{t \cdot z}\\
t_3 := \frac{2}{t} - 2\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{+287}:\\
\;\;\;\;\frac{\frac{2}{z}}{t}\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{+88}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq 500000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 10^{+291}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;\frac{2}{t \cdot z}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -1.0000000000000001e287Initial program 94.7%
Taylor expanded in t around 0
lower-/.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64100.0
Applied rewrites100.0%
Taylor expanded in z around 0
Applied rewrites94.9%
if -1.0000000000000001e287 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -1.99999999999999992e88 or 5e5 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 9.9999999999999996e290Initial program 99.7%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites79.0%
Taylor expanded in z around inf
Applied rewrites47.3%
if -1.99999999999999992e88 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 5e5 or +inf.0 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) Initial program 71.7%
Taylor expanded in t around inf
Applied rewrites94.6%
if 9.9999999999999996e290 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < +inf.0Initial program 89.5%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6494.7
Applied rewrites94.7%
Final simplification76.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ 2.0 (* t z)))
(t_2 (- (/ 2.0 t) 2.0))
(t_3 (/ (- 2.0 (* (- t 1.0) (* z 2.0))) (* t z)))
(t_4 (+ -2.0 (/ x y))))
(if (<= t_3 -1e+287)
t_1
(if (<= t_3 -2e+88)
t_2
(if (<= t_3 500000.0)
t_4
(if (<= t_3 1e+291) t_2 (if (<= t_3 INFINITY) t_1 t_4)))))))
double code(double x, double y, double z, double t) {
double t_1 = 2.0 / (t * z);
double t_2 = (2.0 / t) - 2.0;
double t_3 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_4 = -2.0 + (x / y);
double tmp;
if (t_3 <= -1e+287) {
tmp = t_1;
} else if (t_3 <= -2e+88) {
tmp = t_2;
} else if (t_3 <= 500000.0) {
tmp = t_4;
} else if (t_3 <= 1e+291) {
tmp = t_2;
} else if (t_3 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = t_4;
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double t_1 = 2.0 / (t * z);
double t_2 = (2.0 / t) - 2.0;
double t_3 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_4 = -2.0 + (x / y);
double tmp;
if (t_3 <= -1e+287) {
tmp = t_1;
} else if (t_3 <= -2e+88) {
tmp = t_2;
} else if (t_3 <= 500000.0) {
tmp = t_4;
} else if (t_3 <= 1e+291) {
tmp = t_2;
} else if (t_3 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = t_4;
}
return tmp;
}
def code(x, y, z, t): t_1 = 2.0 / (t * z) t_2 = (2.0 / t) - 2.0 t_3 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z) t_4 = -2.0 + (x / y) tmp = 0 if t_3 <= -1e+287: tmp = t_1 elif t_3 <= -2e+88: tmp = t_2 elif t_3 <= 500000.0: tmp = t_4 elif t_3 <= 1e+291: tmp = t_2 elif t_3 <= math.inf: tmp = t_1 else: tmp = t_4 return tmp
function code(x, y, z, t) t_1 = Float64(2.0 / Float64(t * z)) t_2 = Float64(Float64(2.0 / t) - 2.0) t_3 = Float64(Float64(2.0 - Float64(Float64(t - 1.0) * Float64(z * 2.0))) / Float64(t * z)) t_4 = Float64(-2.0 + Float64(x / y)) tmp = 0.0 if (t_3 <= -1e+287) tmp = t_1; elseif (t_3 <= -2e+88) tmp = t_2; elseif (t_3 <= 500000.0) tmp = t_4; elseif (t_3 <= 1e+291) tmp = t_2; elseif (t_3 <= Inf) tmp = t_1; else tmp = t_4; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 2.0 / (t * z); t_2 = (2.0 / t) - 2.0; t_3 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z); t_4 = -2.0 + (x / y); tmp = 0.0; if (t_3 <= -1e+287) tmp = t_1; elseif (t_3 <= -2e+88) tmp = t_2; elseif (t_3 <= 500000.0) tmp = t_4; elseif (t_3 <= 1e+291) tmp = t_2; elseif (t_3 <= Inf) tmp = t_1; else tmp = t_4; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(2.0 / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision]}, Block[{t$95$3 = N[(N[(2.0 - N[(N[(t - 1.0), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -1e+287], t$95$1, If[LessEqual[t$95$3, -2e+88], t$95$2, If[LessEqual[t$95$3, 500000.0], t$95$4, If[LessEqual[t$95$3, 1e+291], t$95$2, If[LessEqual[t$95$3, Infinity], t$95$1, t$95$4]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{t \cdot z}\\
t_2 := \frac{2}{t} - 2\\
t_3 := \frac{2 - \left(t - 1\right) \cdot \left(z \cdot 2\right)}{t \cdot z}\\
t_4 := -2 + \frac{x}{y}\\
\mathbf{if}\;t\_3 \leq -1 \cdot 10^{+287}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_3 \leq -2 \cdot 10^{+88}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_3 \leq 500000:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_3 \leq 10^{+291}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_3 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_4\\
\end{array}
\end{array}
if (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -1.0000000000000001e287 or 9.9999999999999996e290 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < +inf.0Initial program 92.1%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6494.8
Applied rewrites94.8%
if -1.0000000000000001e287 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -1.99999999999999992e88 or 5e5 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 9.9999999999999996e290Initial program 99.7%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites79.0%
Taylor expanded in z around inf
Applied rewrites47.3%
if -1.99999999999999992e88 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 5e5 or +inf.0 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) Initial program 71.7%
Taylor expanded in t around inf
Applied rewrites94.6%
Final simplification76.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (- (/ 2.0 z) -2.0) t))
(t_2 (/ (- 2.0 (* (- t 1.0) (* z 2.0))) (* t z)))
(t_3 (+ -2.0 (/ x y))))
(if (<= t_2 -4e+86)
t_1
(if (<= t_2 500000.0) t_3 (if (<= t_2 INFINITY) (- t_1 2.0) t_3)))))
double code(double x, double y, double z, double t) {
double t_1 = ((2.0 / z) - -2.0) / t;
double t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -4e+86) {
tmp = t_1;
} else if (t_2 <= 500000.0) {
tmp = t_3;
} else if (t_2 <= ((double) INFINITY)) {
tmp = t_1 - 2.0;
} else {
tmp = t_3;
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double t_1 = ((2.0 / z) - -2.0) / t;
double t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -4e+86) {
tmp = t_1;
} else if (t_2 <= 500000.0) {
tmp = t_3;
} else if (t_2 <= Double.POSITIVE_INFINITY) {
tmp = t_1 - 2.0;
} else {
tmp = t_3;
}
return tmp;
}
def code(x, y, z, t): t_1 = ((2.0 / z) - -2.0) / t t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z) t_3 = -2.0 + (x / y) tmp = 0 if t_2 <= -4e+86: tmp = t_1 elif t_2 <= 500000.0: tmp = t_3 elif t_2 <= math.inf: tmp = t_1 - 2.0 else: tmp = t_3 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(Float64(2.0 / z) - -2.0) / t) t_2 = Float64(Float64(2.0 - Float64(Float64(t - 1.0) * Float64(z * 2.0))) / Float64(t * z)) t_3 = Float64(-2.0 + Float64(x / y)) tmp = 0.0 if (t_2 <= -4e+86) tmp = t_1; elseif (t_2 <= 500000.0) tmp = t_3; elseif (t_2 <= Inf) tmp = Float64(t_1 - 2.0); else tmp = t_3; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = ((2.0 / z) - -2.0) / t; t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z); t_3 = -2.0 + (x / y); tmp = 0.0; if (t_2 <= -4e+86) tmp = t_1; elseif (t_2 <= 500000.0) tmp = t_3; elseif (t_2 <= Inf) tmp = t_1 - 2.0; else tmp = t_3; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision]}, Block[{t$95$2 = N[(N[(2.0 - N[(N[(t - 1.0), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -4e+86], t$95$1, If[LessEqual[t$95$2, 500000.0], t$95$3, If[LessEqual[t$95$2, Infinity], N[(t$95$1 - 2.0), $MachinePrecision], t$95$3]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\frac{2}{z} - -2}{t}\\
t_2 := \frac{2 - \left(t - 1\right) \cdot \left(z \cdot 2\right)}{t \cdot z}\\
t_3 := -2 + \frac{x}{y}\\
\mathbf{if}\;t\_2 \leq -4 \cdot 10^{+86}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 500000:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;t\_1 - 2\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -4.0000000000000001e86Initial program 98.1%
Taylor expanded in t around 0
lower-/.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6488.9
Applied rewrites88.9%
if -4.0000000000000001e86 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 5e5 or +inf.0 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) Initial program 71.4%
Taylor expanded in t around inf
Applied rewrites95.4%
if 5e5 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < +inf.0Initial program 97.2%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites80.6%
Final simplification89.4%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (- (/ 2.0 z) -2.0) t))
(t_2 (/ (- 2.0 (* (- t 1.0) (* z 2.0))) (* t z)))
(t_3 (+ -2.0 (/ x y))))
(if (<= t_2 -4e+86)
t_1
(if (<= t_2 500000.0) t_3 (if (<= t_2 INFINITY) t_1 t_3)))))
double code(double x, double y, double z, double t) {
double t_1 = ((2.0 / z) - -2.0) / t;
double t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -4e+86) {
tmp = t_1;
} else if (t_2 <= 500000.0) {
tmp = t_3;
} else if (t_2 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = t_3;
}
return tmp;
}
public static double code(double x, double y, double z, double t) {
double t_1 = ((2.0 / z) - -2.0) / t;
double t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -4e+86) {
tmp = t_1;
} else if (t_2 <= 500000.0) {
tmp = t_3;
} else if (t_2 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = t_3;
}
return tmp;
}
def code(x, y, z, t): t_1 = ((2.0 / z) - -2.0) / t t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z) t_3 = -2.0 + (x / y) tmp = 0 if t_2 <= -4e+86: tmp = t_1 elif t_2 <= 500000.0: tmp = t_3 elif t_2 <= math.inf: tmp = t_1 else: tmp = t_3 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(Float64(2.0 / z) - -2.0) / t) t_2 = Float64(Float64(2.0 - Float64(Float64(t - 1.0) * Float64(z * 2.0))) / Float64(t * z)) t_3 = Float64(-2.0 + Float64(x / y)) tmp = 0.0 if (t_2 <= -4e+86) tmp = t_1; elseif (t_2 <= 500000.0) tmp = t_3; elseif (t_2 <= Inf) tmp = t_1; else tmp = t_3; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = ((2.0 / z) - -2.0) / t; t_2 = (2.0 - ((t - 1.0) * (z * 2.0))) / (t * z); t_3 = -2.0 + (x / y); tmp = 0.0; if (t_2 <= -4e+86) tmp = t_1; elseif (t_2 <= 500000.0) tmp = t_3; elseif (t_2 <= Inf) tmp = t_1; else tmp = t_3; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision]}, Block[{t$95$2 = N[(N[(2.0 - N[(N[(t - 1.0), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -4e+86], t$95$1, If[LessEqual[t$95$2, 500000.0], t$95$3, If[LessEqual[t$95$2, Infinity], t$95$1, t$95$3]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\frac{2}{z} - -2}{t}\\
t_2 := \frac{2 - \left(t - 1\right) \cdot \left(z \cdot 2\right)}{t \cdot z}\\
t_3 := -2 + \frac{x}{y}\\
\mathbf{if}\;t\_2 \leq -4 \cdot 10^{+86}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 500000:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -4.0000000000000001e86 or 5e5 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < +inf.0Initial program 97.6%
Taylor expanded in t around 0
lower-/.f64N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6484.2
Applied rewrites84.2%
if -4.0000000000000001e86 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 5e5 or +inf.0 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) Initial program 71.4%
Taylor expanded in t around inf
Applied rewrites95.4%
Final simplification89.4%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ (fma 2.0 z 2.0) (* t z)) (/ x y))))
(if (<= (/ x y) -50.0)
t_1
(if (<= (/ x y) 0.0005) (- (/ (- (/ 2.0 z) -2.0) t) 2.0) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (fma(2.0, z, 2.0) / (t * z)) + (x / y);
double tmp;
if ((x / y) <= -50.0) {
tmp = t_1;
} else if ((x / y) <= 0.0005) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t) t_1 = Float64(Float64(fma(2.0, z, 2.0) / Float64(t * z)) + Float64(x / y)) tmp = 0.0 if (Float64(x / y) <= -50.0) tmp = t_1; elseif (Float64(x / y) <= 0.0005) tmp = Float64(Float64(Float64(Float64(2.0 / z) - -2.0) / t) - 2.0); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(N[(2.0 * z + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x / y), $MachinePrecision], -50.0], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 0.0005], N[(N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision] - 2.0), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\mathsf{fma}\left(2, z, 2\right)}{t \cdot z} + \frac{x}{y}\\
\mathbf{if}\;\frac{x}{y} \leq -50:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 0.0005:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 x y) < -50 or 5.0000000000000001e-4 < (/.f64 x y) Initial program 85.7%
Taylor expanded in t around 0
+-commutativeN/A
lower-fma.f6496.3
Applied rewrites96.3%
if -50 < (/.f64 x y) < 5.0000000000000001e-4Initial program 85.1%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites98.7%
Final simplification97.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ 2.0 (* t z)) (/ x y))))
(if (<= (/ x y) -50.0)
t_1
(if (<= (/ x y) 1e+24) (- (/ (- (/ 2.0 z) -2.0) t) 2.0) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (2.0 / (t * z)) + (x / y);
double tmp;
if ((x / y) <= -50.0) {
tmp = t_1;
} else if ((x / y) <= 1e+24) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (2.0d0 / (t * z)) + (x / y)
if ((x / y) <= (-50.0d0)) then
tmp = t_1
else if ((x / y) <= 1d+24) then
tmp = (((2.0d0 / z) - (-2.0d0)) / t) - 2.0d0
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (2.0 / (t * z)) + (x / y);
double tmp;
if ((x / y) <= -50.0) {
tmp = t_1;
} else if ((x / y) <= 1e+24) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (2.0 / (t * z)) + (x / y) tmp = 0 if (x / y) <= -50.0: tmp = t_1 elif (x / y) <= 1e+24: tmp = (((2.0 / z) - -2.0) / t) - 2.0 else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(2.0 / Float64(t * z)) + Float64(x / y)) tmp = 0.0 if (Float64(x / y) <= -50.0) tmp = t_1; elseif (Float64(x / y) <= 1e+24) tmp = Float64(Float64(Float64(Float64(2.0 / z) - -2.0) / t) - 2.0); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (2.0 / (t * z)) + (x / y); tmp = 0.0; if ((x / y) <= -50.0) tmp = t_1; elseif ((x / y) <= 1e+24) tmp = (((2.0 / z) - -2.0) / t) - 2.0; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(2.0 / N[(t * z), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x / y), $MachinePrecision], -50.0], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 1e+24], N[(N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision] - 2.0), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{t \cdot z} + \frac{x}{y}\\
\mathbf{if}\;\frac{x}{y} \leq -50:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 10^{+24}:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 x y) < -50 or 9.9999999999999998e23 < (/.f64 x y) Initial program 84.9%
Taylor expanded in z around 0
Applied rewrites87.5%
if -50 < (/.f64 x y) < 9.9999999999999998e23Initial program 85.9%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites98.0%
Final simplification93.1%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (- -2.0 (/ -2.0 t)) (/ x y))))
(if (<= (/ x y) -4000000.0)
t_1
(if (<= (/ x y) 1e+24) (- (/ (- (/ 2.0 z) -2.0) t) 2.0) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (-2.0 - (-2.0 / t)) + (x / y);
double tmp;
if ((x / y) <= -4000000.0) {
tmp = t_1;
} else if ((x / y) <= 1e+24) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = ((-2.0d0) - ((-2.0d0) / t)) + (x / y)
if ((x / y) <= (-4000000.0d0)) then
tmp = t_1
else if ((x / y) <= 1d+24) then
tmp = (((2.0d0 / z) - (-2.0d0)) / t) - 2.0d0
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (-2.0 - (-2.0 / t)) + (x / y);
double tmp;
if ((x / y) <= -4000000.0) {
tmp = t_1;
} else if ((x / y) <= 1e+24) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (-2.0 - (-2.0 / t)) + (x / y) tmp = 0 if (x / y) <= -4000000.0: tmp = t_1 elif (x / y) <= 1e+24: tmp = (((2.0 / z) - -2.0) / t) - 2.0 else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(-2.0 - Float64(-2.0 / t)) + Float64(x / y)) tmp = 0.0 if (Float64(x / y) <= -4000000.0) tmp = t_1; elseif (Float64(x / y) <= 1e+24) tmp = Float64(Float64(Float64(Float64(2.0 / z) - -2.0) / t) - 2.0); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (-2.0 - (-2.0 / t)) + (x / y); tmp = 0.0; if ((x / y) <= -4000000.0) tmp = t_1; elseif ((x / y) <= 1e+24) tmp = (((2.0 / z) - -2.0) / t) - 2.0; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(-2.0 - N[(-2.0 / t), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x / y), $MachinePrecision], -4000000.0], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 1e+24], N[(N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision] - 2.0), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(-2 - \frac{-2}{t}\right) + \frac{x}{y}\\
\mathbf{if}\;\frac{x}{y} \leq -4000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 10^{+24}:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 x y) < -4e6 or 9.9999999999999998e23 < (/.f64 x y) Initial program 84.8%
Taylor expanded in z around inf
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
metadata-evalN/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6480.9
Applied rewrites80.9%
if -4e6 < (/.f64 x y) < 9.9999999999999998e23Initial program 86.0%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites98.0%
Final simplification90.1%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -230.0) (+ -2.0 (/ x y)) (if (<= (/ x y) 1.65e+24) (- (/ 2.0 t) 2.0) (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -230.0) {
tmp = -2.0 + (x / y);
} else if ((x / y) <= 1.65e+24) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((x / y) <= (-230.0d0)) then
tmp = (-2.0d0) + (x / y)
else if ((x / y) <= 1.65d+24) then
tmp = (2.0d0 / t) - 2.0d0
else
tmp = x / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -230.0) {
tmp = -2.0 + (x / y);
} else if ((x / y) <= 1.65e+24) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -230.0: tmp = -2.0 + (x / y) elif (x / y) <= 1.65e+24: tmp = (2.0 / t) - 2.0 else: tmp = x / y return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -230.0) tmp = Float64(-2.0 + Float64(x / y)); elseif (Float64(x / y) <= 1.65e+24) tmp = Float64(Float64(2.0 / t) - 2.0); else tmp = Float64(x / y); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((x / y) <= -230.0) tmp = -2.0 + (x / y); elseif ((x / y) <= 1.65e+24) tmp = (2.0 / t) - 2.0; else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -230.0], N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 1.65e+24], N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision], N[(x / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -230:\\
\;\;\;\;-2 + \frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 1.65 \cdot 10^{+24}:\\
\;\;\;\;\frac{2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -230Initial program 85.4%
Taylor expanded in t around inf
Applied rewrites64.5%
if -230 < (/.f64 x y) < 1.6499999999999999e24Initial program 86.0%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites98.0%
Taylor expanded in z around inf
Applied rewrites64.1%
if 1.6499999999999999e24 < (/.f64 x y) Initial program 84.2%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites24.2%
Taylor expanded in z around 0
Applied rewrites20.7%
Taylor expanded in x around inf
lower-/.f6479.4
Applied rewrites79.4%
Final simplification67.6%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -19000.0) (/ x y) (if (<= (/ x y) 1.65e+24) (- (/ 2.0 t) 2.0) (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -19000.0) {
tmp = x / y;
} else if ((x / y) <= 1.65e+24) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((x / y) <= (-19000.0d0)) then
tmp = x / y
else if ((x / y) <= 1.65d+24) then
tmp = (2.0d0 / t) - 2.0d0
else
tmp = x / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -19000.0) {
tmp = x / y;
} else if ((x / y) <= 1.65e+24) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -19000.0: tmp = x / y elif (x / y) <= 1.65e+24: tmp = (2.0 / t) - 2.0 else: tmp = x / y return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -19000.0) tmp = Float64(x / y); elseif (Float64(x / y) <= 1.65e+24) tmp = Float64(Float64(2.0 / t) - 2.0); else tmp = Float64(x / y); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((x / y) <= -19000.0) tmp = x / y; elseif ((x / y) <= 1.65e+24) tmp = (2.0 / t) - 2.0; else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -19000.0], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 1.65e+24], N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision], N[(x / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -19000:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 1.65 \cdot 10^{+24}:\\
\;\;\;\;\frac{2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -19000 or 1.6499999999999999e24 < (/.f64 x y) Initial program 84.8%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites31.6%
Taylor expanded in z around 0
Applied rewrites23.2%
Taylor expanded in x around inf
lower-/.f6470.6
Applied rewrites70.6%
if -19000 < (/.f64 x y) < 1.6499999999999999e24Initial program 86.0%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites98.0%
Taylor expanded in z around inf
Applied rewrites64.1%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -85.0) (/ x y) (if (<= (/ x y) 2700000.0) -2.0 (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -85.0) {
tmp = x / y;
} else if ((x / y) <= 2700000.0) {
tmp = -2.0;
} else {
tmp = x / y;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((x / y) <= (-85.0d0)) then
tmp = x / y
else if ((x / y) <= 2700000.0d0) then
tmp = -2.0d0
else
tmp = x / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -85.0) {
tmp = x / y;
} else if ((x / y) <= 2700000.0) {
tmp = -2.0;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -85.0: tmp = x / y elif (x / y) <= 2700000.0: tmp = -2.0 else: tmp = x / y return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -85.0) tmp = Float64(x / y); elseif (Float64(x / y) <= 2700000.0) tmp = -2.0; else tmp = Float64(x / y); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((x / y) <= -85.0) tmp = x / y; elseif ((x / y) <= 2700000.0) tmp = -2.0; else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -85.0], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 2700000.0], -2.0, N[(x / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -85:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 2700000:\\
\;\;\;\;-2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -85 or 2.7e6 < (/.f64 x y) Initial program 85.5%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites34.1%
Taylor expanded in z around 0
Applied rewrites24.3%
Taylor expanded in x around inf
lower-/.f6467.9
Applied rewrites67.9%
if -85 < (/.f64 x y) < 2.7e6Initial program 85.4%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites98.7%
Taylor expanded in t around inf
Applied rewrites38.7%
(FPCore (x y z t) :precision binary64 -2.0)
double code(double x, double y, double z, double t) {
return -2.0;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = -2.0d0
end function
public static double code(double x, double y, double z, double t) {
return -2.0;
}
def code(x, y, z, t): return -2.0
function code(x, y, z, t) return -2.0 end
function tmp = code(x, y, z, t) tmp = -2.0; end
code[x_, y_, z_, t_] := -2.0
\begin{array}{l}
\\
-2
\end{array}
Initial program 85.4%
Taylor expanded in x around 0
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate-*r/N/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
Applied rewrites67.1%
Taylor expanded in t around inf
Applied rewrites21.2%
(FPCore (x y z t) :precision binary64 (- (/ (+ (/ 2.0 z) 2.0) t) (- 2.0 (/ x y))))
double code(double x, double y, double z, double t) {
return (((2.0 / z) + 2.0) / t) - (2.0 - (x / y));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((2.0d0 / z) + 2.0d0) / t) - (2.0d0 - (x / y))
end function
public static double code(double x, double y, double z, double t) {
return (((2.0 / z) + 2.0) / t) - (2.0 - (x / y));
}
def code(x, y, z, t): return (((2.0 / z) + 2.0) / t) - (2.0 - (x / y))
function code(x, y, z, t) return Float64(Float64(Float64(Float64(2.0 / z) + 2.0) / t) - Float64(2.0 - Float64(x / y))) end
function tmp = code(x, y, z, t) tmp = (((2.0 / z) + 2.0) / t) - (2.0 - (x / y)); end
code[x_, y_, z_, t_] := N[(N[(N[(N[(2.0 / z), $MachinePrecision] + 2.0), $MachinePrecision] / t), $MachinePrecision] - N[(2.0 - N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{2}{z} + 2}{t} - \left(2 - \frac{x}{y}\right)
\end{array}
herbie shell --seed 2024298
(FPCore (x y z t)
:name "Data.HashTable.ST.Basic:computeOverhead from hashtables-1.2.0.2"
:precision binary64
:alt
(! :herbie-platform default (- (/ (+ (/ 2 z) 2) t) (- 2 (/ x y))))
(+ (/ x y) (/ (+ 2.0 (* (* z 2.0) (- 1.0 t))) (* t z))))