
(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 24 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 97.7%
Final simplification97.7%
(FPCore (x y z t a b) :precision binary64 (if (or (<= y -1.12e+49) (not (<= y 1.05e-54))) (/ (* x (exp (- (- (* y (log z)) (log a)) b))) y) (/ (* x (exp (- (* (+ t -1.0) (log a)) b))) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((y <= -1.12e+49) || !(y <= 1.05e-54)) {
tmp = (x * exp((((y * log(z)) - log(a)) - b))) / y;
} else {
tmp = (x * exp((((t + -1.0) * log(a)) - 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 ((y <= (-1.12d+49)) .or. (.not. (y <= 1.05d-54))) then
tmp = (x * exp((((y * log(z)) - log(a)) - b))) / y
else
tmp = (x * exp((((t + (-1.0d0)) * log(a)) - 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 ((y <= -1.12e+49) || !(y <= 1.05e-54)) {
tmp = (x * Math.exp((((y * Math.log(z)) - Math.log(a)) - b))) / y;
} else {
tmp = (x * Math.exp((((t + -1.0) * Math.log(a)) - b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (y <= -1.12e+49) or not (y <= 1.05e-54): tmp = (x * math.exp((((y * math.log(z)) - math.log(a)) - b))) / y else: tmp = (x * math.exp((((t + -1.0) * math.log(a)) - b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((y <= -1.12e+49) || !(y <= 1.05e-54)) tmp = Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) - log(a)) - b))) / y); else tmp = Float64(Float64(x * exp(Float64(Float64(Float64(t + -1.0) * log(a)) - b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((y <= -1.12e+49) || ~((y <= 1.05e-54))) tmp = (x * exp((((y * log(z)) - log(a)) - b))) / y; else tmp = (x * exp((((t + -1.0) * log(a)) - b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -1.12e+49], N[Not[LessEqual[y, 1.05e-54]], $MachinePrecision]], N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] - N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x * N[Exp[N[(N[(N[(t + -1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.12 \cdot 10^{+49} \lor \neg \left(y \leq 1.05 \cdot 10^{-54}\right):\\
\;\;\;\;\frac{x \cdot e^{\left(y \cdot \log z - \log a\right) - b}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot e^{\left(t + -1\right) \cdot \log a - b}}{y}\\
\end{array}
\end{array}
if y < -1.12000000000000005e49 or 1.05e-54 < y Initial program 99.2%
Taylor expanded in t around 0 91.8%
+-commutative91.8%
mul-1-neg91.8%
unsub-neg91.8%
Simplified91.8%
if -1.12000000000000005e49 < y < 1.05e-54Initial program 96.2%
Taylor expanded in y around 0 95.7%
Final simplification93.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (* a (* y (exp b)))))
(t_2 (pow a (+ t -1.0)))
(t_3 (/ x (/ y (* (pow z y) t_2)))))
(if (<= b -1.8e+137)
t_1
(if (<= b -1.8e+90)
(/ (* x (/ (pow z y) a)) y)
(if (<= b -1.65e+44)
(/ x (/ y t_2))
(if (<= b -1.05e-85)
t_3
(if (<= b -5.5e-163)
(/ x (/ a (/ (pow z y) y)))
(if (<= b 2.3e-10) t_3 t_1))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / (a * (y * exp(b)));
double t_2 = pow(a, (t + -1.0));
double t_3 = x / (y / (pow(z, y) * t_2));
double tmp;
if (b <= -1.8e+137) {
tmp = t_1;
} else if (b <= -1.8e+90) {
tmp = (x * (pow(z, y) / a)) / y;
} else if (b <= -1.65e+44) {
tmp = x / (y / t_2);
} else if (b <= -1.05e-85) {
tmp = t_3;
} else if (b <= -5.5e-163) {
tmp = x / (a / (pow(z, y) / y));
} else if (b <= 2.3e-10) {
tmp = t_3;
} 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) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = x / (a * (y * exp(b)))
t_2 = a ** (t + (-1.0d0))
t_3 = x / (y / ((z ** y) * t_2))
if (b <= (-1.8d+137)) then
tmp = t_1
else if (b <= (-1.8d+90)) then
tmp = (x * ((z ** y) / a)) / y
else if (b <= (-1.65d+44)) then
tmp = x / (y / t_2)
else if (b <= (-1.05d-85)) then
tmp = t_3
else if (b <= (-5.5d-163)) then
tmp = x / (a / ((z ** y) / y))
else if (b <= 2.3d-10) then
tmp = t_3
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 / (a * (y * Math.exp(b)));
double t_2 = Math.pow(a, (t + -1.0));
double t_3 = x / (y / (Math.pow(z, y) * t_2));
double tmp;
if (b <= -1.8e+137) {
tmp = t_1;
} else if (b <= -1.8e+90) {
tmp = (x * (Math.pow(z, y) / a)) / y;
} else if (b <= -1.65e+44) {
tmp = x / (y / t_2);
} else if (b <= -1.05e-85) {
tmp = t_3;
} else if (b <= -5.5e-163) {
tmp = x / (a / (Math.pow(z, y) / y));
} else if (b <= 2.3e-10) {
tmp = t_3;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x / (a * (y * math.exp(b))) t_2 = math.pow(a, (t + -1.0)) t_3 = x / (y / (math.pow(z, y) * t_2)) tmp = 0 if b <= -1.8e+137: tmp = t_1 elif b <= -1.8e+90: tmp = (x * (math.pow(z, y) / a)) / y elif b <= -1.65e+44: tmp = x / (y / t_2) elif b <= -1.05e-85: tmp = t_3 elif b <= -5.5e-163: tmp = x / (a / (math.pow(z, y) / y)) elif b <= 2.3e-10: tmp = t_3 else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x / Float64(a * Float64(y * exp(b)))) t_2 = a ^ Float64(t + -1.0) t_3 = Float64(x / Float64(y / Float64((z ^ y) * t_2))) tmp = 0.0 if (b <= -1.8e+137) tmp = t_1; elseif (b <= -1.8e+90) tmp = Float64(Float64(x * Float64((z ^ y) / a)) / y); elseif (b <= -1.65e+44) tmp = Float64(x / Float64(y / t_2)); elseif (b <= -1.05e-85) tmp = t_3; elseif (b <= -5.5e-163) tmp = Float64(x / Float64(a / Float64((z ^ y) / y))); elseif (b <= 2.3e-10) tmp = t_3; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x / (a * (y * exp(b))); t_2 = a ^ (t + -1.0); t_3 = x / (y / ((z ^ y) * t_2)); tmp = 0.0; if (b <= -1.8e+137) tmp = t_1; elseif (b <= -1.8e+90) tmp = (x * ((z ^ y) / a)) / y; elseif (b <= -1.65e+44) tmp = x / (y / t_2); elseif (b <= -1.05e-85) tmp = t_3; elseif (b <= -5.5e-163) tmp = x / (a / ((z ^ y) / y)); elseif (b <= 2.3e-10) tmp = t_3; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[(x / N[(y / N[(N[Power[z, y], $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.8e+137], t$95$1, If[LessEqual[b, -1.8e+90], N[(N[(x * N[(N[Power[z, y], $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -1.65e+44], N[(x / N[(y / t$95$2), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.05e-85], t$95$3, If[LessEqual[b, -5.5e-163], N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.3e-10], t$95$3, t$95$1]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{a \cdot \left(y \cdot e^{b}\right)}\\
t_2 := {a}^{\left(t + -1\right)}\\
t_3 := \frac{x}{\frac{y}{{z}^{y} \cdot t_2}}\\
\mathbf{if}\;b \leq -1.8 \cdot 10^{+137}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;b \leq -1.8 \cdot 10^{+90}:\\
\;\;\;\;\frac{x \cdot \frac{{z}^{y}}{a}}{y}\\
\mathbf{elif}\;b \leq -1.65 \cdot 10^{+44}:\\
\;\;\;\;\frac{x}{\frac{y}{t_2}}\\
\mathbf{elif}\;b \leq -1.05 \cdot 10^{-85}:\\
\;\;\;\;t_3\\
\mathbf{elif}\;b \leq -5.5 \cdot 10^{-163}:\\
\;\;\;\;\frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\mathbf{elif}\;b \leq 2.3 \cdot 10^{-10}:\\
\;\;\;\;t_3\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if b < -1.8e137 or 2.30000000000000007e-10 < b Initial program 99.1%
associate-*l/90.3%
*-commutative90.3%
+-commutative90.3%
associate--l+90.3%
exp-sum66.5%
*-commutative66.5%
exp-to-pow66.7%
sub-neg66.7%
metadata-eval66.7%
exp-diff51.4%
*-commutative51.4%
exp-to-pow51.4%
Simplified51.4%
Taylor expanded in t around 0 70.6%
times-frac65.0%
Simplified65.0%
Taylor expanded in y around 0 84.1%
if -1.8e137 < b < -1.8e90Initial program 100.0%
Taylor expanded in t around 0 91.1%
+-commutative91.1%
mul-1-neg91.1%
unsub-neg91.1%
Simplified91.1%
Taylor expanded in b around 0 91.1%
div-exp91.1%
*-commutative91.1%
exp-to-pow91.1%
rem-exp-log91.1%
Simplified91.1%
if -1.8e90 < b < -1.65000000000000007e44Initial program 100.0%
associate-*l/85.7%
*-commutative85.7%
+-commutative85.7%
associate--l+85.7%
exp-sum57.1%
*-commutative57.1%
exp-to-pow57.1%
sub-neg57.1%
metadata-eval57.1%
exp-diff42.9%
*-commutative42.9%
exp-to-pow42.9%
Simplified42.9%
Taylor expanded in b around 0 57.1%
associate-/l*57.1%
*-commutative57.1%
exp-to-pow57.1%
*-commutative57.1%
exp-sum100.0%
exp-sum57.1%
*-commutative57.1%
exp-to-pow57.1%
*-commutative57.1%
exp-to-pow57.1%
sub-neg57.1%
metadata-eval57.1%
Simplified57.1%
Taylor expanded in y around 0 86.2%
exp-to-pow86.2%
Simplified86.2%
if -1.65000000000000007e44 < b < -1.05e-85 or -5.4999999999999998e-163 < b < 2.30000000000000007e-10Initial program 96.8%
associate-*l/89.2%
*-commutative89.2%
+-commutative89.2%
associate--l+89.2%
exp-sum83.9%
*-commutative83.9%
exp-to-pow85.1%
sub-neg85.1%
metadata-eval85.1%
exp-diff84.3%
*-commutative84.3%
exp-to-pow84.3%
Simplified84.3%
Taylor expanded in b around 0 90.8%
associate-/l*90.9%
*-commutative90.9%
exp-to-pow90.9%
*-commutative90.9%
exp-sum96.2%
exp-sum90.9%
*-commutative90.9%
exp-to-pow90.9%
*-commutative90.9%
exp-to-pow92.1%
sub-neg92.1%
metadata-eval92.1%
Simplified92.1%
if -1.05e-85 < b < -5.4999999999999998e-163Initial program 93.5%
associate-*l/82.3%
*-commutative82.3%
+-commutative82.3%
associate--l+82.3%
exp-sum61.3%
*-commutative61.3%
exp-to-pow63.0%
sub-neg63.0%
metadata-eval63.0%
exp-diff63.0%
*-commutative63.0%
exp-to-pow63.0%
Simplified63.0%
Taylor expanded in b around 0 61.9%
associate-/l*66.5%
*-commutative66.5%
exp-to-pow66.5%
*-commutative66.5%
exp-sum98.1%
exp-sum66.5%
*-commutative66.5%
exp-to-pow66.5%
*-commutative66.5%
exp-to-pow68.4%
sub-neg68.4%
metadata-eval68.4%
Simplified68.4%
Taylor expanded in t around 0 74.0%
associate-/l*89.8%
Simplified89.8%
Final simplification88.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (exp (- (* (+ t -1.0) (log a)) b))) y))
(t_2 (/ x (/ a (/ (pow z y) y)))))
(if (<= y -1.6e+103)
t_2
(if (<= y 1.05e-54)
t_1
(if (<= y 9.5e+18)
(/ (* x (pow z y)) (* a (* y (exp b))))
(if (<= y 4.9e+85) t_1 t_2))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * exp((((t + -1.0) * log(a)) - b))) / y;
double t_2 = x / (a / (pow(z, y) / y));
double tmp;
if (y <= -1.6e+103) {
tmp = t_2;
} else if (y <= 1.05e-54) {
tmp = t_1;
} else if (y <= 9.5e+18) {
tmp = (x * pow(z, y)) / (a * (y * exp(b)));
} else if (y <= 4.9e+85) {
tmp = t_1;
} 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 = (x * exp((((t + (-1.0d0)) * log(a)) - b))) / y
t_2 = x / (a / ((z ** y) / y))
if (y <= (-1.6d+103)) then
tmp = t_2
else if (y <= 1.05d-54) then
tmp = t_1
else if (y <= 9.5d+18) then
tmp = (x * (z ** y)) / (a * (y * exp(b)))
else if (y <= 4.9d+85) then
tmp = t_1
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 = (x * Math.exp((((t + -1.0) * Math.log(a)) - b))) / y;
double t_2 = x / (a / (Math.pow(z, y) / y));
double tmp;
if (y <= -1.6e+103) {
tmp = t_2;
} else if (y <= 1.05e-54) {
tmp = t_1;
} else if (y <= 9.5e+18) {
tmp = (x * Math.pow(z, y)) / (a * (y * Math.exp(b)));
} else if (y <= 4.9e+85) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * math.exp((((t + -1.0) * math.log(a)) - b))) / y t_2 = x / (a / (math.pow(z, y) / y)) tmp = 0 if y <= -1.6e+103: tmp = t_2 elif y <= 1.05e-54: tmp = t_1 elif y <= 9.5e+18: tmp = (x * math.pow(z, y)) / (a * (y * math.exp(b))) elif y <= 4.9e+85: tmp = t_1 else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * exp(Float64(Float64(Float64(t + -1.0) * log(a)) - b))) / y) t_2 = Float64(x / Float64(a / Float64((z ^ y) / y))) tmp = 0.0 if (y <= -1.6e+103) tmp = t_2; elseif (y <= 1.05e-54) tmp = t_1; elseif (y <= 9.5e+18) tmp = Float64(Float64(x * (z ^ y)) / Float64(a * Float64(y * exp(b)))); elseif (y <= 4.9e+85) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * exp((((t + -1.0) * log(a)) - b))) / y; t_2 = x / (a / ((z ^ y) / y)); tmp = 0.0; if (y <= -1.6e+103) tmp = t_2; elseif (y <= 1.05e-54) tmp = t_1; elseif (y <= 9.5e+18) tmp = (x * (z ^ y)) / (a * (y * exp(b))); elseif (y <= 4.9e+85) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Exp[N[(N[(N[(t + -1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, Block[{t$95$2 = N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.6e+103], t$95$2, If[LessEqual[y, 1.05e-54], t$95$1, If[LessEqual[y, 9.5e+18], N[(N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision] / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.9e+85], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot e^{\left(t + -1\right) \cdot \log a - b}}{y}\\
t_2 := \frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\mathbf{if}\;y \leq -1.6 \cdot 10^{+103}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;y \leq 1.05 \cdot 10^{-54}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y \leq 9.5 \cdot 10^{+18}:\\
\;\;\;\;\frac{x \cdot {z}^{y}}{a \cdot \left(y \cdot e^{b}\right)}\\
\mathbf{elif}\;y \leq 4.9 \cdot 10^{+85}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
\end{array}
if y < -1.59999999999999996e103 or 4.8999999999999997e85 < y Initial program 100.0%
associate-*l/88.4%
*-commutative88.4%
+-commutative88.4%
associate--l+88.4%
exp-sum62.8%
*-commutative62.8%
exp-to-pow62.8%
sub-neg62.8%
metadata-eval62.8%
exp-diff46.5%
*-commutative46.5%
exp-to-pow46.5%
Simplified46.5%
Taylor expanded in b around 0 66.4%
associate-/l*66.4%
*-commutative66.4%
exp-to-pow66.4%
*-commutative66.4%
exp-sum92.0%
exp-sum66.4%
*-commutative66.4%
exp-to-pow66.4%
*-commutative66.4%
exp-to-pow66.4%
sub-neg66.4%
metadata-eval66.4%
Simplified66.4%
Taylor expanded in t around 0 69.9%
associate-/l*87.4%
Simplified87.4%
if -1.59999999999999996e103 < y < 1.05e-54 or 9.5e18 < y < 4.8999999999999997e85Initial program 96.9%
Taylor expanded in y around 0 92.7%
if 1.05e-54 < y < 9.5e18Initial program 93.9%
associate-*l/98.1%
*-commutative98.1%
+-commutative98.1%
associate--l+98.1%
exp-sum92.8%
*-commutative92.8%
exp-to-pow94.1%
sub-neg94.1%
metadata-eval94.1%
exp-diff94.1%
*-commutative94.1%
exp-to-pow94.1%
Simplified94.1%
Taylor expanded in t around 0 99.8%
Final simplification91.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (pow a (+ t -1.0)))
(t_2 (* (/ t_1 (exp b)) (/ x y)))
(t_3 (/ x (/ a (/ (pow z y) y)))))
(if (<= y -2.8e+49)
t_3
(if (<= y -9.5e-224)
t_2
(if (<= y 1.12e-129) (/ x (/ y t_1)) (if (<= y 160000.0) t_2 t_3))))))
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 = (t_1 / exp(b)) * (x / y);
double t_3 = x / (a / (pow(z, y) / y));
double tmp;
if (y <= -2.8e+49) {
tmp = t_3;
} else if (y <= -9.5e-224) {
tmp = t_2;
} else if (y <= 1.12e-129) {
tmp = x / (y / t_1);
} else if (y <= 160000.0) {
tmp = 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 = a ** (t + (-1.0d0))
t_2 = (t_1 / exp(b)) * (x / y)
t_3 = x / (a / ((z ** y) / y))
if (y <= (-2.8d+49)) then
tmp = t_3
else if (y <= (-9.5d-224)) then
tmp = t_2
else if (y <= 1.12d-129) then
tmp = x / (y / t_1)
else if (y <= 160000.0d0) then
tmp = 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 = Math.pow(a, (t + -1.0));
double t_2 = (t_1 / Math.exp(b)) * (x / y);
double t_3 = x / (a / (Math.pow(z, y) / y));
double tmp;
if (y <= -2.8e+49) {
tmp = t_3;
} else if (y <= -9.5e-224) {
tmp = t_2;
} else if (y <= 1.12e-129) {
tmp = x / (y / t_1);
} else if (y <= 160000.0) {
tmp = t_2;
} else {
tmp = t_3;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = math.pow(a, (t + -1.0)) t_2 = (t_1 / math.exp(b)) * (x / y) t_3 = x / (a / (math.pow(z, y) / y)) tmp = 0 if y <= -2.8e+49: tmp = t_3 elif y <= -9.5e-224: tmp = t_2 elif y <= 1.12e-129: tmp = x / (y / t_1) elif y <= 160000.0: tmp = t_2 else: tmp = t_3 return tmp
function code(x, y, z, t, a, b) t_1 = a ^ Float64(t + -1.0) t_2 = Float64(Float64(t_1 / exp(b)) * Float64(x / y)) t_3 = Float64(x / Float64(a / Float64((z ^ y) / y))) tmp = 0.0 if (y <= -2.8e+49) tmp = t_3; elseif (y <= -9.5e-224) tmp = t_2; elseif (y <= 1.12e-129) tmp = Float64(x / Float64(y / t_1)); elseif (y <= 160000.0) tmp = t_2; else tmp = t_3; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = a ^ (t + -1.0); t_2 = (t_1 / exp(b)) * (x / y); t_3 = x / (a / ((z ^ y) / y)); tmp = 0.0; if (y <= -2.8e+49) tmp = t_3; elseif (y <= -9.5e-224) tmp = t_2; elseif (y <= 1.12e-129) tmp = x / (y / t_1); elseif (y <= 160000.0) tmp = t_2; else tmp = t_3; 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[(t$95$1 / N[Exp[b], $MachinePrecision]), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.8e+49], t$95$3, If[LessEqual[y, -9.5e-224], t$95$2, If[LessEqual[y, 1.12e-129], N[(x / N[(y / t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 160000.0], t$95$2, t$95$3]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {a}^{\left(t + -1\right)}\\
t_2 := \frac{t_1}{e^{b}} \cdot \frac{x}{y}\\
t_3 := \frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\mathbf{if}\;y \leq -2.8 \cdot 10^{+49}:\\
\;\;\;\;t_3\\
\mathbf{elif}\;y \leq -9.5 \cdot 10^{-224}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;y \leq 1.12 \cdot 10^{-129}:\\
\;\;\;\;\frac{x}{\frac{y}{t_1}}\\
\mathbf{elif}\;y \leq 160000:\\
\;\;\;\;t_2\\
\mathbf{else}:\\
\;\;\;\;t_3\\
\end{array}
\end{array}
if y < -2.7999999999999998e49 or 1.6e5 < y Initial program 100.0%
associate-*l/90.0%
*-commutative90.0%
+-commutative90.0%
associate--l+90.0%
exp-sum65.8%
*-commutative65.8%
exp-to-pow65.8%
sub-neg65.8%
metadata-eval65.8%
exp-diff50.0%
*-commutative50.0%
exp-to-pow50.0%
Simplified50.0%
Taylor expanded in b around 0 65.1%
associate-/l*65.1%
*-commutative65.1%
exp-to-pow65.1%
*-commutative65.1%
exp-sum89.3%
exp-sum65.1%
*-commutative65.1%
exp-to-pow65.1%
*-commutative65.1%
exp-to-pow65.1%
sub-neg65.1%
metadata-eval65.1%
Simplified65.1%
Taylor expanded in t around 0 67.7%
associate-/l*82.8%
Simplified82.8%
if -2.7999999999999998e49 < y < -9.5000000000000003e-224 or 1.12000000000000006e-129 < y < 1.6e5Initial program 96.2%
associate-*l/94.6%
*-commutative94.6%
+-commutative94.6%
associate--l+94.6%
exp-sum85.7%
*-commutative85.7%
exp-to-pow87.1%
sub-neg87.1%
metadata-eval87.1%
exp-diff84.5%
*-commutative84.5%
exp-to-pow84.5%
Simplified84.5%
Taylor expanded in y around 0 87.8%
*-commutative87.8%
*-commutative87.8%
times-frac86.2%
exp-to-pow87.5%
sub-neg87.5%
metadata-eval87.5%
Simplified87.5%
if -9.5000000000000003e-224 < y < 1.12000000000000006e-129Initial program 95.1%
associate-*l/79.9%
*-commutative79.9%
+-commutative79.9%
associate--l+79.9%
exp-sum69.5%
*-commutative69.5%
exp-to-pow70.8%
sub-neg70.8%
metadata-eval70.8%
exp-diff70.8%
*-commutative70.8%
exp-to-pow70.8%
Simplified70.8%
Taylor expanded in b around 0 80.4%
associate-/l*77.1%
*-commutative77.1%
exp-to-pow77.1%
*-commutative77.1%
exp-sum77.1%
exp-sum77.1%
*-commutative77.1%
exp-to-pow77.1%
*-commutative77.1%
exp-to-pow78.5%
sub-neg78.5%
metadata-eval78.5%
Simplified78.5%
Taylor expanded in y around 0 77.1%
exp-to-pow78.5%
Simplified78.5%
Final simplification83.2%
(FPCore (x y z t a b) :precision binary64 (if (or (<= (+ t -1.0) -5e+67) (not (<= (+ t -1.0) -0.5))) (/ x (/ y (pow a (+ t -1.0)))) (/ (* x (pow z y)) (* a (* y (exp b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (((t + -1.0) <= -5e+67) || !((t + -1.0) <= -0.5)) {
tmp = x / (y / pow(a, (t + -1.0)));
} else {
tmp = (x * pow(z, y)) / (a * (y * exp(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 + (-1.0d0)) <= (-5d+67)) .or. (.not. ((t + (-1.0d0)) <= (-0.5d0)))) then
tmp = x / (y / (a ** (t + (-1.0d0))))
else
tmp = (x * (z ** y)) / (a * (y * exp(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 + -1.0) <= -5e+67) || !((t + -1.0) <= -0.5)) {
tmp = x / (y / Math.pow(a, (t + -1.0)));
} else {
tmp = (x * Math.pow(z, y)) / (a * (y * Math.exp(b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if ((t + -1.0) <= -5e+67) or not ((t + -1.0) <= -0.5): tmp = x / (y / math.pow(a, (t + -1.0))) else: tmp = (x * math.pow(z, y)) / (a * (y * math.exp(b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((Float64(t + -1.0) <= -5e+67) || !(Float64(t + -1.0) <= -0.5)) tmp = Float64(x / Float64(y / (a ^ Float64(t + -1.0)))); else tmp = Float64(Float64(x * (z ^ y)) / Float64(a * Float64(y * exp(b)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (((t + -1.0) <= -5e+67) || ~(((t + -1.0) <= -0.5))) tmp = x / (y / (a ^ (t + -1.0))); else tmp = (x * (z ^ y)) / (a * (y * exp(b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[N[(t + -1.0), $MachinePrecision], -5e+67], N[Not[LessEqual[N[(t + -1.0), $MachinePrecision], -0.5]], $MachinePrecision]], N[(x / N[(y / N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision] / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t + -1 \leq -5 \cdot 10^{+67} \lor \neg \left(t + -1 \leq -0.5\right):\\
\;\;\;\;\frac{x}{\frac{y}{{a}^{\left(t + -1\right)}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot {z}^{y}}{a \cdot \left(y \cdot e^{b}\right)}\\
\end{array}
\end{array}
if (-.f64 t 1) < -4.99999999999999976e67 or -0.5 < (-.f64 t 1) Initial program 100.0%
associate-*l/90.5%
*-commutative90.5%
+-commutative90.5%
associate--l+90.5%
exp-sum55.2%
*-commutative55.2%
exp-to-pow55.2%
sub-neg55.2%
metadata-eval55.2%
exp-diff50.5%
*-commutative50.5%
exp-to-pow50.5%
Simplified50.5%
Taylor expanded in b around 0 62.9%
associate-/l*62.9%
*-commutative62.9%
exp-to-pow62.9%
*-commutative62.9%
exp-sum89.7%
exp-sum62.9%
*-commutative62.9%
exp-to-pow62.9%
*-commutative62.9%
exp-to-pow62.9%
sub-neg62.9%
metadata-eval62.9%
Simplified62.9%
Taylor expanded in y around 0 78.4%
exp-to-pow78.4%
Simplified78.4%
if -4.99999999999999976e67 < (-.f64 t 1) < -0.5Initial program 96.2%
associate-*l/88.1%
*-commutative88.1%
+-commutative88.1%
associate--l+88.1%
exp-sum84.9%
*-commutative84.9%
exp-to-pow86.1%
sub-neg86.1%
metadata-eval86.1%
exp-diff75.5%
*-commutative75.5%
exp-to-pow75.5%
Simplified75.5%
Taylor expanded in t around 0 84.1%
Final simplification81.8%
(FPCore (x y z t a b) :precision binary64 (if (or (<= b -1.8e+137) (not (<= b 5.2e-10))) (/ x (* a (* y (exp b)))) (/ x (/ a (/ (pow z y) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((b <= -1.8e+137) || !(b <= 5.2e-10)) {
tmp = x / (a * (y * exp(b)));
} else {
tmp = x / (a / (pow(z, y) / 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.8d+137)) .or. (.not. (b <= 5.2d-10))) then
tmp = x / (a * (y * exp(b)))
else
tmp = x / (a / ((z ** y) / 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.8e+137) || !(b <= 5.2e-10)) {
tmp = x / (a * (y * Math.exp(b)));
} else {
tmp = x / (a / (Math.pow(z, y) / y));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (b <= -1.8e+137) or not (b <= 5.2e-10): tmp = x / (a * (y * math.exp(b))) else: tmp = x / (a / (math.pow(z, y) / y)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((b <= -1.8e+137) || !(b <= 5.2e-10)) tmp = Float64(x / Float64(a * Float64(y * exp(b)))); else tmp = Float64(x / Float64(a / Float64((z ^ y) / y))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((b <= -1.8e+137) || ~((b <= 5.2e-10))) tmp = x / (a * (y * exp(b))); else tmp = x / (a / ((z ^ y) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[b, -1.8e+137], N[Not[LessEqual[b, 5.2e-10]], $MachinePrecision]], N[(x / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.8 \cdot 10^{+137} \lor \neg \left(b \leq 5.2 \cdot 10^{-10}\right):\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot e^{b}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\end{array}
\end{array}
if b < -1.8e137 or 5.19999999999999962e-10 < b Initial program 99.1%
associate-*l/90.2%
*-commutative90.2%
+-commutative90.2%
associate--l+90.2%
exp-sum67.2%
*-commutative67.2%
exp-to-pow67.3%
sub-neg67.3%
metadata-eval67.3%
exp-diff51.9%
*-commutative51.9%
exp-to-pow51.9%
Simplified51.9%
Taylor expanded in t around 0 70.3%
times-frac64.6%
Simplified64.6%
Taylor expanded in y around 0 84.8%
if -1.8e137 < b < 5.19999999999999962e-10Initial program 96.8%
associate-*l/88.3%
*-commutative88.3%
+-commutative88.3%
associate--l+88.3%
exp-sum76.5%
*-commutative76.5%
exp-to-pow77.6%
sub-neg77.6%
metadata-eval77.6%
exp-diff74.4%
*-commutative74.4%
exp-to-pow74.4%
Simplified74.4%
Taylor expanded in b around 0 83.1%
associate-/l*83.7%
*-commutative83.7%
exp-to-pow83.7%
*-commutative83.7%
exp-sum96.2%
exp-sum83.7%
*-commutative83.7%
exp-to-pow83.7%
*-commutative83.7%
exp-to-pow84.9%
sub-neg84.9%
metadata-eval84.9%
Simplified84.9%
Taylor expanded in t around 0 69.7%
associate-/l*77.0%
Simplified77.0%
Final simplification80.2%
(FPCore (x y z t a b) :precision binary64 (if (<= y 1.02e+186) (/ x (* a (* y (exp b)))) (/ (* x 2.0) (* y (* a (* b b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 1.02e+186) {
tmp = x / (a * (y * exp(b)));
} else {
tmp = (x * 2.0) / (y * (a * (b * 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 (y <= 1.02d+186) then
tmp = x / (a * (y * exp(b)))
else
tmp = (x * 2.0d0) / (y * (a * (b * 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 (y <= 1.02e+186) {
tmp = x / (a * (y * Math.exp(b)));
} else {
tmp = (x * 2.0) / (y * (a * (b * b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= 1.02e+186: tmp = x / (a * (y * math.exp(b))) else: tmp = (x * 2.0) / (y * (a * (b * b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= 1.02e+186) tmp = Float64(x / Float64(a * Float64(y * exp(b)))); else tmp = Float64(Float64(x * 2.0) / Float64(y * Float64(a * Float64(b * b)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= 1.02e+186) tmp = x / (a * (y * exp(b))); else tmp = (x * 2.0) / (y * (a * (b * b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, 1.02e+186], N[(x / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.02 \cdot 10^{+186}:\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot e^{b}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot 2}{y \cdot \left(a \cdot \left(b \cdot b\right)\right)}\\
\end{array}
\end{array}
if y < 1.01999999999999999e186Initial program 97.5%
associate-*l/90.0%
*-commutative90.0%
+-commutative90.0%
associate--l+90.0%
exp-sum74.7%
*-commutative74.7%
exp-to-pow75.5%
sub-neg75.5%
metadata-eval75.5%
exp-diff67.7%
*-commutative67.7%
exp-to-pow67.7%
Simplified67.7%
Taylor expanded in t around 0 70.3%
times-frac66.0%
Simplified66.0%
Taylor expanded in y around 0 66.6%
if 1.01999999999999999e186 < y Initial program 100.0%
associate-*l/81.5%
*-commutative81.5%
+-commutative81.5%
associate--l+81.5%
exp-sum55.6%
*-commutative55.6%
exp-to-pow55.6%
sub-neg55.6%
metadata-eval55.6%
exp-diff44.4%
*-commutative44.4%
exp-to-pow44.4%
Simplified44.4%
Taylor expanded in t around 0 55.6%
times-frac66.7%
Simplified66.7%
Taylor expanded in y around 0 20.3%
Taylor expanded in b around 0 16.3%
+-commutative16.3%
associate-*r*16.3%
distribute-rgt-out34.9%
unpow234.9%
Simplified34.9%
Taylor expanded in b around inf 45.7%
associate-*r/45.7%
unpow245.7%
associate-*r*42.2%
Simplified42.2%
Final simplification64.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x 2.0) (* y (* a (* b b))))))
(if (<= b -3.8e-146)
(- (/ x (* y a)) (* (/ x a) (/ b y)))
(if (<= b -5.5e-248)
t_1
(if (<= b -5.8e-289)
(/ x (/ y (- (/ 1.0 a) (/ b a))))
(if (<= b 7e-240)
t_1
(/ x (* a (+ y (* y (+ b (* (* b b) 0.5))))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -3.8e-146) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= -5.5e-248) {
tmp = t_1;
} else if (b <= -5.8e-289) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if (b <= 7e-240) {
tmp = t_1;
} else {
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5)))));
}
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 * 2.0d0) / (y * (a * (b * b)))
if (b <= (-3.8d-146)) then
tmp = (x / (y * a)) - ((x / a) * (b / y))
else if (b <= (-5.5d-248)) then
tmp = t_1
else if (b <= (-5.8d-289)) then
tmp = x / (y / ((1.0d0 / a) - (b / a)))
else if (b <= 7d-240) then
tmp = t_1
else
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5d0)))))
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 * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -3.8e-146) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= -5.5e-248) {
tmp = t_1;
} else if (b <= -5.8e-289) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if (b <= 7e-240) {
tmp = t_1;
} else {
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5)))));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * 2.0) / (y * (a * (b * b))) tmp = 0 if b <= -3.8e-146: tmp = (x / (y * a)) - ((x / a) * (b / y)) elif b <= -5.5e-248: tmp = t_1 elif b <= -5.8e-289: tmp = x / (y / ((1.0 / a) - (b / a))) elif b <= 7e-240: tmp = t_1 else: tmp = x / (a * (y + (y * (b + ((b * b) * 0.5))))) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * 2.0) / Float64(y * Float64(a * Float64(b * b)))) tmp = 0.0 if (b <= -3.8e-146) tmp = Float64(Float64(x / Float64(y * a)) - Float64(Float64(x / a) * Float64(b / y))); elseif (b <= -5.5e-248) tmp = t_1; elseif (b <= -5.8e-289) tmp = Float64(x / Float64(y / Float64(Float64(1.0 / a) - Float64(b / a)))); elseif (b <= 7e-240) tmp = t_1; else tmp = Float64(x / Float64(a * Float64(y + Float64(y * Float64(b + Float64(Float64(b * b) * 0.5)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * 2.0) / (y * (a * (b * b))); tmp = 0.0; if (b <= -3.8e-146) tmp = (x / (y * a)) - ((x / a) * (b / y)); elseif (b <= -5.5e-248) tmp = t_1; elseif (b <= -5.8e-289) tmp = x / (y / ((1.0 / a) - (b / a))); elseif (b <= 7e-240) tmp = t_1; else tmp = x / (a * (y + (y * (b + ((b * b) * 0.5))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -3.8e-146], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] - N[(N[(x / a), $MachinePrecision] * N[(b / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -5.5e-248], t$95$1, If[LessEqual[b, -5.8e-289], N[(x / N[(y / N[(N[(1.0 / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 7e-240], t$95$1, N[(x / N[(a * N[(y + N[(y * N[(b + N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot 2}{y \cdot \left(a \cdot \left(b \cdot b\right)\right)}\\
\mathbf{if}\;b \leq -3.8 \cdot 10^{-146}:\\
\;\;\;\;\frac{x}{y \cdot a} - \frac{x}{a} \cdot \frac{b}{y}\\
\mathbf{elif}\;b \leq -5.5 \cdot 10^{-248}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;b \leq -5.8 \cdot 10^{-289}:\\
\;\;\;\;\frac{x}{\frac{y}{\frac{1}{a} - \frac{b}{a}}}\\
\mathbf{elif}\;b \leq 7 \cdot 10^{-240}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot \left(b + \left(b \cdot b\right) \cdot 0.5\right)\right)}\\
\end{array}
\end{array}
if b < -3.79999999999999994e-146Initial program 98.1%
associate-*l/89.5%
*-commutative89.5%
+-commutative89.5%
associate--l+89.5%
exp-sum67.9%
*-commutative67.9%
exp-to-pow68.6%
sub-neg68.6%
metadata-eval68.6%
exp-diff58.9%
*-commutative58.9%
exp-to-pow58.9%
Simplified58.9%
Taylor expanded in t around 0 66.6%
times-frac64.1%
Simplified64.1%
Taylor expanded in y around 0 63.6%
Taylor expanded in b around 0 52.2%
+-commutative52.2%
mul-1-neg52.2%
unsub-neg52.2%
*-commutative52.2%
*-commutative52.2%
times-frac55.7%
Simplified55.7%
if -3.79999999999999994e-146 < b < -5.49999999999999979e-248 or -5.80000000000000012e-289 < b < 7.00000000000000032e-240Initial program 96.7%
associate-*l/90.7%
*-commutative90.7%
+-commutative90.7%
associate--l+90.7%
exp-sum78.9%
*-commutative78.9%
exp-to-pow79.5%
sub-neg79.5%
metadata-eval79.5%
exp-diff79.5%
*-commutative79.5%
exp-to-pow79.5%
Simplified79.5%
Taylor expanded in t around 0 62.4%
times-frac65.3%
Simplified65.3%
Taylor expanded in y around 0 31.3%
Taylor expanded in b around 0 31.3%
+-commutative31.3%
associate-*r*31.3%
distribute-rgt-out31.3%
unpow231.3%
Simplified31.3%
Taylor expanded in b around inf 59.7%
associate-*r/59.7%
unpow259.7%
associate-*r*59.7%
Simplified59.7%
if -5.49999999999999979e-248 < b < -5.80000000000000012e-289Initial program 97.9%
Taylor expanded in t around 0 78.5%
+-commutative78.5%
mul-1-neg78.5%
unsub-neg78.5%
Simplified78.5%
Taylor expanded in y around 0 51.6%
exp-neg51.6%
associate-*r/51.6%
*-rgt-identity51.6%
+-commutative51.6%
exp-sum51.6%
rem-exp-log53.1%
Simplified53.1%
Taylor expanded in b around 0 53.1%
+-commutative53.1%
mul-1-neg53.1%
unsub-neg53.1%
associate-/l*53.1%
Simplified53.1%
Taylor expanded in x around 0 53.2%
associate-/l*57.8%
Simplified57.8%
if 7.00000000000000032e-240 < b Initial program 97.7%
associate-*l/88.4%
*-commutative88.4%
+-commutative88.4%
associate--l+88.4%
exp-sum73.3%
*-commutative73.3%
exp-to-pow74.0%
sub-neg74.0%
metadata-eval74.0%
exp-diff63.1%
*-commutative63.1%
exp-to-pow63.1%
Simplified63.1%
Taylor expanded in t around 0 71.0%
times-frac66.0%
Simplified66.0%
Taylor expanded in y around 0 69.9%
Taylor expanded in b around 0 56.4%
+-commutative56.4%
associate-*r*56.4%
distribute-rgt-out56.4%
unpow256.4%
Simplified56.4%
Final simplification56.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x 2.0) (* y (* a (* b b))))))
(if (<= b -5.8e-144)
(+ (* (/ x (* y a)) (* b (* b 0.5))) (/ (- x (* x b)) (* y a)))
(if (<= b -9e-249)
t_1
(if (<= b -3e-289)
(/ x (/ y (- (/ 1.0 a) (/ b a))))
(if (<= b 1.62e-240)
t_1
(/ x (* a (+ y (* y (+ b (* (* b b) 0.5))))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -5.8e-144) {
tmp = ((x / (y * a)) * (b * (b * 0.5))) + ((x - (x * b)) / (y * a));
} else if (b <= -9e-249) {
tmp = t_1;
} else if (b <= -3e-289) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if (b <= 1.62e-240) {
tmp = t_1;
} else {
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5)))));
}
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 * 2.0d0) / (y * (a * (b * b)))
if (b <= (-5.8d-144)) then
tmp = ((x / (y * a)) * (b * (b * 0.5d0))) + ((x - (x * b)) / (y * a))
else if (b <= (-9d-249)) then
tmp = t_1
else if (b <= (-3d-289)) then
tmp = x / (y / ((1.0d0 / a) - (b / a)))
else if (b <= 1.62d-240) then
tmp = t_1
else
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5d0)))))
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 * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -5.8e-144) {
tmp = ((x / (y * a)) * (b * (b * 0.5))) + ((x - (x * b)) / (y * a));
} else if (b <= -9e-249) {
tmp = t_1;
} else if (b <= -3e-289) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if (b <= 1.62e-240) {
tmp = t_1;
} else {
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5)))));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * 2.0) / (y * (a * (b * b))) tmp = 0 if b <= -5.8e-144: tmp = ((x / (y * a)) * (b * (b * 0.5))) + ((x - (x * b)) / (y * a)) elif b <= -9e-249: tmp = t_1 elif b <= -3e-289: tmp = x / (y / ((1.0 / a) - (b / a))) elif b <= 1.62e-240: tmp = t_1 else: tmp = x / (a * (y + (y * (b + ((b * b) * 0.5))))) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * 2.0) / Float64(y * Float64(a * Float64(b * b)))) tmp = 0.0 if (b <= -5.8e-144) tmp = Float64(Float64(Float64(x / Float64(y * a)) * Float64(b * Float64(b * 0.5))) + Float64(Float64(x - Float64(x * b)) / Float64(y * a))); elseif (b <= -9e-249) tmp = t_1; elseif (b <= -3e-289) tmp = Float64(x / Float64(y / Float64(Float64(1.0 / a) - Float64(b / a)))); elseif (b <= 1.62e-240) tmp = t_1; else tmp = Float64(x / Float64(a * Float64(y + Float64(y * Float64(b + Float64(Float64(b * b) * 0.5)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * 2.0) / (y * (a * (b * b))); tmp = 0.0; if (b <= -5.8e-144) tmp = ((x / (y * a)) * (b * (b * 0.5))) + ((x - (x * b)) / (y * a)); elseif (b <= -9e-249) tmp = t_1; elseif (b <= -3e-289) tmp = x / (y / ((1.0 / a) - (b / a))); elseif (b <= 1.62e-240) tmp = t_1; else tmp = x / (a * (y + (y * (b + ((b * b) * 0.5))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -5.8e-144], N[(N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] * N[(b * N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(x - N[(x * b), $MachinePrecision]), $MachinePrecision] / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -9e-249], t$95$1, If[LessEqual[b, -3e-289], N[(x / N[(y / N[(N[(1.0 / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.62e-240], t$95$1, N[(x / N[(a * N[(y + N[(y * N[(b + N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot 2}{y \cdot \left(a \cdot \left(b \cdot b\right)\right)}\\
\mathbf{if}\;b \leq -5.8 \cdot 10^{-144}:\\
\;\;\;\;\frac{x}{y \cdot a} \cdot \left(b \cdot \left(b \cdot 0.5\right)\right) + \frac{x - x \cdot b}{y \cdot a}\\
\mathbf{elif}\;b \leq -9 \cdot 10^{-249}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;b \leq -3 \cdot 10^{-289}:\\
\;\;\;\;\frac{x}{\frac{y}{\frac{1}{a} - \frac{b}{a}}}\\
\mathbf{elif}\;b \leq 1.62 \cdot 10^{-240}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot \left(b + \left(b \cdot b\right) \cdot 0.5\right)\right)}\\
\end{array}
\end{array}
if b < -5.8000000000000004e-144Initial program 98.1%
associate-*l/89.5%
*-commutative89.5%
+-commutative89.5%
associate--l+89.5%
exp-sum67.9%
*-commutative67.9%
exp-to-pow68.6%
sub-neg68.6%
metadata-eval68.6%
exp-diff58.9%
*-commutative58.9%
exp-to-pow58.9%
Simplified58.9%
Taylor expanded in t around 0 66.6%
times-frac64.1%
Simplified64.1%
Taylor expanded in y around 0 63.6%
Taylor expanded in b around 0 28.1%
+-commutative28.1%
associate-*r*28.1%
distribute-rgt-out36.7%
unpow236.7%
Simplified36.7%
Taylor expanded in b around 0 36.1%
Simplified57.9%
if -5.8000000000000004e-144 < b < -8.99999999999999962e-249 or -2.9999999999999998e-289 < b < 1.61999999999999995e-240Initial program 96.7%
associate-*l/90.7%
*-commutative90.7%
+-commutative90.7%
associate--l+90.7%
exp-sum78.9%
*-commutative78.9%
exp-to-pow79.5%
sub-neg79.5%
metadata-eval79.5%
exp-diff79.5%
*-commutative79.5%
exp-to-pow79.5%
Simplified79.5%
Taylor expanded in t around 0 62.4%
times-frac65.3%
Simplified65.3%
Taylor expanded in y around 0 31.3%
Taylor expanded in b around 0 31.3%
+-commutative31.3%
associate-*r*31.3%
distribute-rgt-out31.3%
unpow231.3%
Simplified31.3%
Taylor expanded in b around inf 59.7%
associate-*r/59.7%
unpow259.7%
associate-*r*59.7%
Simplified59.7%
if -8.99999999999999962e-249 < b < -2.9999999999999998e-289Initial program 97.9%
Taylor expanded in t around 0 78.5%
+-commutative78.5%
mul-1-neg78.5%
unsub-neg78.5%
Simplified78.5%
Taylor expanded in y around 0 51.6%
exp-neg51.6%
associate-*r/51.6%
*-rgt-identity51.6%
+-commutative51.6%
exp-sum51.6%
rem-exp-log53.1%
Simplified53.1%
Taylor expanded in b around 0 53.1%
+-commutative53.1%
mul-1-neg53.1%
unsub-neg53.1%
associate-/l*53.1%
Simplified53.1%
Taylor expanded in x around 0 53.2%
associate-/l*57.8%
Simplified57.8%
if 1.61999999999999995e-240 < b Initial program 97.7%
associate-*l/88.4%
*-commutative88.4%
+-commutative88.4%
associate--l+88.4%
exp-sum73.3%
*-commutative73.3%
exp-to-pow74.0%
sub-neg74.0%
metadata-eval74.0%
exp-diff63.1%
*-commutative63.1%
exp-to-pow63.1%
Simplified63.1%
Taylor expanded in t around 0 71.0%
times-frac66.0%
Simplified66.0%
Taylor expanded in y around 0 69.9%
Taylor expanded in b around 0 56.4%
+-commutative56.4%
associate-*r*56.4%
distribute-rgt-out56.4%
unpow256.4%
Simplified56.4%
Final simplification57.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -6.8e-115)
(/ (+ (- (/ x a) (/ b (/ a x))) (* (* b b) (* (/ x a) 0.5))) y)
(if (<= b -1e-248)
(/ (* x 2.0) (* y (* a (* b b))))
(if (<= b 1e-149)
(/ (* x (/ 1.0 a)) y)
(/ x (* a (+ y (* y (+ b (* (* b b) 0.5))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -6.8e-115) {
tmp = (((x / a) - (b / (a / x))) + ((b * b) * ((x / a) * 0.5))) / y;
} else if (b <= -1e-248) {
tmp = (x * 2.0) / (y * (a * (b * b)));
} else if (b <= 1e-149) {
tmp = (x * (1.0 / a)) / y;
} else {
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5)))));
}
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 <= (-6.8d-115)) then
tmp = (((x / a) - (b / (a / x))) + ((b * b) * ((x / a) * 0.5d0))) / y
else if (b <= (-1d-248)) then
tmp = (x * 2.0d0) / (y * (a * (b * b)))
else if (b <= 1d-149) then
tmp = (x * (1.0d0 / a)) / y
else
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5d0)))))
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 <= -6.8e-115) {
tmp = (((x / a) - (b / (a / x))) + ((b * b) * ((x / a) * 0.5))) / y;
} else if (b <= -1e-248) {
tmp = (x * 2.0) / (y * (a * (b * b)));
} else if (b <= 1e-149) {
tmp = (x * (1.0 / a)) / y;
} else {
tmp = x / (a * (y + (y * (b + ((b * b) * 0.5)))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -6.8e-115: tmp = (((x / a) - (b / (a / x))) + ((b * b) * ((x / a) * 0.5))) / y elif b <= -1e-248: tmp = (x * 2.0) / (y * (a * (b * b))) elif b <= 1e-149: tmp = (x * (1.0 / a)) / y else: tmp = x / (a * (y + (y * (b + ((b * b) * 0.5))))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -6.8e-115) tmp = Float64(Float64(Float64(Float64(x / a) - Float64(b / Float64(a / x))) + Float64(Float64(b * b) * Float64(Float64(x / a) * 0.5))) / y); elseif (b <= -1e-248) tmp = Float64(Float64(x * 2.0) / Float64(y * Float64(a * Float64(b * b)))); elseif (b <= 1e-149) tmp = Float64(Float64(x * Float64(1.0 / a)) / y); else tmp = Float64(x / Float64(a * Float64(y + Float64(y * Float64(b + Float64(Float64(b * b) * 0.5)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -6.8e-115) tmp = (((x / a) - (b / (a / x))) + ((b * b) * ((x / a) * 0.5))) / y; elseif (b <= -1e-248) tmp = (x * 2.0) / (y * (a * (b * b))); elseif (b <= 1e-149) tmp = (x * (1.0 / a)) / y; else tmp = x / (a * (y + (y * (b + ((b * b) * 0.5))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -6.8e-115], N[(N[(N[(N[(x / a), $MachinePrecision] - N[(b / N[(a / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(b * b), $MachinePrecision] * N[(N[(x / a), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -1e-248], N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1e-149], N[(N[(x * N[(1.0 / a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(x / N[(a * N[(y + N[(y * N[(b + N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -6.8 \cdot 10^{-115}:\\
\;\;\;\;\frac{\left(\frac{x}{a} - \frac{b}{\frac{a}{x}}\right) + \left(b \cdot b\right) \cdot \left(\frac{x}{a} \cdot 0.5\right)}{y}\\
\mathbf{elif}\;b \leq -1 \cdot 10^{-248}:\\
\;\;\;\;\frac{x \cdot 2}{y \cdot \left(a \cdot \left(b \cdot b\right)\right)}\\
\mathbf{elif}\;b \leq 10^{-149}:\\
\;\;\;\;\frac{x \cdot \frac{1}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot \left(b + \left(b \cdot b\right) \cdot 0.5\right)\right)}\\
\end{array}
\end{array}
if b < -6.7999999999999996e-115Initial program 99.3%
Taylor expanded in t around 0 84.9%
+-commutative84.9%
mul-1-neg84.9%
unsub-neg84.9%
Simplified84.9%
Taylor expanded in y around 0 67.6%
exp-neg67.6%
associate-*r/67.6%
*-rgt-identity67.6%
+-commutative67.6%
exp-sum67.7%
rem-exp-log68.2%
Simplified68.2%
Taylor expanded in b around 0 49.2%
+-commutative49.2%
+-commutative49.2%
mul-1-neg49.2%
unsub-neg49.2%
associate-/l*49.2%
mul-1-neg49.2%
distribute-rgt-neg-in49.2%
distribute-rgt-out65.4%
metadata-eval65.4%
*-commutative65.4%
distribute-lft-neg-in65.4%
metadata-eval65.4%
unpow265.4%
Simplified65.4%
if -6.7999999999999996e-115 < b < -9.9999999999999998e-249Initial program 91.9%
associate-*l/86.7%
*-commutative86.7%
+-commutative86.7%
associate--l+86.7%
exp-sum70.7%
*-commutative70.7%
exp-to-pow71.8%
sub-neg71.8%
metadata-eval71.8%
exp-diff71.8%
*-commutative71.8%
exp-to-pow71.8%
Simplified71.8%
Taylor expanded in t around 0 68.5%
times-frac60.9%
Simplified60.9%
Taylor expanded in y around 0 33.8%
Taylor expanded in b around 0 33.8%
+-commutative33.8%
associate-*r*33.8%
distribute-rgt-out33.8%
unpow233.8%
Simplified33.8%
Taylor expanded in b around inf 49.1%
associate-*r/49.1%
unpow249.1%
associate-*r*49.1%
Simplified49.1%
if -9.9999999999999998e-249 < b < 9.99999999999999979e-150Initial program 97.9%
Taylor expanded in t around 0 79.9%
+-commutative79.9%
mul-1-neg79.9%
unsub-neg79.9%
Simplified79.9%
Taylor expanded in b around 0 79.9%
div-exp79.9%
*-commutative79.9%
exp-to-pow79.9%
rem-exp-log81.9%
Simplified81.9%
Taylor expanded in y around 0 52.2%
if 9.99999999999999979e-150 < b Initial program 97.9%
associate-*l/87.9%
*-commutative87.9%
+-commutative87.9%
associate--l+87.9%
exp-sum71.5%
*-commutative71.5%
exp-to-pow71.9%
sub-neg71.9%
metadata-eval71.9%
exp-diff59.3%
*-commutative59.3%
exp-to-pow59.3%
Simplified59.3%
Taylor expanded in t around 0 71.1%
times-frac63.6%
Simplified63.6%
Taylor expanded in y around 0 72.6%
Taylor expanded in b around 0 57.0%
+-commutative57.0%
associate-*r*57.0%
distribute-rgt-out57.0%
unpow257.0%
Simplified57.0%
Final simplification57.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x 2.0) (* y (* a (* b b)))))
(t_2 (/ x (/ y (- (/ 1.0 a) (/ b a))))))
(if (<= b -1.8e-144)
t_2
(if (<= b -3.5e-248)
t_1
(if (<= b -5.2e-289)
t_2
(if (or (<= b 9.6e-241) (not (<= b 0.82)))
t_1
(/ (/ x (+ a (* a b))) y)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * 2.0) / (y * (a * (b * b)));
double t_2 = x / (y / ((1.0 / a) - (b / a)));
double tmp;
if (b <= -1.8e-144) {
tmp = t_2;
} else if (b <= -3.5e-248) {
tmp = t_1;
} else if (b <= -5.2e-289) {
tmp = t_2;
} else if ((b <= 9.6e-241) || !(b <= 0.82)) {
tmp = t_1;
} else {
tmp = (x / (a + (a * 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) :: t_2
real(8) :: tmp
t_1 = (x * 2.0d0) / (y * (a * (b * b)))
t_2 = x / (y / ((1.0d0 / a) - (b / a)))
if (b <= (-1.8d-144)) then
tmp = t_2
else if (b <= (-3.5d-248)) then
tmp = t_1
else if (b <= (-5.2d-289)) then
tmp = t_2
else if ((b <= 9.6d-241) .or. (.not. (b <= 0.82d0))) then
tmp = t_1
else
tmp = (x / (a + (a * 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 * 2.0) / (y * (a * (b * b)));
double t_2 = x / (y / ((1.0 / a) - (b / a)));
double tmp;
if (b <= -1.8e-144) {
tmp = t_2;
} else if (b <= -3.5e-248) {
tmp = t_1;
} else if (b <= -5.2e-289) {
tmp = t_2;
} else if ((b <= 9.6e-241) || !(b <= 0.82)) {
tmp = t_1;
} else {
tmp = (x / (a + (a * b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * 2.0) / (y * (a * (b * b))) t_2 = x / (y / ((1.0 / a) - (b / a))) tmp = 0 if b <= -1.8e-144: tmp = t_2 elif b <= -3.5e-248: tmp = t_1 elif b <= -5.2e-289: tmp = t_2 elif (b <= 9.6e-241) or not (b <= 0.82): tmp = t_1 else: tmp = (x / (a + (a * b))) / y return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * 2.0) / Float64(y * Float64(a * Float64(b * b)))) t_2 = Float64(x / Float64(y / Float64(Float64(1.0 / a) - Float64(b / a)))) tmp = 0.0 if (b <= -1.8e-144) tmp = t_2; elseif (b <= -3.5e-248) tmp = t_1; elseif (b <= -5.2e-289) tmp = t_2; elseif ((b <= 9.6e-241) || !(b <= 0.82)) tmp = t_1; else tmp = Float64(Float64(x / Float64(a + Float64(a * b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * 2.0) / (y * (a * (b * b))); t_2 = x / (y / ((1.0 / a) - (b / a))); tmp = 0.0; if (b <= -1.8e-144) tmp = t_2; elseif (b <= -3.5e-248) tmp = t_1; elseif (b <= -5.2e-289) tmp = t_2; elseif ((b <= 9.6e-241) || ~((b <= 0.82))) tmp = t_1; else tmp = (x / (a + (a * b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x / N[(y / N[(N[(1.0 / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.8e-144], t$95$2, If[LessEqual[b, -3.5e-248], t$95$1, If[LessEqual[b, -5.2e-289], t$95$2, If[Or[LessEqual[b, 9.6e-241], N[Not[LessEqual[b, 0.82]], $MachinePrecision]], t$95$1, N[(N[(x / N[(a + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot 2}{y \cdot \left(a \cdot \left(b \cdot b\right)\right)}\\
t_2 := \frac{x}{\frac{y}{\frac{1}{a} - \frac{b}{a}}}\\
\mathbf{if}\;b \leq -1.8 \cdot 10^{-144}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;b \leq -3.5 \cdot 10^{-248}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;b \leq -5.2 \cdot 10^{-289}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;b \leq 9.6 \cdot 10^{-241} \lor \neg \left(b \leq 0.82\right):\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + a \cdot b}}{y}\\
\end{array}
\end{array}
if b < -1.8e-144 or -3.49999999999999983e-248 < b < -5.1999999999999998e-289Initial program 98.1%
Taylor expanded in t around 0 83.9%
+-commutative83.9%
mul-1-neg83.9%
unsub-neg83.9%
Simplified83.9%
Taylor expanded in y around 0 60.6%
exp-neg60.6%
associate-*r/60.6%
*-rgt-identity60.6%
+-commutative60.6%
exp-sum60.7%
rem-exp-log61.5%
Simplified61.5%
Taylor expanded in b around 0 53.3%
+-commutative53.3%
mul-1-neg53.3%
unsub-neg53.3%
associate-/l*51.5%
Simplified51.5%
Taylor expanded in x around 0 52.4%
associate-/l*55.1%
Simplified55.1%
if -1.8e-144 < b < -3.49999999999999983e-248 or -5.1999999999999998e-289 < b < 9.6e-241 or 0.819999999999999951 < b Initial program 99.0%
associate-*l/88.7%
*-commutative88.7%
+-commutative88.7%
associate--l+88.7%
exp-sum70.2%
*-commutative70.2%
exp-to-pow70.4%
sub-neg70.4%
metadata-eval70.4%
exp-diff58.4%
*-commutative58.4%
exp-to-pow58.4%
Simplified58.4%
Taylor expanded in t around 0 66.0%
times-frac62.3%
Simplified62.3%
Taylor expanded in y around 0 66.5%
Taylor expanded in b around 0 52.2%
+-commutative52.2%
associate-*r*52.2%
distribute-rgt-out52.2%
unpow252.2%
Simplified52.2%
Taylor expanded in b around inf 61.2%
associate-*r/61.2%
unpow261.2%
associate-*r*60.3%
Simplified60.3%
if 9.6e-241 < b < 0.819999999999999951Initial program 94.0%
Taylor expanded in t around 0 74.7%
+-commutative74.7%
mul-1-neg74.7%
unsub-neg74.7%
Simplified74.7%
Taylor expanded in y around 0 49.1%
exp-neg49.1%
associate-*r/49.1%
*-rgt-identity49.1%
+-commutative49.1%
exp-sum49.2%
rem-exp-log51.0%
Simplified51.0%
Taylor expanded in b around 0 49.5%
Final simplification56.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x 2.0) (* y (* a (* b b))))))
(if (<= b -2.05e-146)
(- (/ x (* y a)) (* (/ x a) (/ b y)))
(if (<= b -9e-248)
t_1
(if (<= b -6.2e-288)
(/ x (/ y (- (/ 1.0 a) (/ b a))))
(if (or (<= b 7.4e-240) (not (<= b 0.68)))
t_1
(/ (/ x (+ a (* a b))) y)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -2.05e-146) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= -9e-248) {
tmp = t_1;
} else if (b <= -6.2e-288) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if ((b <= 7.4e-240) || !(b <= 0.68)) {
tmp = t_1;
} else {
tmp = (x / (a + (a * 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 * 2.0d0) / (y * (a * (b * b)))
if (b <= (-2.05d-146)) then
tmp = (x / (y * a)) - ((x / a) * (b / y))
else if (b <= (-9d-248)) then
tmp = t_1
else if (b <= (-6.2d-288)) then
tmp = x / (y / ((1.0d0 / a) - (b / a)))
else if ((b <= 7.4d-240) .or. (.not. (b <= 0.68d0))) then
tmp = t_1
else
tmp = (x / (a + (a * 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 * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -2.05e-146) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= -9e-248) {
tmp = t_1;
} else if (b <= -6.2e-288) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if ((b <= 7.4e-240) || !(b <= 0.68)) {
tmp = t_1;
} else {
tmp = (x / (a + (a * b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * 2.0) / (y * (a * (b * b))) tmp = 0 if b <= -2.05e-146: tmp = (x / (y * a)) - ((x / a) * (b / y)) elif b <= -9e-248: tmp = t_1 elif b <= -6.2e-288: tmp = x / (y / ((1.0 / a) - (b / a))) elif (b <= 7.4e-240) or not (b <= 0.68): tmp = t_1 else: tmp = (x / (a + (a * b))) / y return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * 2.0) / Float64(y * Float64(a * Float64(b * b)))) tmp = 0.0 if (b <= -2.05e-146) tmp = Float64(Float64(x / Float64(y * a)) - Float64(Float64(x / a) * Float64(b / y))); elseif (b <= -9e-248) tmp = t_1; elseif (b <= -6.2e-288) tmp = Float64(x / Float64(y / Float64(Float64(1.0 / a) - Float64(b / a)))); elseif ((b <= 7.4e-240) || !(b <= 0.68)) tmp = t_1; else tmp = Float64(Float64(x / Float64(a + Float64(a * b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * 2.0) / (y * (a * (b * b))); tmp = 0.0; if (b <= -2.05e-146) tmp = (x / (y * a)) - ((x / a) * (b / y)); elseif (b <= -9e-248) tmp = t_1; elseif (b <= -6.2e-288) tmp = x / (y / ((1.0 / a) - (b / a))); elseif ((b <= 7.4e-240) || ~((b <= 0.68))) tmp = t_1; else tmp = (x / (a + (a * b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2.05e-146], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] - N[(N[(x / a), $MachinePrecision] * N[(b / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -9e-248], t$95$1, If[LessEqual[b, -6.2e-288], N[(x / N[(y / N[(N[(1.0 / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[b, 7.4e-240], N[Not[LessEqual[b, 0.68]], $MachinePrecision]], t$95$1, N[(N[(x / N[(a + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot 2}{y \cdot \left(a \cdot \left(b \cdot b\right)\right)}\\
\mathbf{if}\;b \leq -2.05 \cdot 10^{-146}:\\
\;\;\;\;\frac{x}{y \cdot a} - \frac{x}{a} \cdot \frac{b}{y}\\
\mathbf{elif}\;b \leq -9 \cdot 10^{-248}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;b \leq -6.2 \cdot 10^{-288}:\\
\;\;\;\;\frac{x}{\frac{y}{\frac{1}{a} - \frac{b}{a}}}\\
\mathbf{elif}\;b \leq 7.4 \cdot 10^{-240} \lor \neg \left(b \leq 0.68\right):\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + a \cdot b}}{y}\\
\end{array}
\end{array}
if b < -2.0499999999999999e-146Initial program 98.1%
associate-*l/89.5%
*-commutative89.5%
+-commutative89.5%
associate--l+89.5%
exp-sum67.9%
*-commutative67.9%
exp-to-pow68.6%
sub-neg68.6%
metadata-eval68.6%
exp-diff58.9%
*-commutative58.9%
exp-to-pow58.9%
Simplified58.9%
Taylor expanded in t around 0 66.6%
times-frac64.1%
Simplified64.1%
Taylor expanded in y around 0 63.6%
Taylor expanded in b around 0 52.2%
+-commutative52.2%
mul-1-neg52.2%
unsub-neg52.2%
*-commutative52.2%
*-commutative52.2%
times-frac55.7%
Simplified55.7%
if -2.0499999999999999e-146 < b < -8.9999999999999992e-248 or -6.19999999999999967e-288 < b < 7.4000000000000003e-240 or 0.680000000000000049 < b Initial program 99.0%
associate-*l/88.7%
*-commutative88.7%
+-commutative88.7%
associate--l+88.7%
exp-sum70.2%
*-commutative70.2%
exp-to-pow70.4%
sub-neg70.4%
metadata-eval70.4%
exp-diff58.4%
*-commutative58.4%
exp-to-pow58.4%
Simplified58.4%
Taylor expanded in t around 0 66.0%
times-frac62.3%
Simplified62.3%
Taylor expanded in y around 0 66.5%
Taylor expanded in b around 0 52.2%
+-commutative52.2%
associate-*r*52.2%
distribute-rgt-out52.2%
unpow252.2%
Simplified52.2%
Taylor expanded in b around inf 61.2%
associate-*r/61.2%
unpow261.2%
associate-*r*60.3%
Simplified60.3%
if -8.9999999999999992e-248 < b < -6.19999999999999967e-288Initial program 97.9%
Taylor expanded in t around 0 78.5%
+-commutative78.5%
mul-1-neg78.5%
unsub-neg78.5%
Simplified78.5%
Taylor expanded in y around 0 51.6%
exp-neg51.6%
associate-*r/51.6%
*-rgt-identity51.6%
+-commutative51.6%
exp-sum51.6%
rem-exp-log53.1%
Simplified53.1%
Taylor expanded in b around 0 53.1%
+-commutative53.1%
mul-1-neg53.1%
unsub-neg53.1%
associate-/l*53.1%
Simplified53.1%
Taylor expanded in x around 0 53.2%
associate-/l*57.8%
Simplified57.8%
if 7.4000000000000003e-240 < b < 0.680000000000000049Initial program 94.0%
Taylor expanded in t around 0 74.7%
+-commutative74.7%
mul-1-neg74.7%
unsub-neg74.7%
Simplified74.7%
Taylor expanded in y around 0 49.1%
exp-neg49.1%
associate-*r/49.1%
*-rgt-identity49.1%
+-commutative49.1%
exp-sum49.2%
rem-exp-log51.0%
Simplified51.0%
Taylor expanded in b around 0 49.5%
Final simplification56.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x 2.0) (* y (* a (* b b))))))
(if (<= b -2.35e-144)
(- (/ x (* y a)) (* (/ x a) (/ b y)))
(if (<= b -9e-249)
t_1
(if (<= b -2e-289)
(/ x (/ y (- (/ 1.0 a) (/ b a))))
(if (<= b 1.7e-241) t_1 (/ x (* a (+ y (* y (* b (* b 0.5))))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -2.35e-144) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= -9e-249) {
tmp = t_1;
} else if (b <= -2e-289) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if (b <= 1.7e-241) {
tmp = t_1;
} else {
tmp = x / (a * (y + (y * (b * (b * 0.5)))));
}
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 * 2.0d0) / (y * (a * (b * b)))
if (b <= (-2.35d-144)) then
tmp = (x / (y * a)) - ((x / a) * (b / y))
else if (b <= (-9d-249)) then
tmp = t_1
else if (b <= (-2d-289)) then
tmp = x / (y / ((1.0d0 / a) - (b / a)))
else if (b <= 1.7d-241) then
tmp = t_1
else
tmp = x / (a * (y + (y * (b * (b * 0.5d0)))))
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 * 2.0) / (y * (a * (b * b)));
double tmp;
if (b <= -2.35e-144) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= -9e-249) {
tmp = t_1;
} else if (b <= -2e-289) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else if (b <= 1.7e-241) {
tmp = t_1;
} else {
tmp = x / (a * (y + (y * (b * (b * 0.5)))));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * 2.0) / (y * (a * (b * b))) tmp = 0 if b <= -2.35e-144: tmp = (x / (y * a)) - ((x / a) * (b / y)) elif b <= -9e-249: tmp = t_1 elif b <= -2e-289: tmp = x / (y / ((1.0 / a) - (b / a))) elif b <= 1.7e-241: tmp = t_1 else: tmp = x / (a * (y + (y * (b * (b * 0.5))))) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * 2.0) / Float64(y * Float64(a * Float64(b * b)))) tmp = 0.0 if (b <= -2.35e-144) tmp = Float64(Float64(x / Float64(y * a)) - Float64(Float64(x / a) * Float64(b / y))); elseif (b <= -9e-249) tmp = t_1; elseif (b <= -2e-289) tmp = Float64(x / Float64(y / Float64(Float64(1.0 / a) - Float64(b / a)))); elseif (b <= 1.7e-241) tmp = t_1; else tmp = Float64(x / Float64(a * Float64(y + Float64(y * Float64(b * Float64(b * 0.5)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * 2.0) / (y * (a * (b * b))); tmp = 0.0; if (b <= -2.35e-144) tmp = (x / (y * a)) - ((x / a) * (b / y)); elseif (b <= -9e-249) tmp = t_1; elseif (b <= -2e-289) tmp = x / (y / ((1.0 / a) - (b / a))); elseif (b <= 1.7e-241) tmp = t_1; else tmp = x / (a * (y + (y * (b * (b * 0.5))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2.35e-144], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] - N[(N[(x / a), $MachinePrecision] * N[(b / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -9e-249], t$95$1, If[LessEqual[b, -2e-289], N[(x / N[(y / N[(N[(1.0 / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.7e-241], t$95$1, N[(x / N[(a * N[(y + N[(y * N[(b * N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot 2}{y \cdot \left(a \cdot \left(b \cdot b\right)\right)}\\
\mathbf{if}\;b \leq -2.35 \cdot 10^{-144}:\\
\;\;\;\;\frac{x}{y \cdot a} - \frac{x}{a} \cdot \frac{b}{y}\\
\mathbf{elif}\;b \leq -9 \cdot 10^{-249}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-289}:\\
\;\;\;\;\frac{x}{\frac{y}{\frac{1}{a} - \frac{b}{a}}}\\
\mathbf{elif}\;b \leq 1.7 \cdot 10^{-241}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot \left(b \cdot \left(b \cdot 0.5\right)\right)\right)}\\
\end{array}
\end{array}
if b < -2.3500000000000001e-144Initial program 98.1%
associate-*l/89.5%
*-commutative89.5%
+-commutative89.5%
associate--l+89.5%
exp-sum67.9%
*-commutative67.9%
exp-to-pow68.6%
sub-neg68.6%
metadata-eval68.6%
exp-diff58.9%
*-commutative58.9%
exp-to-pow58.9%
Simplified58.9%
Taylor expanded in t around 0 66.6%
times-frac64.1%
Simplified64.1%
Taylor expanded in y around 0 63.6%
Taylor expanded in b around 0 52.2%
+-commutative52.2%
mul-1-neg52.2%
unsub-neg52.2%
*-commutative52.2%
*-commutative52.2%
times-frac55.7%
Simplified55.7%
if -2.3500000000000001e-144 < b < -8.99999999999999962e-249 or -2e-289 < b < 1.6999999999999999e-241Initial program 96.7%
associate-*l/90.7%
*-commutative90.7%
+-commutative90.7%
associate--l+90.7%
exp-sum78.9%
*-commutative78.9%
exp-to-pow79.5%
sub-neg79.5%
metadata-eval79.5%
exp-diff79.5%
*-commutative79.5%
exp-to-pow79.5%
Simplified79.5%
Taylor expanded in t around 0 62.4%
times-frac65.3%
Simplified65.3%
Taylor expanded in y around 0 31.3%
Taylor expanded in b around 0 31.3%
+-commutative31.3%
associate-*r*31.3%
distribute-rgt-out31.3%
unpow231.3%
Simplified31.3%
Taylor expanded in b around inf 59.7%
associate-*r/59.7%
unpow259.7%
associate-*r*59.7%
Simplified59.7%
if -8.99999999999999962e-249 < b < -2e-289Initial program 97.9%
Taylor expanded in t around 0 78.5%
+-commutative78.5%
mul-1-neg78.5%
unsub-neg78.5%
Simplified78.5%
Taylor expanded in y around 0 51.6%
exp-neg51.6%
associate-*r/51.6%
*-rgt-identity51.6%
+-commutative51.6%
exp-sum51.6%
rem-exp-log53.1%
Simplified53.1%
Taylor expanded in b around 0 53.1%
+-commutative53.1%
mul-1-neg53.1%
unsub-neg53.1%
associate-/l*53.1%
Simplified53.1%
Taylor expanded in x around 0 53.2%
associate-/l*57.8%
Simplified57.8%
if 1.6999999999999999e-241 < b Initial program 97.7%
associate-*l/88.4%
*-commutative88.4%
+-commutative88.4%
associate--l+88.4%
exp-sum73.3%
*-commutative73.3%
exp-to-pow74.0%
sub-neg74.0%
metadata-eval74.0%
exp-diff63.1%
*-commutative63.1%
exp-to-pow63.1%
Simplified63.1%
Taylor expanded in t around 0 71.0%
times-frac66.0%
Simplified66.0%
Taylor expanded in y around 0 69.9%
Taylor expanded in b around 0 56.4%
+-commutative56.4%
associate-*r*56.4%
distribute-rgt-out56.4%
unpow256.4%
Simplified56.4%
Taylor expanded in b around inf 55.8%
unpow255.8%
associate-*r*55.8%
*-commutative55.8%
*-commutative55.8%
associate-*l*55.8%
*-commutative55.8%
Simplified55.8%
Final simplification56.4%
(FPCore (x y z t a b) :precision binary64 (if (<= b 30.0) (/ (- x (* x b)) (* y a)) (* (/ x y) (/ 2.0 (* b (* a b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 30.0) {
tmp = (x - (x * b)) / (y * a);
} else {
tmp = (x / y) * (2.0 / (b * (a * 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 (b <= 30.0d0) then
tmp = (x - (x * b)) / (y * a)
else
tmp = (x / y) * (2.0d0 / (b * (a * 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 (b <= 30.0) {
tmp = (x - (x * b)) / (y * a);
} else {
tmp = (x / y) * (2.0 / (b * (a * b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 30.0: tmp = (x - (x * b)) / (y * a) else: tmp = (x / y) * (2.0 / (b * (a * b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 30.0) tmp = Float64(Float64(x - Float64(x * b)) / Float64(y * a)); else tmp = Float64(Float64(x / y) * Float64(2.0 / Float64(b * Float64(a * b)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 30.0) tmp = (x - (x * b)) / (y * a); else tmp = (x / y) * (2.0 / (b * (a * b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 30.0], N[(N[(x - N[(x * b), $MachinePrecision]), $MachinePrecision] / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(N[(x / y), $MachinePrecision] * N[(2.0 / N[(b * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 30:\\
\;\;\;\;\frac{x - x \cdot b}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y} \cdot \frac{2}{b \cdot \left(a \cdot b\right)}\\
\end{array}
\end{array}
if b < 30Initial program 96.8%
Taylor expanded in t around 0 79.4%
+-commutative79.4%
mul-1-neg79.4%
unsub-neg79.4%
Simplified79.4%
Taylor expanded in y around 0 51.7%
exp-neg51.7%
associate-*r/51.7%
*-rgt-identity51.7%
+-commutative51.7%
exp-sum51.7%
rem-exp-log52.7%
Simplified52.7%
Taylor expanded in b around 0 47.2%
+-commutative47.2%
mul-1-neg47.2%
unsub-neg47.2%
associate-/l*45.1%
Simplified45.1%
Taylor expanded in x around 0 47.2%
*-commutative47.2%
sub-neg47.2%
mul-1-neg47.2%
+-commutative47.2%
*-commutative47.2%
+-commutative47.2%
distribute-rgt-in46.7%
associate-*l/46.7%
*-lft-identity46.7%
mul-1-neg46.7%
cancel-sign-sub-inv46.7%
associate-*l/47.2%
div-sub47.8%
associate-/r*48.1%
*-commutative48.1%
*-commutative48.1%
Simplified48.1%
if 30 < b Initial program 100.0%
associate-*l/87.5%
*-commutative87.5%
+-commutative87.5%
associate--l+87.5%
exp-sum65.3%
*-commutative65.3%
exp-to-pow65.3%
sub-neg65.3%
metadata-eval65.3%
exp-diff47.2%
*-commutative47.2%
exp-to-pow47.2%
Simplified47.2%
Taylor expanded in t around 0 66.8%
times-frac59.8%
Simplified59.8%
Taylor expanded in y around 0 82.2%
Taylor expanded in b around 0 60.8%
+-commutative60.8%
associate-*r*60.8%
distribute-rgt-out60.8%
unpow260.8%
Simplified60.8%
Taylor expanded in b around inf 60.8%
associate-*r/60.8%
unpow260.8%
associate-*r*59.4%
times-frac53.9%
associate-*r*51.3%
Simplified51.3%
Final simplification49.0%
(FPCore (x y z t a b) :precision binary64 (if (<= b 0.95) (/ (- x (* x b)) (* y a)) (/ x (* (* b (* y b)) (* a 0.5)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 0.95) {
tmp = (x - (x * b)) / (y * a);
} else {
tmp = x / ((b * (y * b)) * (a * 0.5));
}
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 <= 0.95d0) then
tmp = (x - (x * b)) / (y * a)
else
tmp = x / ((b * (y * b)) * (a * 0.5d0))
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 <= 0.95) {
tmp = (x - (x * b)) / (y * a);
} else {
tmp = x / ((b * (y * b)) * (a * 0.5));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 0.95: tmp = (x - (x * b)) / (y * a) else: tmp = x / ((b * (y * b)) * (a * 0.5)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 0.95) tmp = Float64(Float64(x - Float64(x * b)) / Float64(y * a)); else tmp = Float64(x / Float64(Float64(b * Float64(y * b)) * Float64(a * 0.5))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 0.95) tmp = (x - (x * b)) / (y * a); else tmp = x / ((b * (y * b)) * (a * 0.5)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 0.95], N[(N[(x - N[(x * b), $MachinePrecision]), $MachinePrecision] / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(x / N[(N[(b * N[(y * b), $MachinePrecision]), $MachinePrecision] * N[(a * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 0.95:\\
\;\;\;\;\frac{x - x \cdot b}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\left(b \cdot \left(y \cdot b\right)\right) \cdot \left(a \cdot 0.5\right)}\\
\end{array}
\end{array}
if b < 0.94999999999999996Initial program 96.8%
Taylor expanded in t around 0 79.2%
+-commutative79.2%
mul-1-neg79.2%
unsub-neg79.2%
Simplified79.2%
Taylor expanded in y around 0 51.7%
exp-neg51.7%
associate-*r/51.7%
*-rgt-identity51.7%
+-commutative51.7%
exp-sum51.7%
rem-exp-log52.7%
Simplified52.7%
Taylor expanded in b around 0 47.2%
+-commutative47.2%
mul-1-neg47.2%
unsub-neg47.2%
associate-/l*45.0%
Simplified45.0%
Taylor expanded in x around 0 47.2%
*-commutative47.2%
sub-neg47.2%
mul-1-neg47.2%
+-commutative47.2%
*-commutative47.2%
+-commutative47.2%
distribute-rgt-in46.6%
associate-*l/46.6%
*-lft-identity46.6%
mul-1-neg46.6%
cancel-sign-sub-inv46.6%
associate-*l/47.2%
div-sub47.7%
associate-/r*47.5%
*-commutative47.5%
*-commutative47.5%
Simplified47.5%
if 0.94999999999999996 < b Initial program 100.0%
associate-*l/87.8%
*-commutative87.8%
+-commutative87.8%
associate--l+87.8%
exp-sum66.2%
*-commutative66.2%
exp-to-pow66.2%
sub-neg66.2%
metadata-eval66.2%
exp-diff48.6%
*-commutative48.6%
exp-to-pow48.6%
Simplified48.6%
Taylor expanded in t around 0 67.7%
times-frac60.9%
Simplified60.9%
Taylor expanded in y around 0 82.7%
Taylor expanded in b around 0 61.8%
+-commutative61.8%
associate-*r*61.8%
distribute-rgt-out61.8%
unpow261.8%
Simplified61.8%
Taylor expanded in b around inf 61.8%
*-commutative61.8%
*-commutative61.8%
associate-*l*61.8%
unpow261.8%
associate-*l*56.6%
*-commutative56.6%
Simplified56.6%
Final simplification50.1%
(FPCore (x y z t a b) :precision binary64 (if (<= b 0.96) (/ x (/ y (- (/ 1.0 a) (/ b a)))) (/ x (* (* b (* y b)) (* a 0.5)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 0.96) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else {
tmp = x / ((b * (y * b)) * (a * 0.5));
}
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 <= 0.96d0) then
tmp = x / (y / ((1.0d0 / a) - (b / a)))
else
tmp = x / ((b * (y * b)) * (a * 0.5d0))
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 <= 0.96) {
tmp = x / (y / ((1.0 / a) - (b / a)));
} else {
tmp = x / ((b * (y * b)) * (a * 0.5));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 0.96: tmp = x / (y / ((1.0 / a) - (b / a))) else: tmp = x / ((b * (y * b)) * (a * 0.5)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 0.96) tmp = Float64(x / Float64(y / Float64(Float64(1.0 / a) - Float64(b / a)))); else tmp = Float64(x / Float64(Float64(b * Float64(y * b)) * Float64(a * 0.5))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 0.96) tmp = x / (y / ((1.0 / a) - (b / a))); else tmp = x / ((b * (y * b)) * (a * 0.5)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 0.96], N[(x / N[(y / N[(N[(1.0 / a), $MachinePrecision] - N[(b / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(N[(b * N[(y * b), $MachinePrecision]), $MachinePrecision] * N[(a * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 0.96:\\
\;\;\;\;\frac{x}{\frac{y}{\frac{1}{a} - \frac{b}{a}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\left(b \cdot \left(y \cdot b\right)\right) \cdot \left(a \cdot 0.5\right)}\\
\end{array}
\end{array}
if b < 0.95999999999999996Initial program 96.8%
Taylor expanded in t around 0 79.2%
+-commutative79.2%
mul-1-neg79.2%
unsub-neg79.2%
Simplified79.2%
Taylor expanded in y around 0 51.7%
exp-neg51.7%
associate-*r/51.7%
*-rgt-identity51.7%
+-commutative51.7%
exp-sum51.7%
rem-exp-log52.7%
Simplified52.7%
Taylor expanded in b around 0 47.2%
+-commutative47.2%
mul-1-neg47.2%
unsub-neg47.2%
associate-/l*45.0%
Simplified45.0%
Taylor expanded in x around 0 47.2%
associate-/l*48.6%
Simplified48.6%
if 0.95999999999999996 < b Initial program 100.0%
associate-*l/87.8%
*-commutative87.8%
+-commutative87.8%
associate--l+87.8%
exp-sum66.2%
*-commutative66.2%
exp-to-pow66.2%
sub-neg66.2%
metadata-eval66.2%
exp-diff48.6%
*-commutative48.6%
exp-to-pow48.6%
Simplified48.6%
Taylor expanded in t around 0 67.7%
times-frac60.9%
Simplified60.9%
Taylor expanded in y around 0 82.7%
Taylor expanded in b around 0 61.8%
+-commutative61.8%
associate-*r*61.8%
distribute-rgt-out61.8%
unpow261.8%
Simplified61.8%
Taylor expanded in b around inf 61.8%
*-commutative61.8%
*-commutative61.8%
associate-*l*61.8%
unpow261.8%
associate-*l*56.6%
*-commutative56.6%
Simplified56.6%
Final simplification50.9%
(FPCore (x y z t a b) :precision binary64 (if (<= b -3.1e-50) (/ (* b (- (/ x a))) y) (/ x (* a (+ y (* y b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3.1e-50) {
tmp = (b * -(x / a)) / y;
} else {
tmp = x / (a * (y + (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 (b <= (-3.1d-50)) then
tmp = (b * -(x / a)) / y
else
tmp = x / (a * (y + (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 (b <= -3.1e-50) {
tmp = (b * -(x / a)) / y;
} else {
tmp = x / (a * (y + (y * b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -3.1e-50: tmp = (b * -(x / a)) / y else: tmp = x / (a * (y + (y * b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3.1e-50) tmp = Float64(Float64(b * Float64(-Float64(x / a))) / y); else tmp = Float64(x / Float64(a * Float64(y + Float64(y * b)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -3.1e-50) tmp = (b * -(x / a)) / y; else tmp = x / (a * (y + (y * b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3.1e-50], N[(N[(b * (-N[(x / a), $MachinePrecision])), $MachinePrecision] / y), $MachinePrecision], N[(x / N[(a * N[(y + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.1 \cdot 10^{-50}:\\
\;\;\;\;\frac{b \cdot \left(-\frac{x}{a}\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot b\right)}\\
\end{array}
\end{array}
if b < -3.1000000000000002e-50Initial program 99.6%
Taylor expanded in t around 0 84.0%
+-commutative84.0%
mul-1-neg84.0%
unsub-neg84.0%
Simplified84.0%
Taylor expanded in y around 0 69.7%
exp-neg69.7%
associate-*r/69.7%
*-rgt-identity69.7%
+-commutative69.7%
exp-sum69.8%
rem-exp-log70.1%
Simplified70.1%
Taylor expanded in b around 0 56.6%
+-commutative56.6%
mul-1-neg56.6%
unsub-neg56.6%
associate-/l*53.5%
Simplified53.5%
Taylor expanded in b around inf 46.0%
associate-*r/46.0%
associate-/r*53.6%
mul-1-neg53.6%
*-commutative53.6%
distribute-rgt-neg-in53.6%
neg-mul-153.6%
associate-*l/52.1%
neg-mul-152.1%
Simplified52.1%
if -3.1000000000000002e-50 < b Initial program 97.1%
associate-*l/88.3%
*-commutative88.3%
+-commutative88.3%
associate--l+88.3%
exp-sum75.4%
*-commutative75.4%
exp-to-pow76.3%
sub-neg76.3%
metadata-eval76.3%
exp-diff69.6%
*-commutative69.6%
exp-to-pow69.6%
Simplified69.6%
Taylor expanded in t around 0 71.5%
times-frac68.5%
Simplified68.5%
Taylor expanded in y around 0 60.1%
Taylor expanded in b around 0 40.7%
distribute-lft-out42.7%
*-commutative42.7%
Simplified42.7%
Final simplification45.0%
(FPCore (x y z t a b) :precision binary64 (if (<= b -1.7e-50) (/ (- (/ (* x b) a)) y) (/ x (* a (+ y (* y b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.7e-50) {
tmp = -((x * b) / a) / y;
} else {
tmp = x / (a * (y + (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 (b <= (-1.7d-50)) then
tmp = -((x * b) / a) / y
else
tmp = x / (a * (y + (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 (b <= -1.7e-50) {
tmp = -((x * b) / a) / y;
} else {
tmp = x / (a * (y + (y * b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1.7e-50: tmp = -((x * b) / a) / y else: tmp = x / (a * (y + (y * b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.7e-50) tmp = Float64(Float64(-Float64(Float64(x * b) / a)) / y); else tmp = Float64(x / Float64(a * Float64(y + Float64(y * b)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -1.7e-50) tmp = -((x * b) / a) / y; else tmp = x / (a * (y + (y * b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.7e-50], N[((-N[(N[(x * b), $MachinePrecision] / a), $MachinePrecision]) / y), $MachinePrecision], N[(x / N[(a * N[(y + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.7 \cdot 10^{-50}:\\
\;\;\;\;\frac{-\frac{x \cdot b}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot b\right)}\\
\end{array}
\end{array}
if b < -1.70000000000000007e-50Initial program 99.6%
Taylor expanded in t around 0 84.0%
+-commutative84.0%
mul-1-neg84.0%
unsub-neg84.0%
Simplified84.0%
Taylor expanded in y around 0 69.7%
exp-neg69.7%
associate-*r/69.7%
*-rgt-identity69.7%
+-commutative69.7%
exp-sum69.8%
rem-exp-log70.1%
Simplified70.1%
Taylor expanded in b around 0 56.6%
+-commutative56.6%
mul-1-neg56.6%
unsub-neg56.6%
associate-/l*53.5%
Simplified53.5%
Taylor expanded in b around inf 53.6%
if -1.70000000000000007e-50 < b Initial program 97.1%
associate-*l/88.3%
*-commutative88.3%
+-commutative88.3%
associate--l+88.3%
exp-sum75.4%
*-commutative75.4%
exp-to-pow76.3%
sub-neg76.3%
metadata-eval76.3%
exp-diff69.6%
*-commutative69.6%
exp-to-pow69.6%
Simplified69.6%
Taylor expanded in t around 0 71.5%
times-frac68.5%
Simplified68.5%
Taylor expanded in y around 0 60.1%
Taylor expanded in b around 0 40.7%
distribute-lft-out42.7%
*-commutative42.7%
Simplified42.7%
Final simplification45.4%
(FPCore (x y z t a b) :precision binary64 (if (<= b -4e-292) (/ (- x (* x b)) (* y a)) (/ (/ x (+ a (* a b))) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -4e-292) {
tmp = (x - (x * b)) / (y * a);
} else {
tmp = (x / (a + (a * 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 <= (-4d-292)) then
tmp = (x - (x * b)) / (y * a)
else
tmp = (x / (a + (a * 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 <= -4e-292) {
tmp = (x - (x * b)) / (y * a);
} else {
tmp = (x / (a + (a * b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -4e-292: tmp = (x - (x * b)) / (y * a) else: tmp = (x / (a + (a * b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -4e-292) tmp = Float64(Float64(x - Float64(x * b)) / Float64(y * a)); else tmp = Float64(Float64(x / Float64(a + Float64(a * b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -4e-292) tmp = (x - (x * b)) / (y * a); else tmp = (x / (a + (a * b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -4e-292], N[(N[(x - N[(x * b), $MachinePrecision]), $MachinePrecision] / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(N[(x / N[(a + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -4 \cdot 10^{-292}:\\
\;\;\;\;\frac{x - x \cdot b}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + a \cdot b}}{y}\\
\end{array}
\end{array}
if b < -4.0000000000000002e-292Initial program 97.5%
Taylor expanded in t around 0 80.3%
+-commutative80.3%
mul-1-neg80.3%
unsub-neg80.3%
Simplified80.3%
Taylor expanded in y around 0 55.3%
exp-neg55.3%
associate-*r/55.3%
*-rgt-identity55.3%
+-commutative55.3%
exp-sum55.3%
rem-exp-log56.0%
Simplified56.0%
Taylor expanded in b around 0 49.0%
+-commutative49.0%
mul-1-neg49.0%
unsub-neg49.0%
associate-/l*47.4%
Simplified47.4%
Taylor expanded in x around 0 48.2%
*-commutative48.2%
sub-neg48.2%
mul-1-neg48.2%
+-commutative48.2%
*-commutative48.2%
+-commutative48.2%
distribute-rgt-in48.2%
associate-*l/48.2%
*-lft-identity48.2%
mul-1-neg48.2%
cancel-sign-sub-inv48.2%
associate-*l/49.0%
div-sub49.0%
associate-/r*49.7%
*-commutative49.7%
*-commutative49.7%
Simplified49.7%
if -4.0000000000000002e-292 < b Initial program 97.9%
Taylor expanded in t around 0 84.4%
+-commutative84.4%
mul-1-neg84.4%
unsub-neg84.4%
Simplified84.4%
Taylor expanded in y around 0 64.6%
exp-neg64.6%
associate-*r/64.6%
*-rgt-identity64.6%
+-commutative64.6%
exp-sum64.6%
rem-exp-log65.4%
Simplified65.4%
Taylor expanded in b around 0 41.6%
Final simplification45.4%
(FPCore (x y z t a b) :precision binary64 (if (<= b -0.45) (/ (- (* x b)) (* y a)) (/ (/ x a) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -0.45) {
tmp = -(x * b) / (y * a);
} else {
tmp = (x / a) / 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 <= (-0.45d0)) then
tmp = -(x * b) / (y * a)
else
tmp = (x / a) / 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 <= -0.45) {
tmp = -(x * b) / (y * a);
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -0.45: tmp = -(x * b) / (y * a) else: tmp = (x / a) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -0.45) tmp = Float64(Float64(-Float64(x * b)) / Float64(y * a)); else tmp = Float64(Float64(x / a) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -0.45) tmp = -(x * b) / (y * a); else tmp = (x / a) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -0.45], N[((-N[(x * b), $MachinePrecision]) / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.45:\\
\;\;\;\;\frac{-x \cdot b}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if b < -0.450000000000000011Initial program 99.9%
Taylor expanded in t around 0 88.1%
+-commutative88.1%
mul-1-neg88.1%
unsub-neg88.1%
Simplified88.1%
Taylor expanded in y around 0 74.4%
exp-neg74.4%
associate-*r/74.4%
*-rgt-identity74.4%
+-commutative74.4%
exp-sum74.4%
rem-exp-log74.4%
Simplified74.4%
Taylor expanded in b around 0 57.7%
+-commutative57.7%
mul-1-neg57.7%
unsub-neg57.7%
associate-/l*53.9%
Simplified53.9%
Taylor expanded in b around inf 55.6%
associate-*r/55.6%
mul-1-neg55.6%
*-commutative55.6%
*-commutative55.6%
Simplified55.6%
if -0.450000000000000011 < b Initial program 97.2%
Taylor expanded in t around 0 81.2%
+-commutative81.2%
mul-1-neg81.2%
unsub-neg81.2%
Simplified81.2%
Taylor expanded in y around 0 56.9%
exp-neg56.9%
associate-*r/56.9%
*-rgt-identity56.9%
+-commutative56.9%
exp-sum56.9%
rem-exp-log57.8%
Simplified57.8%
Taylor expanded in b around 0 35.2%
Final simplification39.2%
(FPCore (x y z t a b) :precision binary64 (if (<= t 3.3e-118) (/ (/ x a) y) (* (/ x y) (/ 1.0 a))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 3.3e-118) {
tmp = (x / a) / y;
} else {
tmp = (x / y) * (1.0 / a);
}
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 <= 3.3d-118) then
tmp = (x / a) / y
else
tmp = (x / y) * (1.0d0 / a)
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 <= 3.3e-118) {
tmp = (x / a) / y;
} else {
tmp = (x / y) * (1.0 / a);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= 3.3e-118: tmp = (x / a) / y else: tmp = (x / y) * (1.0 / a) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= 3.3e-118) tmp = Float64(Float64(x / a) / y); else tmp = Float64(Float64(x / y) * Float64(1.0 / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= 3.3e-118) tmp = (x / a) / y; else tmp = (x / y) * (1.0 / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, 3.3e-118], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / y), $MachinePrecision] * N[(1.0 / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 3.3 \cdot 10^{-118}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y} \cdot \frac{1}{a}\\
\end{array}
\end{array}
if t < 3.3e-118Initial program 98.4%
Taylor expanded in t around 0 88.5%
+-commutative88.5%
mul-1-neg88.5%
unsub-neg88.5%
Simplified88.5%
Taylor expanded in y around 0 66.2%
exp-neg66.2%
associate-*r/66.2%
*-rgt-identity66.2%
+-commutative66.2%
exp-sum66.3%
rem-exp-log67.2%
Simplified67.2%
Taylor expanded in b around 0 39.6%
if 3.3e-118 < t Initial program 96.2%
associate-*l/87.9%
*-commutative87.9%
+-commutative87.9%
associate--l+87.9%
exp-sum66.4%
*-commutative66.4%
exp-to-pow67.0%
sub-neg67.0%
metadata-eval67.0%
exp-diff62.0%
*-commutative62.0%
exp-to-pow62.0%
Simplified62.0%
Taylor expanded in t around 0 59.1%
times-frac50.4%
Simplified50.4%
Taylor expanded in y around 0 50.8%
Taylor expanded in b around 0 25.8%
*-commutative25.8%
Simplified25.8%
div-inv25.8%
add-exp-log12.8%
rec-exp12.8%
log-prod12.7%
pow112.7%
metadata-eval12.7%
pow-flip12.7%
inv-pow12.7%
log-prod12.8%
div-inv12.8%
rec-exp12.8%
add-exp-log25.8%
div-inv25.8%
associate-/r/32.9%
Applied egg-rr32.9%
Final simplification37.5%
(FPCore (x y z t a b) :precision binary64 (if (<= z 9.8e-82) (/ x (* y a)) (/ (/ x a) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= 9.8e-82) {
tmp = x / (y * a);
} else {
tmp = (x / a) / 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 (z <= 9.8d-82) then
tmp = x / (y * a)
else
tmp = (x / a) / 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 (z <= 9.8e-82) {
tmp = x / (y * a);
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= 9.8e-82: tmp = x / (y * a) else: tmp = (x / a) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= 9.8e-82) tmp = Float64(x / Float64(y * a)); else tmp = Float64(Float64(x / a) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= 9.8e-82) tmp = x / (y * a); else tmp = (x / a) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, 9.8e-82], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq 9.8 \cdot 10^{-82}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if z < 9.8000000000000006e-82Initial program 97.4%
associate-*l/85.6%
*-commutative85.6%
+-commutative85.6%
associate--l+85.6%
exp-sum68.2%
*-commutative68.2%
exp-to-pow68.9%
sub-neg68.9%
metadata-eval68.9%
exp-diff64.1%
*-commutative64.1%
exp-to-pow64.1%
Simplified64.1%
Taylor expanded in t around 0 73.1%
times-frac66.5%
Simplified66.5%
Taylor expanded in y around 0 62.2%
Taylor expanded in b around 0 36.8%
*-commutative36.8%
Simplified36.8%
if 9.8000000000000006e-82 < z Initial program 97.9%
Taylor expanded in t around 0 82.7%
+-commutative82.7%
mul-1-neg82.7%
unsub-neg82.7%
Simplified82.7%
Taylor expanded in y around 0 61.3%
exp-neg61.3%
associate-*r/61.3%
*-rgt-identity61.3%
+-commutative61.3%
exp-sum61.3%
rem-exp-log62.0%
Simplified62.0%
Taylor expanded in b around 0 36.7%
Final simplification36.8%
(FPCore (x y z t a b) :precision binary64 (/ x (* y a)))
double code(double x, double y, double z, double t, double a, double b) {
return x / (y * a);
}
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 * a)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / (y * a);
}
def code(x, y, z, t, a, b): return x / (y * a)
function code(x, y, z, t, a, b) return Float64(x / Float64(y * a)) end
function tmp = code(x, y, z, t, a, b) tmp = x / (y * a); end
code[x_, y_, z_, t_, a_, b_] := N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y \cdot a}
\end{array}
Initial program 97.7%
associate-*l/89.1%
*-commutative89.1%
+-commutative89.1%
associate--l+89.1%
exp-sum72.7%
*-commutative72.7%
exp-to-pow73.4%
sub-neg73.4%
metadata-eval73.4%
exp-diff65.2%
*-commutative65.2%
exp-to-pow65.2%
Simplified65.2%
Taylor expanded in t around 0 68.8%
times-frac66.1%
Simplified66.1%
Taylor expanded in y around 0 61.8%
Taylor expanded in b around 0 35.0%
*-commutative35.0%
Simplified35.0%
Final simplification35.0%
(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 2023290
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:precision binary64
:herbie-target
(if (< t -0.8845848504127471) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z)))) (if (< t 852031.2288374073) (/ (* (/ x y) (pow a (- t 1.0))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))