
(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 (+ (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z)) (/ x y)))) (if (<= t_1 1e+286) t_1 (/ (fma (- (/ (- (/ 2.0 z) -2.0) t) 2.0) y x) y))))
double code(double x, double y, double z, double t) {
double t_1 = ((((1.0 - t) * (z * 2.0)) + 2.0) / (t * z)) + (x / y);
double tmp;
if (t_1 <= 1e+286) {
tmp = t_1;
} else {
tmp = fma(((((2.0 / z) - -2.0) / t) - 2.0), y, x) / y;
}
return tmp;
}
function code(x, y, z, t) t_1 = Float64(Float64(Float64(Float64(Float64(1.0 - t) * Float64(z * 2.0)) + 2.0) / Float64(t * z)) + Float64(x / y)) tmp = 0.0 if (t_1 <= 1e+286) tmp = t_1; else tmp = Float64(fma(Float64(Float64(Float64(Float64(2.0 / z) - -2.0) / t) - 2.0), y, x) / y); end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = 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]}, If[LessEqual[t$95$1, 1e+286], t$95$1, N[(N[(N[(N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision] - 2.0), $MachinePrecision] * y + x), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\left(1 - t\right) \cdot \left(z \cdot 2\right) + 2}{t \cdot z} + \frac{x}{y}\\
\mathbf{if}\;t\_1 \leq 10^{+286}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\frac{\frac{2}{z} - -2}{t} - 2, y, x\right)}{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))) < 1.00000000000000003e286Initial program 99.8%
if 1.00000000000000003e286 < (+.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 54.7%
Taylor expanded in y around 0
Applied rewrites100.0%
Final simplification99.9%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ (fma z 2.0 2.0) (* t z)) (/ x y)))
(t_2 (+ (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z)) (/ x y)))
(t_3 (/ (fma (- (/ (- (/ 2.0 z) -2.0) t) 2.0) y x) y)))
(if (<= t_2 -1e+41)
t_1
(if (<= t_2 10000000000.0) t_3 (if (<= t_2 1e+286) t_1 t_3)))))
double code(double x, double y, double z, double t) {
double t_1 = (fma(z, 2.0, 2.0) / (t * z)) + (x / y);
double t_2 = ((((1.0 - t) * (z * 2.0)) + 2.0) / (t * z)) + (x / y);
double t_3 = fma(((((2.0 / z) - -2.0) / t) - 2.0), y, x) / y;
double tmp;
if (t_2 <= -1e+41) {
tmp = t_1;
} else if (t_2 <= 10000000000.0) {
tmp = t_3;
} else if (t_2 <= 1e+286) {
tmp = t_1;
} else {
tmp = t_3;
}
return tmp;
}
function code(x, y, z, t) t_1 = Float64(Float64(fma(z, 2.0, 2.0) / Float64(t * z)) + Float64(x / y)) t_2 = Float64(Float64(Float64(Float64(Float64(1.0 - t) * Float64(z * 2.0)) + 2.0) / Float64(t * z)) + Float64(x / y)) t_3 = Float64(fma(Float64(Float64(Float64(Float64(2.0 / z) - -2.0) / t) - 2.0), y, x) / y) tmp = 0.0 if (t_2 <= -1e+41) tmp = t_1; elseif (t_2 <= 10000000000.0) tmp = t_3; elseif (t_2 <= 1e+286) tmp = t_1; else tmp = t_3; end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(N[(z * 2.0 + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = 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]}, Block[{t$95$3 = N[(N[(N[(N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision] - 2.0), $MachinePrecision] * y + x), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[t$95$2, -1e+41], t$95$1, If[LessEqual[t$95$2, 10000000000.0], t$95$3, If[LessEqual[t$95$2, 1e+286], t$95$1, t$95$3]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\mathsf{fma}\left(z, 2, 2\right)}{t \cdot z} + \frac{x}{y}\\
t_2 := \frac{\left(1 - t\right) \cdot \left(z \cdot 2\right) + 2}{t \cdot z} + \frac{x}{y}\\
t_3 := \frac{\mathsf{fma}\left(\frac{\frac{2}{z} - -2}{t} - 2, y, x\right)}{y}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{+41}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 10000000000:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq 10^{+286}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\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))) < -1.00000000000000001e41 or 1e10 < (+.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))) < 1.00000000000000003e286Initial program 99.8%
Taylor expanded in t around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.7
Applied rewrites99.7%
if -1.00000000000000001e41 < (+.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))) < 1e10 or 1.00000000000000003e286 < (+.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 70.6%
Taylor expanded in y around 0
Applied rewrites99.9%
Final simplification99.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ 2.0 (* t z)))
(t_2 (/ (+ (* (- 1.0 t) (* z 2.0)) 2.0) (* t z)))
(t_3 (+ -2.0 (/ x y))))
(if (<= t_2 (- INFINITY))
t_1
(if (<= t_2 1e+181) 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 / (t * z);
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 <= -((double) INFINITY)) {
tmp = t_1;
} else if (t_2 <= 1e+181) {
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 / (t * z);
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 <= -Double.POSITIVE_INFINITY) {
tmp = t_1;
} else if (t_2 <= 1e+181) {
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 / (t * z) t_2 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z) t_3 = -2.0 + (x / y) tmp = 0 if t_2 <= -math.inf: tmp = t_1 elif t_2 <= 1e+181: 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(2.0 / Float64(t * z)) 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 <= Float64(-Inf)) tmp = t_1; elseif (t_2 <= 1e+181) 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 / (t * z); t_2 = (((1.0 - t) * (z * 2.0)) + 2.0) / (t * z); t_3 = -2.0 + (x / y); tmp = 0.0; if (t_2 <= -Inf) tmp = t_1; elseif (t_2 <= 1e+181) 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[(2.0 / N[(t * z), $MachinePrecision]), $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, (-Infinity)], t$95$1, If[LessEqual[t$95$2, 1e+181], t$95$3, If[LessEqual[t$95$2, Infinity], t$95$1, t$95$3]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{t \cdot z}\\
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 -\infty:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 10^{+181}:\\
\;\;\;\;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)) < -inf.0 or 9.9999999999999992e180 < (/.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 95.2%
Taylor expanded in z around 0
associate-/r*N/A
metadata-evalN/A
associate-*r/N/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6482.3
Applied rewrites82.3%
Applied rewrites82.4%
if -inf.0 < (/.f64 (+.f64 #s(literal 2 binary64) (*.f64 (*.f64 z #s(literal 2 binary64)) (-.f64 #s(literal 1 binary64) t))) (*.f64 t z)) < 9.9999999999999992e180 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 86.3%
Taylor expanded in t around inf
Applied rewrites68.0%
Final simplification71.7%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ (fma z 2.0 2.0) (* t z)) (/ x y))))
(if (<= (/ x y) -5e-8)
t_1
(if (<= (/ x y) 1e-8) (- (/ (- (/ 2.0 z) -2.0) t) 2.0) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (fma(z, 2.0, 2.0) / (t * z)) + (x / y);
double tmp;
if ((x / y) <= -5e-8) {
tmp = t_1;
} else if ((x / y) <= 1e-8) {
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(z, 2.0, 2.0) / Float64(t * z)) + Float64(x / y)) tmp = 0.0 if (Float64(x / y) <= -5e-8) tmp = t_1; elseif (Float64(x / y) <= 1e-8) 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[(z * 2.0 + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x / y), $MachinePrecision], -5e-8], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 1e-8], 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(z, 2, 2\right)}{t \cdot z} + \frac{x}{y}\\
\mathbf{if}\;\frac{x}{y} \leq -5 \cdot 10^{-8}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 10^{-8}:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 x y) < -4.9999999999999998e-8 or 1e-8 < (/.f64 x y) Initial program 89.0%
Taylor expanded in t around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f6496.9
Applied rewrites96.9%
if -4.9999999999999998e-8 < (/.f64 x y) < 1e-8Initial program 88.0%
Taylor expanded in y around inf
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 rewrites99.5%
Final simplification98.1%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ 2.0 (* t z)) (/ x y))))
(if (<= (/ x y) -5e+85)
t_1
(if (<= (/ x y) 1e+37) (- (/ (- (/ 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) <= -5e+85) {
tmp = t_1;
} else if ((x / y) <= 1e+37) {
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) <= (-5d+85)) then
tmp = t_1
else if ((x / y) <= 1d+37) 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) <= -5e+85) {
tmp = t_1;
} else if ((x / y) <= 1e+37) {
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) <= -5e+85: tmp = t_1 elif (x / y) <= 1e+37: 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) <= -5e+85) tmp = t_1; elseif (Float64(x / y) <= 1e+37) 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) <= -5e+85) tmp = t_1; elseif ((x / y) <= 1e+37) 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], -5e+85], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 1e+37], 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 -5 \cdot 10^{+85}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 10^{+37}:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 x y) < -5.0000000000000001e85 or 9.99999999999999954e36 < (/.f64 x y) Initial program 86.2%
Taylor expanded in z around 0
Applied rewrites93.8%
if -5.0000000000000001e85 < (/.f64 x y) < 9.99999999999999954e36Initial program 90.5%
Taylor expanded in y around inf
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 rewrites96.1%
Final simplification95.1%
(FPCore (x y z t)
:precision binary64
(if (<= (/ x y) -1.8e+93)
(+ (/ 2.0 t) (/ x y))
(if (<= (/ x y) 102000000.0)
(- (/ (- (/ 2.0 z) -2.0) t) 2.0)
(+ (- (/ 2.0 t) 2.0) (/ x y)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -1.8e+93) {
tmp = (2.0 / t) + (x / y);
} else if ((x / y) <= 102000000.0) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else {
tmp = ((2.0 / t) - 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) <= (-1.8d+93)) then
tmp = (2.0d0 / t) + (x / y)
else if ((x / y) <= 102000000.0d0) then
tmp = (((2.0d0 / z) - (-2.0d0)) / t) - 2.0d0
else
tmp = ((2.0d0 / t) - 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) <= -1.8e+93) {
tmp = (2.0 / t) + (x / y);
} else if ((x / y) <= 102000000.0) {
tmp = (((2.0 / z) - -2.0) / t) - 2.0;
} else {
tmp = ((2.0 / t) - 2.0) + (x / y);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -1.8e+93: tmp = (2.0 / t) + (x / y) elif (x / y) <= 102000000.0: tmp = (((2.0 / z) - -2.0) / t) - 2.0 else: tmp = ((2.0 / t) - 2.0) + (x / y) return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -1.8e+93) tmp = Float64(Float64(2.0 / t) + Float64(x / y)); elseif (Float64(x / y) <= 102000000.0) tmp = Float64(Float64(Float64(Float64(2.0 / z) - -2.0) / t) - 2.0); else tmp = Float64(Float64(Float64(2.0 / t) - 2.0) + Float64(x / y)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((x / y) <= -1.8e+93) tmp = (2.0 / t) + (x / y); elseif ((x / y) <= 102000000.0) tmp = (((2.0 / z) - -2.0) / t) - 2.0; else tmp = ((2.0 / t) - 2.0) + (x / y); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -1.8e+93], N[(N[(2.0 / t), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 102000000.0], N[(N[(N[(N[(2.0 / z), $MachinePrecision] - -2.0), $MachinePrecision] / t), $MachinePrecision] - 2.0), $MachinePrecision], N[(N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -1.8 \cdot 10^{+93}:\\
\;\;\;\;\frac{2}{t} + \frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 102000000:\\
\;\;\;\;\frac{\frac{2}{z} - -2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{2}{t} - 2\right) + \frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -1.8e93Initial program 86.2%
Taylor expanded in z around inf
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6490.6
Applied rewrites90.6%
Taylor expanded in t around 0
Applied rewrites90.6%
if -1.8e93 < (/.f64 x y) < 1.02e8Initial program 90.4%
Taylor expanded in y around inf
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 rewrites96.8%
if 1.02e8 < (/.f64 x y) Initial program 86.3%
Taylor expanded in z around inf
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6479.2
Applied rewrites79.2%
Final simplification91.0%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ (/ 2.0 t) (/ x y))))
(if (<= (/ x y) -1.8e+93)
t_1
(if (<= (/ x y) 6.2e+70) (- (/ (fma z 2.0 2.0) (* t z)) 2.0) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (2.0 / t) + (x / y);
double tmp;
if ((x / y) <= -1.8e+93) {
tmp = t_1;
} else if ((x / y) <= 6.2e+70) {
tmp = (fma(z, 2.0, 2.0) / (t * z)) - 2.0;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t) t_1 = Float64(Float64(2.0 / t) + Float64(x / y)) tmp = 0.0 if (Float64(x / y) <= -1.8e+93) tmp = t_1; elseif (Float64(x / y) <= 6.2e+70) tmp = Float64(Float64(fma(z, 2.0, 2.0) / Float64(t * z)) - 2.0); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(2.0 / t), $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x / y), $MachinePrecision], -1.8e+93], t$95$1, If[LessEqual[N[(x / y), $MachinePrecision], 6.2e+70], N[(N[(N[(z * 2.0 + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{t} + \frac{x}{y}\\
\mathbf{if}\;\frac{x}{y} \leq -1.8 \cdot 10^{+93}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{x}{y} \leq 6.2 \cdot 10^{+70}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, 2, 2\right)}{t \cdot z} - 2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 x y) < -1.8e93 or 6.2000000000000006e70 < (/.f64 x y) Initial program 86.4%
Taylor expanded in z around inf
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6485.1
Applied rewrites85.1%
Taylor expanded in t around 0
Applied rewrites85.1%
if -1.8e93 < (/.f64 x y) < 6.2000000000000006e70Initial program 90.1%
Taylor expanded in y around inf
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%
Taylor expanded in z around 0
Applied rewrites95.6%
Final simplification91.0%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -7.2e+93) (/ x y) (if (<= (/ x y) 2.5e+79) (- (/ (fma z 2.0 2.0) (* t z)) 2.0) (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -7.2e+93) {
tmp = x / y;
} else if ((x / y) <= 2.5e+79) {
tmp = (fma(z, 2.0, 2.0) / (t * z)) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -7.2e+93) tmp = Float64(x / y); elseif (Float64(x / y) <= 2.5e+79) tmp = Float64(Float64(fma(z, 2.0, 2.0) / Float64(t * z)) - 2.0); else tmp = Float64(x / y); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -7.2e+93], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 2.5e+79], N[(N[(N[(z * 2.0 + 2.0), $MachinePrecision] / N[(t * z), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision], N[(x / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -7.2 \cdot 10^{+93}:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 2.5 \cdot 10^{+79}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z, 2, 2\right)}{t \cdot z} - 2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -7.1999999999999998e93 or 2.5e79 < (/.f64 x y) Initial program 86.3%
Taylor expanded in y around 0
lower-/.f6482.0
Applied rewrites82.0%
if -7.1999999999999998e93 < (/.f64 x y) < 2.5e79Initial program 90.2%
Taylor expanded in y around inf
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.6%
Taylor expanded in z around 0
Applied rewrites95.6%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -1.8e+93) (/ x y) (if (<= (/ x y) 7e+77) (- (/ 2.0 (* t z)) 2.0) (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -1.8e+93) {
tmp = x / y;
} else if ((x / y) <= 7e+77) {
tmp = (2.0 / (t * z)) - 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) <= (-1.8d+93)) then
tmp = x / y
else if ((x / y) <= 7d+77) then
tmp = (2.0d0 / (t * z)) - 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) <= -1.8e+93) {
tmp = x / y;
} else if ((x / y) <= 7e+77) {
tmp = (2.0 / (t * z)) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -1.8e+93: tmp = x / y elif (x / y) <= 7e+77: tmp = (2.0 / (t * z)) - 2.0 else: tmp = x / y return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -1.8e+93) tmp = Float64(x / y); elseif (Float64(x / y) <= 7e+77) tmp = Float64(Float64(2.0 / Float64(t * z)) - 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) <= -1.8e+93) tmp = x / y; elseif ((x / y) <= 7e+77) tmp = (2.0 / (t * z)) - 2.0; else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -1.8e+93], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 7e+77], N[(N[(2.0 / N[(t * z), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision], N[(x / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -1.8 \cdot 10^{+93}:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 7 \cdot 10^{+77}:\\
\;\;\;\;\frac{2}{t \cdot z} - 2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -1.8e93 or 7.0000000000000003e77 < (/.f64 x y) Initial program 86.3%
Taylor expanded in y around 0
lower-/.f6482.0
Applied rewrites82.0%
if -1.8e93 < (/.f64 x y) < 7.0000000000000003e77Initial program 90.2%
Taylor expanded in y around inf
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.6%
Taylor expanded in z around 0
Applied rewrites95.6%
Taylor expanded in z around 0
Applied rewrites65.3%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -1.05e+90) (/ x y) (if (<= (/ x y) 1.4e+45) (- (/ 2.0 t) 2.0) (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -1.05e+90) {
tmp = x / y;
} else if ((x / y) <= 1.4e+45) {
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) <= (-1.05d+90)) then
tmp = x / y
else if ((x / y) <= 1.4d+45) 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) <= -1.05e+90) {
tmp = x / y;
} else if ((x / y) <= 1.4e+45) {
tmp = (2.0 / t) - 2.0;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -1.05e+90: tmp = x / y elif (x / y) <= 1.4e+45: 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) <= -1.05e+90) tmp = Float64(x / y); elseif (Float64(x / y) <= 1.4e+45) 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) <= -1.05e+90) tmp = x / y; elseif ((x / y) <= 1.4e+45) 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], -1.05e+90], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 1.4e+45], 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 -1.05 \cdot 10^{+90}:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 1.4 \cdot 10^{+45}:\\
\;\;\;\;\frac{2}{t} - 2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -1.0499999999999999e90 or 1.4e45 < (/.f64 x y) Initial program 85.9%
Taylor expanded in y around 0
lower-/.f6480.1
Applied rewrites80.1%
if -1.0499999999999999e90 < (/.f64 x y) < 1.4e45Initial program 90.6%
Taylor expanded in y around inf
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 rewrites96.2%
Taylor expanded in z around inf
Applied rewrites56.6%
(FPCore (x y z t) :precision binary64 (if (<= (/ x y) -1.05) (/ x y) (if (<= (/ x y) 9e-9) -2.0 (/ x y))))
double code(double x, double y, double z, double t) {
double tmp;
if ((x / y) <= -1.05) {
tmp = x / y;
} else if ((x / y) <= 9e-9) {
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) <= (-1.05d0)) then
tmp = x / y
else if ((x / y) <= 9d-9) 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) <= -1.05) {
tmp = x / y;
} else if ((x / y) <= 9e-9) {
tmp = -2.0;
} else {
tmp = x / y;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (x / y) <= -1.05: tmp = x / y elif (x / y) <= 9e-9: tmp = -2.0 else: tmp = x / y return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(x / y) <= -1.05) tmp = Float64(x / y); elseif (Float64(x / y) <= 9e-9) 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) <= -1.05) tmp = x / y; elseif ((x / y) <= 9e-9) tmp = -2.0; else tmp = x / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(x / y), $MachinePrecision], -1.05], N[(x / y), $MachinePrecision], If[LessEqual[N[(x / y), $MachinePrecision], 9e-9], -2.0, N[(x / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{y} \leq -1.05:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{elif}\;\frac{x}{y} \leq 9 \cdot 10^{-9}:\\
\;\;\;\;-2\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y}\\
\end{array}
\end{array}
if (/.f64 x y) < -1.05000000000000004 or 8.99999999999999953e-9 < (/.f64 x y) Initial program 88.9%
Taylor expanded in y around 0
lower-/.f6465.9
Applied rewrites65.9%
if -1.05000000000000004 < (/.f64 x y) < 8.99999999999999953e-9Initial program 88.1%
Taylor expanded in y around inf
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 rewrites99.5%
Taylor expanded in t around inf
Applied rewrites33.9%
(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 88.5%
Taylor expanded in y around inf
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 rewrites63.6%
Taylor expanded in t around inf
Applied rewrites16.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 2024277
(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))))