
(FPCore (x y z t a b) :precision binary64 (+ (- (+ (+ x y) z) (* z (log t))) (* (- a 0.5) b)))
double code(double x, double y, double z, double t, double a, double b) {
return (((x + y) + z) - (z * log(t))) + ((a - 0.5) * b);
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (((x + y) + z) - (z * log(t))) + ((a - 0.5d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (((x + y) + z) - (z * Math.log(t))) + ((a - 0.5) * b);
}
def code(x, y, z, t, a, b): return (((x + y) + z) - (z * math.log(t))) + ((a - 0.5) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(Float64(x + y) + z) - Float64(z * log(t))) + Float64(Float64(a - 0.5) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = (((x + y) + z) - (z * log(t))) + ((a - 0.5) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(N[(x + y), $MachinePrecision] + z), $MachinePrecision] - N[(z * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(a - 0.5), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(x + y\right) + z\right) - z \cdot \log t\right) + \left(a - 0.5\right) \cdot b
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (+ (- (+ (+ x y) z) (* z (log t))) (* (- a 0.5) b)))
double code(double x, double y, double z, double t, double a, double b) {
return (((x + y) + z) - (z * log(t))) + ((a - 0.5) * b);
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (((x + y) + z) - (z * log(t))) + ((a - 0.5d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (((x + y) + z) - (z * Math.log(t))) + ((a - 0.5) * b);
}
def code(x, y, z, t, a, b): return (((x + y) + z) - (z * math.log(t))) + ((a - 0.5) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(Float64(x + y) + z) - Float64(z * log(t))) + Float64(Float64(a - 0.5) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = (((x + y) + z) - (z * log(t))) + ((a - 0.5) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(N[(x + y), $MachinePrecision] + z), $MachinePrecision] - N[(z * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(a - 0.5), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(x + y\right) + z\right) - z \cdot \log t\right) + \left(a - 0.5\right) \cdot b
\end{array}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (+ (* z (- 1.0 (log t))) (+ (+ (* b (+ a -0.5)) x) y)))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
return (z * (1.0 - log(t))) + (((b * (a + -0.5)) + x) + y);
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (z * (1.0d0 - log(t))) + (((b * (a + (-0.5d0))) + x) + y)
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
return (z * (1.0 - Math.log(t))) + (((b * (a + -0.5)) + x) + y);
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): return (z * (1.0 - math.log(t))) + (((b * (a + -0.5)) + x) + y)
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) return Float64(Float64(z * Float64(1.0 - log(t))) + Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y)) end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp = code(x, y, z, t, a, b)
tmp = (z * (1.0 - log(t))) + (((b * (a + -0.5)) + x) + y);
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := N[(N[(z * N[(1.0 - N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
z \cdot \left(1 - \log t\right) + \left(\left(b \cdot \left(a + -0.5\right) + x\right) + y\right)
\end{array}
Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
sub-negN/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6499.9%
Applied egg-rr99.9%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (if (<= x -1.05e+14) (+ (+ (* b (+ a -0.5)) x) y) (+ (- (+ z y) (* z (log t))) (* b (- a 0.5)))))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= -1.05e+14) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = ((z + y) - (z * log(t))) + (b * (a - 0.5));
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (x <= (-1.05d+14)) then
tmp = ((b * (a + (-0.5d0))) + x) + y
else
tmp = ((z + y) - (z * log(t))) + (b * (a - 0.5d0))
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= -1.05e+14) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = ((z + y) - (z * Math.log(t))) + (b * (a - 0.5));
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): tmp = 0 if x <= -1.05e+14: tmp = ((b * (a + -0.5)) + x) + y else: tmp = ((z + y) - (z * math.log(t))) + (b * (a - 0.5)) return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) tmp = 0.0 if (x <= -1.05e+14) tmp = Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y); else tmp = Float64(Float64(Float64(z + y) - Float64(z * log(t))) + Float64(b * Float64(a - 0.5))); end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
tmp = 0.0;
if (x <= -1.05e+14)
tmp = ((b * (a + -0.5)) + x) + y;
else
tmp = ((z + y) - (z * log(t))) + (b * (a - 0.5));
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, -1.05e+14], N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision], N[(N[(N[(z + y), $MachinePrecision] - N[(z * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.05 \cdot 10^{+14}:\\
\;\;\;\;\left(b \cdot \left(a + -0.5\right) + x\right) + y\\
\mathbf{else}:\\
\;\;\;\;\left(\left(z + y\right) - z \cdot \log t\right) + b \cdot \left(a - 0.5\right)\\
\end{array}
\end{array}
if x < -1.05e14Initial program 99.8%
Taylor expanded in z around 0
+-lowering-+.f6487.7%
Simplified87.7%
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6487.7%
Applied egg-rr87.7%
if -1.05e14 < x Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f6486.5%
Simplified86.5%
Final simplification86.8%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
(FPCore (x y z t a b)
:precision binary64
(if (<= z -1.25e+119)
(+ (* b (- a 0.5)) (- z (* z (log t))))
(if (<= z 3.4e+226)
(+ (+ (* b (+ a -0.5)) x) y)
(+ (* z (- 1.0 (log t))) y))))assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -1.25e+119) {
tmp = (b * (a - 0.5)) + (z - (z * log(t)));
} else if (z <= 3.4e+226) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = (z * (1.0 - log(t))) + y;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (z <= (-1.25d+119)) then
tmp = (b * (a - 0.5d0)) + (z - (z * log(t)))
else if (z <= 3.4d+226) then
tmp = ((b * (a + (-0.5d0))) + x) + y
else
tmp = (z * (1.0d0 - log(t))) + y
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -1.25e+119) {
tmp = (b * (a - 0.5)) + (z - (z * Math.log(t)));
} else if (z <= 3.4e+226) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = (z * (1.0 - Math.log(t))) + y;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): tmp = 0 if z <= -1.25e+119: tmp = (b * (a - 0.5)) + (z - (z * math.log(t))) elif z <= 3.4e+226: tmp = ((b * (a + -0.5)) + x) + y else: tmp = (z * (1.0 - math.log(t))) + y return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -1.25e+119) tmp = Float64(Float64(b * Float64(a - 0.5)) + Float64(z - Float64(z * log(t)))); elseif (z <= 3.4e+226) tmp = Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y); else tmp = Float64(Float64(z * Float64(1.0 - log(t))) + y); end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
tmp = 0.0;
if (z <= -1.25e+119)
tmp = (b * (a - 0.5)) + (z - (z * log(t)));
elseif (z <= 3.4e+226)
tmp = ((b * (a + -0.5)) + x) + y;
else
tmp = (z * (1.0 - log(t))) + y;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -1.25e+119], N[(N[(b * N[(a - 0.5), $MachinePrecision]), $MachinePrecision] + N[(z - N[(z * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.4e+226], N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision], N[(N[(z * N[(1.0 - N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.25 \cdot 10^{+119}:\\
\;\;\;\;b \cdot \left(a - 0.5\right) + \left(z - z \cdot \log t\right)\\
\mathbf{elif}\;z \leq 3.4 \cdot 10^{+226}:\\
\;\;\;\;\left(b \cdot \left(a + -0.5\right) + x\right) + y\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(1 - \log t\right) + y\\
\end{array}
\end{array}
if z < -1.25e119Initial program 99.6%
Taylor expanded in z around inf
Simplified89.2%
if -1.25e119 < z < 3.39999999999999979e226Initial program 100.0%
Taylor expanded in z around 0
+-lowering-+.f6493.1%
Simplified93.1%
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6493.1%
Applied egg-rr93.1%
if 3.39999999999999979e226 < z Initial program 99.4%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around inf
Simplified81.4%
Final simplification91.7%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* z (- 1.0 (log t)))))
(if (<= z -2.9e+167)
(+ t_1 (* b a))
(if (<= z 3.4e+226) (+ (+ (* b (+ a -0.5)) x) y) (+ t_1 y)))))assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - log(t));
double tmp;
if (z <= -2.9e+167) {
tmp = t_1 + (b * a);
} else if (z <= 3.4e+226) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1 + y;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = z * (1.0d0 - log(t))
if (z <= (-2.9d+167)) then
tmp = t_1 + (b * a)
else if (z <= 3.4d+226) then
tmp = ((b * (a + (-0.5d0))) + x) + y
else
tmp = t_1 + y
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - Math.log(t));
double tmp;
if (z <= -2.9e+167) {
tmp = t_1 + (b * a);
} else if (z <= 3.4e+226) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1 + y;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): t_1 = z * (1.0 - math.log(t)) tmp = 0 if z <= -2.9e+167: tmp = t_1 + (b * a) elif z <= 3.4e+226: tmp = ((b * (a + -0.5)) + x) + y else: tmp = t_1 + y return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) t_1 = Float64(z * Float64(1.0 - log(t))) tmp = 0.0 if (z <= -2.9e+167) tmp = Float64(t_1 + Float64(b * a)); elseif (z <= 3.4e+226) tmp = Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y); else tmp = Float64(t_1 + y); end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
t_1 = z * (1.0 - log(t));
tmp = 0.0;
if (z <= -2.9e+167)
tmp = t_1 + (b * a);
elseif (z <= 3.4e+226)
tmp = ((b * (a + -0.5)) + x) + y;
else
tmp = t_1 + y;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z * N[(1.0 - N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -2.9e+167], N[(t$95$1 + N[(b * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.4e+226], N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision], N[(t$95$1 + y), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - \log t\right)\\
\mathbf{if}\;z \leq -2.9 \cdot 10^{+167}:\\
\;\;\;\;t\_1 + b \cdot a\\
\mathbf{elif}\;z \leq 3.4 \cdot 10^{+226}:\\
\;\;\;\;\left(b \cdot \left(a + -0.5\right) + x\right) + y\\
\mathbf{else}:\\
\;\;\;\;t\_1 + y\\
\end{array}
\end{array}
if z < -2.89999999999999975e167Initial program 99.6%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.7%
Simplified99.7%
Taylor expanded in a around inf
*-commutativeN/A
*-lowering-*.f6485.4%
Simplified85.4%
if -2.89999999999999975e167 < z < 3.39999999999999979e226Initial program 99.9%
Taylor expanded in z around 0
+-lowering-+.f6491.2%
Simplified91.2%
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6491.2%
Applied egg-rr91.2%
if 3.39999999999999979e226 < z Initial program 99.4%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around inf
Simplified81.4%
Final simplification89.9%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ (* z (- 1.0 (log t))) y)))
(if (<= z -1.85e+124)
t_1
(if (<= z 5.6e+226) (+ (+ (* b (+ a -0.5)) x) y) t_1))))assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z * (1.0 - log(t))) + y;
double tmp;
if (z <= -1.85e+124) {
tmp = t_1;
} else if (z <= 5.6e+226) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = (z * (1.0d0 - log(t))) + y
if (z <= (-1.85d+124)) then
tmp = t_1
else if (z <= 5.6d+226) then
tmp = ((b * (a + (-0.5d0))) + x) + y
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (z * (1.0 - Math.log(t))) + y;
double tmp;
if (z <= -1.85e+124) {
tmp = t_1;
} else if (z <= 5.6e+226) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): t_1 = (z * (1.0 - math.log(t))) + y tmp = 0 if z <= -1.85e+124: tmp = t_1 elif z <= 5.6e+226: tmp = ((b * (a + -0.5)) + x) + y else: tmp = t_1 return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) t_1 = Float64(Float64(z * Float64(1.0 - log(t))) + y) tmp = 0.0 if (z <= -1.85e+124) tmp = t_1; elseif (z <= 5.6e+226) tmp = Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y); else tmp = t_1; end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
t_1 = (z * (1.0 - log(t))) + y;
tmp = 0.0;
if (z <= -1.85e+124)
tmp = t_1;
elseif (z <= 5.6e+226)
tmp = ((b * (a + -0.5)) + x) + y;
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(z * N[(1.0 - N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]}, If[LessEqual[z, -1.85e+124], t$95$1, If[LessEqual[z, 5.6e+226], N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - \log t\right) + y\\
\mathbf{if}\;z \leq -1.85 \cdot 10^{+124}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 5.6 \cdot 10^{+226}:\\
\;\;\;\;\left(b \cdot \left(a + -0.5\right) + x\right) + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.85000000000000004e124 or 5.6000000000000005e226 < z Initial program 99.6%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around inf
Simplified70.5%
if -1.85000000000000004e124 < z < 5.6000000000000005e226Initial program 100.0%
Taylor expanded in z around 0
+-lowering-+.f6492.7%
Simplified92.7%
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6492.7%
Applied egg-rr92.7%
Final simplification87.9%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* z (- 1.0 (log t)))))
(if (<= z -1.34e+225)
t_1
(if (<= z 1.65e+254) (+ (+ (* b (+ a -0.5)) x) y) (+ t_1 x)))))assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - log(t));
double tmp;
if (z <= -1.34e+225) {
tmp = t_1;
} else if (z <= 1.65e+254) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1 + x;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = z * (1.0d0 - log(t))
if (z <= (-1.34d+225)) then
tmp = t_1
else if (z <= 1.65d+254) then
tmp = ((b * (a + (-0.5d0))) + x) + y
else
tmp = t_1 + x
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - Math.log(t));
double tmp;
if (z <= -1.34e+225) {
tmp = t_1;
} else if (z <= 1.65e+254) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1 + x;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): t_1 = z * (1.0 - math.log(t)) tmp = 0 if z <= -1.34e+225: tmp = t_1 elif z <= 1.65e+254: tmp = ((b * (a + -0.5)) + x) + y else: tmp = t_1 + x return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) t_1 = Float64(z * Float64(1.0 - log(t))) tmp = 0.0 if (z <= -1.34e+225) tmp = t_1; elseif (z <= 1.65e+254) tmp = Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y); else tmp = Float64(t_1 + x); end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
t_1 = z * (1.0 - log(t));
tmp = 0.0;
if (z <= -1.34e+225)
tmp = t_1;
elseif (z <= 1.65e+254)
tmp = ((b * (a + -0.5)) + x) + y;
else
tmp = t_1 + x;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z * N[(1.0 - N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.34e+225], t$95$1, If[LessEqual[z, 1.65e+254], N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision], N[(t$95$1 + x), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - \log t\right)\\
\mathbf{if}\;z \leq -1.34 \cdot 10^{+225}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.65 \cdot 10^{+254}:\\
\;\;\;\;\left(b \cdot \left(a + -0.5\right) + x\right) + y\\
\mathbf{else}:\\
\;\;\;\;t\_1 + x\\
\end{array}
\end{array}
if z < -1.33999999999999995e225Initial program 99.6%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6484.9%
Simplified84.9%
if -1.33999999999999995e225 < z < 1.64999999999999996e254Initial program 99.9%
Taylor expanded in z around 0
+-lowering-+.f6488.1%
Simplified88.1%
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6488.1%
Applied egg-rr88.1%
if 1.64999999999999996e254 < z Initial program 99.4%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf
Simplified88.8%
Final simplification87.9%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* z (- 1.0 (log t)))))
(if (<= z -5.9e+225)
t_1
(if (<= z 9.2e+254) (+ (+ (* b (+ a -0.5)) x) y) t_1))))assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - log(t));
double tmp;
if (z <= -5.9e+225) {
tmp = t_1;
} else if (z <= 9.2e+254) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = z * (1.0d0 - log(t))
if (z <= (-5.9d+225)) then
tmp = t_1
else if (z <= 9.2d+254) then
tmp = ((b * (a + (-0.5d0))) + x) + y
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = z * (1.0 - Math.log(t));
double tmp;
if (z <= -5.9e+225) {
tmp = t_1;
} else if (z <= 9.2e+254) {
tmp = ((b * (a + -0.5)) + x) + y;
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): t_1 = z * (1.0 - math.log(t)) tmp = 0 if z <= -5.9e+225: tmp = t_1 elif z <= 9.2e+254: tmp = ((b * (a + -0.5)) + x) + y else: tmp = t_1 return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) t_1 = Float64(z * Float64(1.0 - log(t))) tmp = 0.0 if (z <= -5.9e+225) tmp = t_1; elseif (z <= 9.2e+254) tmp = Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y); else tmp = t_1; end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
t_1 = z * (1.0 - log(t));
tmp = 0.0;
if (z <= -5.9e+225)
tmp = t_1;
elseif (z <= 9.2e+254)
tmp = ((b * (a + -0.5)) + x) + y;
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(z * N[(1.0 - N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -5.9e+225], t$95$1, If[LessEqual[z, 9.2e+254], N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
t_1 := z \cdot \left(1 - \log t\right)\\
\mathbf{if}\;z \leq -5.9 \cdot 10^{+225}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 9.2 \cdot 10^{+254}:\\
\;\;\;\;\left(b \cdot \left(a + -0.5\right) + x\right) + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -5.8999999999999998e225 or 9.19999999999999994e254 < z Initial program 99.6%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6486.4%
Simplified86.4%
if -5.8999999999999998e225 < z < 9.19999999999999994e254Initial program 99.9%
Taylor expanded in z around 0
+-lowering-+.f6488.1%
Simplified88.1%
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6488.1%
Applied egg-rr88.1%
Final simplification87.9%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (+ (* z (- 1.0 (log t))) (+ (* b (+ a -0.5)) (+ x y))))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
return (z * (1.0 - log(t))) + ((b * (a + -0.5)) + (x + y));
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (z * (1.0d0 - log(t))) + ((b * (a + (-0.5d0))) + (x + y))
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
return (z * (1.0 - Math.log(t))) + ((b * (a + -0.5)) + (x + y));
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): return (z * (1.0 - math.log(t))) + ((b * (a + -0.5)) + (x + y))
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) return Float64(Float64(z * Float64(1.0 - log(t))) + Float64(Float64(b * Float64(a + -0.5)) + Float64(x + y))) end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp = code(x, y, z, t, a, b)
tmp = (z * (1.0 - log(t))) + ((b * (a + -0.5)) + (x + y));
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := N[(N[(z * N[(1.0 - N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + N[(x + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
z \cdot \left(1 - \log t\right) + \left(b \cdot \left(a + -0.5\right) + \left(x + y\right)\right)
\end{array}
Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
Final simplification99.9%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* b (+ a -0.5)))) (if (<= b -8e+39) t_1 (if (<= b 2e+22) (+ x y) t_1))))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (a + -0.5);
double tmp;
if (b <= -8e+39) {
tmp = t_1;
} else if (b <= 2e+22) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = b * (a + (-0.5d0))
if (b <= (-8d+39)) then
tmp = t_1
else if (b <= 2d+22) then
tmp = x + y
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (a + -0.5);
double tmp;
if (b <= -8e+39) {
tmp = t_1;
} else if (b <= 2e+22) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): t_1 = b * (a + -0.5) tmp = 0 if b <= -8e+39: tmp = t_1 elif b <= 2e+22: tmp = x + y else: tmp = t_1 return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) t_1 = Float64(b * Float64(a + -0.5)) tmp = 0.0 if (b <= -8e+39) tmp = t_1; elseif (b <= 2e+22) tmp = Float64(x + y); else tmp = t_1; end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
t_1 = b * (a + -0.5);
tmp = 0.0;
if (b <= -8e+39)
tmp = t_1;
elseif (b <= 2e+22)
tmp = x + y;
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -8e+39], t$95$1, If[LessEqual[b, 2e+22], N[(x + y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
t_1 := b \cdot \left(a + -0.5\right)\\
\mathbf{if}\;b \leq -8 \cdot 10^{+39}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2 \cdot 10^{+22}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -7.99999999999999952e39 or 2e22 < b Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6471.8%
Simplified71.8%
if -7.99999999999999952e39 < b < 2e22Initial program 99.8%
Taylor expanded in z around 0
+-lowering-+.f6473.2%
Simplified73.2%
Taylor expanded in b around 0
+-commutativeN/A
+-lowering-+.f6462.1%
Simplified62.1%
Final simplification66.4%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* b (- a 0.5)))) (if (<= (+ x y) -4e-85) (+ x t_1) (+ y t_1))))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (a - 0.5);
double tmp;
if ((x + y) <= -4e-85) {
tmp = x + t_1;
} else {
tmp = y + t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = b * (a - 0.5d0)
if ((x + y) <= (-4d-85)) then
tmp = x + t_1
else
tmp = y + t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = b * (a - 0.5);
double tmp;
if ((x + y) <= -4e-85) {
tmp = x + t_1;
} else {
tmp = y + t_1;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): t_1 = b * (a - 0.5) tmp = 0 if (x + y) <= -4e-85: tmp = x + t_1 else: tmp = y + t_1 return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) t_1 = Float64(b * Float64(a - 0.5)) tmp = 0.0 if (Float64(x + y) <= -4e-85) tmp = Float64(x + t_1); else tmp = Float64(y + t_1); end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
t_1 = b * (a - 0.5);
tmp = 0.0;
if ((x + y) <= -4e-85)
tmp = x + t_1;
else
tmp = y + t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(b * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x + y), $MachinePrecision], -4e-85], N[(x + t$95$1), $MachinePrecision], N[(y + t$95$1), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
t_1 := b \cdot \left(a - 0.5\right)\\
\mathbf{if}\;x + y \leq -4 \cdot 10^{-85}:\\
\;\;\;\;x + t\_1\\
\mathbf{else}:\\
\;\;\;\;y + t\_1\\
\end{array}
\end{array}
if (+.f64 x y) < -3.9999999999999999e-85Initial program 99.9%
Taylor expanded in x around inf
Simplified57.9%
if -3.9999999999999999e-85 < (+.f64 x y) Initial program 99.9%
Taylor expanded in y around inf
Simplified61.2%
Final simplification59.7%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (if (<= a -7.4e+110) (* b a) (if (<= a 1.36e+23) (+ x y) (* b a))))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (a <= -7.4e+110) {
tmp = b * a;
} else if (a <= 1.36e+23) {
tmp = x + y;
} else {
tmp = b * a;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (a <= (-7.4d+110)) then
tmp = b * a
else if (a <= 1.36d+23) then
tmp = x + y
else
tmp = b * a
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (a <= -7.4e+110) {
tmp = b * a;
} else if (a <= 1.36e+23) {
tmp = x + y;
} else {
tmp = b * a;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): tmp = 0 if a <= -7.4e+110: tmp = b * a elif a <= 1.36e+23: tmp = x + y else: tmp = b * a return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) tmp = 0.0 if (a <= -7.4e+110) tmp = Float64(b * a); elseif (a <= 1.36e+23) tmp = Float64(x + y); else tmp = Float64(b * a); end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
tmp = 0.0;
if (a <= -7.4e+110)
tmp = b * a;
elseif (a <= 1.36e+23)
tmp = x + y;
else
tmp = b * a;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := If[LessEqual[a, -7.4e+110], N[(b * a), $MachinePrecision], If[LessEqual[a, 1.36e+23], N[(x + y), $MachinePrecision], N[(b * a), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
\mathbf{if}\;a \leq -7.4 \cdot 10^{+110}:\\
\;\;\;\;b \cdot a\\
\mathbf{elif}\;a \leq 1.36 \cdot 10^{+23}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;b \cdot a\\
\end{array}
\end{array}
if a < -7.40000000000000024e110 or 1.36e23 < a Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
Taylor expanded in a around inf
*-commutativeN/A
*-lowering-*.f6463.1%
Simplified63.1%
if -7.40000000000000024e110 < a < 1.36e23Initial program 99.9%
Taylor expanded in z around 0
+-lowering-+.f6474.6%
Simplified74.6%
Taylor expanded in b around 0
+-commutativeN/A
+-lowering-+.f6450.6%
Simplified50.6%
Final simplification55.1%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (if (<= y 1.72e-279) x (if (<= y 3.2e+56) (* b a) y)))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 1.72e-279) {
tmp = x;
} else if (y <= 3.2e+56) {
tmp = b * a;
} else {
tmp = y;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= 1.72d-279) then
tmp = x
else if (y <= 3.2d+56) then
tmp = b * a
else
tmp = y
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 1.72e-279) {
tmp = x;
} else if (y <= 3.2e+56) {
tmp = b * a;
} else {
tmp = y;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): tmp = 0 if y <= 1.72e-279: tmp = x elif y <= 3.2e+56: tmp = b * a else: tmp = y return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) tmp = 0.0 if (y <= 1.72e-279) tmp = x; elseif (y <= 3.2e+56) tmp = Float64(b * a); else tmp = y; end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
tmp = 0.0;
if (y <= 1.72e-279)
tmp = x;
elseif (y <= 3.2e+56)
tmp = b * a;
else
tmp = y;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, 1.72e-279], x, If[LessEqual[y, 3.2e+56], N[(b * a), $MachinePrecision], y]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.72 \cdot 10^{-279}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 3.2 \cdot 10^{+56}:\\
\;\;\;\;b \cdot a\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if y < 1.7199999999999999e-279Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around inf
Simplified21.4%
if 1.7199999999999999e-279 < y < 3.20000000000000003e56Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
Taylor expanded in a around inf
*-commutativeN/A
*-lowering-*.f6434.2%
Simplified34.2%
if 3.20000000000000003e56 < y Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf
Simplified49.3%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (if (<= y 1.15e+66) (+ x (* b (- a 0.5))) (+ y (* b a))))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 1.15e+66) {
tmp = x + (b * (a - 0.5));
} else {
tmp = y + (b * a);
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= 1.15d+66) then
tmp = x + (b * (a - 0.5d0))
else
tmp = y + (b * a)
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 1.15e+66) {
tmp = x + (b * (a - 0.5));
} else {
tmp = y + (b * a);
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): tmp = 0 if y <= 1.15e+66: tmp = x + (b * (a - 0.5)) else: tmp = y + (b * a) return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) tmp = 0.0 if (y <= 1.15e+66) tmp = Float64(x + Float64(b * Float64(a - 0.5))); else tmp = Float64(y + Float64(b * a)); end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
tmp = 0.0;
if (y <= 1.15e+66)
tmp = x + (b * (a - 0.5));
else
tmp = y + (b * a);
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, 1.15e+66], N[(x + N[(b * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y + N[(b * a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.15 \cdot 10^{+66}:\\
\;\;\;\;x + b \cdot \left(a - 0.5\right)\\
\mathbf{else}:\\
\;\;\;\;y + b \cdot a\\
\end{array}
\end{array}
if y < 1.15e66Initial program 99.9%
Taylor expanded in x around inf
Simplified63.2%
if 1.15e66 < y Initial program 99.9%
Taylor expanded in y around inf
Simplified73.9%
Taylor expanded in a around inf
*-commutativeN/A
*-lowering-*.f6463.1%
Simplified63.1%
Final simplification63.2%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (+ (+ (* b (+ a -0.5)) x) y))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
return ((b * (a + -0.5)) + x) + y;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = ((b * (a + (-0.5d0))) + x) + y
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
return ((b * (a + -0.5)) + x) + y;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): return ((b * (a + -0.5)) + x) + y
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) return Float64(Float64(Float64(b * Float64(a + -0.5)) + x) + y) end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp = code(x, y, z, t, a, b)
tmp = ((b * (a + -0.5)) + x) + y;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(b * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + y), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\left(b \cdot \left(a + -0.5\right) + x\right) + y
\end{array}
Initial program 99.9%
Taylor expanded in z around 0
+-lowering-+.f6479.8%
Simplified79.8%
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6479.8%
Applied egg-rr79.8%
Final simplification79.8%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (+ (* b (- a 0.5)) (+ x y)))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
return (b * (a - 0.5)) + (x + y);
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = (b * (a - 0.5d0)) + (x + y)
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
return (b * (a - 0.5)) + (x + y);
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): return (b * (a - 0.5)) + (x + y)
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) return Float64(Float64(b * Float64(a - 0.5)) + Float64(x + y)) end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp = code(x, y, z, t, a, b)
tmp = (b * (a - 0.5)) + (x + y);
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := N[(N[(b * N[(a - 0.5), $MachinePrecision]), $MachinePrecision] + N[(x + y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
b \cdot \left(a - 0.5\right) + \left(x + y\right)
\end{array}
Initial program 99.9%
Taylor expanded in z around 0
+-lowering-+.f6479.8%
Simplified79.8%
Final simplification79.8%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 (if (<= y 2.9e+43) x y))
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 2.9e+43) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= 2.9d+43) then
tmp = x
else
tmp = y
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 2.9e+43) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): tmp = 0 if y <= 2.9e+43: tmp = x else: tmp = y return tmp
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) tmp = 0.0 if (y <= 2.9e+43) tmp = x; else tmp = y; end return tmp end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp_2 = code(x, y, z, t, a, b)
tmp = 0.0;
if (y <= 2.9e+43)
tmp = x;
else
tmp = y;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, 2.9e+43], x, y]
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2.9 \cdot 10^{+43}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if y < 2.9000000000000002e43Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around inf
Simplified23.2%
if 2.9000000000000002e43 < y Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf
Simplified47.6%
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. (FPCore (x y z t a b) :precision binary64 x)
assert(x < y && y < z && z < t && t < a && a < b);
double code(double x, double y, double z, double t, double a, double b) {
return x;
}
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x
end function
assert x < y && y < z && z < t && t < a && a < b;
public static double code(double x, double y, double z, double t, double a, double b) {
return x;
}
[x, y, z, t, a, b] = sort([x, y, z, t, a, b]) def code(x, y, z, t, a, b): return x
x, y, z, t, a, b = sort([x, y, z, t, a, b]) function code(x, y, z, t, a, b) return x end
x, y, z, t, a, b = num2cell(sort([x, y, z, t, a, b])){:}
function tmp = code(x, y, z, t, a, b)
tmp = x;
end
NOTE: x, y, z, t, a, and b should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_] := x
\begin{array}{l}
[x, y, z, t, a, b] = \mathsf{sort}([x, y, z, t, a, b])\\
\\
x
\end{array}
Initial program 99.9%
+-commutativeN/A
associate--l+N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
log-lowering-log.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around inf
Simplified20.8%
(FPCore (x y z t a b) :precision binary64 (+ (+ (+ x y) (/ (* (- 1.0 (pow (log t) 2.0)) z) (+ 1.0 (log t)))) (* (- a 0.5) b)))
double code(double x, double y, double z, double t, double a, double b) {
return ((x + y) + (((1.0 - pow(log(t), 2.0)) * z) / (1.0 + log(t)))) + ((a - 0.5) * b);
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = ((x + y) + (((1.0d0 - (log(t) ** 2.0d0)) * z) / (1.0d0 + log(t)))) + ((a - 0.5d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x + y) + (((1.0 - Math.pow(Math.log(t), 2.0)) * z) / (1.0 + Math.log(t)))) + ((a - 0.5) * b);
}
def code(x, y, z, t, a, b): return ((x + y) + (((1.0 - math.pow(math.log(t), 2.0)) * z) / (1.0 + math.log(t)))) + ((a - 0.5) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x + y) + Float64(Float64(Float64(1.0 - (log(t) ^ 2.0)) * z) / Float64(1.0 + log(t)))) + Float64(Float64(a - 0.5) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = ((x + y) + (((1.0 - (log(t) ^ 2.0)) * z) / (1.0 + log(t)))) + ((a - 0.5) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x + y), $MachinePrecision] + N[(N[(N[(1.0 - N[Power[N[Log[t], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision] / N[(1.0 + N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(a - 0.5), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x + y\right) + \frac{\left(1 - {\log t}^{2}\right) \cdot z}{1 + \log t}\right) + \left(a - 0.5\right) \cdot b
\end{array}
herbie shell --seed 2024140
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:logBeta from math-functions-0.1.5.2, A"
:precision binary64
:alt
(! :herbie-platform default (+ (+ (+ x y) (/ (* (- 1 (pow (log t) 2)) z) (+ 1 (log t)))) (* (- a 1/2) b)))
(+ (- (+ (+ x y) z) (* z (log t))) (* (- a 0.5) b)))