
(FPCore (x y z t a b) :precision binary64 (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))
double code(double x, double y, double z, double t, double a, double b) {
return (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y;
}
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 * exp((((y * log(z)) + ((t - 1.0d0) * log(a))) - b))) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.exp((((y * Math.log(z)) + ((t - 1.0) * Math.log(a))) - b))) / y;
}
def code(x, y, z, t, a, b): return (x * math.exp((((y * math.log(z)) + ((t - 1.0) * math.log(a))) - b))) / y
function code(x, y, z, t, a, b) return Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(Float64(t - 1.0) * log(a))) - b))) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[(t - 1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))
double code(double x, double y, double z, double t, double a, double b) {
return (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y;
}
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 * exp((((y * log(z)) + ((t - 1.0d0) * log(a))) - b))) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.exp((((y * Math.log(z)) + ((t - 1.0) * Math.log(a))) - b))) / y;
}
def code(x, y, z, t, a, b): return (x * math.exp((((y * math.log(z)) + ((t - 1.0) * math.log(a))) - b))) / y
function code(x, y, z, t, a, b) return Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(Float64(t - 1.0) * log(a))) - b))) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[(t - 1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}
\end{array}
(FPCore (x y z t a b) :precision binary64 (/ (* x (exp (- (+ (* y (log z)) (* (+ t -1.0) (log a))) b))) y))
double code(double x, double y, double z, double t, double a, double b) {
return (x * exp((((y * log(z)) + ((t + -1.0) * log(a))) - b))) / y;
}
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 * exp((((y * log(z)) + ((t + (-1.0d0)) * log(a))) - b))) / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return (x * Math.exp((((y * Math.log(z)) + ((t + -1.0) * Math.log(a))) - b))) / y;
}
def code(x, y, z, t, a, b): return (x * math.exp((((y * math.log(z)) + ((t + -1.0) * math.log(a))) - b))) / y
function code(x, y, z, t, a, b) return Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) + Float64(Float64(t + -1.0) * log(a))) - b))) / y) end
function tmp = code(x, y, z, t, a, b) tmp = (x * exp((((y * log(z)) + ((t + -1.0) * log(a))) - b))) / y; end
code[x_, y_, z_, t_, a_, b_] := N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] + N[(N[(t + -1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot e^{\left(y \cdot \log z + \left(t + -1\right) \cdot \log a\right) - b}}{y}
\end{array}
Initial program 99.3%
Final simplification99.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (pow z y) y))))
(if (<= y -1.2e+55)
t_1
(if (<= y 8.5e+17) (/ (* (pow a (+ t -1.0)) (/ x (exp b))) y) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (pow(z, y) / y);
double tmp;
if (y <= -1.2e+55) {
tmp = t_1;
} else if (y <= 8.5e+17) {
tmp = (pow(a, (t + -1.0)) * (x / exp(b))) / y;
} else {
tmp = t_1;
}
return tmp;
}
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 = x * ((z ** y) / y)
if (y <= (-1.2d+55)) then
tmp = t_1
else if (y <= 8.5d+17) then
tmp = ((a ** (t + (-1.0d0))) * (x / exp(b))) / y
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (Math.pow(z, y) / y);
double tmp;
if (y <= -1.2e+55) {
tmp = t_1;
} else if (y <= 8.5e+17) {
tmp = (Math.pow(a, (t + -1.0)) * (x / Math.exp(b))) / y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (math.pow(z, y) / y) tmp = 0 if y <= -1.2e+55: tmp = t_1 elif y <= 8.5e+17: tmp = (math.pow(a, (t + -1.0)) * (x / math.exp(b))) / y else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64((z ^ y) / y)) tmp = 0.0 if (y <= -1.2e+55) tmp = t_1; elseif (y <= 8.5e+17) tmp = Float64(Float64((a ^ Float64(t + -1.0)) * Float64(x / exp(b))) / y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * ((z ^ y) / y); tmp = 0.0; if (y <= -1.2e+55) tmp = t_1; elseif (y <= 8.5e+17) tmp = ((a ^ (t + -1.0)) * (x / exp(b))) / y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.2e+55], t$95$1, If[LessEqual[y, 8.5e+17], N[(N[(N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision] * N[(x / N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{{z}^{y}}{y}\\
\mathbf{if}\;y \leq -1.2 \cdot 10^{+55}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 8.5 \cdot 10^{+17}:\\
\;\;\;\;\frac{{a}^{\left(t + -1\right)} \cdot \frac{x}{e^{b}}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.2e55 or 8.5e17 < y Initial program 100.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
log-lowering-log.f6482.4%
Simplified82.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
pow-to-expN/A
/-lowering-/.f64N/A
pow-lowering-pow.f6482.4%
Applied egg-rr82.4%
if -1.2e55 < y < 8.5e17Initial program 98.7%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6479.9%
Simplified79.9%
Final simplification81.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (pow a t) y))))
(if (<= t -1.45e+56)
t_1
(if (<= t 0.00055) (/ (* x (pow z y)) (* a (* y (exp b)))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (pow(a, t) / y);
double tmp;
if (t <= -1.45e+56) {
tmp = t_1;
} else if (t <= 0.00055) {
tmp = (x * pow(z, y)) / (a * (y * exp(b)));
} else {
tmp = t_1;
}
return tmp;
}
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 = x * ((a ** t) / y)
if (t <= (-1.45d+56)) then
tmp = t_1
else if (t <= 0.00055d0) then
tmp = (x * (z ** y)) / (a * (y * exp(b)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (Math.pow(a, t) / y);
double tmp;
if (t <= -1.45e+56) {
tmp = t_1;
} else if (t <= 0.00055) {
tmp = (x * Math.pow(z, y)) / (a * (y * Math.exp(b)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (math.pow(a, t) / y) tmp = 0 if t <= -1.45e+56: tmp = t_1 elif t <= 0.00055: tmp = (x * math.pow(z, y)) / (a * (y * math.exp(b))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64((a ^ t) / y)) tmp = 0.0 if (t <= -1.45e+56) tmp = t_1; elseif (t <= 0.00055) tmp = Float64(Float64(x * (z ^ y)) / Float64(a * Float64(y * exp(b)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * ((a ^ t) / y); tmp = 0.0; if (t <= -1.45e+56) tmp = t_1; elseif (t <= 0.00055) tmp = (x * (z ^ y)) / (a * (y * exp(b))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[Power[a, t], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.45e+56], t$95$1, If[LessEqual[t, 0.00055], N[(N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision] / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{{a}^{t}}{y}\\
\mathbf{if}\;t \leq -1.45 \cdot 10^{+56}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 0.00055:\\
\;\;\;\;\frac{x \cdot {z}^{y}}{a \cdot \left(y \cdot e^{b}\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.45000000000000004e56 or 5.50000000000000033e-4 < t Initial program 100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6483.7%
Simplified83.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
exp-to-powN/A
pow-lowering-pow.f6483.7%
Applied egg-rr83.7%
if -1.45000000000000004e56 < t < 5.50000000000000033e-4Initial program 98.5%
associate-*l/N/A
exp-diffN/A
associate-*r/N/A
associate-*l/N/A
associate-/r/N/A
exp-diffN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
exp-diffN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
exp-sumN/A
*-lowering-*.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
Simplified68.4%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6475.4%
Simplified75.4%
Final simplification79.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (* b (+ 0.5 (* b -0.16666666666666666)))))
(t_2 (+ x (* b (- x t_1))))
(t_3 (/ (/ x (exp b)) y)))
(if (<= b -1.55e-45)
t_3
(if (<= b -1.52e-251)
(/ x (* y a))
(if (<= b 7.3e-211)
(*
b
(*
b
(*
b
(+
(/ (+ -0.16666666666666666 (/ 0.5 b)) (/ y x))
(/ (- (/ x b) x) (* y (* b b)))))))
(if (<= b 1.75)
(/ (* (/ 1.0 y) (* (+ x (* b (- t_1 x))) t_2)) t_2)
t_3))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (b * (0.5 + (b * -0.16666666666666666)));
double t_2 = x + (b * (x - t_1));
double t_3 = (x / exp(b)) / y;
double tmp;
if (b <= -1.55e-45) {
tmp = t_3;
} else if (b <= -1.52e-251) {
tmp = x / (y * a);
} else if (b <= 7.3e-211) {
tmp = b * (b * (b * (((-0.16666666666666666 + (0.5 / b)) / (y / x)) + (((x / b) - x) / (y * (b * b))))));
} else if (b <= 1.75) {
tmp = ((1.0 / y) * ((x + (b * (t_1 - x))) * t_2)) / t_2;
} else {
tmp = t_3;
}
return tmp;
}
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) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = x * (b * (0.5d0 + (b * (-0.16666666666666666d0))))
t_2 = x + (b * (x - t_1))
t_3 = (x / exp(b)) / y
if (b <= (-1.55d-45)) then
tmp = t_3
else if (b <= (-1.52d-251)) then
tmp = x / (y * a)
else if (b <= 7.3d-211) then
tmp = b * (b * (b * ((((-0.16666666666666666d0) + (0.5d0 / b)) / (y / x)) + (((x / b) - x) / (y * (b * b))))))
else if (b <= 1.75d0) then
tmp = ((1.0d0 / y) * ((x + (b * (t_1 - x))) * t_2)) / t_2
else
tmp = t_3
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (b * (0.5 + (b * -0.16666666666666666)));
double t_2 = x + (b * (x - t_1));
double t_3 = (x / Math.exp(b)) / y;
double tmp;
if (b <= -1.55e-45) {
tmp = t_3;
} else if (b <= -1.52e-251) {
tmp = x / (y * a);
} else if (b <= 7.3e-211) {
tmp = b * (b * (b * (((-0.16666666666666666 + (0.5 / b)) / (y / x)) + (((x / b) - x) / (y * (b * b))))));
} else if (b <= 1.75) {
tmp = ((1.0 / y) * ((x + (b * (t_1 - x))) * t_2)) / t_2;
} else {
tmp = t_3;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (b * (0.5 + (b * -0.16666666666666666))) t_2 = x + (b * (x - t_1)) t_3 = (x / math.exp(b)) / y tmp = 0 if b <= -1.55e-45: tmp = t_3 elif b <= -1.52e-251: tmp = x / (y * a) elif b <= 7.3e-211: tmp = b * (b * (b * (((-0.16666666666666666 + (0.5 / b)) / (y / x)) + (((x / b) - x) / (y * (b * b)))))) elif b <= 1.75: tmp = ((1.0 / y) * ((x + (b * (t_1 - x))) * t_2)) / t_2 else: tmp = t_3 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(b * Float64(0.5 + Float64(b * -0.16666666666666666)))) t_2 = Float64(x + Float64(b * Float64(x - t_1))) t_3 = Float64(Float64(x / exp(b)) / y) tmp = 0.0 if (b <= -1.55e-45) tmp = t_3; elseif (b <= -1.52e-251) tmp = Float64(x / Float64(y * a)); elseif (b <= 7.3e-211) tmp = Float64(b * Float64(b * Float64(b * Float64(Float64(Float64(-0.16666666666666666 + Float64(0.5 / b)) / Float64(y / x)) + Float64(Float64(Float64(x / b) - x) / Float64(y * Float64(b * b))))))); elseif (b <= 1.75) tmp = Float64(Float64(Float64(1.0 / y) * Float64(Float64(x + Float64(b * Float64(t_1 - x))) * t_2)) / t_2); else tmp = t_3; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * (b * (0.5 + (b * -0.16666666666666666))); t_2 = x + (b * (x - t_1)); t_3 = (x / exp(b)) / y; tmp = 0.0; if (b <= -1.55e-45) tmp = t_3; elseif (b <= -1.52e-251) tmp = x / (y * a); elseif (b <= 7.3e-211) tmp = b * (b * (b * (((-0.16666666666666666 + (0.5 / b)) / (y / x)) + (((x / b) - x) / (y * (b * b)))))); elseif (b <= 1.75) tmp = ((1.0 / y) * ((x + (b * (t_1 - x))) * t_2)) / t_2; else tmp = t_3; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(b * N[(0.5 + N[(b * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(b * N[(x - t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(x / N[Exp[b], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[b, -1.55e-45], t$95$3, If[LessEqual[b, -1.52e-251], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 7.3e-211], N[(b * N[(b * N[(b * N[(N[(N[(-0.16666666666666666 + N[(0.5 / b), $MachinePrecision]), $MachinePrecision] / N[(y / x), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(x / b), $MachinePrecision] - x), $MachinePrecision] / N[(y * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.75], N[(N[(N[(1.0 / y), $MachinePrecision] * N[(N[(x + N[(b * N[(t$95$1 - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision], t$95$3]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(b \cdot \left(0.5 + b \cdot -0.16666666666666666\right)\right)\\
t_2 := x + b \cdot \left(x - t\_1\right)\\
t_3 := \frac{\frac{x}{e^{b}}}{y}\\
\mathbf{if}\;b \leq -1.55 \cdot 10^{-45}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;b \leq -1.52 \cdot 10^{-251}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{elif}\;b \leq 7.3 \cdot 10^{-211}:\\
\;\;\;\;b \cdot \left(b \cdot \left(b \cdot \left(\frac{-0.16666666666666666 + \frac{0.5}{b}}{\frac{y}{x}} + \frac{\frac{x}{b} - x}{y \cdot \left(b \cdot b\right)}\right)\right)\right)\\
\mathbf{elif}\;b \leq 1.75:\\
\;\;\;\;\frac{\frac{1}{y} \cdot \left(\left(x + b \cdot \left(t\_1 - x\right)\right) \cdot t\_2\right)}{t\_2}\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if b < -1.55e-45 or 1.75 < b Initial program 99.9%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6472.2%
Simplified72.2%
sub0-negN/A
rec-expN/A
div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6472.2%
Applied egg-rr72.2%
if -1.55e-45 < b < -1.52e-251Initial program 98.2%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6469.8%
Simplified69.8%
Taylor expanded in b around 0
Simplified69.8%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f6448.3%
Simplified48.3%
if -1.52e-251 < b < 7.29999999999999968e-211Initial program 98.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6410.3%
Simplified10.3%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6410.3%
Simplified10.3%
Taylor expanded in b around inf
cube-multN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
Simplified0.0%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
Applied egg-rr54.5%
if 7.29999999999999968e-211 < b < 1.75Initial program 98.5%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6425.2%
Simplified25.2%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6425.2%
Simplified25.2%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6425.2%
Applied egg-rr25.2%
associate-/r/N/A
flip-+N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr51.6%
Final simplification62.7%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -1.65e+56)
(* x (/ (pow a t) y))
(if (<= t 1.65e-159)
(/ (/ x (* a (exp b))) y)
(if (<= t 5.9e-21)
(* x (/ (pow z y) y))
(/ (* x (pow a (+ t -1.0))) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.65e+56) {
tmp = x * (pow(a, t) / y);
} else if (t <= 1.65e-159) {
tmp = (x / (a * exp(b))) / y;
} else if (t <= 5.9e-21) {
tmp = x * (pow(z, y) / y);
} else {
tmp = (x * pow(a, (t + -1.0))) / y;
}
return tmp;
}
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 (t <= (-1.65d+56)) then
tmp = x * ((a ** t) / y)
else if (t <= 1.65d-159) then
tmp = (x / (a * exp(b))) / y
else if (t <= 5.9d-21) then
tmp = x * ((z ** y) / y)
else
tmp = (x * (a ** (t + (-1.0d0)))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.65e+56) {
tmp = x * (Math.pow(a, t) / y);
} else if (t <= 1.65e-159) {
tmp = (x / (a * Math.exp(b))) / y;
} else if (t <= 5.9e-21) {
tmp = x * (Math.pow(z, y) / y);
} else {
tmp = (x * Math.pow(a, (t + -1.0))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.65e+56: tmp = x * (math.pow(a, t) / y) elif t <= 1.65e-159: tmp = (x / (a * math.exp(b))) / y elif t <= 5.9e-21: tmp = x * (math.pow(z, y) / y) else: tmp = (x * math.pow(a, (t + -1.0))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.65e+56) tmp = Float64(x * Float64((a ^ t) / y)); elseif (t <= 1.65e-159) tmp = Float64(Float64(x / Float64(a * exp(b))) / y); elseif (t <= 5.9e-21) tmp = Float64(x * Float64((z ^ y) / y)); else tmp = Float64(Float64(x * (a ^ Float64(t + -1.0))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.65e+56) tmp = x * ((a ^ t) / y); elseif (t <= 1.65e-159) tmp = (x / (a * exp(b))) / y; elseif (t <= 5.9e-21) tmp = x * ((z ^ y) / y); else tmp = (x * (a ^ (t + -1.0))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.65e+56], N[(x * N[(N[Power[a, t], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.65e-159], N[(N[(x / N[(a * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[t, 5.9e-21], N[(x * N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.65 \cdot 10^{+56}:\\
\;\;\;\;x \cdot \frac{{a}^{t}}{y}\\
\mathbf{elif}\;t \leq 1.65 \cdot 10^{-159}:\\
\;\;\;\;\frac{\frac{x}{a \cdot e^{b}}}{y}\\
\mathbf{elif}\;t \leq 5.9 \cdot 10^{-21}:\\
\;\;\;\;x \cdot \frac{{z}^{y}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot {a}^{\left(t + -1\right)}}{y}\\
\end{array}
\end{array}
if t < -1.65000000000000001e56Initial program 100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6487.7%
Simplified87.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
exp-to-powN/A
pow-lowering-pow.f6487.7%
Applied egg-rr87.7%
if -1.65000000000000001e56 < t < 1.6500000000000001e-159Initial program 98.5%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6474.7%
Simplified74.7%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6475.9%
Simplified75.9%
if 1.6500000000000001e-159 < t < 5.9000000000000003e-21Initial program 98.7%
Taylor expanded in y around inf
*-lowering-*.f64N/A
log-lowering-log.f6473.8%
Simplified73.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
pow-to-expN/A
/-lowering-/.f64N/A
pow-lowering-pow.f6473.8%
Applied egg-rr73.8%
if 5.9000000000000003e-21 < t Initial program 99.9%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.3%
Simplified66.3%
Taylor expanded in b around 0
Simplified78.8%
Final simplification79.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (pow a t) y))))
(if (<= t -1.75e+56)
t_1
(if (<= t 2.4e-153)
(/ (/ x (* a (exp b))) y)
(if (<= t 0.00055) (* x (/ (pow z y) y)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (pow(a, t) / y);
double tmp;
if (t <= -1.75e+56) {
tmp = t_1;
} else if (t <= 2.4e-153) {
tmp = (x / (a * exp(b))) / y;
} else if (t <= 0.00055) {
tmp = x * (pow(z, y) / y);
} else {
tmp = t_1;
}
return tmp;
}
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 = x * ((a ** t) / y)
if (t <= (-1.75d+56)) then
tmp = t_1
else if (t <= 2.4d-153) then
tmp = (x / (a * exp(b))) / y
else if (t <= 0.00055d0) then
tmp = x * ((z ** y) / y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (Math.pow(a, t) / y);
double tmp;
if (t <= -1.75e+56) {
tmp = t_1;
} else if (t <= 2.4e-153) {
tmp = (x / (a * Math.exp(b))) / y;
} else if (t <= 0.00055) {
tmp = x * (Math.pow(z, y) / y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (math.pow(a, t) / y) tmp = 0 if t <= -1.75e+56: tmp = t_1 elif t <= 2.4e-153: tmp = (x / (a * math.exp(b))) / y elif t <= 0.00055: tmp = x * (math.pow(z, y) / y) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64((a ^ t) / y)) tmp = 0.0 if (t <= -1.75e+56) tmp = t_1; elseif (t <= 2.4e-153) tmp = Float64(Float64(x / Float64(a * exp(b))) / y); elseif (t <= 0.00055) tmp = Float64(x * Float64((z ^ y) / y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * ((a ^ t) / y); tmp = 0.0; if (t <= -1.75e+56) tmp = t_1; elseif (t <= 2.4e-153) tmp = (x / (a * exp(b))) / y; elseif (t <= 0.00055) tmp = x * ((z ^ y) / y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[Power[a, t], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.75e+56], t$95$1, If[LessEqual[t, 2.4e-153], N[(N[(x / N[(a * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[t, 0.00055], N[(x * N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{{a}^{t}}{y}\\
\mathbf{if}\;t \leq -1.75 \cdot 10^{+56}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-153}:\\
\;\;\;\;\frac{\frac{x}{a \cdot e^{b}}}{y}\\
\mathbf{elif}\;t \leq 0.00055:\\
\;\;\;\;x \cdot \frac{{z}^{y}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.75e56 or 5.50000000000000033e-4 < t Initial program 100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6483.7%
Simplified83.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
exp-to-powN/A
pow-lowering-pow.f6483.7%
Applied egg-rr83.7%
if -1.75e56 < t < 2.4000000000000002e-153Initial program 98.5%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6474.7%
Simplified74.7%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6475.9%
Simplified75.9%
if 2.4000000000000002e-153 < t < 5.50000000000000033e-4Initial program 98.6%
Taylor expanded in y around inf
*-lowering-*.f64N/A
log-lowering-log.f6468.2%
Simplified68.2%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
pow-to-expN/A
/-lowering-/.f64N/A
pow-lowering-pow.f6468.2%
Applied egg-rr68.2%
Final simplification78.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (pow a t) y))))
(if (<= t -1.35e+56)
t_1
(if (<= t 4.1e-154)
(/ x (* a (* y (exp b))))
(if (<= t 0.00055) (* x (/ (pow z y) y)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (pow(a, t) / y);
double tmp;
if (t <= -1.35e+56) {
tmp = t_1;
} else if (t <= 4.1e-154) {
tmp = x / (a * (y * exp(b)));
} else if (t <= 0.00055) {
tmp = x * (pow(z, y) / y);
} else {
tmp = t_1;
}
return tmp;
}
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 = x * ((a ** t) / y)
if (t <= (-1.35d+56)) then
tmp = t_1
else if (t <= 4.1d-154) then
tmp = x / (a * (y * exp(b)))
else if (t <= 0.00055d0) then
tmp = x * ((z ** y) / y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (Math.pow(a, t) / y);
double tmp;
if (t <= -1.35e+56) {
tmp = t_1;
} else if (t <= 4.1e-154) {
tmp = x / (a * (y * Math.exp(b)));
} else if (t <= 0.00055) {
tmp = x * (Math.pow(z, y) / y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (math.pow(a, t) / y) tmp = 0 if t <= -1.35e+56: tmp = t_1 elif t <= 4.1e-154: tmp = x / (a * (y * math.exp(b))) elif t <= 0.00055: tmp = x * (math.pow(z, y) / y) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64((a ^ t) / y)) tmp = 0.0 if (t <= -1.35e+56) tmp = t_1; elseif (t <= 4.1e-154) tmp = Float64(x / Float64(a * Float64(y * exp(b)))); elseif (t <= 0.00055) tmp = Float64(x * Float64((z ^ y) / y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * ((a ^ t) / y); tmp = 0.0; if (t <= -1.35e+56) tmp = t_1; elseif (t <= 4.1e-154) tmp = x / (a * (y * exp(b))); elseif (t <= 0.00055) tmp = x * ((z ^ y) / y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[Power[a, t], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.35e+56], t$95$1, If[LessEqual[t, 4.1e-154], N[(x / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.00055], N[(x * N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{{a}^{t}}{y}\\
\mathbf{if}\;t \leq -1.35 \cdot 10^{+56}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 4.1 \cdot 10^{-154}:\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot e^{b}\right)}\\
\mathbf{elif}\;t \leq 0.00055:\\
\;\;\;\;x \cdot \frac{{z}^{y}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.35000000000000005e56 or 5.50000000000000033e-4 < t Initial program 100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6483.7%
Simplified83.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
exp-to-powN/A
pow-lowering-pow.f6483.7%
Applied egg-rr83.7%
if -1.35000000000000005e56 < t < 4.1e-154Initial program 98.5%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6474.7%
Simplified74.7%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6475.0%
Simplified75.0%
if 4.1e-154 < t < 5.50000000000000033e-4Initial program 98.6%
Taylor expanded in y around inf
*-lowering-*.f64N/A
log-lowering-log.f6468.2%
Simplified68.2%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
pow-to-expN/A
/-lowering-/.f64N/A
pow-lowering-pow.f6468.2%
Applied egg-rr68.2%
Final simplification78.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (pow a t) y))))
(if (<= t -1.3e+56)
t_1
(if (<= t 2.7e-155)
(/ (/ x (exp b)) y)
(if (<= t 0.00055) (* x (/ (pow z y) y)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (pow(a, t) / y);
double tmp;
if (t <= -1.3e+56) {
tmp = t_1;
} else if (t <= 2.7e-155) {
tmp = (x / exp(b)) / y;
} else if (t <= 0.00055) {
tmp = x * (pow(z, y) / y);
} else {
tmp = t_1;
}
return tmp;
}
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 = x * ((a ** t) / y)
if (t <= (-1.3d+56)) then
tmp = t_1
else if (t <= 2.7d-155) then
tmp = (x / exp(b)) / y
else if (t <= 0.00055d0) then
tmp = x * ((z ** y) / y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (Math.pow(a, t) / y);
double tmp;
if (t <= -1.3e+56) {
tmp = t_1;
} else if (t <= 2.7e-155) {
tmp = (x / Math.exp(b)) / y;
} else if (t <= 0.00055) {
tmp = x * (Math.pow(z, y) / y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (math.pow(a, t) / y) tmp = 0 if t <= -1.3e+56: tmp = t_1 elif t <= 2.7e-155: tmp = (x / math.exp(b)) / y elif t <= 0.00055: tmp = x * (math.pow(z, y) / y) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64((a ^ t) / y)) tmp = 0.0 if (t <= -1.3e+56) tmp = t_1; elseif (t <= 2.7e-155) tmp = Float64(Float64(x / exp(b)) / y); elseif (t <= 0.00055) tmp = Float64(x * Float64((z ^ y) / y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * ((a ^ t) / y); tmp = 0.0; if (t <= -1.3e+56) tmp = t_1; elseif (t <= 2.7e-155) tmp = (x / exp(b)) / y; elseif (t <= 0.00055) tmp = x * ((z ^ y) / y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[Power[a, t], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.3e+56], t$95$1, If[LessEqual[t, 2.7e-155], N[(N[(x / N[Exp[b], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[t, 0.00055], N[(x * N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{{a}^{t}}{y}\\
\mathbf{if}\;t \leq -1.3 \cdot 10^{+56}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.7 \cdot 10^{-155}:\\
\;\;\;\;\frac{\frac{x}{e^{b}}}{y}\\
\mathbf{elif}\;t \leq 0.00055:\\
\;\;\;\;x \cdot \frac{{z}^{y}}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.30000000000000005e56 or 5.50000000000000033e-4 < t Initial program 100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6483.7%
Simplified83.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
exp-to-powN/A
pow-lowering-pow.f6483.7%
Applied egg-rr83.7%
if -1.30000000000000005e56 < t < 2.69999999999999981e-155Initial program 98.5%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6457.8%
Simplified57.8%
sub0-negN/A
rec-expN/A
div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6457.8%
Applied egg-rr57.8%
if 2.69999999999999981e-155 < t < 5.50000000000000033e-4Initial program 98.6%
Taylor expanded in y around inf
*-lowering-*.f64N/A
log-lowering-log.f6468.2%
Simplified68.2%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
pow-to-expN/A
/-lowering-/.f64N/A
pow-lowering-pow.f6468.2%
Applied egg-rr68.2%
Final simplification72.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (/ (pow a t) y))))
(if (<= t -1.75e+56)
t_1
(if (<= t 8.2e-120)
(/ (/ x (exp b)) y)
(if (<= t 3.4e-14) (/ 1.0 (/ (/ (* y -6.0) x) (* b (* b b)))) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (pow(a, t) / y);
double tmp;
if (t <= -1.75e+56) {
tmp = t_1;
} else if (t <= 8.2e-120) {
tmp = (x / exp(b)) / y;
} else if (t <= 3.4e-14) {
tmp = 1.0 / (((y * -6.0) / x) / (b * (b * b)));
} else {
tmp = t_1;
}
return tmp;
}
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 = x * ((a ** t) / y)
if (t <= (-1.75d+56)) then
tmp = t_1
else if (t <= 8.2d-120) then
tmp = (x / exp(b)) / y
else if (t <= 3.4d-14) then
tmp = 1.0d0 / (((y * (-6.0d0)) / x) / (b * (b * b)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (Math.pow(a, t) / y);
double tmp;
if (t <= -1.75e+56) {
tmp = t_1;
} else if (t <= 8.2e-120) {
tmp = (x / Math.exp(b)) / y;
} else if (t <= 3.4e-14) {
tmp = 1.0 / (((y * -6.0) / x) / (b * (b * b)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (math.pow(a, t) / y) tmp = 0 if t <= -1.75e+56: tmp = t_1 elif t <= 8.2e-120: tmp = (x / math.exp(b)) / y elif t <= 3.4e-14: tmp = 1.0 / (((y * -6.0) / x) / (b * (b * b))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64((a ^ t) / y)) tmp = 0.0 if (t <= -1.75e+56) tmp = t_1; elseif (t <= 8.2e-120) tmp = Float64(Float64(x / exp(b)) / y); elseif (t <= 3.4e-14) tmp = Float64(1.0 / Float64(Float64(Float64(y * -6.0) / x) / Float64(b * Float64(b * b)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * ((a ^ t) / y); tmp = 0.0; if (t <= -1.75e+56) tmp = t_1; elseif (t <= 8.2e-120) tmp = (x / exp(b)) / y; elseif (t <= 3.4e-14) tmp = 1.0 / (((y * -6.0) / x) / (b * (b * b))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[Power[a, t], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.75e+56], t$95$1, If[LessEqual[t, 8.2e-120], N[(N[(x / N[Exp[b], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[t, 3.4e-14], N[(1.0 / N[(N[(N[(y * -6.0), $MachinePrecision] / x), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \frac{{a}^{t}}{y}\\
\mathbf{if}\;t \leq -1.75 \cdot 10^{+56}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 8.2 \cdot 10^{-120}:\\
\;\;\;\;\frac{\frac{x}{e^{b}}}{y}\\
\mathbf{elif}\;t \leq 3.4 \cdot 10^{-14}:\\
\;\;\;\;\frac{1}{\frac{\frac{y \cdot -6}{x}}{b \cdot \left(b \cdot b\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.75e56 or 3.40000000000000003e-14 < t Initial program 100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6483.1%
Simplified83.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
exp-to-powN/A
pow-lowering-pow.f6483.1%
Applied egg-rr83.1%
if -1.75e56 < t < 8.20000000000000068e-120Initial program 98.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6455.5%
Simplified55.5%
sub0-negN/A
rec-expN/A
div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6455.5%
Applied egg-rr55.5%
if 8.20000000000000068e-120 < t < 3.40000000000000003e-14Initial program 98.3%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6436.2%
Simplified36.2%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6426.8%
Simplified26.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6428.6%
Applied egg-rr28.6%
Taylor expanded in b around inf
associate-*r/N/A
*-commutativeN/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.0%
Simplified53.0%
Final simplification69.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (* (* b b) 0.5)) y)))
(if (<= b -4.6e+92)
t_1
(if (<= b -3.8e-231)
(/ x (* y a))
(if (<= b 1.4e-247)
t_1
(if (<= b 8e+135) (/ (* x (/ 1.0 a)) y) (* b (/ x (* y b)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * ((b * b) * 0.5)) / y;
double tmp;
if (b <= -4.6e+92) {
tmp = t_1;
} else if (b <= -3.8e-231) {
tmp = x / (y * a);
} else if (b <= 1.4e-247) {
tmp = t_1;
} else if (b <= 8e+135) {
tmp = (x * (1.0 / a)) / y;
} else {
tmp = b * (x / (y * b));
}
return tmp;
}
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 = (x * ((b * b) * 0.5d0)) / y
if (b <= (-4.6d+92)) then
tmp = t_1
else if (b <= (-3.8d-231)) then
tmp = x / (y * a)
else if (b <= 1.4d-247) then
tmp = t_1
else if (b <= 8d+135) then
tmp = (x * (1.0d0 / a)) / y
else
tmp = b * (x / (y * b))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * ((b * b) * 0.5)) / y;
double tmp;
if (b <= -4.6e+92) {
tmp = t_1;
} else if (b <= -3.8e-231) {
tmp = x / (y * a);
} else if (b <= 1.4e-247) {
tmp = t_1;
} else if (b <= 8e+135) {
tmp = (x * (1.0 / a)) / y;
} else {
tmp = b * (x / (y * b));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * ((b * b) * 0.5)) / y tmp = 0 if b <= -4.6e+92: tmp = t_1 elif b <= -3.8e-231: tmp = x / (y * a) elif b <= 1.4e-247: tmp = t_1 elif b <= 8e+135: tmp = (x * (1.0 / a)) / y else: tmp = b * (x / (y * b)) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * Float64(Float64(b * b) * 0.5)) / y) tmp = 0.0 if (b <= -4.6e+92) tmp = t_1; elseif (b <= -3.8e-231) tmp = Float64(x / Float64(y * a)); elseif (b <= 1.4e-247) tmp = t_1; elseif (b <= 8e+135) tmp = Float64(Float64(x * Float64(1.0 / a)) / y); else tmp = Float64(b * Float64(x / Float64(y * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * ((b * b) * 0.5)) / y; tmp = 0.0; if (b <= -4.6e+92) tmp = t_1; elseif (b <= -3.8e-231) tmp = x / (y * a); elseif (b <= 1.4e-247) tmp = t_1; elseif (b <= 8e+135) tmp = (x * (1.0 / a)) / y; else tmp = b * (x / (y * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[b, -4.6e+92], t$95$1, If[LessEqual[b, -3.8e-231], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.4e-247], t$95$1, If[LessEqual[b, 8e+135], N[(N[(x * N[(1.0 / a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(b * N[(x / N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot \left(\left(b \cdot b\right) \cdot 0.5\right)}{y}\\
\mathbf{if}\;b \leq -4.6 \cdot 10^{+92}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -3.8 \cdot 10^{-231}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{elif}\;b \leq 1.4 \cdot 10^{-247}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 8 \cdot 10^{+135}:\\
\;\;\;\;\frac{x \cdot \frac{1}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;b \cdot \frac{x}{y \cdot b}\\
\end{array}
\end{array}
if b < -4.59999999999999997e92 or -3.80000000000000013e-231 < b < 1.39999999999999993e-247Initial program 99.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6457.1%
Simplified57.1%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6450.9%
Simplified50.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.2%
Simplified63.2%
if -4.59999999999999997e92 < b < -3.80000000000000013e-231Initial program 98.8%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6464.1%
Simplified64.1%
Taylor expanded in b around 0
Simplified68.7%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f6438.8%
Simplified38.8%
if 1.39999999999999993e-247 < b < 7.99999999999999969e135Initial program 98.6%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6467.3%
Simplified67.3%
Taylor expanded in b around 0
Simplified60.3%
Taylor expanded in t around 0
/-lowering-/.f6435.5%
Simplified35.5%
if 7.99999999999999969e135 < b Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6485.9%
Simplified85.9%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f641.3%
Simplified1.3%
Taylor expanded in b around inf
cube-multN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
Simplified1.3%
Taylor expanded in b around 0
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6442.0%
Simplified42.0%
Final simplification45.8%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -1.95e-31)
(* (* b (* (* b b) -0.16666666666666666)) (/ x y))
(if (<= b -4.2e-226)
(/ x (* y a))
(if (<= b 1.8e-247)
(/ (* x (* (* b b) 0.5)) y)
(/ (/ x (* a (+ 1.0 b))) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.95e-31) {
tmp = (b * ((b * b) * -0.16666666666666666)) * (x / y);
} else if (b <= -4.2e-226) {
tmp = x / (y * a);
} else if (b <= 1.8e-247) {
tmp = (x * ((b * b) * 0.5)) / y;
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
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 (b <= (-1.95d-31)) then
tmp = (b * ((b * b) * (-0.16666666666666666d0))) * (x / y)
else if (b <= (-4.2d-226)) then
tmp = x / (y * a)
else if (b <= 1.8d-247) then
tmp = (x * ((b * b) * 0.5d0)) / y
else
tmp = (x / (a * (1.0d0 + b))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.95e-31) {
tmp = (b * ((b * b) * -0.16666666666666666)) * (x / y);
} else if (b <= -4.2e-226) {
tmp = x / (y * a);
} else if (b <= 1.8e-247) {
tmp = (x * ((b * b) * 0.5)) / y;
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1.95e-31: tmp = (b * ((b * b) * -0.16666666666666666)) * (x / y) elif b <= -4.2e-226: tmp = x / (y * a) elif b <= 1.8e-247: tmp = (x * ((b * b) * 0.5)) / y else: tmp = (x / (a * (1.0 + b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.95e-31) tmp = Float64(Float64(b * Float64(Float64(b * b) * -0.16666666666666666)) * Float64(x / y)); elseif (b <= -4.2e-226) tmp = Float64(x / Float64(y * a)); elseif (b <= 1.8e-247) tmp = Float64(Float64(x * Float64(Float64(b * b) * 0.5)) / y); else tmp = Float64(Float64(x / Float64(a * Float64(1.0 + b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -1.95e-31) tmp = (b * ((b * b) * -0.16666666666666666)) * (x / y); elseif (b <= -4.2e-226) tmp = x / (y * a); elseif (b <= 1.8e-247) tmp = (x * ((b * b) * 0.5)) / y; else tmp = (x / (a * (1.0 + b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.95e-31], N[(N[(b * N[(N[(b * b), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -4.2e-226], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.8e-247], N[(N[(x * N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a * N[(1.0 + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.95 \cdot 10^{-31}:\\
\;\;\;\;\left(b \cdot \left(\left(b \cdot b\right) \cdot -0.16666666666666666\right)\right) \cdot \frac{x}{y}\\
\mathbf{elif}\;b \leq -4.2 \cdot 10^{-226}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{elif}\;b \leq 1.8 \cdot 10^{-247}:\\
\;\;\;\;\frac{x \cdot \left(\left(b \cdot b\right) \cdot 0.5\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a \cdot \left(1 + b\right)}}{y}\\
\end{array}
\end{array}
if b < -1.9500000000000001e-31Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6466.5%
Simplified66.5%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6454.1%
Simplified54.1%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.3%
Simplified56.3%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6453.8%
Applied egg-rr53.8%
if -1.9500000000000001e-31 < b < -4.2000000000000003e-226Initial program 98.0%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6474.5%
Simplified74.5%
Taylor expanded in b around 0
Simplified74.5%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f6449.2%
Simplified49.2%
if -4.2000000000000003e-226 < b < 1.7999999999999998e-247Initial program 99.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6411.5%
Simplified11.5%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6411.5%
Simplified11.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.9%
Simplified50.9%
if 1.7999999999999998e-247 < b Initial program 99.1%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.1%
Simplified66.1%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6463.9%
Simplified63.9%
Taylor expanded in b around 0
*-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6441.8%
Simplified41.8%
Final simplification47.7%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -2.5e+58)
(* b (* b (/ (* x (* b -0.16666666666666666)) y)))
(if (<= b -5.5e-229)
(/ x (* y a))
(if (<= b 1.35e-247)
(/ (* x (* (* b b) 0.5)) y)
(/ (/ x (* a (+ 1.0 b))) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -2.5e+58) {
tmp = b * (b * ((x * (b * -0.16666666666666666)) / y));
} else if (b <= -5.5e-229) {
tmp = x / (y * a);
} else if (b <= 1.35e-247) {
tmp = (x * ((b * b) * 0.5)) / y;
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
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 (b <= (-2.5d+58)) then
tmp = b * (b * ((x * (b * (-0.16666666666666666d0))) / y))
else if (b <= (-5.5d-229)) then
tmp = x / (y * a)
else if (b <= 1.35d-247) then
tmp = (x * ((b * b) * 0.5d0)) / y
else
tmp = (x / (a * (1.0d0 + b))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -2.5e+58) {
tmp = b * (b * ((x * (b * -0.16666666666666666)) / y));
} else if (b <= -5.5e-229) {
tmp = x / (y * a);
} else if (b <= 1.35e-247) {
tmp = (x * ((b * b) * 0.5)) / y;
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -2.5e+58: tmp = b * (b * ((x * (b * -0.16666666666666666)) / y)) elif b <= -5.5e-229: tmp = x / (y * a) elif b <= 1.35e-247: tmp = (x * ((b * b) * 0.5)) / y else: tmp = (x / (a * (1.0 + b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -2.5e+58) tmp = Float64(b * Float64(b * Float64(Float64(x * Float64(b * -0.16666666666666666)) / y))); elseif (b <= -5.5e-229) tmp = Float64(x / Float64(y * a)); elseif (b <= 1.35e-247) tmp = Float64(Float64(x * Float64(Float64(b * b) * 0.5)) / y); else tmp = Float64(Float64(x / Float64(a * Float64(1.0 + b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -2.5e+58) tmp = b * (b * ((x * (b * -0.16666666666666666)) / y)); elseif (b <= -5.5e-229) tmp = x / (y * a); elseif (b <= 1.35e-247) tmp = (x * ((b * b) * 0.5)) / y; else tmp = (x / (a * (1.0 + b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -2.5e+58], N[(b * N[(b * N[(N[(x * N[(b * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -5.5e-229], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.35e-247], N[(N[(x * N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a * N[(1.0 + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -2.5 \cdot 10^{+58}:\\
\;\;\;\;b \cdot \left(b \cdot \frac{x \cdot \left(b \cdot -0.16666666666666666\right)}{y}\right)\\
\mathbf{elif}\;b \leq -5.5 \cdot 10^{-229}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{elif}\;b \leq 1.35 \cdot 10^{-247}:\\
\;\;\;\;\frac{x \cdot \left(\left(b \cdot b\right) \cdot 0.5\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a \cdot \left(1 + b\right)}}{y}\\
\end{array}
\end{array}
if b < -2.49999999999999993e58Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6475.0%
Simplified75.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6467.0%
Simplified67.0%
Taylor expanded in b around inf
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
associate-*r/N/A
associate-*r/N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r/N/A
associate-*r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
Simplified63.8%
if -2.49999999999999993e58 < b < -5.5000000000000001e-229Initial program 98.7%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6465.4%
Simplified65.4%
Taylor expanded in b around 0
Simplified70.3%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f6440.5%
Simplified40.5%
if -5.5000000000000001e-229 < b < 1.35000000000000004e-247Initial program 99.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6411.5%
Simplified11.5%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6411.5%
Simplified11.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.9%
Simplified50.9%
if 1.35000000000000004e-247 < b Initial program 99.1%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.1%
Simplified66.1%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6463.9%
Simplified63.9%
Taylor expanded in b around 0
*-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6441.8%
Simplified41.8%
Final simplification47.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (* (* b b) 0.5)) y)))
(if (<= b -3.4e+92)
t_1
(if (<= b -1.35e-224)
(/ x (* y a))
(if (<= b 1.8e-247) t_1 (/ (/ x (* a (+ 1.0 b))) y))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * ((b * b) * 0.5)) / y;
double tmp;
if (b <= -3.4e+92) {
tmp = t_1;
} else if (b <= -1.35e-224) {
tmp = x / (y * a);
} else if (b <= 1.8e-247) {
tmp = t_1;
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
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 = (x * ((b * b) * 0.5d0)) / y
if (b <= (-3.4d+92)) then
tmp = t_1
else if (b <= (-1.35d-224)) then
tmp = x / (y * a)
else if (b <= 1.8d-247) then
tmp = t_1
else
tmp = (x / (a * (1.0d0 + b))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * ((b * b) * 0.5)) / y;
double tmp;
if (b <= -3.4e+92) {
tmp = t_1;
} else if (b <= -1.35e-224) {
tmp = x / (y * a);
} else if (b <= 1.8e-247) {
tmp = t_1;
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * ((b * b) * 0.5)) / y tmp = 0 if b <= -3.4e+92: tmp = t_1 elif b <= -1.35e-224: tmp = x / (y * a) elif b <= 1.8e-247: tmp = t_1 else: tmp = (x / (a * (1.0 + b))) / y return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * Float64(Float64(b * b) * 0.5)) / y) tmp = 0.0 if (b <= -3.4e+92) tmp = t_1; elseif (b <= -1.35e-224) tmp = Float64(x / Float64(y * a)); elseif (b <= 1.8e-247) tmp = t_1; else tmp = Float64(Float64(x / Float64(a * Float64(1.0 + b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * ((b * b) * 0.5)) / y; tmp = 0.0; if (b <= -3.4e+92) tmp = t_1; elseif (b <= -1.35e-224) tmp = x / (y * a); elseif (b <= 1.8e-247) tmp = t_1; else tmp = (x / (a * (1.0 + b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[b, -3.4e+92], t$95$1, If[LessEqual[b, -1.35e-224], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.8e-247], t$95$1, N[(N[(x / N[(a * N[(1.0 + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot \left(\left(b \cdot b\right) \cdot 0.5\right)}{y}\\
\mathbf{if}\;b \leq -3.4 \cdot 10^{+92}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -1.35 \cdot 10^{-224}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{elif}\;b \leq 1.8 \cdot 10^{-247}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a \cdot \left(1 + b\right)}}{y}\\
\end{array}
\end{array}
if b < -3.3999999999999998e92 or -1.34999999999999999e-224 < b < 1.7999999999999998e-247Initial program 99.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6457.1%
Simplified57.1%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6450.9%
Simplified50.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.2%
Simplified63.2%
if -3.3999999999999998e92 < b < -1.34999999999999999e-224Initial program 98.8%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6464.1%
Simplified64.1%
Taylor expanded in b around 0
Simplified68.7%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f6438.8%
Simplified38.8%
if 1.7999999999999998e-247 < b Initial program 99.1%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.1%
Simplified66.1%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6463.9%
Simplified63.9%
Taylor expanded in b around 0
*-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6441.8%
Simplified41.8%
Final simplification47.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= b 5e-122)
(* x (/ (* b (* (* b b) -0.16666666666666666)) y))
(/
(/ x (* a (+ 1.0 (* b (+ 1.0 (* b (+ 0.5 (* b 0.16666666666666666))))))))
y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 5e-122) {
tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y);
} else {
tmp = (x / (a * (1.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666)))))))) / y;
}
return tmp;
}
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 (b <= 5d-122) then
tmp = x * ((b * ((b * b) * (-0.16666666666666666d0))) / y)
else
tmp = (x / (a * (1.0d0 + (b * (1.0d0 + (b * (0.5d0 + (b * 0.16666666666666666d0)))))))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 5e-122) {
tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y);
} else {
tmp = (x / (a * (1.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666)))))))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 5e-122: tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y) else: tmp = (x / (a * (1.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666)))))))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 5e-122) tmp = Float64(x * Float64(Float64(b * Float64(Float64(b * b) * -0.16666666666666666)) / y)); else tmp = Float64(Float64(x / Float64(a * Float64(1.0 + Float64(b * Float64(1.0 + Float64(b * Float64(0.5 + Float64(b * 0.16666666666666666)))))))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 5e-122) tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y); else tmp = (x / (a * (1.0 + (b * (1.0 + (b * (0.5 + (b * 0.16666666666666666)))))))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 5e-122], N[(x * N[(N[(b * N[(N[(b * b), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(a * N[(1.0 + N[(b * N[(1.0 + N[(b * N[(0.5 + N[(b * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 5 \cdot 10^{-122}:\\
\;\;\;\;x \cdot \frac{b \cdot \left(\left(b \cdot b\right) \cdot -0.16666666666666666\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a \cdot \left(1 + b \cdot \left(1 + b \cdot \left(0.5 + b \cdot 0.16666666666666666\right)\right)\right)}}{y}\\
\end{array}
\end{array}
if b < 4.9999999999999999e-122Initial program 99.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6441.1%
Simplified41.1%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6435.4%
Simplified35.4%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.1%
Simplified46.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6449.0%
Applied egg-rr49.0%
if 4.9999999999999999e-122 < b Initial program 99.6%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.8%
Simplified66.8%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6473.7%
Simplified73.7%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6466.5%
Simplified66.5%
Final simplification54.3%
(FPCore (x y z t a b) :precision binary64 (if (<= b 5.6e-121) (* x (/ (* b (* (* b b) -0.16666666666666666)) y)) (/ (/ x (* a (+ 1.0 (* b (+ 1.0 (* b 0.5)))))) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 5.6e-121) {
tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y);
} else {
tmp = (x / (a * (1.0 + (b * (1.0 + (b * 0.5)))))) / y;
}
return tmp;
}
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 (b <= 5.6d-121) then
tmp = x * ((b * ((b * b) * (-0.16666666666666666d0))) / y)
else
tmp = (x / (a * (1.0d0 + (b * (1.0d0 + (b * 0.5d0)))))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 5.6e-121) {
tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y);
} else {
tmp = (x / (a * (1.0 + (b * (1.0 + (b * 0.5)))))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 5.6e-121: tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y) else: tmp = (x / (a * (1.0 + (b * (1.0 + (b * 0.5)))))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 5.6e-121) tmp = Float64(x * Float64(Float64(b * Float64(Float64(b * b) * -0.16666666666666666)) / y)); else tmp = Float64(Float64(x / Float64(a * Float64(1.0 + Float64(b * Float64(1.0 + Float64(b * 0.5)))))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 5.6e-121) tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y); else tmp = (x / (a * (1.0 + (b * (1.0 + (b * 0.5)))))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 5.6e-121], N[(x * N[(N[(b * N[(N[(b * b), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(a * N[(1.0 + N[(b * N[(1.0 + N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 5.6 \cdot 10^{-121}:\\
\;\;\;\;x \cdot \frac{b \cdot \left(\left(b \cdot b\right) \cdot -0.16666666666666666\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a \cdot \left(1 + b \cdot \left(1 + b \cdot 0.5\right)\right)}}{y}\\
\end{array}
\end{array}
if b < 5.6000000000000002e-121Initial program 99.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6441.1%
Simplified41.1%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6435.4%
Simplified35.4%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.1%
Simplified46.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6449.0%
Applied egg-rr49.0%
if 5.6000000000000002e-121 < b Initial program 99.6%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.8%
Simplified66.8%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6473.7%
Simplified73.7%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6462.8%
Simplified62.8%
Final simplification53.2%
(FPCore (x y z t a b) :precision binary64 (if (<= b 5.5e-121) (* x (/ (* b (* (* b b) -0.16666666666666666)) y)) (/ (/ x (* a (+ 1.0 b))) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 5.5e-121) {
tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y);
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
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 (b <= 5.5d-121) then
tmp = x * ((b * ((b * b) * (-0.16666666666666666d0))) / y)
else
tmp = (x / (a * (1.0d0 + b))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 5.5e-121) {
tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y);
} else {
tmp = (x / (a * (1.0 + b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 5.5e-121: tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y) else: tmp = (x / (a * (1.0 + b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 5.5e-121) tmp = Float64(x * Float64(Float64(b * Float64(Float64(b * b) * -0.16666666666666666)) / y)); else tmp = Float64(Float64(x / Float64(a * Float64(1.0 + b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 5.5e-121) tmp = x * ((b * ((b * b) * -0.16666666666666666)) / y); else tmp = (x / (a * (1.0 + b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 5.5e-121], N[(x * N[(N[(b * N[(N[(b * b), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(a * N[(1.0 + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 5.5 \cdot 10^{-121}:\\
\;\;\;\;x \cdot \frac{b \cdot \left(\left(b \cdot b\right) \cdot -0.16666666666666666\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a \cdot \left(1 + b\right)}}{y}\\
\end{array}
\end{array}
if b < 5.50000000000000031e-121Initial program 99.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6441.1%
Simplified41.1%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6435.4%
Simplified35.4%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.1%
Simplified46.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6449.0%
Applied egg-rr49.0%
if 5.50000000000000031e-121 < b Initial program 99.6%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.8%
Simplified66.8%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6473.7%
Simplified73.7%
Taylor expanded in b around 0
*-commutativeN/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f6443.1%
Simplified43.1%
Final simplification47.2%
(FPCore (x y z t a b) :precision binary64 (if (<= t 4.8e-18) (/ x (* y a)) (* b (/ x (* y b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 4.8e-18) {
tmp = x / (y * a);
} else {
tmp = b * (x / (y * b));
}
return tmp;
}
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 (t <= 4.8d-18) then
tmp = x / (y * a)
else
tmp = b * (x / (y * b))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 4.8e-18) {
tmp = x / (y * a);
} else {
tmp = b * (x / (y * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= 4.8e-18: tmp = x / (y * a) else: tmp = b * (x / (y * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= 4.8e-18) tmp = Float64(x / Float64(y * a)); else tmp = Float64(b * Float64(x / Float64(y * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= 4.8e-18) tmp = x / (y * a); else tmp = b * (x / (y * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, 4.8e-18], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(b * N[(x / N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 4.8 \cdot 10^{-18}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;b \cdot \frac{x}{y \cdot b}\\
\end{array}
\end{array}
if t < 4.79999999999999988e-18Initial program 99.0%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6465.6%
Simplified65.6%
Taylor expanded in b around 0
Simplified52.3%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f6434.4%
Simplified34.4%
if 4.79999999999999988e-18 < t Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6448.3%
Simplified48.3%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6429.5%
Simplified29.5%
Taylor expanded in b around inf
cube-multN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
Simplified24.5%
Taylor expanded in b around 0
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6437.5%
Simplified37.5%
Final simplification35.3%
(FPCore (x y z t a b) :precision binary64 (if (<= t 1.28e-5) (/ x (* y a)) (/ 1.0 (/ y x))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 1.28e-5) {
tmp = x / (y * a);
} else {
tmp = 1.0 / (y / x);
}
return tmp;
}
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 (t <= 1.28d-5) then
tmp = x / (y * a)
else
tmp = 1.0d0 / (y / x)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 1.28e-5) {
tmp = x / (y * a);
} else {
tmp = 1.0 / (y / x);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= 1.28e-5: tmp = x / (y * a) else: tmp = 1.0 / (y / x) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= 1.28e-5) tmp = Float64(x / Float64(y * a)); else tmp = Float64(1.0 / Float64(y / x)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= 1.28e-5) tmp = x / (y * a); else tmp = 1.0 / (y / x); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, 1.28e-5], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(y / x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 1.28 \cdot 10^{-5}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{y}{x}}\\
\end{array}
\end{array}
if t < 1.2799999999999999e-5Initial program 99.0%
Taylor expanded in y around 0
exp-diffN/A
associate-/l*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
exp-to-powN/A
pow-lowering-pow.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
exp-lowering-exp.f6466.0%
Simplified66.0%
Taylor expanded in b around 0
Simplified52.3%
Taylor expanded in t around 0
/-lowering-/.f64N/A
*-lowering-*.f6434.6%
Simplified34.6%
if 1.2799999999999999e-5 < t Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6448.1%
Simplified48.1%
Taylor expanded in b around 0
/-lowering-/.f6422.3%
Simplified22.3%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6422.5%
Applied egg-rr22.5%
Final simplification31.3%
(FPCore (x y z t a b) :precision binary64 (/ 1.0 (/ y x)))
double code(double x, double y, double z, double t, double a, double b) {
return 1.0 / (y / x);
}
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 = 1.0d0 / (y / x)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return 1.0 / (y / x);
}
def code(x, y, z, t, a, b): return 1.0 / (y / x)
function code(x, y, z, t, a, b) return Float64(1.0 / Float64(y / x)) end
function tmp = code(x, y, z, t, a, b) tmp = 1.0 / (y / x); end
code[x_, y_, z_, t_, a_, b_] := N[(1.0 / N[(y / x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{y}{x}}
\end{array}
Initial program 99.3%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6448.5%
Simplified48.5%
Taylor expanded in b around 0
/-lowering-/.f6416.1%
Simplified16.1%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6416.5%
Applied egg-rr16.5%
(FPCore (x y z t a b) :precision binary64 (/ x y))
double code(double x, double y, double z, double t, double a, double b) {
return x / y;
}
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
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / y;
}
def code(x, y, z, t, a, b): return x / y
function code(x, y, z, t, a, b) return Float64(x / y) end
function tmp = code(x, y, z, t, a, b) tmp = x / y; end
code[x_, y_, z_, t_, a_, b_] := N[(x / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y}
\end{array}
Initial program 99.3%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6448.5%
Simplified48.5%
Taylor expanded in b around 0
/-lowering-/.f6416.1%
Simplified16.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (pow a (- t 1.0)))
(t_2 (/ (* x (/ t_1 y)) (- (+ b 1.0) (* y (log z))))))
(if (< t -0.8845848504127471)
t_2
(if (< t 852031.2288374073)
(/ (* (/ x y) t_1) (exp (- b (* (log z) y))))
t_2))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = pow(a, (t - 1.0));
double t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * log(z)));
double tmp;
if (t < -0.8845848504127471) {
tmp = t_2;
} else if (t < 852031.2288374073) {
tmp = ((x / y) * t_1) / exp((b - (log(z) * y)));
} else {
tmp = t_2;
}
return tmp;
}
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) :: t_2
real(8) :: tmp
t_1 = a ** (t - 1.0d0)
t_2 = (x * (t_1 / y)) / ((b + 1.0d0) - (y * log(z)))
if (t < (-0.8845848504127471d0)) then
tmp = t_2
else if (t < 852031.2288374073d0) then
tmp = ((x / y) * t_1) / exp((b - (log(z) * y)))
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = Math.pow(a, (t - 1.0));
double t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * Math.log(z)));
double tmp;
if (t < -0.8845848504127471) {
tmp = t_2;
} else if (t < 852031.2288374073) {
tmp = ((x / y) * t_1) / Math.exp((b - (Math.log(z) * y)));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = math.pow(a, (t - 1.0)) t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * math.log(z))) tmp = 0 if t < -0.8845848504127471: tmp = t_2 elif t < 852031.2288374073: tmp = ((x / y) * t_1) / math.exp((b - (math.log(z) * y))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = a ^ Float64(t - 1.0) t_2 = Float64(Float64(x * Float64(t_1 / y)) / Float64(Float64(b + 1.0) - Float64(y * log(z)))) tmp = 0.0 if (t < -0.8845848504127471) tmp = t_2; elseif (t < 852031.2288374073) tmp = Float64(Float64(Float64(x / y) * t_1) / exp(Float64(b - Float64(log(z) * y)))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = a ^ (t - 1.0); t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * log(z))); tmp = 0.0; if (t < -0.8845848504127471) tmp = t_2; elseif (t < 852031.2288374073) tmp = ((x / y) * t_1) / exp((b - (log(z) * y))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * N[(t$95$1 / y), $MachinePrecision]), $MachinePrecision] / N[(N[(b + 1.0), $MachinePrecision] - N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -0.8845848504127471], t$95$2, If[Less[t, 852031.2288374073], N[(N[(N[(x / y), $MachinePrecision] * t$95$1), $MachinePrecision] / N[Exp[N[(b - N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {a}^{\left(t - 1\right)}\\
t_2 := \frac{x \cdot \frac{t\_1}{y}}{\left(b + 1\right) - y \cdot \log z}\\
\mathbf{if}\;t < -0.8845848504127471:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t < 852031.2288374073:\\
\;\;\;\;\frac{\frac{x}{y} \cdot t\_1}{e^{b - \log z \cdot y}}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024152
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:precision binary64
:alt
(! :herbie-platform default (if (< t -8845848504127471/10000000000000000) (/ (* x (/ (pow a (- t 1)) y)) (- (+ b 1) (* y (log z)))) (if (< t 8520312288374073/10000000000) (/ (* (/ x y) (pow a (- t 1))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1)) y)) (- (+ b 1) (* y (log z)))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))