
(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 98.3%
Final simplification98.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= (+ t -1.0) -5e+67)
(/ (* x (exp (- (* (+ t -1.0) (log a)) b))) y)
(if (<= (+ t -1.0) 4e+105)
(/ (* x (exp (- (- (* y (log z)) (log a)) b))) y)
(* x (/ (pow a (+ t -1.0)) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((t + -1.0) <= -5e+67) {
tmp = (x * exp((((t + -1.0) * log(a)) - b))) / y;
} else if ((t + -1.0) <= 4e+105) {
tmp = (x * exp((((y * log(z)) - log(a)) - b))) / y;
} else {
tmp = x * (pow(a, (t + -1.0)) / y);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((t + (-1.0d0)) <= (-5d+67)) then
tmp = (x * exp((((t + (-1.0d0)) * log(a)) - b))) / y
else if ((t + (-1.0d0)) <= 4d+105) then
tmp = (x * exp((((y * log(z)) - log(a)) - b))) / y
else
tmp = x * ((a ** (t + (-1.0d0))) / y)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((t + -1.0) <= -5e+67) {
tmp = (x * Math.exp((((t + -1.0) * Math.log(a)) - b))) / y;
} else if ((t + -1.0) <= 4e+105) {
tmp = (x * Math.exp((((y * Math.log(z)) - Math.log(a)) - b))) / y;
} else {
tmp = x * (Math.pow(a, (t + -1.0)) / y);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (t + -1.0) <= -5e+67: tmp = (x * math.exp((((t + -1.0) * math.log(a)) - b))) / y elif (t + -1.0) <= 4e+105: tmp = (x * math.exp((((y * math.log(z)) - math.log(a)) - b))) / y else: tmp = x * (math.pow(a, (t + -1.0)) / y) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(t + -1.0) <= -5e+67) tmp = Float64(Float64(x * exp(Float64(Float64(Float64(t + -1.0) * log(a)) - b))) / y); elseif (Float64(t + -1.0) <= 4e+105) tmp = Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) - log(a)) - b))) / y); else tmp = Float64(x * Float64((a ^ Float64(t + -1.0)) / y)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((t + -1.0) <= -5e+67) tmp = (x * exp((((t + -1.0) * log(a)) - b))) / y; elseif ((t + -1.0) <= 4e+105) tmp = (x * exp((((y * log(z)) - log(a)) - b))) / y; else tmp = x * ((a ^ (t + -1.0)) / y); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(t + -1.0), $MachinePrecision], -5e+67], N[(N[(x * N[Exp[N[(N[(N[(t + -1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[N[(t + -1.0), $MachinePrecision], 4e+105], 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[(x * N[(N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t + -1 \leq -5 \cdot 10^{+67}:\\
\;\;\;\;\frac{x \cdot e^{\left(t + -1\right) \cdot \log a - b}}{y}\\
\mathbf{elif}\;t + -1 \leq 4 \cdot 10^{+105}:\\
\;\;\;\;\frac{x \cdot e^{\left(y \cdot \log z - \log a\right) - b}}{y}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{{a}^{\left(t + -1\right)}}{y}\\
\end{array}
\end{array}
if (-.f64 t #s(literal 1 binary64)) < -4.99999999999999976e67Initial program 100.0%
Taylor expanded in y around 0 95.7%
if -4.99999999999999976e67 < (-.f64 t #s(literal 1 binary64)) < 3.9999999999999998e105Initial program 97.5%
Taylor expanded in t around 0 95.6%
+-commutative95.6%
mul-1-neg95.6%
unsub-neg95.6%
Simplified95.6%
if 3.9999999999999998e105 < (-.f64 t #s(literal 1 binary64)) Initial program 100.0%
associate-/l*100.0%
exp-diff92.7%
associate-/l/92.7%
exp-sum75.6%
associate-/l*75.6%
*-commutative75.6%
exp-to-pow75.6%
*-commutative75.6%
exp-to-pow75.6%
sub-neg75.6%
metadata-eval75.6%
Simplified75.6%
Taylor expanded in y around 0 92.7%
exp-to-pow92.7%
sub-neg92.7%
metadata-eval92.7%
associate-/r*92.7%
Simplified92.7%
Taylor expanded in b around 0 97.6%
Simplified97.6%
Final simplification95.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= (+ t -1.0) -5e+55)
(* x (/ (pow a (+ t -1.0)) y))
(if (<= (+ t -1.0) 2e+79)
(* (/ (pow z y) a) (/ x (* y (exp b))))
(/ (* x (/ (pow a t) a)) y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((t + -1.0) <= -5e+55) {
tmp = x * (pow(a, (t + -1.0)) / y);
} else if ((t + -1.0) <= 2e+79) {
tmp = (pow(z, y) / a) * (x / (y * exp(b)));
} else {
tmp = (x * (pow(a, t) / 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 ((t + (-1.0d0)) <= (-5d+55)) then
tmp = x * ((a ** (t + (-1.0d0))) / y)
else if ((t + (-1.0d0)) <= 2d+79) then
tmp = ((z ** y) / a) * (x / (y * exp(b)))
else
tmp = (x * ((a ** t) / 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 ((t + -1.0) <= -5e+55) {
tmp = x * (Math.pow(a, (t + -1.0)) / y);
} else if ((t + -1.0) <= 2e+79) {
tmp = (Math.pow(z, y) / a) * (x / (y * Math.exp(b)));
} else {
tmp = (x * (Math.pow(a, t) / a)) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (t + -1.0) <= -5e+55: tmp = x * (math.pow(a, (t + -1.0)) / y) elif (t + -1.0) <= 2e+79: tmp = (math.pow(z, y) / a) * (x / (y * math.exp(b))) else: tmp = (x * (math.pow(a, t) / a)) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(t + -1.0) <= -5e+55) tmp = Float64(x * Float64((a ^ Float64(t + -1.0)) / y)); elseif (Float64(t + -1.0) <= 2e+79) tmp = Float64(Float64((z ^ y) / a) * Float64(x / Float64(y * exp(b)))); else tmp = Float64(Float64(x * Float64((a ^ t) / a)) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((t + -1.0) <= -5e+55) tmp = x * ((a ^ (t + -1.0)) / y); elseif ((t + -1.0) <= 2e+79) tmp = ((z ^ y) / a) * (x / (y * exp(b))); else tmp = (x * ((a ^ t) / a)) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(t + -1.0), $MachinePrecision], -5e+55], N[(x * N[(N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(t + -1.0), $MachinePrecision], 2e+79], N[(N[(N[Power[z, y], $MachinePrecision] / a), $MachinePrecision] * N[(x / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x * N[(N[Power[a, t], $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t + -1 \leq -5 \cdot 10^{+55}:\\
\;\;\;\;x \cdot \frac{{a}^{\left(t + -1\right)}}{y}\\
\mathbf{elif}\;t + -1 \leq 2 \cdot 10^{+79}:\\
\;\;\;\;\frac{{z}^{y}}{a} \cdot \frac{x}{y \cdot e^{b}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \frac{{a}^{t}}{a}}{y}\\
\end{array}
\end{array}
if (-.f64 t #s(literal 1 binary64)) < -5.00000000000000046e55Initial program 100.0%
associate-/l*100.0%
exp-diff82.4%
associate-/l/82.4%
exp-sum70.6%
associate-/l*70.6%
*-commutative70.6%
exp-to-pow70.6%
*-commutative70.6%
exp-to-pow70.6%
sub-neg70.6%
metadata-eval70.6%
Simplified70.6%
Taylor expanded in y around 0 82.4%
exp-to-pow82.4%
sub-neg82.4%
metadata-eval82.4%
associate-/r*82.4%
Simplified82.4%
Taylor expanded in b around 0 92.3%
Simplified92.3%
if -5.00000000000000046e55 < (-.f64 t #s(literal 1 binary64)) < 1.99999999999999993e79Initial program 97.3%
associate-/l*97.9%
exp-diff82.3%
associate-/l/82.3%
exp-sum78.5%
associate-/l*76.7%
*-commutative76.7%
exp-to-pow76.7%
*-commutative76.7%
exp-to-pow77.4%
sub-neg77.4%
metadata-eval77.4%
Simplified77.4%
Taylor expanded in t around 0 82.5%
*-commutative82.5%
times-frac82.5%
Simplified82.5%
if 1.99999999999999993e79 < (-.f64 t #s(literal 1 binary64)) Initial program 100.0%
Taylor expanded in y around 0 95.6%
div-exp89.0%
exp-to-pow89.0%
sub-neg89.0%
metadata-eval89.0%
Simplified89.0%
unpow-prod-up89.0%
unpow-189.0%
Applied egg-rr89.0%
associate-*r/89.0%
*-rgt-identity89.0%
Simplified89.0%
Taylor expanded in b around 0 93.4%
associate-/l*93.4%
Simplified93.4%
Final simplification86.4%
(FPCore (x y z t a b) :precision binary64 (if (or (<= y -8e+16) (not (<= y 2.6e+95))) (/ (/ (* x (pow z y)) a) 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 <= -8e+16) || !(y <= 2.6e+95)) {
tmp = ((x * pow(z, y)) / a) / 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 <= (-8d+16)) .or. (.not. (y <= 2.6d+95))) then
tmp = ((x * (z ** y)) / a) / 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 <= -8e+16) || !(y <= 2.6e+95)) {
tmp = ((x * Math.pow(z, y)) / a) / 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 <= -8e+16) or not (y <= 2.6e+95): tmp = ((x * math.pow(z, y)) / a) / 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 <= -8e+16) || !(y <= 2.6e+95)) tmp = Float64(Float64(Float64(x * (z ^ y)) / a) / 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 <= -8e+16) || ~((y <= 2.6e+95))) tmp = ((x * (z ^ y)) / a) / 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, -8e+16], N[Not[LessEqual[y, 2.6e+95]], $MachinePrecision]], N[(N[(N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision] / a), $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 -8 \cdot 10^{+16} \lor \neg \left(y \leq 2.6 \cdot 10^{+95}\right):\\
\;\;\;\;\frac{\frac{x \cdot {z}^{y}}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot e^{\left(t + -1\right) \cdot \log a - b}}{y}\\
\end{array}
\end{array}
if y < -8e16 or 2.5999999999999999e95 < y Initial program 100.0%
Taylor expanded in b around 0 93.2%
exp-sum68.4%
*-commutative68.4%
exp-to-pow68.4%
exp-to-pow68.4%
sub-neg68.4%
metadata-eval68.4%
Simplified68.4%
Taylor expanded in t around 0 81.5%
if -8e16 < y < 2.5999999999999999e95Initial program 97.2%
Taylor expanded in y around 0 95.4%
Final simplification89.9%
(FPCore (x y z t a b) :precision binary64 (if (or (<= y -5.5e+16) (not (<= y 3.9e+87))) (/ (/ (* x (pow z y)) a) y) (/ (* x (/ (/ (pow a t) a) (exp b))) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((y <= -5.5e+16) || !(y <= 3.9e+87)) {
tmp = ((x * pow(z, y)) / a) / y;
} else {
tmp = (x * ((pow(a, t) / a) / exp(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 <= (-5.5d+16)) .or. (.not. (y <= 3.9d+87))) then
tmp = ((x * (z ** y)) / a) / y
else
tmp = (x * (((a ** t) / a) / exp(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 <= -5.5e+16) || !(y <= 3.9e+87)) {
tmp = ((x * Math.pow(z, y)) / a) / y;
} else {
tmp = (x * ((Math.pow(a, t) / a) / Math.exp(b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (y <= -5.5e+16) or not (y <= 3.9e+87): tmp = ((x * math.pow(z, y)) / a) / y else: tmp = (x * ((math.pow(a, t) / a) / math.exp(b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((y <= -5.5e+16) || !(y <= 3.9e+87)) tmp = Float64(Float64(Float64(x * (z ^ y)) / a) / y); else tmp = Float64(Float64(x * Float64(Float64((a ^ t) / a) / exp(b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((y <= -5.5e+16) || ~((y <= 3.9e+87))) tmp = ((x * (z ^ y)) / a) / y; else tmp = (x * (((a ^ t) / a) / exp(b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -5.5e+16], N[Not[LessEqual[y, 3.9e+87]], $MachinePrecision]], N[(N[(N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision] / y), $MachinePrecision], N[(N[(x * N[(N[(N[Power[a, t], $MachinePrecision] / a), $MachinePrecision] / N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.5 \cdot 10^{+16} \lor \neg \left(y \leq 3.9 \cdot 10^{+87}\right):\\
\;\;\;\;\frac{\frac{x \cdot {z}^{y}}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \frac{\frac{{a}^{t}}{a}}{e^{b}}}{y}\\
\end{array}
\end{array}
if y < -5.5e16 or 3.9000000000000002e87 < y Initial program 100.0%
Taylor expanded in b around 0 92.3%
exp-sum67.8%
*-commutative67.8%
exp-to-pow67.8%
exp-to-pow67.8%
sub-neg67.8%
metadata-eval67.8%
Simplified67.8%
Taylor expanded in t around 0 80.7%
if -5.5e16 < y < 3.9000000000000002e87Initial program 97.2%
Taylor expanded in y around 0 95.3%
div-exp90.2%
exp-to-pow91.0%
sub-neg91.0%
metadata-eval91.0%
Simplified91.0%
unpow-prod-up91.1%
unpow-191.1%
Applied egg-rr91.1%
associate-*r/91.1%
*-rgt-identity91.1%
Simplified91.1%
Final simplification86.9%
(FPCore (x y z t a b) :precision binary64 (if (or (<= y -3.6e+16) (not (<= y 1.15e+92))) (/ (/ (* x (pow z y)) a) y) (* x (/ (/ (pow a (+ t -1.0)) y) (exp b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((y <= -3.6e+16) || !(y <= 1.15e+92)) {
tmp = ((x * pow(z, y)) / a) / y;
} else {
tmp = x * ((pow(a, (t + -1.0)) / 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 ((y <= (-3.6d+16)) .or. (.not. (y <= 1.15d+92))) then
tmp = ((x * (z ** y)) / a) / y
else
tmp = x * (((a ** (t + (-1.0d0))) / 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 ((y <= -3.6e+16) || !(y <= 1.15e+92)) {
tmp = ((x * Math.pow(z, y)) / a) / y;
} else {
tmp = x * ((Math.pow(a, (t + -1.0)) / y) / Math.exp(b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (y <= -3.6e+16) or not (y <= 1.15e+92): tmp = ((x * math.pow(z, y)) / a) / y else: tmp = x * ((math.pow(a, (t + -1.0)) / y) / math.exp(b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((y <= -3.6e+16) || !(y <= 1.15e+92)) tmp = Float64(Float64(Float64(x * (z ^ y)) / a) / y); else tmp = Float64(x * Float64(Float64((a ^ Float64(t + -1.0)) / y) / exp(b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((y <= -3.6e+16) || ~((y <= 1.15e+92))) tmp = ((x * (z ^ y)) / a) / y; else tmp = x * (((a ^ (t + -1.0)) / y) / exp(b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -3.6e+16], N[Not[LessEqual[y, 1.15e+92]], $MachinePrecision]], N[(N[(N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision] / y), $MachinePrecision], N[(x * N[(N[(N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision] / y), $MachinePrecision] / N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.6 \cdot 10^{+16} \lor \neg \left(y \leq 1.15 \cdot 10^{+92}\right):\\
\;\;\;\;\frac{\frac{x \cdot {z}^{y}}{a}}{y}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{\frac{{a}^{\left(t + -1\right)}}{y}}{e^{b}}\\
\end{array}
\end{array}
if y < -3.6e16 or 1.14999999999999999e92 < y Initial program 100.0%
Taylor expanded in b around 0 92.3%
exp-sum67.8%
*-commutative67.8%
exp-to-pow67.8%
exp-to-pow67.8%
sub-neg67.8%
metadata-eval67.8%
Simplified67.8%
Taylor expanded in t around 0 80.7%
if -3.6e16 < y < 1.14999999999999999e92Initial program 97.2%
associate-/l*97.8%
exp-diff89.4%
associate-/l/89.4%
exp-sum86.8%
associate-/l*86.8%
*-commutative86.8%
exp-to-pow86.8%
*-commutative86.8%
exp-to-pow87.6%
sub-neg87.6%
metadata-eval87.6%
Simplified87.6%
Taylor expanded in y around 0 90.8%
exp-to-pow91.6%
sub-neg91.6%
metadata-eval91.6%
associate-/r*89.0%
Simplified89.0%
Final simplification85.7%
(FPCore (x y z t a b) :precision binary64 (if (or (<= (+ t -1.0) -2e+62) (not (<= (+ t -1.0) 4e+105))) (* x (/ (pow a (+ t -1.0)) y)) (/ x (* a (* y (exp b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (((t + -1.0) <= -2e+62) || !((t + -1.0) <= 4e+105)) {
tmp = x * (pow(a, (t + -1.0)) / y);
} else {
tmp = x / (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)) <= (-2d+62)) .or. (.not. ((t + (-1.0d0)) <= 4d+105))) then
tmp = x * ((a ** (t + (-1.0d0))) / y)
else
tmp = x / (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) <= -2e+62) || !((t + -1.0) <= 4e+105)) {
tmp = x * (Math.pow(a, (t + -1.0)) / y);
} else {
tmp = x / (a * (y * Math.exp(b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if ((t + -1.0) <= -2e+62) or not ((t + -1.0) <= 4e+105): tmp = x * (math.pow(a, (t + -1.0)) / y) else: tmp = x / (a * (y * math.exp(b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((Float64(t + -1.0) <= -2e+62) || !(Float64(t + -1.0) <= 4e+105)) tmp = Float64(x * Float64((a ^ Float64(t + -1.0)) / y)); else tmp = Float64(x / 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) <= -2e+62) || ~(((t + -1.0) <= 4e+105))) tmp = x * ((a ^ (t + -1.0)) / y); else tmp = x / (a * (y * exp(b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[N[(t + -1.0), $MachinePrecision], -2e+62], N[Not[LessEqual[N[(t + -1.0), $MachinePrecision], 4e+105]], $MachinePrecision]], N[(x * N[(N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(x / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t + -1 \leq -2 \cdot 10^{+62} \lor \neg \left(t + -1 \leq 4 \cdot 10^{+105}\right):\\
\;\;\;\;x \cdot \frac{{a}^{\left(t + -1\right)}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot e^{b}\right)}\\
\end{array}
\end{array}
if (-.f64 t #s(literal 1 binary64)) < -2.00000000000000007e62 or 3.9999999999999998e105 < (-.f64 t #s(literal 1 binary64)) Initial program 100.0%
associate-/l*100.0%
exp-diff87.8%
associate-/l/87.8%
exp-sum73.3%
associate-/l*73.3%
*-commutative73.3%
exp-to-pow73.3%
*-commutative73.3%
exp-to-pow73.3%
sub-neg73.3%
metadata-eval73.3%
Simplified73.3%
Taylor expanded in y around 0 87.8%
exp-to-pow87.8%
sub-neg87.8%
metadata-eval87.8%
associate-/r*87.8%
Simplified87.8%
Taylor expanded in b around 0 95.6%
Simplified95.6%
if -2.00000000000000007e62 < (-.f64 t #s(literal 1 binary64)) < 3.9999999999999998e105Initial program 97.4%
associate-/l*98.0%
exp-diff81.7%
associate-/l/81.7%
exp-sum76.9%
associate-/l*75.1%
*-commutative75.1%
exp-to-pow75.1%
*-commutative75.1%
exp-to-pow75.8%
sub-neg75.8%
metadata-eval75.8%
Simplified75.8%
Taylor expanded in y around 0 67.3%
exp-to-pow68.1%
sub-neg68.1%
metadata-eval68.1%
associate-/r*63.9%
Simplified63.9%
Taylor expanded in t around 0 71.8%
Final simplification80.2%
(FPCore (x y z t a b) :precision binary64 (/ x (* a (* y (exp b)))))
double code(double x, double y, double z, double t, double a, double b) {
return x / (a * (y * exp(b)));
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x / (a * (y * exp(b)))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / (a * (y * Math.exp(b)));
}
def code(x, y, z, t, a, b): return x / (a * (y * math.exp(b)))
function code(x, y, z, t, a, b) return Float64(x / Float64(a * Float64(y * exp(b)))) end
function tmp = code(x, y, z, t, a, b) tmp = x / (a * (y * exp(b))); end
code[x_, y_, z_, t_, a_, b_] := N[(x / N[(a * N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{a \cdot \left(y \cdot e^{b}\right)}
\end{array}
Initial program 98.3%
associate-/l*98.7%
exp-diff83.9%
associate-/l/83.9%
exp-sum75.7%
associate-/l*74.5%
*-commutative74.5%
exp-to-pow74.5%
*-commutative74.5%
exp-to-pow75.0%
sub-neg75.0%
metadata-eval75.0%
Simplified75.0%
Taylor expanded in y around 0 74.5%
exp-to-pow75.0%
sub-neg75.0%
metadata-eval75.0%
associate-/r*72.3%
Simplified72.3%
Taylor expanded in t around 0 62.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -7.4e+128)
(/ (* x (+ (/ 1.0 a) (* b (+ (* 0.5 (/ b a)) (/ -1.0 a))))) y)
(if (<= b -5.2e+90)
(/
x
(*
a
(+ y (* b (+ y (* b (+ (* 0.16666666666666666 (* y b)) (* y 0.5))))))))
(if (<= b -1.1e-143)
(* x (/ 1.0 (* y a)))
(if (<= b 2e-184)
(/ (* b (- (/ x (* a b)) (/ x a))) y)
(/
(/
x
(+
a
(* b (+ a (* b (+ (* 0.16666666666666666 (* a b)) (* a 0.5)))))))
y))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -7.4e+128) {
tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y;
} else if (b <= -5.2e+90) {
tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5)))))));
} else if (b <= -1.1e-143) {
tmp = x * (1.0 / (y * a));
} else if (b <= 2e-184) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (b * ((0.16666666666666666 * (a * b)) + (a * 0.5))))))) / y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-7.4d+128)) then
tmp = (x * ((1.0d0 / a) + (b * ((0.5d0 * (b / a)) + ((-1.0d0) / a))))) / y
else if (b <= (-5.2d+90)) then
tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666d0 * (y * b)) + (y * 0.5d0)))))))
else if (b <= (-1.1d-143)) then
tmp = x * (1.0d0 / (y * a))
else if (b <= 2d-184) then
tmp = (b * ((x / (a * b)) - (x / a))) / y
else
tmp = (x / (a + (b * (a + (b * ((0.16666666666666666d0 * (a * b)) + (a * 0.5d0))))))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -7.4e+128) {
tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y;
} else if (b <= -5.2e+90) {
tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5)))))));
} else if (b <= -1.1e-143) {
tmp = x * (1.0 / (y * a));
} else if (b <= 2e-184) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (b * ((0.16666666666666666 * (a * b)) + (a * 0.5))))))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -7.4e+128: tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y elif b <= -5.2e+90: tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5))))))) elif b <= -1.1e-143: tmp = x * (1.0 / (y * a)) elif b <= 2e-184: tmp = (b * ((x / (a * b)) - (x / a))) / y else: tmp = (x / (a + (b * (a + (b * ((0.16666666666666666 * (a * b)) + (a * 0.5))))))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -7.4e+128) tmp = Float64(Float64(x * Float64(Float64(1.0 / a) + Float64(b * Float64(Float64(0.5 * Float64(b / a)) + Float64(-1.0 / a))))) / y); elseif (b <= -5.2e+90) tmp = Float64(x / Float64(a * Float64(y + Float64(b * Float64(y + Float64(b * Float64(Float64(0.16666666666666666 * Float64(y * b)) + Float64(y * 0.5)))))))); elseif (b <= -1.1e-143) tmp = Float64(x * Float64(1.0 / Float64(y * a))); elseif (b <= 2e-184) tmp = Float64(Float64(b * Float64(Float64(x / Float64(a * b)) - Float64(x / a))) / y); else tmp = Float64(Float64(x / Float64(a + Float64(b * Float64(a + Float64(b * Float64(Float64(0.16666666666666666 * Float64(a * b)) + Float64(a * 0.5))))))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -7.4e+128) tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y; elseif (b <= -5.2e+90) tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5))))))); elseif (b <= -1.1e-143) tmp = x * (1.0 / (y * a)); elseif (b <= 2e-184) tmp = (b * ((x / (a * b)) - (x / a))) / y; else tmp = (x / (a + (b * (a + (b * ((0.16666666666666666 * (a * b)) + (a * 0.5))))))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -7.4e+128], N[(N[(x * N[(N[(1.0 / a), $MachinePrecision] + N[(b * N[(N[(0.5 * N[(b / a), $MachinePrecision]), $MachinePrecision] + N[(-1.0 / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -5.2e+90], N[(x / N[(a * N[(y + N[(b * N[(y + N[(b * N[(N[(0.16666666666666666 * N[(y * b), $MachinePrecision]), $MachinePrecision] + N[(y * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.1e-143], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2e-184], N[(N[(b * N[(N[(x / N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a + N[(b * N[(a + N[(b * N[(N[(0.16666666666666666 * N[(a * b), $MachinePrecision]), $MachinePrecision] + N[(a * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -7.4 \cdot 10^{+128}:\\
\;\;\;\;\frac{x \cdot \left(\frac{1}{a} + b \cdot \left(0.5 \cdot \frac{b}{a} + \frac{-1}{a}\right)\right)}{y}\\
\mathbf{elif}\;b \leq -5.2 \cdot 10^{+90}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + b \cdot \left(y + b \cdot \left(0.16666666666666666 \cdot \left(y \cdot b\right) + y \cdot 0.5\right)\right)\right)}\\
\mathbf{elif}\;b \leq -1.1 \cdot 10^{-143}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{elif}\;b \leq 2 \cdot 10^{-184}:\\
\;\;\;\;\frac{b \cdot \left(\frac{x}{a \cdot b} - \frac{x}{a}\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + b \cdot \left(a + b \cdot \left(0.16666666666666666 \cdot \left(a \cdot b\right) + a \cdot 0.5\right)\right)}}{y}\\
\end{array}
\end{array}
if b < -7.4000000000000002e128Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 70.0%
Taylor expanded in x around 0 83.8%
if -7.4000000000000002e128 < b < -5.1999999999999997e90Initial program 100.0%
associate-/l*100.0%
exp-diff27.3%
associate-/l/27.3%
exp-sum27.3%
associate-/l*27.3%
*-commutative27.3%
exp-to-pow27.3%
*-commutative27.3%
exp-to-pow27.3%
sub-neg27.3%
metadata-eval27.3%
Simplified27.3%
Taylor expanded in y around 0 36.8%
exp-to-pow36.8%
sub-neg36.8%
metadata-eval36.8%
associate-/r*36.8%
Simplified36.8%
Taylor expanded in t around 0 37.4%
Taylor expanded in b around 0 64.2%
if -5.1999999999999997e90 < b < -1.09999999999999995e-143Initial program 95.4%
associate-/l*97.4%
exp-diff86.7%
associate-/l/86.7%
exp-sum80.4%
associate-/l*78.2%
*-commutative78.2%
exp-to-pow78.2%
*-commutative78.2%
exp-to-pow78.6%
sub-neg78.6%
metadata-eval78.6%
Simplified78.6%
Taylor expanded in y around 0 70.4%
exp-to-pow70.8%
sub-neg70.8%
metadata-eval70.8%
associate-/r*68.7%
Simplified68.7%
Taylor expanded in t around 0 52.8%
Taylor expanded in b around 0 32.2%
div-inv34.3%
Applied egg-rr34.3%
if -1.09999999999999995e-143 < b < 2.0000000000000001e-184Initial program 97.9%
Taylor expanded in y around 0 80.6%
div-exp80.6%
exp-to-pow81.1%
sub-neg81.1%
metadata-eval81.1%
Simplified81.1%
Taylor expanded in t around 0 39.7%
Taylor expanded in b around 0 37.9%
+-commutative37.9%
mul-1-neg37.9%
unsub-neg37.9%
associate-/l*34.3%
Simplified34.3%
Taylor expanded in b around inf 51.5%
if 2.0000000000000001e-184 < b Initial program 99.1%
Taylor expanded in y around 0 80.7%
div-exp76.9%
exp-to-pow77.6%
sub-neg77.6%
metadata-eval77.6%
Simplified77.6%
Taylor expanded in t around 0 69.3%
Taylor expanded in b around 0 63.8%
Final simplification58.5%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -3.5e+128)
(/ (* x (+ (/ 1.0 a) (* b (+ (* 0.5 (/ b a)) (/ -1.0 a))))) y)
(if (<= b -3.5e+91)
(/
x
(*
a
(+ y (* b (+ y (* b (+ (* 0.16666666666666666 (* y b)) (* y 0.5))))))))
(if (<= b -1.2e-142)
(* x (/ 1.0 (* y a)))
(if (<= b 2.4e-181)
(/ (* b (- (/ x (* a b)) (/ x a))) y)
(/ (/ x (+ a (* b (+ a (* 0.5 (* a b)))))) y))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3.5e+128) {
tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y;
} else if (b <= -3.5e+91) {
tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5)))))));
} else if (b <= -1.2e-142) {
tmp = x * (1.0 / (y * a));
} else if (b <= 2.4e-181) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (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 <= (-3.5d+128)) then
tmp = (x * ((1.0d0 / a) + (b * ((0.5d0 * (b / a)) + ((-1.0d0) / a))))) / y
else if (b <= (-3.5d+91)) then
tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666d0 * (y * b)) + (y * 0.5d0)))))))
else if (b <= (-1.2d-142)) then
tmp = x * (1.0d0 / (y * a))
else if (b <= 2.4d-181) then
tmp = (b * ((x / (a * b)) - (x / a))) / y
else
tmp = (x / (a + (b * (a + (0.5d0 * (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 <= -3.5e+128) {
tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y;
} else if (b <= -3.5e+91) {
tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5)))))));
} else if (b <= -1.2e-142) {
tmp = x * (1.0 / (y * a));
} else if (b <= 2.4e-181) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -3.5e+128: tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y elif b <= -3.5e+91: tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5))))))) elif b <= -1.2e-142: tmp = x * (1.0 / (y * a)) elif b <= 2.4e-181: tmp = (b * ((x / (a * b)) - (x / a))) / y else: tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3.5e+128) tmp = Float64(Float64(x * Float64(Float64(1.0 / a) + Float64(b * Float64(Float64(0.5 * Float64(b / a)) + Float64(-1.0 / a))))) / y); elseif (b <= -3.5e+91) tmp = Float64(x / Float64(a * Float64(y + Float64(b * Float64(y + Float64(b * Float64(Float64(0.16666666666666666 * Float64(y * b)) + Float64(y * 0.5)))))))); elseif (b <= -1.2e-142) tmp = Float64(x * Float64(1.0 / Float64(y * a))); elseif (b <= 2.4e-181) tmp = Float64(Float64(b * Float64(Float64(x / Float64(a * b)) - Float64(x / a))) / y); else tmp = Float64(Float64(x / Float64(a + Float64(b * Float64(a + Float64(0.5 * Float64(a * b)))))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -3.5e+128) tmp = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y; elseif (b <= -3.5e+91) tmp = x / (a * (y + (b * (y + (b * ((0.16666666666666666 * (y * b)) + (y * 0.5))))))); elseif (b <= -1.2e-142) tmp = x * (1.0 / (y * a)); elseif (b <= 2.4e-181) tmp = (b * ((x / (a * b)) - (x / a))) / y; else tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3.5e+128], N[(N[(x * N[(N[(1.0 / a), $MachinePrecision] + N[(b * N[(N[(0.5 * N[(b / a), $MachinePrecision]), $MachinePrecision] + N[(-1.0 / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -3.5e+91], N[(x / N[(a * N[(y + N[(b * N[(y + N[(b * N[(N[(0.16666666666666666 * N[(y * b), $MachinePrecision]), $MachinePrecision] + N[(y * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.2e-142], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.4e-181], N[(N[(b * N[(N[(x / N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a + N[(b * N[(a + N[(0.5 * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.5 \cdot 10^{+128}:\\
\;\;\;\;\frac{x \cdot \left(\frac{1}{a} + b \cdot \left(0.5 \cdot \frac{b}{a} + \frac{-1}{a}\right)\right)}{y}\\
\mathbf{elif}\;b \leq -3.5 \cdot 10^{+91}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + b \cdot \left(y + b \cdot \left(0.16666666666666666 \cdot \left(y \cdot b\right) + y \cdot 0.5\right)\right)\right)}\\
\mathbf{elif}\;b \leq -1.2 \cdot 10^{-142}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{elif}\;b \leq 2.4 \cdot 10^{-181}:\\
\;\;\;\;\frac{b \cdot \left(\frac{x}{a \cdot b} - \frac{x}{a}\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + b \cdot \left(a + 0.5 \cdot \left(a \cdot b\right)\right)}}{y}\\
\end{array}
\end{array}
if b < -3.49999999999999969e128Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 70.0%
Taylor expanded in x around 0 83.8%
if -3.49999999999999969e128 < b < -3.50000000000000001e91Initial program 100.0%
associate-/l*100.0%
exp-diff27.3%
associate-/l/27.3%
exp-sum27.3%
associate-/l*27.3%
*-commutative27.3%
exp-to-pow27.3%
*-commutative27.3%
exp-to-pow27.3%
sub-neg27.3%
metadata-eval27.3%
Simplified27.3%
Taylor expanded in y around 0 36.8%
exp-to-pow36.8%
sub-neg36.8%
metadata-eval36.8%
associate-/r*36.8%
Simplified36.8%
Taylor expanded in t around 0 37.4%
Taylor expanded in b around 0 64.2%
if -3.50000000000000001e91 < b < -1.19999999999999994e-142Initial program 95.4%
associate-/l*97.4%
exp-diff86.7%
associate-/l/86.7%
exp-sum80.4%
associate-/l*78.2%
*-commutative78.2%
exp-to-pow78.2%
*-commutative78.2%
exp-to-pow78.6%
sub-neg78.6%
metadata-eval78.6%
Simplified78.6%
Taylor expanded in y around 0 70.4%
exp-to-pow70.8%
sub-neg70.8%
metadata-eval70.8%
associate-/r*68.7%
Simplified68.7%
Taylor expanded in t around 0 52.8%
Taylor expanded in b around 0 32.2%
div-inv34.3%
Applied egg-rr34.3%
if -1.19999999999999994e-142 < b < 2.4000000000000001e-181Initial program 97.9%
Taylor expanded in y around 0 80.6%
div-exp80.6%
exp-to-pow81.1%
sub-neg81.1%
metadata-eval81.1%
Simplified81.1%
Taylor expanded in t around 0 39.7%
Taylor expanded in b around 0 37.9%
+-commutative37.9%
mul-1-neg37.9%
unsub-neg37.9%
associate-/l*34.3%
Simplified34.3%
Taylor expanded in b around inf 51.5%
if 2.4000000000000001e-181 < b Initial program 99.1%
Taylor expanded in y around 0 80.7%
div-exp76.9%
exp-to-pow77.6%
sub-neg77.6%
metadata-eval77.6%
Simplified77.6%
Taylor expanded in t around 0 69.3%
Taylor expanded in b around 0 60.3%
Final simplification57.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (+ (/ 1.0 a) (* b (+ (* 0.5 (/ b a)) (/ -1.0 a))))) y)))
(if (<= b -3.5e+128)
t_1
(if (<= b -1.35e+111)
(/ x (* a (* y (+ 1.0 (* b (+ 1.0 (* b 0.5)))))))
(if (<= b -1.35e-100)
t_1
(if (<= b 4.5e-181)
(/ (* b (- (/ x (* a b)) (/ x a))) y)
(/ (/ x (+ a (* b (+ a (* 0.5 (* a b)))))) y)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y;
double tmp;
if (b <= -3.5e+128) {
tmp = t_1;
} else if (b <= -1.35e+111) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else if (b <= -1.35e-100) {
tmp = t_1;
} else if (b <= 4.5e-181) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (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 * ((1.0d0 / a) + (b * ((0.5d0 * (b / a)) + ((-1.0d0) / a))))) / y
if (b <= (-3.5d+128)) then
tmp = t_1
else if (b <= (-1.35d+111)) then
tmp = x / (a * (y * (1.0d0 + (b * (1.0d0 + (b * 0.5d0))))))
else if (b <= (-1.35d-100)) then
tmp = t_1
else if (b <= 4.5d-181) then
tmp = (b * ((x / (a * b)) - (x / a))) / y
else
tmp = (x / (a + (b * (a + (0.5d0 * (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 * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y;
double tmp;
if (b <= -3.5e+128) {
tmp = t_1;
} else if (b <= -1.35e+111) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else if (b <= -1.35e-100) {
tmp = t_1;
} else if (b <= 4.5e-181) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y tmp = 0 if b <= -3.5e+128: tmp = t_1 elif b <= -1.35e+111: tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))) elif b <= -1.35e-100: tmp = t_1 elif b <= 4.5e-181: tmp = (b * ((x / (a * b)) - (x / a))) / y else: tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * Float64(Float64(1.0 / a) + Float64(b * Float64(Float64(0.5 * Float64(b / a)) + Float64(-1.0 / a))))) / y) tmp = 0.0 if (b <= -3.5e+128) tmp = t_1; elseif (b <= -1.35e+111) tmp = Float64(x / Float64(a * Float64(y * Float64(1.0 + Float64(b * Float64(1.0 + Float64(b * 0.5))))))); elseif (b <= -1.35e-100) tmp = t_1; elseif (b <= 4.5e-181) tmp = Float64(Float64(b * Float64(Float64(x / Float64(a * b)) - Float64(x / a))) / y); else tmp = Float64(Float64(x / Float64(a + Float64(b * Float64(a + Float64(0.5 * Float64(a * b)))))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * ((1.0 / a) + (b * ((0.5 * (b / a)) + (-1.0 / a))))) / y; tmp = 0.0; if (b <= -3.5e+128) tmp = t_1; elseif (b <= -1.35e+111) tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))); elseif (b <= -1.35e-100) tmp = t_1; elseif (b <= 4.5e-181) tmp = (b * ((x / (a * b)) - (x / a))) / y; else tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[(N[(1.0 / a), $MachinePrecision] + N[(b * N[(N[(0.5 * N[(b / a), $MachinePrecision]), $MachinePrecision] + N[(-1.0 / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, If[LessEqual[b, -3.5e+128], t$95$1, If[LessEqual[b, -1.35e+111], N[(x / N[(a * N[(y * N[(1.0 + N[(b * N[(1.0 + N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.35e-100], t$95$1, If[LessEqual[b, 4.5e-181], N[(N[(b * N[(N[(x / N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a + N[(b * N[(a + N[(0.5 * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot \left(\frac{1}{a} + b \cdot \left(0.5 \cdot \frac{b}{a} + \frac{-1}{a}\right)\right)}{y}\\
\mathbf{if}\;b \leq -3.5 \cdot 10^{+128}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -1.35 \cdot 10^{+111}:\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot \left(1 + b \cdot \left(1 + b \cdot 0.5\right)\right)\right)}\\
\mathbf{elif}\;b \leq -1.35 \cdot 10^{-100}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 4.5 \cdot 10^{-181}:\\
\;\;\;\;\frac{b \cdot \left(\frac{x}{a \cdot b} - \frac{x}{a}\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + b \cdot \left(a + 0.5 \cdot \left(a \cdot b\right)\right)}}{y}\\
\end{array}
\end{array}
if b < -3.49999999999999969e128 or -1.3499999999999999e111 < b < -1.35000000000000008e-100Initial program 99.8%
Taylor expanded in y around 0 87.2%
div-exp72.8%
exp-to-pow73.1%
sub-neg73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in t around 0 74.7%
Taylor expanded in b around 0 46.7%
Taylor expanded in x around 0 59.6%
if -3.49999999999999969e128 < b < -1.3499999999999999e111Initial program 100.0%
associate-/l*100.0%
exp-diff14.3%
associate-/l/14.3%
exp-sum14.3%
associate-/l*14.3%
*-commutative14.3%
exp-to-pow14.3%
*-commutative14.3%
exp-to-pow14.3%
sub-neg14.3%
metadata-eval14.3%
Simplified14.3%
Taylor expanded in y around 0 29.0%
exp-to-pow29.0%
sub-neg29.0%
metadata-eval29.0%
associate-/r*29.0%
Simplified29.0%
Taylor expanded in t around 0 29.7%
Taylor expanded in b around 0 71.9%
Taylor expanded in y around 0 71.9%
if -1.35000000000000008e-100 < b < 4.4999999999999999e-181Initial program 95.1%
Taylor expanded in y around 0 76.2%
div-exp76.2%
exp-to-pow76.7%
sub-neg76.7%
metadata-eval76.7%
Simplified76.7%
Taylor expanded in t around 0 37.2%
Taylor expanded in b around 0 35.7%
+-commutative35.7%
mul-1-neg35.7%
unsub-neg35.7%
associate-/l*32.7%
Simplified32.7%
Taylor expanded in b around inf 47.2%
if 4.4999999999999999e-181 < b Initial program 99.1%
Taylor expanded in y around 0 80.7%
div-exp76.9%
exp-to-pow77.6%
sub-neg77.6%
metadata-eval77.6%
Simplified77.6%
Taylor expanded in t around 0 69.3%
Taylor expanded in b around 0 60.3%
Final simplification57.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -3.5e+128)
(/ (+ (/ x a) (* b (/ (- (* x (* b 0.5)) x) a))) y)
(if (<= b -1.46e+109)
(/ x (* a (* y (+ 1.0 (* b (+ 1.0 (* b 0.5)))))))
(if (<= b -3.45e-158)
(* b (- (/ x (* a (* y b))) (/ x (* y a))))
(if (<= b 1.7e-188)
(/ (* b (- (/ x (* a b)) (/ x a))) y)
(/ (/ x (+ a (* b (+ a (* 0.5 (* a b)))))) y))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3.5e+128) {
tmp = ((x / a) + (b * (((x * (b * 0.5)) - x) / a))) / y;
} else if (b <= -1.46e+109) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else if (b <= -3.45e-158) {
tmp = b * ((x / (a * (y * b))) - (x / (y * a)));
} else if (b <= 1.7e-188) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (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 <= (-3.5d+128)) then
tmp = ((x / a) + (b * (((x * (b * 0.5d0)) - x) / a))) / y
else if (b <= (-1.46d+109)) then
tmp = x / (a * (y * (1.0d0 + (b * (1.0d0 + (b * 0.5d0))))))
else if (b <= (-3.45d-158)) then
tmp = b * ((x / (a * (y * b))) - (x / (y * a)))
else if (b <= 1.7d-188) then
tmp = (b * ((x / (a * b)) - (x / a))) / y
else
tmp = (x / (a + (b * (a + (0.5d0 * (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 <= -3.5e+128) {
tmp = ((x / a) + (b * (((x * (b * 0.5)) - x) / a))) / y;
} else if (b <= -1.46e+109) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else if (b <= -3.45e-158) {
tmp = b * ((x / (a * (y * b))) - (x / (y * a)));
} else if (b <= 1.7e-188) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -3.5e+128: tmp = ((x / a) + (b * (((x * (b * 0.5)) - x) / a))) / y elif b <= -1.46e+109: tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))) elif b <= -3.45e-158: tmp = b * ((x / (a * (y * b))) - (x / (y * a))) elif b <= 1.7e-188: tmp = (b * ((x / (a * b)) - (x / a))) / y else: tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3.5e+128) tmp = Float64(Float64(Float64(x / a) + Float64(b * Float64(Float64(Float64(x * Float64(b * 0.5)) - x) / a))) / y); elseif (b <= -1.46e+109) tmp = Float64(x / Float64(a * Float64(y * Float64(1.0 + Float64(b * Float64(1.0 + Float64(b * 0.5))))))); elseif (b <= -3.45e-158) tmp = Float64(b * Float64(Float64(x / Float64(a * Float64(y * b))) - Float64(x / Float64(y * a)))); elseif (b <= 1.7e-188) tmp = Float64(Float64(b * Float64(Float64(x / Float64(a * b)) - Float64(x / a))) / y); else tmp = Float64(Float64(x / Float64(a + Float64(b * Float64(a + Float64(0.5 * Float64(a * b)))))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -3.5e+128) tmp = ((x / a) + (b * (((x * (b * 0.5)) - x) / a))) / y; elseif (b <= -1.46e+109) tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))); elseif (b <= -3.45e-158) tmp = b * ((x / (a * (y * b))) - (x / (y * a))); elseif (b <= 1.7e-188) tmp = (b * ((x / (a * b)) - (x / a))) / y; else tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3.5e+128], N[(N[(N[(x / a), $MachinePrecision] + N[(b * N[(N[(N[(x * N[(b * 0.5), $MachinePrecision]), $MachinePrecision] - x), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -1.46e+109], N[(x / N[(a * N[(y * N[(1.0 + N[(b * N[(1.0 + N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -3.45e-158], N[(b * N[(N[(x / N[(a * N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.7e-188], N[(N[(b * N[(N[(x / N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a + N[(b * N[(a + N[(0.5 * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.5 \cdot 10^{+128}:\\
\;\;\;\;\frac{\frac{x}{a} + b \cdot \frac{x \cdot \left(b \cdot 0.5\right) - x}{a}}{y}\\
\mathbf{elif}\;b \leq -1.46 \cdot 10^{+109}:\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot \left(1 + b \cdot \left(1 + b \cdot 0.5\right)\right)\right)}\\
\mathbf{elif}\;b \leq -3.45 \cdot 10^{-158}:\\
\;\;\;\;b \cdot \left(\frac{x}{a \cdot \left(y \cdot b\right)} - \frac{x}{y \cdot a}\right)\\
\mathbf{elif}\;b \leq 1.7 \cdot 10^{-188}:\\
\;\;\;\;\frac{b \cdot \left(\frac{x}{a \cdot b} - \frac{x}{a}\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + b \cdot \left(a + 0.5 \cdot \left(a \cdot b\right)\right)}}{y}\\
\end{array}
\end{array}
if b < -3.49999999999999969e128Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 70.0%
Taylor expanded in b around 0 70.0%
Simplified75.9%
if -3.49999999999999969e128 < b < -1.46e109Initial program 100.0%
associate-/l*100.0%
exp-diff14.3%
associate-/l/14.3%
exp-sum14.3%
associate-/l*14.3%
*-commutative14.3%
exp-to-pow14.3%
*-commutative14.3%
exp-to-pow14.3%
sub-neg14.3%
metadata-eval14.3%
Simplified14.3%
Taylor expanded in y around 0 29.0%
exp-to-pow29.0%
sub-neg29.0%
metadata-eval29.0%
associate-/r*29.0%
Simplified29.0%
Taylor expanded in t around 0 29.7%
Taylor expanded in b around 0 71.9%
Taylor expanded in y around 0 71.9%
if -1.46e109 < b < -3.4499999999999998e-158Initial program 96.0%
Taylor expanded in y around 0 73.3%
div-exp66.0%
exp-to-pow66.5%
sub-neg66.5%
metadata-eval66.5%
Simplified66.5%
Taylor expanded in t around 0 47.3%
Taylor expanded in b around 0 29.8%
+-commutative29.8%
mul-1-neg29.8%
unsub-neg29.8%
associate-/l*29.8%
Simplified29.8%
Taylor expanded in b around inf 34.3%
if -3.4499999999999998e-158 < b < 1.70000000000000014e-188Initial program 97.7%
Taylor expanded in y around 0 81.0%
div-exp81.0%
exp-to-pow81.4%
sub-neg81.4%
metadata-eval81.4%
Simplified81.4%
Taylor expanded in t around 0 42.3%
Taylor expanded in b around 0 40.4%
+-commutative40.4%
mul-1-neg40.4%
unsub-neg40.4%
associate-/l*36.6%
Simplified36.6%
Taylor expanded in b around inf 53.1%
if 1.70000000000000014e-188 < b Initial program 99.1%
Taylor expanded in y around 0 80.7%
div-exp76.9%
exp-to-pow77.6%
sub-neg77.6%
metadata-eval77.6%
Simplified77.6%
Taylor expanded in t around 0 69.3%
Taylor expanded in b around 0 60.3%
Final simplification55.8%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -4.4e+128)
(/ (/ (* x b) a) (- y))
(if (or (<= b -1.2e+111) (not (<= b 1.95e-180)))
(/ x (* a (* y (+ 1.0 (* b (+ 1.0 (* b 0.5)))))))
(/ (* b (- (/ x (* a b)) (/ x a))) y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -4.4e+128) {
tmp = ((x * b) / a) / -y;
} else if ((b <= -1.2e+111) || !(b <= 1.95e-180)) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else {
tmp = (b * ((x / (a * b)) - (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 <= (-4.4d+128)) then
tmp = ((x * b) / a) / -y
else if ((b <= (-1.2d+111)) .or. (.not. (b <= 1.95d-180))) then
tmp = x / (a * (y * (1.0d0 + (b * (1.0d0 + (b * 0.5d0))))))
else
tmp = (b * ((x / (a * b)) - (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 <= -4.4e+128) {
tmp = ((x * b) / a) / -y;
} else if ((b <= -1.2e+111) || !(b <= 1.95e-180)) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -4.4e+128: tmp = ((x * b) / a) / -y elif (b <= -1.2e+111) or not (b <= 1.95e-180): tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))) else: tmp = (b * ((x / (a * b)) - (x / a))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -4.4e+128) tmp = Float64(Float64(Float64(x * b) / a) / Float64(-y)); elseif ((b <= -1.2e+111) || !(b <= 1.95e-180)) tmp = Float64(x / Float64(a * Float64(y * Float64(1.0 + Float64(b * Float64(1.0 + Float64(b * 0.5))))))); else tmp = Float64(Float64(b * Float64(Float64(x / Float64(a * b)) - Float64(x / a))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -4.4e+128) tmp = ((x * b) / a) / -y; elseif ((b <= -1.2e+111) || ~((b <= 1.95e-180))) tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))); else tmp = (b * ((x / (a * b)) - (x / a))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -4.4e+128], N[(N[(N[(x * b), $MachinePrecision] / a), $MachinePrecision] / (-y)), $MachinePrecision], If[Or[LessEqual[b, -1.2e+111], N[Not[LessEqual[b, 1.95e-180]], $MachinePrecision]], N[(x / N[(a * N[(y * N[(1.0 + N[(b * N[(1.0 + N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(b * N[(N[(x / N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -4.4 \cdot 10^{+128}:\\
\;\;\;\;\frac{\frac{x \cdot b}{a}}{-y}\\
\mathbf{elif}\;b \leq -1.2 \cdot 10^{+111} \lor \neg \left(b \leq 1.95 \cdot 10^{-180}\right):\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot \left(1 + b \cdot \left(1 + b \cdot 0.5\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{b \cdot \left(\frac{x}{a \cdot b} - \frac{x}{a}\right)}{y}\\
\end{array}
\end{array}
if b < -4.40000000000000033e128Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 57.4%
+-commutative57.4%
mul-1-neg57.4%
unsub-neg57.4%
associate-/l*46.7%
Simplified46.7%
Taylor expanded in b around inf 57.4%
if -4.40000000000000033e128 < b < -1.20000000000000003e111 or 1.9500000000000001e-180 < b Initial program 99.2%
associate-/l*98.5%
exp-diff79.9%
associate-/l/79.9%
exp-sum75.5%
associate-/l*73.7%
*-commutative73.7%
exp-to-pow73.7%
*-commutative73.7%
exp-to-pow74.3%
sub-neg74.3%
metadata-eval74.3%
Simplified74.3%
Taylor expanded in y around 0 73.2%
exp-to-pow73.8%
sub-neg73.8%
metadata-eval73.8%
associate-/r*70.3%
Simplified70.3%
Taylor expanded in t around 0 66.0%
Taylor expanded in b around 0 53.7%
Taylor expanded in y around 0 56.2%
if -1.20000000000000003e111 < b < 1.9500000000000001e-180Initial program 96.9%
Taylor expanded in y around 0 77.0%
div-exp73.3%
exp-to-pow73.7%
sub-neg73.7%
metadata-eval73.7%
Simplified73.7%
Taylor expanded in t around 0 44.9%
Taylor expanded in b around 0 35.0%
+-commutative35.0%
mul-1-neg35.0%
unsub-neg35.0%
associate-/l*33.1%
Simplified33.1%
Taylor expanded in b around inf 42.0%
Final simplification50.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -4.5e+133)
(/ (/ (* x b) a) (- y))
(if (<= b -5.5e+109)
(/ x (* a (* y (+ 1.0 (* b (+ 1.0 (* b 0.5)))))))
(if (<= b 1.65e-172)
(/ (* b (- (/ x (* a b)) (/ x a))) y)
(/ (/ x (+ a (* b (+ a (* 0.5 (* a b)))))) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -4.5e+133) {
tmp = ((x * b) / a) / -y;
} else if (b <= -5.5e+109) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else if (b <= 1.65e-172) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (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 <= (-4.5d+133)) then
tmp = ((x * b) / a) / -y
else if (b <= (-5.5d+109)) then
tmp = x / (a * (y * (1.0d0 + (b * (1.0d0 + (b * 0.5d0))))))
else if (b <= 1.65d-172) then
tmp = (b * ((x / (a * b)) - (x / a))) / y
else
tmp = (x / (a + (b * (a + (0.5d0 * (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 <= -4.5e+133) {
tmp = ((x * b) / a) / -y;
} else if (b <= -5.5e+109) {
tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5))))));
} else if (b <= 1.65e-172) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -4.5e+133: tmp = ((x * b) / a) / -y elif b <= -5.5e+109: tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))) elif b <= 1.65e-172: tmp = (b * ((x / (a * b)) - (x / a))) / y else: tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -4.5e+133) tmp = Float64(Float64(Float64(x * b) / a) / Float64(-y)); elseif (b <= -5.5e+109) tmp = Float64(x / Float64(a * Float64(y * Float64(1.0 + Float64(b * Float64(1.0 + Float64(b * 0.5))))))); elseif (b <= 1.65e-172) tmp = Float64(Float64(b * Float64(Float64(x / Float64(a * b)) - Float64(x / a))) / y); else tmp = Float64(Float64(x / Float64(a + Float64(b * Float64(a + Float64(0.5 * Float64(a * b)))))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -4.5e+133) tmp = ((x * b) / a) / -y; elseif (b <= -5.5e+109) tmp = x / (a * (y * (1.0 + (b * (1.0 + (b * 0.5)))))); elseif (b <= 1.65e-172) tmp = (b * ((x / (a * b)) - (x / a))) / y; else tmp = (x / (a + (b * (a + (0.5 * (a * b)))))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -4.5e+133], N[(N[(N[(x * b), $MachinePrecision] / a), $MachinePrecision] / (-y)), $MachinePrecision], If[LessEqual[b, -5.5e+109], N[(x / N[(a * N[(y * N[(1.0 + N[(b * N[(1.0 + N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.65e-172], N[(N[(b * N[(N[(x / N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a + N[(b * N[(a + N[(0.5 * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -4.5 \cdot 10^{+133}:\\
\;\;\;\;\frac{\frac{x \cdot b}{a}}{-y}\\
\mathbf{elif}\;b \leq -5.5 \cdot 10^{+109}:\\
\;\;\;\;\frac{x}{a \cdot \left(y \cdot \left(1 + b \cdot \left(1 + b \cdot 0.5\right)\right)\right)}\\
\mathbf{elif}\;b \leq 1.65 \cdot 10^{-172}:\\
\;\;\;\;\frac{b \cdot \left(\frac{x}{a \cdot b} - \frac{x}{a}\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + b \cdot \left(a + 0.5 \cdot \left(a \cdot b\right)\right)}}{y}\\
\end{array}
\end{array}
if b < -4.49999999999999985e133Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 57.4%
+-commutative57.4%
mul-1-neg57.4%
unsub-neg57.4%
associate-/l*46.7%
Simplified46.7%
Taylor expanded in b around inf 57.4%
if -4.49999999999999985e133 < b < -5.4999999999999998e109Initial program 100.0%
associate-/l*100.0%
exp-diff14.3%
associate-/l/14.3%
exp-sum14.3%
associate-/l*14.3%
*-commutative14.3%
exp-to-pow14.3%
*-commutative14.3%
exp-to-pow14.3%
sub-neg14.3%
metadata-eval14.3%
Simplified14.3%
Taylor expanded in y around 0 29.0%
exp-to-pow29.0%
sub-neg29.0%
metadata-eval29.0%
associate-/r*29.0%
Simplified29.0%
Taylor expanded in t around 0 29.7%
Taylor expanded in b around 0 71.9%
Taylor expanded in y around 0 71.9%
if -5.4999999999999998e109 < b < 1.65e-172Initial program 96.9%
Taylor expanded in y around 0 77.0%
div-exp73.3%
exp-to-pow73.7%
sub-neg73.7%
metadata-eval73.7%
Simplified73.7%
Taylor expanded in t around 0 44.9%
Taylor expanded in b around 0 35.0%
+-commutative35.0%
mul-1-neg35.0%
unsub-neg35.0%
associate-/l*33.1%
Simplified33.1%
Taylor expanded in b around inf 42.0%
if 1.65e-172 < b Initial program 99.1%
Taylor expanded in y around 0 80.7%
div-exp76.9%
exp-to-pow77.6%
sub-neg77.6%
metadata-eval77.6%
Simplified77.6%
Taylor expanded in t around 0 69.3%
Taylor expanded in b around 0 60.3%
Final simplification52.5%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -5.8e+131)
(/ (/ (* x b) a) (- y))
(if (<= b -9.5e-283)
(* x (/ 1.0 (* y a)))
(if (<= b 4.8e-190)
(/ (* b (- (/ x (* a b)) (/ x a))) y)
(/ (/ x (+ a (* a b))) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -5.8e+131) {
tmp = ((x * b) / a) / -y;
} else if (b <= -9.5e-283) {
tmp = x * (1.0 / (y * a));
} else if (b <= 4.8e-190) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} 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 <= (-5.8d+131)) then
tmp = ((x * b) / a) / -y
else if (b <= (-9.5d-283)) then
tmp = x * (1.0d0 / (y * a))
else if (b <= 4.8d-190) then
tmp = (b * ((x / (a * b)) - (x / a))) / y
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 <= -5.8e+131) {
tmp = ((x * b) / a) / -y;
} else if (b <= -9.5e-283) {
tmp = x * (1.0 / (y * a));
} else if (b <= 4.8e-190) {
tmp = (b * ((x / (a * b)) - (x / a))) / y;
} else {
tmp = (x / (a + (a * b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -5.8e+131: tmp = ((x * b) / a) / -y elif b <= -9.5e-283: tmp = x * (1.0 / (y * a)) elif b <= 4.8e-190: tmp = (b * ((x / (a * b)) - (x / a))) / y else: tmp = (x / (a + (a * b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -5.8e+131) tmp = Float64(Float64(Float64(x * b) / a) / Float64(-y)); elseif (b <= -9.5e-283) tmp = Float64(x * Float64(1.0 / Float64(y * a))); elseif (b <= 4.8e-190) tmp = Float64(Float64(b * Float64(Float64(x / Float64(a * b)) - Float64(x / a))) / y); 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 <= -5.8e+131) tmp = ((x * b) / a) / -y; elseif (b <= -9.5e-283) tmp = x * (1.0 / (y * a)); elseif (b <= 4.8e-190) tmp = (b * ((x / (a * b)) - (x / a))) / y; else tmp = (x / (a + (a * b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -5.8e+131], N[(N[(N[(x * b), $MachinePrecision] / a), $MachinePrecision] / (-y)), $MachinePrecision], If[LessEqual[b, -9.5e-283], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 4.8e-190], N[(N[(b * N[(N[(x / N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(x / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x / N[(a + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -5.8 \cdot 10^{+131}:\\
\;\;\;\;\frac{\frac{x \cdot b}{a}}{-y}\\
\mathbf{elif}\;b \leq -9.5 \cdot 10^{-283}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{elif}\;b \leq 4.8 \cdot 10^{-190}:\\
\;\;\;\;\frac{b \cdot \left(\frac{x}{a \cdot b} - \frac{x}{a}\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + a \cdot b}}{y}\\
\end{array}
\end{array}
if b < -5.8000000000000002e131Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 57.4%
+-commutative57.4%
mul-1-neg57.4%
unsub-neg57.4%
associate-/l*46.7%
Simplified46.7%
Taylor expanded in b around inf 57.4%
if -5.8000000000000002e131 < b < -9.49999999999999979e-283Initial program 96.2%
associate-/l*98.3%
exp-diff82.7%
associate-/l/82.7%
exp-sum75.4%
associate-/l*74.2%
*-commutative74.2%
exp-to-pow74.2%
*-commutative74.2%
exp-to-pow74.6%
sub-neg74.6%
metadata-eval74.6%
Simplified74.6%
Taylor expanded in y around 0 72.4%
exp-to-pow72.8%
sub-neg72.8%
metadata-eval72.8%
associate-/r*71.6%
Simplified71.6%
Taylor expanded in t around 0 48.7%
Taylor expanded in b around 0 34.9%
div-inv37.3%
Applied egg-rr37.3%
if -9.49999999999999979e-283 < b < 4.8000000000000001e-190Initial program 99.3%
Taylor expanded in y around 0 74.3%
div-exp74.3%
exp-to-pow75.0%
sub-neg75.0%
metadata-eval75.0%
Simplified75.0%
Taylor expanded in t around 0 40.7%
Taylor expanded in b around 0 37.4%
+-commutative37.4%
mul-1-neg37.4%
unsub-neg37.4%
associate-/l*31.0%
Simplified31.0%
Taylor expanded in b around inf 52.6%
if 4.8000000000000001e-190 < b Initial program 99.1%
Taylor expanded in y around 0 80.7%
div-exp76.9%
exp-to-pow77.6%
sub-neg77.6%
metadata-eval77.6%
Simplified77.6%
Taylor expanded in t around 0 69.3%
Taylor expanded in b around 0 49.0%
Final simplification46.8%
(FPCore (x y z t a b) :precision binary64 (if (<= b -2e+129) (/ (/ (* x b) a) (- y)) (if (<= b 4.8e-251) (* x (/ 1.0 (* 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 <= -2e+129) {
tmp = ((x * b) / a) / -y;
} else if (b <= 4.8e-251) {
tmp = x * (1.0 / (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 <= (-2d+129)) then
tmp = ((x * b) / a) / -y
else if (b <= 4.8d-251) then
tmp = x * (1.0d0 / (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 <= -2e+129) {
tmp = ((x * b) / a) / -y;
} else if (b <= 4.8e-251) {
tmp = x * (1.0 / (y * a));
} else {
tmp = (x / (a + (a * b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -2e+129: tmp = ((x * b) / a) / -y elif b <= 4.8e-251: tmp = x * (1.0 / (y * a)) else: tmp = (x / (a + (a * b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -2e+129) tmp = Float64(Float64(Float64(x * b) / a) / Float64(-y)); elseif (b <= 4.8e-251) tmp = Float64(x * Float64(1.0 / 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 <= -2e+129) tmp = ((x * b) / a) / -y; elseif (b <= 4.8e-251) tmp = x * (1.0 / (y * a)); else tmp = (x / (a + (a * b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -2e+129], N[(N[(N[(x * b), $MachinePrecision] / a), $MachinePrecision] / (-y)), $MachinePrecision], If[LessEqual[b, 4.8e-251], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $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 -2 \cdot 10^{+129}:\\
\;\;\;\;\frac{\frac{x \cdot b}{a}}{-y}\\
\mathbf{elif}\;b \leq 4.8 \cdot 10^{-251}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a + a \cdot b}}{y}\\
\end{array}
\end{array}
if b < -2e129Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 57.4%
+-commutative57.4%
mul-1-neg57.4%
unsub-neg57.4%
associate-/l*46.7%
Simplified46.7%
Taylor expanded in b around inf 57.4%
if -2e129 < b < 4.79999999999999992e-251Initial program 96.7%
associate-/l*98.4%
exp-diff85.7%
associate-/l/85.7%
exp-sum76.8%
associate-/l*75.9%
*-commutative75.9%
exp-to-pow75.9%
*-commutative75.9%
exp-to-pow76.4%
sub-neg76.4%
metadata-eval76.4%
Simplified76.4%
Taylor expanded in y around 0 72.5%
exp-to-pow73.1%
sub-neg73.1%
metadata-eval73.1%
associate-/r*72.1%
Simplified72.1%
Taylor expanded in t around 0 50.7%
Taylor expanded in b around 0 39.4%
div-inv41.3%
Applied egg-rr41.3%
if 4.79999999999999992e-251 < b Initial program 99.2%
Taylor expanded in y around 0 80.2%
div-exp76.8%
exp-to-pow77.5%
sub-neg77.5%
metadata-eval77.5%
Simplified77.5%
Taylor expanded in t around 0 64.1%
Taylor expanded in b around 0 45.9%
Final simplification45.7%
(FPCore (x y z t a b) :precision binary64 (if (<= b -5.8e+138) (/ (* x b) (* y (- a))) (if (<= b 4.3e-249) (* x (/ 1.0 (* y a))) (/ (/ x a) y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -5.8e+138) {
tmp = (x * b) / (y * -a);
} else if (b <= 4.3e-249) {
tmp = x * (1.0 / (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 <= (-5.8d+138)) then
tmp = (x * b) / (y * -a)
else if (b <= 4.3d-249) then
tmp = x * (1.0d0 / (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 <= -5.8e+138) {
tmp = (x * b) / (y * -a);
} else if (b <= 4.3e-249) {
tmp = x * (1.0 / (y * a));
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -5.8e+138: tmp = (x * b) / (y * -a) elif b <= 4.3e-249: tmp = x * (1.0 / (y * a)) else: tmp = (x / a) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -5.8e+138) tmp = Float64(Float64(x * b) / Float64(y * Float64(-a))); elseif (b <= 4.3e-249) tmp = Float64(x * Float64(1.0 / 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 <= -5.8e+138) tmp = (x * b) / (y * -a); elseif (b <= 4.3e-249) tmp = x * (1.0 / (y * a)); else tmp = (x / a) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -5.8e+138], N[(N[(x * b), $MachinePrecision] / N[(y * (-a)), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 4.3e-249], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -5.8 \cdot 10^{+138}:\\
\;\;\;\;\frac{x \cdot b}{y \cdot \left(-a\right)}\\
\mathbf{elif}\;b \leq 4.3 \cdot 10^{-249}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if b < -5.80000000000000019e138Initial program 100.0%
Taylor expanded in y around 0 97.2%
div-exp77.1%
exp-to-pow77.1%
sub-neg77.1%
metadata-eval77.1%
Simplified77.1%
Taylor expanded in t around 0 94.4%
Taylor expanded in b around 0 56.2%
+-commutative56.2%
mul-1-neg56.2%
unsub-neg56.2%
associate-/l*47.9%
Simplified47.9%
Taylor expanded in b around inf 50.7%
if -5.80000000000000019e138 < b < 4.3000000000000002e-249Initial program 96.7%
associate-/l*98.4%
exp-diff85.8%
associate-/l/85.8%
exp-sum77.1%
associate-/l*76.1%
*-commutative76.1%
exp-to-pow76.1%
*-commutative76.1%
exp-to-pow76.6%
sub-neg76.6%
metadata-eval76.6%
Simplified76.6%
Taylor expanded in y around 0 72.8%
exp-to-pow73.3%
sub-neg73.3%
metadata-eval73.3%
associate-/r*71.4%
Simplified71.4%
Taylor expanded in t around 0 51.1%
Taylor expanded in b around 0 39.1%
div-inv41.0%
Applied egg-rr41.0%
if 4.3000000000000002e-249 < b Initial program 99.2%
Taylor expanded in y around 0 80.2%
div-exp76.8%
exp-to-pow77.5%
sub-neg77.5%
metadata-eval77.5%
Simplified77.5%
Taylor expanded in t around 0 64.1%
Taylor expanded in b around 0 39.5%
Final simplification41.6%
(FPCore (x y z t a b) :precision binary64 (if (<= b -8e+136) (/ (* (/ x a) (- b)) y) (if (<= b 2.9e-244) (* x (/ 1.0 (* y a))) (/ (/ x a) y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -8e+136) {
tmp = ((x / a) * -b) / y;
} else if (b <= 2.9e-244) {
tmp = x * (1.0 / (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 <= (-8d+136)) then
tmp = ((x / a) * -b) / y
else if (b <= 2.9d-244) then
tmp = x * (1.0d0 / (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 <= -8e+136) {
tmp = ((x / a) * -b) / y;
} else if (b <= 2.9e-244) {
tmp = x * (1.0 / (y * a));
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -8e+136: tmp = ((x / a) * -b) / y elif b <= 2.9e-244: tmp = x * (1.0 / (y * a)) else: tmp = (x / a) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -8e+136) tmp = Float64(Float64(Float64(x / a) * Float64(-b)) / y); elseif (b <= 2.9e-244) tmp = Float64(x * Float64(1.0 / 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 <= -8e+136) tmp = ((x / a) * -b) / y; elseif (b <= 2.9e-244) tmp = x * (1.0 / (y * a)); else tmp = (x / a) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -8e+136], N[(N[(N[(x / a), $MachinePrecision] * (-b)), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, 2.9e-244], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -8 \cdot 10^{+136}:\\
\;\;\;\;\frac{\frac{x}{a} \cdot \left(-b\right)}{y}\\
\mathbf{elif}\;b \leq 2.9 \cdot 10^{-244}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if b < -8.00000000000000047e136Initial program 100.0%
Taylor expanded in y around 0 97.3%
div-exp77.8%
exp-to-pow77.8%
sub-neg77.8%
metadata-eval77.8%
Simplified77.8%
Taylor expanded in t around 0 94.5%
Taylor expanded in b around 0 57.4%
+-commutative57.4%
mul-1-neg57.4%
unsub-neg57.4%
associate-/l*46.7%
Simplified46.7%
Taylor expanded in b around inf 57.4%
neg-mul-157.4%
associate-*r/46.7%
*-commutative46.7%
distribute-rgt-neg-in46.7%
Simplified46.7%
if -8.00000000000000047e136 < b < 2.89999999999999996e-244Initial program 96.7%
associate-/l*98.4%
exp-diff85.7%
associate-/l/85.7%
exp-sum76.8%
associate-/l*75.9%
*-commutative75.9%
exp-to-pow75.9%
*-commutative75.9%
exp-to-pow76.4%
sub-neg76.4%
metadata-eval76.4%
Simplified76.4%
Taylor expanded in y around 0 72.5%
exp-to-pow73.1%
sub-neg73.1%
metadata-eval73.1%
associate-/r*72.1%
Simplified72.1%
Taylor expanded in t around 0 50.7%
Taylor expanded in b around 0 39.4%
div-inv41.3%
Applied egg-rr41.3%
if 2.89999999999999996e-244 < b Initial program 99.2%
Taylor expanded in y around 0 80.2%
div-exp76.8%
exp-to-pow77.5%
sub-neg77.5%
metadata-eval77.5%
Simplified77.5%
Taylor expanded in t around 0 64.1%
Taylor expanded in b around 0 39.5%
Final simplification41.2%
(FPCore (x y z t a b) :precision binary64 (if (<= b -6.4e+106) (* (/ x (* y a)) (- b)) (if (<= b 2.1e-245) (* x (/ 1.0 (* y a))) (/ (/ x a) y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -6.4e+106) {
tmp = (x / (y * a)) * -b;
} else if (b <= 2.1e-245) {
tmp = x * (1.0 / (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 <= (-6.4d+106)) then
tmp = (x / (y * a)) * -b
else if (b <= 2.1d-245) then
tmp = x * (1.0d0 / (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 <= -6.4e+106) {
tmp = (x / (y * a)) * -b;
} else if (b <= 2.1e-245) {
tmp = x * (1.0 / (y * a));
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -6.4e+106: tmp = (x / (y * a)) * -b elif b <= 2.1e-245: tmp = x * (1.0 / (y * a)) else: tmp = (x / a) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -6.4e+106) tmp = Float64(Float64(x / Float64(y * a)) * Float64(-b)); elseif (b <= 2.1e-245) tmp = Float64(x * Float64(1.0 / 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 <= -6.4e+106) tmp = (x / (y * a)) * -b; elseif (b <= 2.1e-245) tmp = x * (1.0 / (y * a)); else tmp = (x / a) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -6.4e+106], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] * (-b)), $MachinePrecision], If[LessEqual[b, 2.1e-245], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -6.4 \cdot 10^{+106}:\\
\;\;\;\;\frac{x}{y \cdot a} \cdot \left(-b\right)\\
\mathbf{elif}\;b \leq 2.1 \cdot 10^{-245}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if b < -6.3999999999999996e106Initial program 100.0%
Taylor expanded in y around 0 91.1%
div-exp70.5%
exp-to-pow70.5%
sub-neg70.5%
metadata-eval70.5%
Simplified70.5%
Taylor expanded in t around 0 84.3%
Taylor expanded in b around 0 49.6%
+-commutative49.6%
mul-1-neg49.6%
unsub-neg49.6%
associate-/l*40.9%
Simplified40.9%
Taylor expanded in b around inf 43.0%
mul-1-neg43.0%
associate-/l*38.9%
distribute-rgt-neg-in38.9%
distribute-frac-neg238.9%
distribute-rgt-neg-in38.9%
Simplified38.9%
if -6.3999999999999996e106 < b < 2.1000000000000001e-245Initial program 96.4%
associate-/l*98.3%
exp-diff90.8%
associate-/l/90.8%
exp-sum81.3%
associate-/l*80.2%
*-commutative80.2%
exp-to-pow80.2%
*-commutative80.2%
exp-to-pow80.8%
sub-neg80.8%
metadata-eval80.8%
Simplified80.8%
Taylor expanded in y around 0 75.5%
exp-to-pow76.1%
sub-neg76.1%
metadata-eval76.1%
associate-/r*75.0%
Simplified75.0%
Taylor expanded in t around 0 51.7%
Taylor expanded in b around 0 41.5%
div-inv43.5%
Applied egg-rr43.5%
if 2.1000000000000001e-245 < b Initial program 99.2%
Taylor expanded in y around 0 80.2%
div-exp76.8%
exp-to-pow77.5%
sub-neg77.5%
metadata-eval77.5%
Simplified77.5%
Taylor expanded in t around 0 64.1%
Taylor expanded in b around 0 39.5%
Final simplification40.9%
(FPCore (x y z t a b) :precision binary64 (if (<= b -3.9e-5) (/ (/ (* 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 <= -3.9e-5) {
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 <= (-3.9d-5)) 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 <= -3.9e-5) {
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 <= -3.9e-5: 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 <= -3.9e-5) tmp = Float64(Float64(Float64(x * b) / a) / Float64(-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.9e-5) 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, -3.9e-5], 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 -3.9 \cdot 10^{-5}:\\
\;\;\;\;\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 < -3.8999999999999999e-5Initial program 100.0%
Taylor expanded in y around 0 86.9%
div-exp67.4%
exp-to-pow67.4%
sub-neg67.4%
metadata-eval67.4%
Simplified67.4%
Taylor expanded in t around 0 78.1%
Taylor expanded in b around 0 40.9%
+-commutative40.9%
mul-1-neg40.9%
unsub-neg40.9%
associate-/l*35.2%
Simplified35.2%
Taylor expanded in b around inf 40.9%
if -3.8999999999999999e-5 < b Initial program 97.7%
associate-/l*98.2%
exp-diff90.3%
associate-/l/90.3%
exp-sum82.9%
associate-/l*81.3%
*-commutative81.3%
exp-to-pow81.3%
*-commutative81.3%
exp-to-pow82.0%
sub-neg82.0%
metadata-eval82.0%
Simplified82.0%
Taylor expanded in y around 0 77.1%
exp-to-pow77.7%
sub-neg77.7%
metadata-eval77.7%
associate-/r*75.6%
Simplified75.6%
Taylor expanded in t around 0 56.7%
Taylor expanded in b around 0 41.8%
distribute-lft-out44.5%
*-commutative44.5%
Simplified44.5%
Final simplification43.6%
(FPCore (x y z t a b) :precision binary64 (if (<= b -5.5e-5) (/ (* x b) (* y (- a))) (/ x (* a (+ y (* y b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -5.5e-5) {
tmp = (x * b) / (y * -a);
} 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 <= (-5.5d-5)) then
tmp = (x * b) / (y * -a)
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 <= -5.5e-5) {
tmp = (x * b) / (y * -a);
} else {
tmp = x / (a * (y + (y * b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -5.5e-5: tmp = (x * b) / (y * -a) else: tmp = x / (a * (y + (y * b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -5.5e-5) tmp = Float64(Float64(x * b) / Float64(y * Float64(-a))); 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 <= -5.5e-5) tmp = (x * b) / (y * -a); else tmp = x / (a * (y + (y * b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -5.5e-5], N[(N[(x * b), $MachinePrecision] / N[(y * (-a)), $MachinePrecision]), $MachinePrecision], N[(x / N[(a * N[(y + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -5.5 \cdot 10^{-5}:\\
\;\;\;\;\frac{x \cdot b}{y \cdot \left(-a\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot b\right)}\\
\end{array}
\end{array}
if b < -5.5000000000000002e-5Initial program 100.0%
Taylor expanded in y around 0 86.9%
div-exp67.4%
exp-to-pow67.4%
sub-neg67.4%
metadata-eval67.4%
Simplified67.4%
Taylor expanded in t around 0 78.1%
Taylor expanded in b around 0 40.9%
+-commutative40.9%
mul-1-neg40.9%
unsub-neg40.9%
associate-/l*35.2%
Simplified35.2%
Taylor expanded in b around inf 36.6%
if -5.5000000000000002e-5 < b Initial program 97.7%
associate-/l*98.2%
exp-diff90.3%
associate-/l/90.3%
exp-sum82.9%
associate-/l*81.3%
*-commutative81.3%
exp-to-pow81.3%
*-commutative81.3%
exp-to-pow82.0%
sub-neg82.0%
metadata-eval82.0%
Simplified82.0%
Taylor expanded in y around 0 77.1%
exp-to-pow77.7%
sub-neg77.7%
metadata-eval77.7%
associate-/r*75.6%
Simplified75.6%
Taylor expanded in t around 0 56.7%
Taylor expanded in b around 0 41.8%
distribute-lft-out44.5%
*-commutative44.5%
Simplified44.5%
Final simplification42.4%
(FPCore (x y z t a b) :precision binary64 (if (<= b 1.75e-245) (* x (/ 1.0 (* y a))) (/ (/ x a) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 1.75e-245) {
tmp = x * (1.0 / (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 <= 1.75d-245) then
tmp = x * (1.0d0 / (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 <= 1.75e-245) {
tmp = x * (1.0 / (y * a));
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 1.75e-245: tmp = x * (1.0 / (y * a)) else: tmp = (x / a) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 1.75e-245) tmp = Float64(x * Float64(1.0 / 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 <= 1.75e-245) tmp = x * (1.0 / (y * a)); else tmp = (x / a) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 1.75e-245], N[(x * N[(1.0 / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 1.75 \cdot 10^{-245}:\\
\;\;\;\;x \cdot \frac{1}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if b < 1.75000000000000008e-245Initial program 97.6%
associate-/l*98.8%
exp-diff82.2%
associate-/l/82.2%
exp-sum72.0%
associate-/l*71.3%
*-commutative71.3%
exp-to-pow71.3%
*-commutative71.3%
exp-to-pow71.7%
sub-neg71.7%
metadata-eval71.7%
Simplified71.7%
Taylor expanded in y around 0 73.9%
exp-to-pow74.3%
sub-neg74.3%
metadata-eval74.3%
associate-/r*72.1%
Simplified72.1%
Taylor expanded in t around 0 62.1%
Taylor expanded in b around 0 35.6%
div-inv37.0%
Applied egg-rr37.0%
if 1.75000000000000008e-245 < b Initial program 99.2%
Taylor expanded in y around 0 80.2%
div-exp76.8%
exp-to-pow77.5%
sub-neg77.5%
metadata-eval77.5%
Simplified77.5%
Taylor expanded in t around 0 64.1%
Taylor expanded in b around 0 39.5%
Final simplification38.2%
(FPCore (x y z t a b) :precision binary64 (if (<= b 3.75e+49) (/ x (* y a)) (/ (/ x a) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 3.75e+49) {
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 (b <= 3.75d+49) 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 (b <= 3.75e+49) {
tmp = x / (y * a);
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 3.75e+49: tmp = x / (y * a) else: tmp = (x / a) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 3.75e+49) 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 (b <= 3.75e+49) tmp = x / (y * a); else tmp = (x / a) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 3.75e+49], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(N[(x / a), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 3.75 \cdot 10^{+49}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if b < 3.7499999999999998e49Initial program 97.9%
associate-/l*98.4%
exp-diff86.7%
associate-/l/86.7%
exp-sum77.0%
associate-/l*75.6%
*-commutative75.6%
exp-to-pow75.6%
*-commutative75.6%
exp-to-pow76.2%
sub-neg76.2%
metadata-eval76.2%
Simplified76.2%
Taylor expanded in y around 0 73.2%
exp-to-pow73.8%
sub-neg73.8%
metadata-eval73.8%
associate-/r*72.3%
Simplified72.3%
Taylor expanded in t around 0 56.5%
Taylor expanded in b around 0 37.3%
if 3.7499999999999998e49 < b Initial program 100.0%
Taylor expanded in y around 0 88.2%
div-exp80.2%
exp-to-pow80.2%
sub-neg80.2%
metadata-eval80.2%
Simplified80.2%
Taylor expanded in t around 0 86.2%
Taylor expanded in b around 0 37.7%
Final simplification37.4%
(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 98.3%
associate-/l*98.7%
exp-diff83.9%
associate-/l/83.9%
exp-sum75.7%
associate-/l*74.5%
*-commutative74.5%
exp-to-pow74.5%
*-commutative74.5%
exp-to-pow75.0%
sub-neg75.0%
metadata-eval75.0%
Simplified75.0%
Taylor expanded in y around 0 74.5%
exp-to-pow75.0%
sub-neg75.0%
metadata-eval75.0%
associate-/r*72.3%
Simplified72.3%
Taylor expanded in t around 0 62.3%
Taylor expanded in b around 0 34.8%
Final simplification34.8%
(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 2024096
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:precision binary64
:alt
(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))