
(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 (if (<= (+ (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z)) (/ x y)) INFINITY) (fma (fma (fma 2.0 t -2.0) z -2.0) (/ -1.0 (* t z)) (/ x y)) (+ -2.0 (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((((((1.0 - t) * (z * 2.0)) + 2.0) / (t * z)) + (x / y)) <= ((double) INFINITY)) {
tmp = fma(fma(fma(2.0, t, -2.0), z, -2.0), (-1.0 / (t * z)), (x / y));
} else {
tmp = -2.0 + (x / y);
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if (Float64(Float64(Float64(Float64(Float64(1.0 - t) * Float64(z * 2.0)) + 2.0) / Float64(t * z)) + Float64(x / y)) <= Inf) tmp = fma(fma(fma(2.0, t, -2.0), z, -2.0), Float64(-1.0 / Float64(t * z)), Float64(x / y)); else tmp = Float64(-2.0 + Float64(x / y)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[N[(N[(N[(N[(N[(1.0 - t), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(2.0 * t + -2.0), $MachinePrecision] * z + -2.0), $MachinePrecision] * N[(-1.0 / N[(t * z), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision], N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\left(1 - t\right) \cdot \left(z \cdot 2\right) + 2}{t \cdot z} + \frac{x}{y} \leq \infty:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, t, -2\right), z, -2\right), \frac{-1}{t \cdot z}, \frac{x}{y}\right)\\
\mathbf{else}:\\
\;\;\;\;-2 + \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.4%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites99.4%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6499.4
Applied rewrites99.4%
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 t around inf
Applied rewrites100.0%
Final simplification99.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ 2.0 (* t z)))
(t_2 (+ -2.0 (/ x y)))
(t_3 (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z))))
(if (<= t_3 -2e+148)
t_1
(if (<= t_3 -1.0)
t_2
(if (<= t_3 5e+294)
(- (/ 2.0 t) 2.0)
(if (<= t_3 INFINITY) t_1 t_2))))))
double code(double x, double y, double z, double t) {
double t_1 = 2.0 / (t * z);
double t_2 = -2.0 + (x / y);
double t_3 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z);
double tmp;
if (t_3 <= -2e+148) {
tmp = t_1;
} else if (t_3 <= -1.0) {
tmp = t_2;
} else if (t_3 <= 5e+294) {
tmp = (2.0 / t) - 2.0;
} else if (t_3 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = t_2;
}
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 + (x / y);
double t_3 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z);
double tmp;
if (t_3 <= -2e+148) {
tmp = t_1;
} else if (t_3 <= -1.0) {
tmp = t_2;
} else if (t_3 <= 5e+294) {
tmp = (2.0 / t) - 2.0;
} else if (t_3 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t): t_1 = 2.0 / (t * z) t_2 = -2.0 + (x / y) t_3 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z) tmp = 0 if t_3 <= -2e+148: tmp = t_1 elif t_3 <= -1.0: tmp = t_2 elif t_3 <= 5e+294: tmp = (2.0 / t) - 2.0 elif t_3 <= math.inf: tmp = t_1 else: tmp = t_2 return tmp
function code(x, y, z, t) t_1 = Float64(2.0 / Float64(t * z)) t_2 = Float64(-2.0 + Float64(x / y)) t_3 = Float64(Float64(Float64(Float64(1.0 - t) * Float64(z * 2.0)) + 2.0) / Float64(t * z)) tmp = 0.0 if (t_3 <= -2e+148) tmp = t_1; elseif (t_3 <= -1.0) tmp = t_2; elseif (t_3 <= 5e+294) tmp = Float64(Float64(2.0 / t) - 2.0); elseif (t_3 <= Inf) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 2.0 / (t * z); t_2 = -2.0 + (x / y); t_3 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z); tmp = 0.0; if (t_3 <= -2e+148) tmp = t_1; elseif (t_3 <= -1.0) tmp = t_2; elseif (t_3 <= 5e+294) tmp = (2.0 / t) - 2.0; elseif (t_3 <= Inf) tmp = t_1; else tmp = t_2; 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[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(N[(1.0 - t), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -2e+148], t$95$1, If[LessEqual[t$95$3, -1.0], t$95$2, If[LessEqual[t$95$3, 5e+294], N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision], If[LessEqual[t$95$3, Infinity], t$95$1, t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{t \cdot z}\\
t_2 := -2 + \frac{x}{y}\\
t_3 := \frac{\left(1 - t\right) \cdot \left(z \cdot 2\right) + 2}{t \cdot z}\\
\mathbf{if}\;t\_3 \leq -2 \cdot 10^{+148}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_3 \leq -1:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_3 \leq 5 \cdot 10^{+294}:\\
\;\;\;\;\frac{2}{t} - 2\\
\mathbf{elif}\;t\_3 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\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)) < -2.0000000000000001e148 or 4.9999999999999999e294 < (/.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.7%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6479.6
Applied rewrites79.6%
if -2.0000000000000001e148 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -1 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 66.4%
Taylor expanded in t around inf
Applied rewrites89.1%
if -1 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 4.9999999999999999e294Initial program 99.6%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6428.1
Applied rewrites28.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
+-commutativeN/A
associate-*r/N/A
metadata-evalN/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
distribute-lft-outN/A
Applied rewrites82.5%
Taylor expanded in z around inf
Applied rewrites54.5%
Final simplification78.9%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (- (/ 2.0 z) -2.0) t))
(t_2 (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z)))
(t_3 (+ -2.0 (/ x y))))
(if (<= t_2 -1e+26)
t_1
(if (<= t_2 -2.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -1e+26) {
tmp = t_1;
} else if (t_2 <= -2.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -1e+26) {
tmp = t_1;
} else if (t_2 <= -2.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z) t_3 = -2.0 + (x / y) tmp = 0 if t_2 <= -1e+26: tmp = t_1 elif t_2 <= -2.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(Float64(Float64(1.0 - t) * Float64(z * 2.0)) + 2.0) / Float64(t * z)) t_3 = Float64(-2.0 + Float64(x / y)) tmp = 0.0 if (t_2 <= -1e+26) tmp = t_1; elseif (t_2 <= -2.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z); t_3 = -2.0 + (x / y); tmp = 0.0; if (t_2 <= -1e+26) tmp = t_1; elseif (t_2 <= -2.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[(N[(N[(1.0 - t), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e+26], t$95$1, If[LessEqual[t$95$2, -2.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{\left(1 - t\right) \cdot \left(z \cdot 2\right) + 2}{t \cdot z}\\
t_3 := -2 + \frac{x}{y}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq -2:\\
\;\;\;\;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)) < -1.00000000000000005e26Initial program 94.5%
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-/.f6477.7
Applied rewrites77.7%
if -1.00000000000000005e26 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -2 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 57.7%
Taylor expanded in t around inf
Applied rewrites99.1%
if -2 < (/.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 98.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
associate-+r+N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
associate-*r/N/A
metadata-evalN/A
associate--l+N/A
lower--.f64N/A
Applied rewrites87.2%
Final simplification89.3%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z)))
(t_2 (+ -2.0 (/ x y))))
(if (<= t_1 -1e+26)
(/ (- (/ 2.0 z) -2.0) t)
(if (<= t_1 -2.0)
t_2
(if (<= t_1 INFINITY) (/ (fma z (fma -2.0 t 2.0) 2.0) (* t z)) t_2)))))
double code(double x, double y, double z, double t) {
double t_1 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z);
double t_2 = -2.0 + (x / y);
double tmp;
if (t_1 <= -1e+26) {
tmp = ((2.0 / z) - -2.0) / t;
} else if (t_1 <= -2.0) {
tmp = t_2;
} else if (t_1 <= ((double) INFINITY)) {
tmp = fma(z, fma(-2.0, t, 2.0), 2.0) / (t * z);
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t) t_1 = Float64(Float64(Float64(Float64(1.0 - t) * Float64(z * 2.0)) + 2.0) / Float64(t * z)) t_2 = Float64(-2.0 + Float64(x / y)) tmp = 0.0 if (t_1 <= -1e+26) tmp = Float64(Float64(Float64(2.0 / z) - -2.0) / t); elseif (t_1 <= -2.0) tmp = t_2; elseif (t_1 <= Inf) tmp = Float64(fma(z, fma(-2.0, t, 2.0), 2.0) / Float64(t * z)); else tmp = t_2; end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(N[(N[(1.0 - t), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e+26], N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision], If[LessEqual[t$95$1, -2.0], t$95$2, If[LessEqual[t$95$1, Infinity], N[(N[(z * N[(-2.0 * t + 2.0), $MachinePrecision] + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\left(1 - t\right) \cdot \left(z \cdot 2\right) + 2}{t \cdot z}\\
t_2 := -2 + \frac{x}{y}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{+26}:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t}\\
\mathbf{elif}\;t\_1 \leq -2:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, \mathsf{fma}\left(-2, t, 2\right), 2\right)}{t \cdot z}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\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.00000000000000005e26Initial program 94.5%
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-/.f6477.7
Applied rewrites77.7%
if -1.00000000000000005e26 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -2 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 57.7%
Taylor expanded in t around inf
Applied rewrites99.1%
if -2 < (/.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 98.5%
Taylor expanded in z around 0
Applied rewrites57.9%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
clear-numN/A
frac-addN/A
lower-/.f64N/A
lower-fma.f64N/A
lower-/.f64N/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6443.8
Applied rewrites43.8%
lift-/.f64N/A
lift-*.f64N/A
*-rgt-identityN/A
lift-*.f64N/A
associate-/r*N/A
div-invN/A
lift-/.f64N/A
clear-numN/A
lift-/.f64N/A
lower-*.f64N/A
Applied rewrites43.7%
Taylor expanded in x around 0
lower-/.f64N/A
+-commutativeN/A
associate-*r*N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-inN/A
mul-1-negN/A
distribute-lft-neg-inN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-neg-inN/A
metadata-evalN/A
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-outN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6485.8
Applied rewrites85.8%
Final simplification88.9%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ (- (/ 2.0 z) -2.0) t))
(t_2 (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z)))
(t_3 (+ -2.0 (/ x y))))
(if (<= t_2 -1e+26)
t_1
(if (<= t_2 -1.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -1e+26) {
tmp = t_1;
} else if (t_2 <= -1.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z);
double t_3 = -2.0 + (x / y);
double tmp;
if (t_2 <= -1e+26) {
tmp = t_1;
} else if (t_2 <= -1.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z) t_3 = -2.0 + (x / y) tmp = 0 if t_2 <= -1e+26: tmp = t_1 elif t_2 <= -1.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(Float64(Float64(1.0 - t) * Float64(z * 2.0)) + 2.0) / Float64(t * z)) t_3 = Float64(-2.0 + Float64(x / y)) tmp = 0.0 if (t_2 <= -1e+26) tmp = t_1; elseif (t_2 <= -1.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 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z); t_3 = -2.0 + (x / y); tmp = 0.0; if (t_2 <= -1e+26) tmp = t_1; elseif (t_2 <= -1.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[(N[(N[(1.0 - t), $MachinePrecision] * N[(z * 2.0), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e+26], t$95$1, If[LessEqual[t$95$2, -1.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{\left(1 - t\right) \cdot \left(z \cdot 2\right) + 2}{t \cdot z}\\
t_3 := -2 + \frac{x}{y}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq -1:\\
\;\;\;\;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)) < -1.00000000000000005e26 or -1 < (/.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 96.5%
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-/.f6481.8
Applied rewrites81.8%
if -1.00000000000000005e26 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < -1 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 58.1%
Taylor expanded in t around inf
Applied rewrites99.0%
Final simplification88.9%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ (fma 2.0 z 2.0) (* t z)) (/ x y))))
(if (<= (/ x y) -5e+23)
t_1
(if (<= (/ x y) 0.001)
(fma (/ 1.0 t) (- (/ 2.0 z) -2.0) -2.0)
(if (<= (/ x y) 1e+307) t_1 (/ x y))))))
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) <= -5e+23) {
tmp = t_1;
} else if ((x / y) <= 0.001) {
tmp = fma((1.0 / t), ((2.0 / z) - -2.0), -2.0);
} else if ((x / y) <= 1e+307) {
tmp = t_1;
} else {
tmp = x / y;
}
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) <= -5e+23) tmp = t_1; elseif (Float64(x / y) <= 0.001) tmp = fma(Float64(1.0 / t), Float64(Float64(2.0 / z) - -2.0), -2.0); elseif (Float64(x / y) <= 1e+307) tmp = t_1; else tmp = Float64(x / y); 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], -5e+23], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 0.001], N[(N[(1.0 / t), $MachinePrecision] * N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] + -2.0), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 1e+307], t$95$1, N[(x / y), $MachinePrecision]]]]]
\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 -5 \cdot 10^{+23}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 0.001:\\
\;\;\;\;\mathsf{fma}\left(\frac{1}{t}, \frac{2}{z} - -2, -2\right)\\
\mathbf{elif}\;\frac{x}{y} \leq 10^{+307}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -4.9999999999999999e23 or 1e-3 < (/.f64 x y) < 9.99999999999999986e306Initial program 81.3%
Taylor expanded in t around 0
+-commutativeN/A
lower-fma.f6499.7
Applied rewrites99.7%
if -4.9999999999999999e23 < (/.f64 x y) < 1e-3Initial program 82.3%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6437.2
Applied rewrites37.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
+-commutativeN/A
associate-*r/N/A
metadata-evalN/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
distribute-lft-outN/A
Applied rewrites98.6%
if 9.99999999999999986e306 < (/.f64 x y) Initial program 61.5%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites61.5%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6461.5
Applied rewrites61.5%
Taylor expanded in x around inf
lower-/.f64100.0
Applied rewrites100.0%
Final simplification99.2%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ 2.0 (* t z)) (/ x y))))
(if (<= (/ x y) -4.1e+93)
t_1
(if (<= (/ x y) 780000000.0)
(fma (/ 1.0 t) (- (/ 2.0 z) -2.0) -2.0)
(if (<= (/ x y) INFINITY) t_1 (/ x y))))))
double code(double x, double y, double z, double t) {
double t_1 = (2.0 / (t * z)) + (x / y);
double tmp;
if ((x / y) <= -4.1e+93) {
tmp = t_1;
} else if ((x / y) <= 780000000.0) {
tmp = fma((1.0 / t), ((2.0 / z) - -2.0), -2.0);
} else if ((x / y) <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = x / y;
}
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) <= -4.1e+93) tmp = t_1; elseif (Float64(x / y) <= 780000000.0) tmp = fma(Float64(1.0 / t), Float64(Float64(2.0 / z) - -2.0), -2.0); elseif (Float64(x / y) <= Inf) tmp = t_1; else tmp = Float64(x / y); end return 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], -4.1e+93], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 780000000.0], N[(N[(1.0 / t), $MachinePrecision] * N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] + -2.0), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], Infinity], t$95$1, N[(x / y), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{t \cdot z} + \frac{x}{y}\\
\mathbf{if}\;\frac{x}{y} \leq -4.1 \cdot 10^{+93}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 780000000:\\
\;\;\;\;\mathsf{fma}\left(\frac{1}{t}, \frac{2}{z} - -2, -2\right)\\
\mathbf{elif}\;\frac{x}{y} \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -4.1000000000000001e93 or 7.8e8 < (/.f64 x y) < +inf.0Initial program 78.7%
Taylor expanded in z around 0
Applied rewrites87.9%
if -4.1000000000000001e93 < (/.f64 x y) < 7.8e8Initial program 82.4%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6435.6
Applied rewrites35.6%
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
+-commutativeN/A
associate-*r/N/A
metadata-evalN/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
distribute-lft-outN/A
Applied rewrites95.6%
if +inf.0 < (/.f64 x y) Initial program 80.7%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites80.7%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6480.7
Applied rewrites80.7%
Taylor expanded in x around inf
lower-/.f6435.8
Applied rewrites35.8%
Final simplification92.2%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ 2.0 (* t z)) (/ x y))))
(if (<= (/ x y) -4.1e+93)
t_1
(if (<= (/ x y) 780000000.0)
(- (/ (- (/ 2.0 z) -2.0) t) 2.0)
(if (<= (/ x y) INFINITY) t_1 (/ x y))))))
double code(double x, double y, double z, double t) {
double t_1 = (2.0 / (t * z)) + (x / y);
double tmp;
if ((x / y) <= -4.1e+93) {
tmp = t_1;
} else if ((x / y) <= 780000000.0) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else if ((x / y) <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = x / y;
}
return tmp;
}
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) <= -4.1e+93) {
tmp = t_1;
} else if ((x / y) <= 780000000.0) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else if ((x / y) <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): t_1 = (2.0 / (t * z)) + (x / y) tmp = 0 if (x / y) <= -4.1e+93: tmp = t_1 elif (x / y) <= 780000000.0: tmp = (((2.0 / z) - -2.0) / t) - 2.0 elif (x / y) <= math.inf: tmp = t_1 else: tmp = x / y 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) <= -4.1e+93) tmp = t_1; elseif (Float64(x / y) <= 780000000.0) tmp = Float64(Float64(Float64(Float64(2.0 / z) - -2.0) / t) - 2.0); elseif (Float64(x / y) <= Inf) tmp = t_1; else tmp = Float64(x / y); 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) <= -4.1e+93) tmp = t_1; elseif ((x / y) <= 780000000.0) tmp = (((2.0 / z) - -2.0) / t) - 2.0; elseif ((x / y) <= Inf) tmp = t_1; else tmp = x / y; 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], -4.1e+93], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 780000000.0], N[(N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision] - 2.0), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], Infinity], t$95$1, N[(x / y), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{t \cdot z} + \frac{x}{y}\\
\mathbf{if}\;\frac{x}{y} \leq -4.1 \cdot 10^{+93}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 780000000:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t} - 2\\
\mathbf{elif}\;\frac{x}{y} \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -4.1000000000000001e93 or 7.8e8 < (/.f64 x y) < +inf.0Initial program 78.7%
Taylor expanded in z around 0
Applied rewrites87.9%
if -4.1000000000000001e93 < (/.f64 x y) < 7.8e8Initial program 82.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
associate-+r+N/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
associate-*r/N/A
metadata-evalN/A
associate--l+N/A
lower--.f64N/A
Applied rewrites95.5%
if +inf.0 < (/.f64 x y) Initial program 80.7%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites80.7%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6480.7
Applied rewrites80.7%
Taylor expanded in x around inf
lower-/.f6435.8
Applied rewrites35.8%
Final simplification92.2%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -4.1e+93) (/ x y) (if (<= (/ x y) 16500000.0) (- (/ 2.0 t) 2.0) (+ -2.0 (/ x y)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -4.1e+93) {
tmp = x / y;
} else if ((x / y) <= 16500000.0) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = -2.0 + (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) <= (-4.1d+93)) then
tmp = x / y
else if ((x / y) <= 16500000.0d0) then
tmp = (2.0d0 / t) - 2.0d0
else
tmp = (-2.0d0) + (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) <= -4.1e+93) {
tmp = x / y;
} else if ((x / y) <= 16500000.0) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = -2.0 + (x / y);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -4.1e+93: tmp = x / y elif (x / y) <= 16500000.0: tmp = (2.0 / t) - 2.0 else: tmp = -2.0 + (x / y) return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -4.1e+93) tmp = Float64(x / y); elseif (Float64(x / y) <= 16500000.0) tmp = Float64(Float64(2.0 / t) - 2.0); else tmp = Float64(-2.0 + Float64(x / y)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((x / y) <= -4.1e+93) tmp = x / y; elseif ((x / y) <= 16500000.0) tmp = (2.0 / t) - 2.0; else tmp = -2.0 + (x / y); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -4.1e+93], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 16500000.0], N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision], N[(-2.0 + N[(x / y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -4.1 \cdot 10^{+93}:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 16500000:\\
\;\;\;\;\frac{2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;-2 + \frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -4.1000000000000001e93Initial program 79.6%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites79.6%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6479.6
Applied rewrites79.6%
Taylor expanded in x around inf
lower-/.f6476.2
Applied rewrites76.2%
if -4.1000000000000001e93 < (/.f64 x y) < 1.65e7Initial program 82.4%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6435.6
Applied rewrites35.6%
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
+-commutativeN/A
associate-*r/N/A
metadata-evalN/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
distribute-lft-outN/A
Applied rewrites95.6%
Taylor expanded in z around inf
Applied rewrites62.3%
if 1.65e7 < (/.f64 x y) Initial program 78.0%
Taylor expanded in t around inf
Applied rewrites71.2%
Final simplification67.2%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -4.1e+93) (/ x y) (if (<= (/ x y) 780000000.0) (- (/ 2.0 t) 2.0) (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -4.1e+93) {
tmp = x / y;
} else if ((x / y) <= 780000000.0) {
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) <= (-4.1d+93)) then
tmp = x / y
else if ((x / y) <= 780000000.0d0) 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) <= -4.1e+93) {
tmp = x / y;
} else if ((x / y) <= 780000000.0) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -4.1e+93: tmp = x / y elif (x / y) <= 780000000.0: 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) <= -4.1e+93) tmp = Float64(x / y); elseif (Float64(x / y) <= 780000000.0) 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) <= -4.1e+93) tmp = x / y; elseif ((x / y) <= 780000000.0) 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], -4.1e+93], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 780000000.0], 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 -4.1 \cdot 10^{+93}:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 780000000:\\
\;\;\;\;\frac{2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -4.1000000000000001e93 or 7.8e8 < (/.f64 x y) Initial program 78.7%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites78.7%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6478.7
Applied rewrites78.7%
Taylor expanded in x around inf
lower-/.f6473.1
Applied rewrites73.1%
if -4.1000000000000001e93 < (/.f64 x y) < 7.8e8Initial program 82.4%
Taylor expanded in z around 0
lower-/.f64N/A
lower-*.f6435.6
Applied rewrites35.6%
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
+-commutativeN/A
associate-*r/N/A
metadata-evalN/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
distribute-lft-outN/A
Applied rewrites95.6%
Taylor expanded in z around inf
Applied rewrites62.3%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -4.1e+93) (/ x y) (if (<= (/ x y) 780000000.0) (/ 2.0 t) (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -4.1e+93) {
tmp = x / y;
} else if ((x / y) <= 780000000.0) {
tmp = 2.0 / t;
} 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) <= (-4.1d+93)) then
tmp = x / y
else if ((x / y) <= 780000000.0d0) then
tmp = 2.0d0 / t
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) <= -4.1e+93) {
tmp = x / y;
} else if ((x / y) <= 780000000.0) {
tmp = 2.0 / t;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -4.1e+93: tmp = x / y elif (x / y) <= 780000000.0: tmp = 2.0 / t else: tmp = x / y return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -4.1e+93) tmp = Float64(x / y); elseif (Float64(x / y) <= 780000000.0) tmp = Float64(2.0 / t); else tmp = Float64(x / y); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((x / y) <= -4.1e+93) tmp = x / y; elseif ((x / y) <= 780000000.0) tmp = 2.0 / t; else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -4.1e+93], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 780000000.0], N[(2.0 / t), $MachinePrecision], N[(x / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -4.1 \cdot 10^{+93}:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 780000000:\\
\;\;\;\;\frac{2}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -4.1000000000000001e93 or 7.8e8 < (/.f64 x y) Initial program 78.7%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites78.7%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6478.7
Applied rewrites78.7%
Taylor expanded in x around inf
lower-/.f6473.1
Applied rewrites73.1%
if -4.1000000000000001e93 < (/.f64 x y) < 7.8e8Initial program 82.4%
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-/.f6463.3
Applied rewrites63.3%
Taylor expanded in z around 0
Applied rewrites35.6%
Taylor expanded in z around inf
Applied rewrites30.4%
(FPCore (x y z t) :precision binary64 (/ x y))
double code(double x, double y, double z, double t) {
return 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 = x / y
end function
public static double code(double x, double y, double z, double t) {
return x / y;
}
def code(x, y, z, t): return x / y
function code(x, y, z, t) return Float64(x / y) end
function tmp = code(x, y, z, t) tmp = x / y; end
code[x_, y_, z_, t_] := N[(x / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y}
\end{array}
Initial program 80.7%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
lower-fma.f64N/A
Applied rewrites80.7%
Taylor expanded in z around 0
sub-negN/A
metadata-evalN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6480.7
Applied rewrites80.7%
Taylor expanded in x around inf
lower-/.f6435.8
Applied rewrites35.8%
(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 2024296
(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))))