
(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 22 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.7%
Final simplification98.7%
(FPCore (x y z t a b) :precision binary64 (if (or (<= (+ t -1.0) -2e+100) (not (<= (+ t -1.0) -1.0))) (/ (* x (exp (- (* (+ t -1.0) (log a)) b))) y) (/ (* x (exp (- (- (* y (log z)) (log a)) b))) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (((t + -1.0) <= -2e+100) || !((t + -1.0) <= -1.0)) {
tmp = (x * exp((((t + -1.0) * log(a)) - b))) / y;
} else {
tmp = (x * exp((((y * log(z)) - 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 (((t + (-1.0d0)) <= (-2d+100)) .or. (.not. ((t + (-1.0d0)) <= (-1.0d0)))) then
tmp = (x * exp((((t + (-1.0d0)) * log(a)) - b))) / y
else
tmp = (x * exp((((y * log(z)) - 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 (((t + -1.0) <= -2e+100) || !((t + -1.0) <= -1.0)) {
tmp = (x * Math.exp((((t + -1.0) * Math.log(a)) - b))) / y;
} else {
tmp = (x * Math.exp((((y * Math.log(z)) - Math.log(a)) - b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if ((t + -1.0) <= -2e+100) or not ((t + -1.0) <= -1.0): tmp = (x * math.exp((((t + -1.0) * math.log(a)) - b))) / y else: tmp = (x * math.exp((((y * math.log(z)) - math.log(a)) - b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((Float64(t + -1.0) <= -2e+100) || !(Float64(t + -1.0) <= -1.0)) tmp = Float64(Float64(x * exp(Float64(Float64(Float64(t + -1.0) * log(a)) - b))) / y); else tmp = Float64(Float64(x * exp(Float64(Float64(Float64(y * log(z)) - log(a)) - b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (((t + -1.0) <= -2e+100) || ~(((t + -1.0) <= -1.0))) tmp = (x * exp((((t + -1.0) * log(a)) - b))) / y; else tmp = (x * exp((((y * log(z)) - log(a)) - b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[N[(t + -1.0), $MachinePrecision], -2e+100], N[Not[LessEqual[N[(t + -1.0), $MachinePrecision], -1.0]], $MachinePrecision]], N[(N[(x * N[Exp[N[(N[(N[(t + -1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(N[(x * N[Exp[N[(N[(N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision] - N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t + -1 \leq -2 \cdot 10^{+100} \lor \neg \left(t + -1 \leq -1\right):\\
\;\;\;\;\frac{x \cdot e^{\left(t + -1\right) \cdot \log a - b}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot e^{\left(y \cdot \log z - \log a\right) - b}}{y}\\
\end{array}
\end{array}
if (-.f64 t 1) < -2.00000000000000003e100 or -1 < (-.f64 t 1) Initial program 100.0%
Taylor expanded in y around 0 96.2%
if -2.00000000000000003e100 < (-.f64 t 1) < -1Initial program 97.8%
Taylor expanded in t around 0 95.9%
+-commutative95.9%
mul-1-neg95.9%
unsub-neg95.9%
Simplified95.9%
Final simplification96.0%
(FPCore (x y z t a b) :precision binary64 (if (or (<= y -1.3e+229) (not (<= y 9.4e+138))) (/ x (/ a (/ (pow z y) y))) (/ (* x (exp (- (* (+ t -1.0) (log a)) b))) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((y <= -1.3e+229) || !(y <= 9.4e+138)) {
tmp = x / (a / (pow(z, y) / y));
} else {
tmp = (x * exp((((t + -1.0) * log(a)) - b))) / y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((y <= (-1.3d+229)) .or. (.not. (y <= 9.4d+138))) then
tmp = x / (a / ((z ** y) / y))
else
tmp = (x * exp((((t + (-1.0d0)) * log(a)) - b))) / y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((y <= -1.3e+229) || !(y <= 9.4e+138)) {
tmp = x / (a / (Math.pow(z, y) / y));
} else {
tmp = (x * Math.exp((((t + -1.0) * Math.log(a)) - b))) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (y <= -1.3e+229) or not (y <= 9.4e+138): tmp = x / (a / (math.pow(z, y) / y)) else: tmp = (x * math.exp((((t + -1.0) * math.log(a)) - b))) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((y <= -1.3e+229) || !(y <= 9.4e+138)) tmp = Float64(x / Float64(a / Float64((z ^ y) / y))); else tmp = Float64(Float64(x * exp(Float64(Float64(Float64(t + -1.0) * log(a)) - b))) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((y <= -1.3e+229) || ~((y <= 9.4e+138))) tmp = x / (a / ((z ^ y) / y)); else tmp = (x * exp((((t + -1.0) * log(a)) - b))) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -1.3e+229], N[Not[LessEqual[y, 9.4e+138]], $MachinePrecision]], N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x * N[Exp[N[(N[(N[(t + -1.0), $MachinePrecision] * N[Log[a], $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.3 \cdot 10^{+229} \lor \neg \left(y \leq 9.4 \cdot 10^{+138}\right):\\
\;\;\;\;\frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot e^{\left(t + -1\right) \cdot \log a - b}}{y}\\
\end{array}
\end{array}
if y < -1.3e229 or 9.3999999999999996e138 < y Initial program 100.0%
associate-*l/78.3%
*-commutative78.3%
exp-diff56.5%
+-commutative56.5%
exp-sum50.0%
associate-/l*50.0%
*-commutative50.0%
exp-to-pow50.0%
sub-neg50.0%
metadata-eval50.0%
*-commutative50.0%
exp-to-pow50.0%
Simplified50.0%
Taylor expanded in b around 0 84.9%
associate-/l*84.9%
*-commutative84.9%
exp-to-pow84.9%
*-commutative84.9%
exp-sum95.7%
exp-sum84.9%
*-commutative84.9%
exp-to-pow84.9%
exp-to-pow84.9%
sub-neg84.9%
metadata-eval84.9%
Simplified84.9%
Taylor expanded in t around 0 80.5%
associate-/l*93.6%
Simplified93.6%
if -1.3e229 < y < 9.3999999999999996e138Initial program 98.4%
Taylor expanded in y around 0 91.2%
Final simplification91.7%
(FPCore (x y z t a b) :precision binary64 (if (or (<= y -1.3e+229) (not (<= y 5.4e+123))) (/ x (/ a (/ (pow z y) y))) (/ x (/ y (/ (pow a (+ t -1.0)) (exp b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((y <= -1.3e+229) || !(y <= 5.4e+123)) {
tmp = x / (a / (pow(z, y) / y));
} else {
tmp = x / (y / (pow(a, (t + -1.0)) / 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 <= (-1.3d+229)) .or. (.not. (y <= 5.4d+123))) then
tmp = x / (a / ((z ** y) / y))
else
tmp = x / (y / ((a ** (t + (-1.0d0))) / 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 <= -1.3e+229) || !(y <= 5.4e+123)) {
tmp = x / (a / (Math.pow(z, y) / y));
} else {
tmp = x / (y / (Math.pow(a, (t + -1.0)) / Math.exp(b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (y <= -1.3e+229) or not (y <= 5.4e+123): tmp = x / (a / (math.pow(z, y) / y)) else: tmp = x / (y / (math.pow(a, (t + -1.0)) / math.exp(b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((y <= -1.3e+229) || !(y <= 5.4e+123)) tmp = Float64(x / Float64(a / Float64((z ^ y) / y))); else tmp = Float64(x / Float64(y / Float64((a ^ Float64(t + -1.0)) / exp(b)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((y <= -1.3e+229) || ~((y <= 5.4e+123))) tmp = x / (a / ((z ^ y) / y)); else tmp = x / (y / ((a ^ (t + -1.0)) / exp(b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -1.3e+229], N[Not[LessEqual[y, 5.4e+123]], $MachinePrecision]], N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(y / N[(N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision] / N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.3 \cdot 10^{+229} \lor \neg \left(y \leq 5.4 \cdot 10^{+123}\right):\\
\;\;\;\;\frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{y}{\frac{{a}^{\left(t + -1\right)}}{e^{b}}}}\\
\end{array}
\end{array}
if y < -1.3e229 or 5.40000000000000026e123 < y Initial program 100.0%
associate-*l/79.2%
*-commutative79.2%
exp-diff56.3%
+-commutative56.3%
exp-sum50.0%
associate-/l*50.0%
*-commutative50.0%
exp-to-pow50.0%
sub-neg50.0%
metadata-eval50.0%
*-commutative50.0%
exp-to-pow50.0%
Simplified50.0%
Taylor expanded in b around 0 83.4%
associate-/l*83.4%
*-commutative83.4%
exp-to-pow83.4%
*-commutative83.4%
exp-sum93.8%
exp-sum83.4%
*-commutative83.4%
exp-to-pow83.4%
exp-to-pow83.4%
sub-neg83.4%
metadata-eval83.4%
Simplified83.4%
Taylor expanded in t around 0 79.3%
associate-/l*91.8%
Simplified91.8%
if -1.3e229 < y < 5.40000000000000026e123Initial program 98.4%
Taylor expanded in y around 0 91.6%
associate-/l*92.0%
div-exp78.6%
exp-to-pow79.2%
sub-neg79.2%
metadata-eval79.2%
Simplified79.2%
Final simplification81.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (pow a (+ t -1.0))) y)) (t_2 (* y (exp b))))
(if (<= t -2.8e+122)
t_1
(if (<= t 3.6e-72)
(* (/ x a) (/ (pow z y) t_2))
(if (<= t 1.05e-31)
(/ x (/ a (/ (pow z y) y)))
(if (<= t 3.1e-17) (/ x t_2) t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * pow(a, (t + -1.0))) / y;
double t_2 = y * exp(b);
double tmp;
if (t <= -2.8e+122) {
tmp = t_1;
} else if (t <= 3.6e-72) {
tmp = (x / a) * (pow(z, y) / t_2);
} else if (t <= 1.05e-31) {
tmp = x / (a / (pow(z, y) / y));
} else if (t <= 3.1e-17) {
tmp = x / t_2;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (x * (a ** (t + (-1.0d0)))) / y
t_2 = y * exp(b)
if (t <= (-2.8d+122)) then
tmp = t_1
else if (t <= 3.6d-72) then
tmp = (x / a) * ((z ** y) / t_2)
else if (t <= 1.05d-31) then
tmp = x / (a / ((z ** y) / y))
else if (t <= 3.1d-17) then
tmp = x / t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * Math.pow(a, (t + -1.0))) / y;
double t_2 = y * Math.exp(b);
double tmp;
if (t <= -2.8e+122) {
tmp = t_1;
} else if (t <= 3.6e-72) {
tmp = (x / a) * (Math.pow(z, y) / t_2);
} else if (t <= 1.05e-31) {
tmp = x / (a / (Math.pow(z, y) / y));
} else if (t <= 3.1e-17) {
tmp = x / t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * math.pow(a, (t + -1.0))) / y t_2 = y * math.exp(b) tmp = 0 if t <= -2.8e+122: tmp = t_1 elif t <= 3.6e-72: tmp = (x / a) * (math.pow(z, y) / t_2) elif t <= 1.05e-31: tmp = x / (a / (math.pow(z, y) / y)) elif t <= 3.1e-17: tmp = x / t_2 else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * (a ^ Float64(t + -1.0))) / y) t_2 = Float64(y * exp(b)) tmp = 0.0 if (t <= -2.8e+122) tmp = t_1; elseif (t <= 3.6e-72) tmp = Float64(Float64(x / a) * Float64((z ^ y) / t_2)); elseif (t <= 1.05e-31) tmp = Float64(x / Float64(a / Float64((z ^ y) / y))); elseif (t <= 3.1e-17) tmp = Float64(x / t_2); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * (a ^ (t + -1.0))) / y; t_2 = y * exp(b); tmp = 0.0; if (t <= -2.8e+122) tmp = t_1; elseif (t <= 3.6e-72) tmp = (x / a) * ((z ^ y) / t_2); elseif (t <= 1.05e-31) tmp = x / (a / ((z ^ y) / y)); elseif (t <= 3.1e-17) tmp = x / t_2; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, Block[{t$95$2 = N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.8e+122], t$95$1, If[LessEqual[t, 3.6e-72], N[(N[(x / a), $MachinePrecision] * N[(N[Power[z, y], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.05e-31], N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.1e-17], N[(x / t$95$2), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot {a}^{\left(t + -1\right)}}{y}\\
t_2 := y \cdot e^{b}\\
\mathbf{if}\;t \leq -2.8 \cdot 10^{+122}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 3.6 \cdot 10^{-72}:\\
\;\;\;\;\frac{x}{a} \cdot \frac{{z}^{y}}{t\_2}\\
\mathbf{elif}\;t \leq 1.05 \cdot 10^{-31}:\\
\;\;\;\;\frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\mathbf{elif}\;t \leq 3.1 \cdot 10^{-17}:\\
\;\;\;\;\frac{x}{t\_2}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -2.8e122 or 3.0999999999999998e-17 < t Initial program 100.0%
Taylor expanded in y around 0 95.9%
Taylor expanded in b around 0 87.8%
*-commutative87.8%
exp-to-pow87.8%
sub-neg87.8%
metadata-eval87.8%
+-commutative87.8%
Simplified87.8%
if -2.8e122 < t < 3.6e-72Initial program 97.8%
associate-*l/87.5%
*-commutative87.5%
exp-diff69.7%
+-commutative69.7%
exp-sum66.3%
associate-/l*66.3%
*-commutative66.3%
exp-to-pow67.1%
sub-neg67.1%
metadata-eval67.1%
*-commutative67.1%
exp-to-pow67.1%
Simplified67.1%
Taylor expanded in t around 0 78.2%
times-frac78.2%
Simplified78.2%
if 3.6e-72 < t < 1.04999999999999996e-31Initial program 98.6%
associate-*l/88.6%
*-commutative88.6%
exp-diff58.6%
+-commutative58.6%
exp-sum58.6%
associate-/l*58.6%
*-commutative58.6%
exp-to-pow59.8%
sub-neg59.8%
metadata-eval59.8%
*-commutative59.8%
exp-to-pow59.8%
Simplified59.8%
Taylor expanded in b around 0 69.3%
associate-/l*78.9%
*-commutative78.9%
exp-to-pow78.9%
*-commutative78.9%
exp-sum78.9%
exp-sum78.9%
*-commutative78.9%
exp-to-pow78.9%
exp-to-pow80.3%
sub-neg80.3%
metadata-eval80.3%
Simplified80.3%
Taylor expanded in t around 0 60.3%
associate-/l*80.3%
Simplified80.3%
if 1.04999999999999996e-31 < t < 3.0999999999999998e-17Initial program 100.0%
Taylor expanded in y around 0 100.0%
Taylor expanded in b around inf 100.0%
neg-mul-1100.0%
Simplified100.0%
associate-/l*100.0%
div-inv100.0%
div-inv100.0%
add-sqr-sqrt66.7%
sqrt-unprod67.2%
sqr-neg67.2%
sqrt-unprod0.5%
add-sqr-sqrt1.6%
exp-neg1.6%
add-sqr-sqrt1.0%
sqrt-unprod34.4%
sqr-neg34.4%
sqrt-unprod33.3%
add-sqr-sqrt100.0%
Applied egg-rr100.0%
associate-*r/100.0%
*-rgt-identity100.0%
Simplified100.0%
Final simplification82.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ (* x (pow a (+ t -1.0))) y)) (t_2 (/ x (* y (exp b)))))
(if (<= b -2.45e+182)
t_2
(if (<= b -5e+106)
t_1
(if (<= b -220000.0)
t_2
(if (<= b -6e-255)
t_1
(if (<= b 4.1e-165)
(/ x (/ a (/ (pow z y) y)))
(if (<= b 2.1e+28) t_1 t_2))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * pow(a, (t + -1.0))) / y;
double t_2 = x / (y * exp(b));
double tmp;
if (b <= -2.45e+182) {
tmp = t_2;
} else if (b <= -5e+106) {
tmp = t_1;
} else if (b <= -220000.0) {
tmp = t_2;
} else if (b <= -6e-255) {
tmp = t_1;
} else if (b <= 4.1e-165) {
tmp = x / (a / (pow(z, y) / y));
} else if (b <= 2.1e+28) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (x * (a ** (t + (-1.0d0)))) / y
t_2 = x / (y * exp(b))
if (b <= (-2.45d+182)) then
tmp = t_2
else if (b <= (-5d+106)) then
tmp = t_1
else if (b <= (-220000.0d0)) then
tmp = t_2
else if (b <= (-6d-255)) then
tmp = t_1
else if (b <= 4.1d-165) then
tmp = x / (a / ((z ** y) / y))
else if (b <= 2.1d+28) then
tmp = t_1
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (x * Math.pow(a, (t + -1.0))) / y;
double t_2 = x / (y * Math.exp(b));
double tmp;
if (b <= -2.45e+182) {
tmp = t_2;
} else if (b <= -5e+106) {
tmp = t_1;
} else if (b <= -220000.0) {
tmp = t_2;
} else if (b <= -6e-255) {
tmp = t_1;
} else if (b <= 4.1e-165) {
tmp = x / (a / (Math.pow(z, y) / y));
} else if (b <= 2.1e+28) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (x * math.pow(a, (t + -1.0))) / y t_2 = x / (y * math.exp(b)) tmp = 0 if b <= -2.45e+182: tmp = t_2 elif b <= -5e+106: tmp = t_1 elif b <= -220000.0: tmp = t_2 elif b <= -6e-255: tmp = t_1 elif b <= 4.1e-165: tmp = x / (a / (math.pow(z, y) / y)) elif b <= 2.1e+28: tmp = t_1 else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(x * (a ^ Float64(t + -1.0))) / y) t_2 = Float64(x / Float64(y * exp(b))) tmp = 0.0 if (b <= -2.45e+182) tmp = t_2; elseif (b <= -5e+106) tmp = t_1; elseif (b <= -220000.0) tmp = t_2; elseif (b <= -6e-255) tmp = t_1; elseif (b <= 4.1e-165) tmp = Float64(x / Float64(a / Float64((z ^ y) / y))); elseif (b <= 2.1e+28) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (x * (a ^ (t + -1.0))) / y; t_2 = x / (y * exp(b)); tmp = 0.0; if (b <= -2.45e+182) tmp = t_2; elseif (b <= -5e+106) tmp = t_1; elseif (b <= -220000.0) tmp = t_2; elseif (b <= -6e-255) tmp = t_1; elseif (b <= 4.1e-165) tmp = x / (a / ((z ^ y) / y)); elseif (b <= 2.1e+28) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(x * N[Power[a, N[(t + -1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]}, Block[{t$95$2 = N[(x / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2.45e+182], t$95$2, If[LessEqual[b, -5e+106], t$95$1, If[LessEqual[b, -220000.0], t$95$2, If[LessEqual[b, -6e-255], t$95$1, If[LessEqual[b, 4.1e-165], N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.1e+28], t$95$1, t$95$2]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x \cdot {a}^{\left(t + -1\right)}}{y}\\
t_2 := \frac{x}{y \cdot e^{b}}\\
\mathbf{if}\;b \leq -2.45 \cdot 10^{+182}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;b \leq -5 \cdot 10^{+106}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -220000:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;b \leq -6 \cdot 10^{-255}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 4.1 \cdot 10^{-165}:\\
\;\;\;\;\frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\mathbf{elif}\;b \leq 2.1 \cdot 10^{+28}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if b < -2.45e182 or -4.9999999999999998e106 < b < -2.2e5 or 2.09999999999999989e28 < b Initial program 100.0%
Taylor expanded in y around 0 93.6%
Taylor expanded in b around inf 89.1%
neg-mul-189.1%
Simplified89.1%
associate-/l*89.1%
div-inv89.1%
div-inv89.1%
add-sqr-sqrt44.5%
sqrt-unprod49.9%
sqr-neg49.9%
sqrt-unprod5.3%
add-sqr-sqrt12.5%
exp-neg12.5%
add-sqr-sqrt7.2%
sqrt-unprod51.7%
sqr-neg51.7%
sqrt-unprod44.5%
add-sqr-sqrt89.1%
Applied egg-rr89.1%
associate-*r/89.1%
*-rgt-identity89.1%
Simplified89.1%
if -2.45e182 < b < -4.9999999999999998e106 or -2.2e5 < b < -6.00000000000000004e-255 or 4.1000000000000002e-165 < b < 2.09999999999999989e28Initial program 98.9%
Taylor expanded in y around 0 79.6%
Taylor expanded in b around 0 77.5%
*-commutative77.5%
exp-to-pow78.2%
sub-neg78.2%
metadata-eval78.2%
+-commutative78.2%
Simplified78.2%
if -6.00000000000000004e-255 < b < 4.1000000000000002e-165Initial program 94.5%
associate-*l/88.7%
*-commutative88.7%
exp-diff88.7%
+-commutative88.7%
exp-sum81.2%
associate-/l*81.2%
*-commutative81.2%
exp-to-pow82.4%
sub-neg82.4%
metadata-eval82.4%
*-commutative82.4%
exp-to-pow82.4%
Simplified82.4%
Taylor expanded in b around 0 82.0%
associate-/l*86.2%
*-commutative86.2%
exp-to-pow86.2%
*-commutative86.2%
exp-sum98.7%
exp-sum86.2%
*-commutative86.2%
exp-to-pow86.2%
exp-to-pow87.5%
sub-neg87.5%
metadata-eval87.5%
Simplified87.5%
Taylor expanded in t around 0 82.8%
associate-/l*85.3%
Simplified85.3%
Final simplification83.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (* y (exp b)))))
(if (<= b -450.0)
t_1
(if (<= b -3.4e-285)
(- (/ x (* y a)) (* (/ x a) (/ b y)))
(if (<= b 3.3e-219)
(/ (- x) (/ y b))
(if (<= b 230000.0) (* (/ x a) (/ 1.0 y)) t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / (y * exp(b));
double tmp;
if (b <= -450.0) {
tmp = t_1;
} else if (b <= -3.4e-285) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= 3.3e-219) {
tmp = -x / (y / b);
} else if (b <= 230000.0) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = x / (y * exp(b))
if (b <= (-450.0d0)) then
tmp = t_1
else if (b <= (-3.4d-285)) then
tmp = (x / (y * a)) - ((x / a) * (b / y))
else if (b <= 3.3d-219) then
tmp = -x / (y / b)
else if (b <= 230000.0d0) then
tmp = (x / a) * (1.0d0 / y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / (y * Math.exp(b));
double tmp;
if (b <= -450.0) {
tmp = t_1;
} else if (b <= -3.4e-285) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= 3.3e-219) {
tmp = -x / (y / b);
} else if (b <= 230000.0) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x / (y * math.exp(b)) tmp = 0 if b <= -450.0: tmp = t_1 elif b <= -3.4e-285: tmp = (x / (y * a)) - ((x / a) * (b / y)) elif b <= 3.3e-219: tmp = -x / (y / b) elif b <= 230000.0: tmp = (x / a) * (1.0 / y) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x / Float64(y * exp(b))) tmp = 0.0 if (b <= -450.0) tmp = t_1; elseif (b <= -3.4e-285) tmp = Float64(Float64(x / Float64(y * a)) - Float64(Float64(x / a) * Float64(b / y))); elseif (b <= 3.3e-219) tmp = Float64(Float64(-x) / Float64(y / b)); elseif (b <= 230000.0) tmp = Float64(Float64(x / a) * Float64(1.0 / y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x / (y * exp(b)); tmp = 0.0; if (b <= -450.0) tmp = t_1; elseif (b <= -3.4e-285) tmp = (x / (y * a)) - ((x / a) * (b / y)); elseif (b <= 3.3e-219) tmp = -x / (y / b); elseif (b <= 230000.0) tmp = (x / a) * (1.0 / y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -450.0], t$95$1, If[LessEqual[b, -3.4e-285], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] - N[(N[(x / a), $MachinePrecision] * N[(b / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.3e-219], N[((-x) / N[(y / b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 230000.0], N[(N[(x / a), $MachinePrecision] * N[(1.0 / y), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{y \cdot e^{b}}\\
\mathbf{if}\;b \leq -450:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -3.4 \cdot 10^{-285}:\\
\;\;\;\;\frac{x}{y \cdot a} - \frac{x}{a} \cdot \frac{b}{y}\\
\mathbf{elif}\;b \leq 3.3 \cdot 10^{-219}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\mathbf{elif}\;b \leq 230000:\\
\;\;\;\;\frac{x}{a} \cdot \frac{1}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -450 or 2.3e5 < b Initial program 100.0%
Taylor expanded in y around 0 91.3%
Taylor expanded in b around inf 81.9%
neg-mul-181.9%
Simplified81.9%
associate-/l*81.9%
div-inv81.9%
div-inv81.9%
add-sqr-sqrt45.1%
sqrt-unprod50.1%
sqr-neg50.1%
sqrt-unprod5.0%
add-sqr-sqrt19.7%
exp-neg19.7%
add-sqr-sqrt14.7%
sqrt-unprod51.5%
sqr-neg51.5%
sqrt-unprod36.8%
add-sqr-sqrt81.9%
Applied egg-rr81.9%
associate-*r/81.9%
*-rgt-identity81.9%
Simplified81.9%
if -450 < b < -3.3999999999999999e-285Initial program 97.0%
associate-*l/94.6%
*-commutative94.6%
exp-diff94.7%
+-commutative94.7%
exp-sum75.5%
associate-/l*75.5%
*-commutative75.5%
exp-to-pow76.5%
sub-neg76.5%
metadata-eval76.5%
*-commutative76.5%
exp-to-pow76.5%
Simplified76.5%
Taylor expanded in t around 0 64.6%
times-frac66.8%
Simplified66.8%
Taylor expanded in y around 0 38.6%
associate-*r*38.6%
Simplified38.6%
Taylor expanded in b around 0 38.6%
+-commutative38.6%
mul-1-neg38.6%
unsub-neg38.6%
*-commutative38.6%
times-frac42.6%
Simplified42.6%
if -3.3999999999999999e-285 < b < 3.3000000000000002e-219Initial program 99.1%
Taylor expanded in y around 0 70.9%
Taylor expanded in b around inf 20.7%
neg-mul-120.7%
Simplified20.7%
Taylor expanded in b around 0 20.7%
Taylor expanded in b around inf 40.5%
mul-1-neg40.5%
associate-*r/40.5%
*-commutative40.5%
associate-/r/54.9%
distribute-neg-frac54.9%
Simplified54.9%
if 3.3000000000000002e-219 < b < 2.3e5Initial program 96.4%
associate-*l/86.5%
*-commutative86.5%
exp-diff84.5%
+-commutative84.5%
exp-sum74.3%
associate-/l*74.3%
*-commutative74.3%
exp-to-pow75.6%
sub-neg75.6%
metadata-eval75.6%
*-commutative75.6%
exp-to-pow75.6%
Simplified75.6%
Taylor expanded in t around 0 64.0%
times-frac68.1%
Simplified68.1%
Taylor expanded in b around 0 63.0%
times-frac67.0%
Simplified67.0%
Taylor expanded in y around 0 36.3%
Final simplification63.4%
(FPCore (x y z t a b) :precision binary64 (if (or (<= b -6600000000.0) (not (<= b 1400000000.0))) (/ x (* y (exp b))) (* (/ (pow z y) y) (/ x a))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((b <= -6600000000.0) || !(b <= 1400000000.0)) {
tmp = x / (y * exp(b));
} else {
tmp = (pow(z, y) / y) * (x / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b <= (-6600000000.0d0)) .or. (.not. (b <= 1400000000.0d0))) then
tmp = x / (y * exp(b))
else
tmp = ((z ** y) / y) * (x / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((b <= -6600000000.0) || !(b <= 1400000000.0)) {
tmp = x / (y * Math.exp(b));
} else {
tmp = (Math.pow(z, y) / y) * (x / a);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (b <= -6600000000.0) or not (b <= 1400000000.0): tmp = x / (y * math.exp(b)) else: tmp = (math.pow(z, y) / y) * (x / a) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((b <= -6600000000.0) || !(b <= 1400000000.0)) tmp = Float64(x / Float64(y * exp(b))); else tmp = Float64(Float64((z ^ y) / y) * Float64(x / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((b <= -6600000000.0) || ~((b <= 1400000000.0))) tmp = x / (y * exp(b)); else tmp = ((z ^ y) / y) * (x / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[b, -6600000000.0], N[Not[LessEqual[b, 1400000000.0]], $MachinePrecision]], N[(x / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision] * N[(x / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -6600000000 \lor \neg \left(b \leq 1400000000\right):\\
\;\;\;\;\frac{x}{y \cdot e^{b}}\\
\mathbf{else}:\\
\;\;\;\;\frac{{z}^{y}}{y} \cdot \frac{x}{a}\\
\end{array}
\end{array}
if b < -6.6e9 or 1.4e9 < b Initial program 100.0%
Taylor expanded in y around 0 91.1%
Taylor expanded in b around inf 82.2%
neg-mul-182.2%
Simplified82.2%
associate-/l*82.2%
div-inv82.2%
div-inv82.2%
add-sqr-sqrt44.6%
sqrt-unprod49.7%
sqr-neg49.7%
sqrt-unprod5.1%
add-sqr-sqrt19.3%
exp-neg19.3%
add-sqr-sqrt14.2%
sqrt-unprod51.9%
sqr-neg51.9%
sqrt-unprod37.7%
add-sqr-sqrt82.2%
Applied egg-rr82.2%
associate-*r/82.2%
*-rgt-identity82.2%
Simplified82.2%
if -6.6e9 < b < 1.4e9Initial program 97.3%
associate-*l/89.9%
*-commutative89.9%
exp-diff88.4%
+-commutative88.4%
exp-sum75.4%
associate-/l*75.4%
*-commutative75.4%
exp-to-pow76.4%
sub-neg76.4%
metadata-eval76.4%
*-commutative76.4%
exp-to-pow76.4%
Simplified76.4%
Taylor expanded in t around 0 66.5%
times-frac69.0%
Simplified69.0%
Taylor expanded in b around 0 65.9%
times-frac68.4%
Simplified68.4%
Final simplification75.6%
(FPCore (x y z t a b) :precision binary64 (if (or (<= b -1550000000.0) (not (<= b 14000000.0))) (/ x (* y (exp b))) (/ x (/ a (/ (pow z y) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((b <= -1550000000.0) || !(b <= 14000000.0)) {
tmp = x / (y * exp(b));
} else {
tmp = x / (a / (pow(z, y) / y));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if ((b <= (-1550000000.0d0)) .or. (.not. (b <= 14000000.0d0))) then
tmp = x / (y * exp(b))
else
tmp = x / (a / ((z ** y) / y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((b <= -1550000000.0) || !(b <= 14000000.0)) {
tmp = x / (y * Math.exp(b));
} else {
tmp = x / (a / (Math.pow(z, y) / y));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (b <= -1550000000.0) or not (b <= 14000000.0): tmp = x / (y * math.exp(b)) else: tmp = x / (a / (math.pow(z, y) / y)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((b <= -1550000000.0) || !(b <= 14000000.0)) tmp = Float64(x / Float64(y * exp(b))); else tmp = Float64(x / Float64(a / Float64((z ^ y) / y))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((b <= -1550000000.0) || ~((b <= 14000000.0))) tmp = x / (y * exp(b)); else tmp = x / (a / ((z ^ y) / y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[b, -1550000000.0], N[Not[LessEqual[b, 14000000.0]], $MachinePrecision]], N[(x / N[(y * N[Exp[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(a / N[(N[Power[z, y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1550000000 \lor \neg \left(b \leq 14000000\right):\\
\;\;\;\;\frac{x}{y \cdot e^{b}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{a}{\frac{{z}^{y}}{y}}}\\
\end{array}
\end{array}
if b < -1.55e9 or 1.4e7 < b Initial program 100.0%
Taylor expanded in y around 0 91.1%
Taylor expanded in b around inf 82.2%
neg-mul-182.2%
Simplified82.2%
associate-/l*82.2%
div-inv82.2%
div-inv82.2%
add-sqr-sqrt44.6%
sqrt-unprod49.7%
sqr-neg49.7%
sqrt-unprod5.1%
add-sqr-sqrt19.3%
exp-neg19.3%
add-sqr-sqrt14.2%
sqrt-unprod51.9%
sqr-neg51.9%
sqrt-unprod37.7%
add-sqr-sqrt82.2%
Applied egg-rr82.2%
associate-*r/82.2%
*-rgt-identity82.2%
Simplified82.2%
if -1.55e9 < b < 1.4e7Initial program 97.3%
associate-*l/89.9%
*-commutative89.9%
exp-diff88.4%
+-commutative88.4%
exp-sum75.4%
associate-/l*75.4%
*-commutative75.4%
exp-to-pow76.4%
sub-neg76.4%
metadata-eval76.4%
*-commutative76.4%
exp-to-pow76.4%
Simplified76.4%
Taylor expanded in b around 0 80.6%
associate-/l*81.3%
*-commutative81.3%
exp-to-pow81.3%
*-commutative81.3%
exp-sum96.7%
exp-sum81.3%
*-commutative81.3%
exp-to-pow81.3%
exp-to-pow82.3%
sub-neg82.3%
metadata-eval82.3%
Simplified82.3%
Taylor expanded in t around 0 65.9%
associate-/l*70.8%
Simplified70.8%
Final simplification76.7%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -1.45e-284)
(- (/ x (* y a)) (* (/ x a) (/ b y)))
(if (<= b 9.5e-164)
(/ (- x) (/ y b))
(if (<= b 0.038)
(* (/ x a) (/ 1.0 y))
(/ 1.0 (* a (+ (/ y x) (/ b (/ x y)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.45e-284) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= 9.5e-164) {
tmp = -x / (y / b);
} else if (b <= 0.038) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = 1.0 / (a * ((y / x) + (b / (x / 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.45d-284)) then
tmp = (x / (y * a)) - ((x / a) * (b / y))
else if (b <= 9.5d-164) then
tmp = -x / (y / b)
else if (b <= 0.038d0) then
tmp = (x / a) * (1.0d0 / y)
else
tmp = 1.0d0 / (a * ((y / x) + (b / (x / 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.45e-284) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= 9.5e-164) {
tmp = -x / (y / b);
} else if (b <= 0.038) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = 1.0 / (a * ((y / x) + (b / (x / y))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1.45e-284: tmp = (x / (y * a)) - ((x / a) * (b / y)) elif b <= 9.5e-164: tmp = -x / (y / b) elif b <= 0.038: tmp = (x / a) * (1.0 / y) else: tmp = 1.0 / (a * ((y / x) + (b / (x / y)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.45e-284) tmp = Float64(Float64(x / Float64(y * a)) - Float64(Float64(x / a) * Float64(b / y))); elseif (b <= 9.5e-164) tmp = Float64(Float64(-x) / Float64(y / b)); elseif (b <= 0.038) tmp = Float64(Float64(x / a) * Float64(1.0 / y)); else tmp = Float64(1.0 / Float64(a * Float64(Float64(y / x) + Float64(b / Float64(x / y))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -1.45e-284) tmp = (x / (y * a)) - ((x / a) * (b / y)); elseif (b <= 9.5e-164) tmp = -x / (y / b); elseif (b <= 0.038) tmp = (x / a) * (1.0 / y); else tmp = 1.0 / (a * ((y / x) + (b / (x / y)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.45e-284], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] - N[(N[(x / a), $MachinePrecision] * N[(b / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 9.5e-164], N[((-x) / N[(y / b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 0.038], N[(N[(x / a), $MachinePrecision] * N[(1.0 / y), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(a * N[(N[(y / x), $MachinePrecision] + N[(b / N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.45 \cdot 10^{-284}:\\
\;\;\;\;\frac{x}{y \cdot a} - \frac{x}{a} \cdot \frac{b}{y}\\
\mathbf{elif}\;b \leq 9.5 \cdot 10^{-164}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\mathbf{elif}\;b \leq 0.038:\\
\;\;\;\;\frac{x}{a} \cdot \frac{1}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{a \cdot \left(\frac{y}{x} + \frac{b}{\frac{x}{y}}\right)}\\
\end{array}
\end{array}
if b < -1.4500000000000001e-284Initial program 98.9%
associate-*l/94.1%
*-commutative94.1%
exp-diff66.5%
+-commutative66.5%
exp-sum58.7%
associate-/l*58.7%
*-commutative58.7%
exp-to-pow59.0%
sub-neg59.0%
metadata-eval59.0%
*-commutative59.0%
exp-to-pow59.0%
Simplified59.0%
Taylor expanded in t around 0 61.0%
times-frac58.7%
Simplified58.7%
Taylor expanded in y around 0 62.5%
associate-*r*55.4%
Simplified55.4%
Taylor expanded in b around 0 39.1%
+-commutative39.1%
mul-1-neg39.1%
unsub-neg39.1%
*-commutative39.1%
times-frac39.7%
Simplified39.7%
if -1.4500000000000001e-284 < b < 9.5000000000000001e-164Initial program 96.6%
Taylor expanded in y around 0 67.2%
Taylor expanded in b around inf 22.0%
neg-mul-122.0%
Simplified22.0%
Taylor expanded in b around 0 22.0%
Taylor expanded in b around inf 39.1%
mul-1-neg39.1%
associate-*r/36.2%
*-commutative36.2%
associate-/r/49.6%
distribute-neg-frac49.6%
Simplified49.6%
if 9.5000000000000001e-164 < b < 0.0379999999999999991Initial program 98.2%
associate-*l/84.9%
*-commutative84.9%
exp-diff84.9%
+-commutative84.9%
exp-sum73.4%
associate-/l*73.4%
*-commutative73.4%
exp-to-pow74.5%
sub-neg74.5%
metadata-eval74.5%
*-commutative74.5%
exp-to-pow74.5%
Simplified74.5%
Taylor expanded in t around 0 58.1%
times-frac66.5%
Simplified66.5%
Taylor expanded in b around 0 58.1%
times-frac66.5%
Simplified66.5%
Taylor expanded in y around 0 42.6%
if 0.0379999999999999991 < b Initial program 99.6%
associate-*l/86.5%
*-commutative86.5%
exp-diff58.7%
+-commutative58.7%
exp-sum57.1%
associate-/l*57.1%
*-commutative57.1%
exp-to-pow57.3%
sub-neg57.3%
metadata-eval57.3%
*-commutative57.3%
exp-to-pow57.3%
Simplified57.3%
Taylor expanded in t around 0 80.4%
times-frac73.8%
Simplified73.8%
Taylor expanded in y around 0 85.5%
associate-*r*74.0%
Simplified74.0%
clear-num74.0%
inv-pow74.0%
Applied egg-rr74.0%
unpow-174.0%
associate-*r*85.5%
*-lft-identity85.5%
times-frac85.5%
/-rgt-identity85.5%
*-commutative85.5%
associate-*r/75.6%
Simplified75.6%
Taylor expanded in b around 0 34.6%
associate-/l*34.6%
Simplified34.6%
Final simplification40.2%
(FPCore (x y z t a b)
:precision binary64
(if (<= b 1.08e-296)
(- (/ x (* y a)) (/ (* x b) (* y a)))
(if (<= b 4.2e-165)
(/ (- x) (/ y b))
(if (<= b 0.038)
(* (/ x a) (/ 1.0 y))
(/ 1.0 (* a (+ (/ y x) (/ b (/ x y)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 1.08e-296) {
tmp = (x / (y * a)) - ((x * b) / (y * a));
} else if (b <= 4.2e-165) {
tmp = -x / (y / b);
} else if (b <= 0.038) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = 1.0 / (a * ((y / x) + (b / (x / 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.08d-296) then
tmp = (x / (y * a)) - ((x * b) / (y * a))
else if (b <= 4.2d-165) then
tmp = -x / (y / b)
else if (b <= 0.038d0) then
tmp = (x / a) * (1.0d0 / y)
else
tmp = 1.0d0 / (a * ((y / x) + (b / (x / 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.08e-296) {
tmp = (x / (y * a)) - ((x * b) / (y * a));
} else if (b <= 4.2e-165) {
tmp = -x / (y / b);
} else if (b <= 0.038) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = 1.0 / (a * ((y / x) + (b / (x / y))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 1.08e-296: tmp = (x / (y * a)) - ((x * b) / (y * a)) elif b <= 4.2e-165: tmp = -x / (y / b) elif b <= 0.038: tmp = (x / a) * (1.0 / y) else: tmp = 1.0 / (a * ((y / x) + (b / (x / y)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 1.08e-296) tmp = Float64(Float64(x / Float64(y * a)) - Float64(Float64(x * b) / Float64(y * a))); elseif (b <= 4.2e-165) tmp = Float64(Float64(-x) / Float64(y / b)); elseif (b <= 0.038) tmp = Float64(Float64(x / a) * Float64(1.0 / y)); else tmp = Float64(1.0 / Float64(a * Float64(Float64(y / x) + Float64(b / Float64(x / y))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 1.08e-296) tmp = (x / (y * a)) - ((x * b) / (y * a)); elseif (b <= 4.2e-165) tmp = -x / (y / b); elseif (b <= 0.038) tmp = (x / a) * (1.0 / y); else tmp = 1.0 / (a * ((y / x) + (b / (x / y)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 1.08e-296], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] - N[(N[(x * b), $MachinePrecision] / N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 4.2e-165], N[((-x) / N[(y / b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 0.038], N[(N[(x / a), $MachinePrecision] * N[(1.0 / y), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(a * N[(N[(y / x), $MachinePrecision] + N[(b / N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 1.08 \cdot 10^{-296}:\\
\;\;\;\;\frac{x}{y \cdot a} - \frac{x \cdot b}{y \cdot a}\\
\mathbf{elif}\;b \leq 4.2 \cdot 10^{-165}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\mathbf{elif}\;b \leq 0.038:\\
\;\;\;\;\frac{x}{a} \cdot \frac{1}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{a \cdot \left(\frac{y}{x} + \frac{b}{\frac{x}{y}}\right)}\\
\end{array}
\end{array}
if b < 1.08e-296Initial program 98.9%
associate-*l/93.0%
*-commutative93.0%
exp-diff67.5%
+-commutative67.5%
exp-sum60.2%
associate-/l*60.2%
*-commutative60.2%
exp-to-pow60.5%
sub-neg60.5%
metadata-eval60.5%
*-commutative60.5%
exp-to-pow60.5%
Simplified60.5%
Taylor expanded in t around 0 60.2%
times-frac58.1%
Simplified58.1%
Taylor expanded in y around 0 61.0%
associate-*r*54.3%
Simplified54.3%
Taylor expanded in b around 0 39.2%
if 1.08e-296 < b < 4.1999999999999999e-165Initial program 95.3%
Taylor expanded in y around 0 61.7%
Taylor expanded in b around inf 25.6%
neg-mul-125.6%
Simplified25.6%
Taylor expanded in b around 0 25.6%
Taylor expanded in b around inf 49.9%
mul-1-neg49.9%
associate-*r/45.8%
*-commutative45.8%
associate-/r/57.5%
distribute-neg-frac57.5%
Simplified57.5%
if 4.1999999999999999e-165 < b < 0.0379999999999999991Initial program 98.2%
associate-*l/84.9%
*-commutative84.9%
exp-diff84.9%
+-commutative84.9%
exp-sum73.4%
associate-/l*73.4%
*-commutative73.4%
exp-to-pow74.5%
sub-neg74.5%
metadata-eval74.5%
*-commutative74.5%
exp-to-pow74.5%
Simplified74.5%
Taylor expanded in t around 0 58.1%
times-frac66.5%
Simplified66.5%
Taylor expanded in b around 0 58.1%
times-frac66.5%
Simplified66.5%
Taylor expanded in y around 0 42.6%
if 0.0379999999999999991 < b Initial program 99.6%
associate-*l/86.5%
*-commutative86.5%
exp-diff58.7%
+-commutative58.7%
exp-sum57.1%
associate-/l*57.1%
*-commutative57.1%
exp-to-pow57.3%
sub-neg57.3%
metadata-eval57.3%
*-commutative57.3%
exp-to-pow57.3%
Simplified57.3%
Taylor expanded in t around 0 80.4%
times-frac73.8%
Simplified73.8%
Taylor expanded in y around 0 85.5%
associate-*r*74.0%
Simplified74.0%
clear-num74.0%
inv-pow74.0%
Applied egg-rr74.0%
unpow-174.0%
associate-*r*85.5%
*-lft-identity85.5%
times-frac85.5%
/-rgt-identity85.5%
*-commutative85.5%
associate-*r/75.6%
Simplified75.6%
Taylor expanded in b around 0 34.6%
associate-/l*34.6%
Simplified34.6%
Final simplification40.2%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -5000.0)
(/ (* x (- b)) y)
(if (<= b -5e-277)
(/ 1.0 (* a (/ y x)))
(if (<= b 1.7e-222)
(/ (- x) (/ y b))
(if (<= b 7.5e+91) (* (/ x a) (/ 1.0 y)) (/ x (+ y (* y b))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -5000.0) {
tmp = (x * -b) / y;
} else if (b <= -5e-277) {
tmp = 1.0 / (a * (y / x));
} else if (b <= 1.7e-222) {
tmp = -x / (y / b);
} else if (b <= 7.5e+91) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = x / (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 <= (-5000.0d0)) then
tmp = (x * -b) / y
else if (b <= (-5d-277)) then
tmp = 1.0d0 / (a * (y / x))
else if (b <= 1.7d-222) then
tmp = -x / (y / b)
else if (b <= 7.5d+91) then
tmp = (x / a) * (1.0d0 / y)
else
tmp = x / (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 <= -5000.0) {
tmp = (x * -b) / y;
} else if (b <= -5e-277) {
tmp = 1.0 / (a * (y / x));
} else if (b <= 1.7e-222) {
tmp = -x / (y / b);
} else if (b <= 7.5e+91) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = x / (y + (y * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -5000.0: tmp = (x * -b) / y elif b <= -5e-277: tmp = 1.0 / (a * (y / x)) elif b <= 1.7e-222: tmp = -x / (y / b) elif b <= 7.5e+91: tmp = (x / a) * (1.0 / y) else: tmp = x / (y + (y * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -5000.0) tmp = Float64(Float64(x * Float64(-b)) / y); elseif (b <= -5e-277) tmp = Float64(1.0 / Float64(a * Float64(y / x))); elseif (b <= 1.7e-222) tmp = Float64(Float64(-x) / Float64(y / b)); elseif (b <= 7.5e+91) tmp = Float64(Float64(x / a) * Float64(1.0 / y)); else tmp = Float64(x / Float64(y + Float64(y * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -5000.0) tmp = (x * -b) / y; elseif (b <= -5e-277) tmp = 1.0 / (a * (y / x)); elseif (b <= 1.7e-222) tmp = -x / (y / b); elseif (b <= 7.5e+91) tmp = (x / a) * (1.0 / y); else tmp = x / (y + (y * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -5000.0], N[(N[(x * (-b)), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -5e-277], N[(1.0 / N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.7e-222], N[((-x) / N[(y / b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 7.5e+91], N[(N[(x / a), $MachinePrecision] * N[(1.0 / y), $MachinePrecision]), $MachinePrecision], N[(x / N[(y + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -5000:\\
\;\;\;\;\frac{x \cdot \left(-b\right)}{y}\\
\mathbf{elif}\;b \leq -5 \cdot 10^{-277}:\\
\;\;\;\;\frac{1}{a \cdot \frac{y}{x}}\\
\mathbf{elif}\;b \leq 1.7 \cdot 10^{-222}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\mathbf{elif}\;b \leq 7.5 \cdot 10^{+91}:\\
\;\;\;\;\frac{x}{a} \cdot \frac{1}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y + y \cdot b}\\
\end{array}
\end{array}
if b < -5e3Initial program 100.0%
Taylor expanded in y around 0 88.9%
Taylor expanded in b around inf 76.6%
neg-mul-176.6%
Simplified76.6%
Taylor expanded in b around 0 37.4%
Taylor expanded in b around inf 37.4%
if -5e3 < b < -5e-277Initial program 99.0%
associate-*l/94.6%
*-commutative94.6%
exp-diff94.7%
+-commutative94.7%
exp-sum74.7%
associate-/l*74.7%
*-commutative74.7%
exp-to-pow75.5%
sub-neg75.5%
metadata-eval75.5%
*-commutative75.5%
exp-to-pow75.5%
Simplified75.5%
Taylor expanded in t around 0 63.0%
times-frac67.4%
Simplified67.4%
Taylor expanded in y around 0 35.9%
associate-*r*35.9%
Simplified35.9%
clear-num35.8%
inv-pow35.8%
Applied egg-rr35.8%
unpow-135.8%
associate-*r*35.8%
*-lft-identity35.8%
times-frac42.3%
/-rgt-identity42.3%
*-commutative42.3%
associate-*r/42.3%
Simplified42.3%
Taylor expanded in b around 0 35.4%
associate-*r/41.7%
Simplified41.7%
if -5e-277 < b < 1.7000000000000001e-222Initial program 95.6%
Taylor expanded in y around 0 65.8%
Taylor expanded in b around inf 19.4%
neg-mul-119.4%
Simplified19.4%
Taylor expanded in b around 0 19.4%
Taylor expanded in b around inf 38.1%
mul-1-neg38.1%
associate-*r/38.1%
*-commutative38.1%
associate-/r/54.9%
distribute-neg-frac54.9%
Simplified54.9%
if 1.7000000000000001e-222 < b < 7.50000000000000033e91Initial program 97.2%
associate-*l/89.6%
*-commutative89.6%
exp-diff80.3%
+-commutative80.3%
exp-sum72.5%
associate-/l*72.5%
*-commutative72.5%
exp-to-pow73.5%
sub-neg73.5%
metadata-eval73.5%
*-commutative73.5%
exp-to-pow73.5%
Simplified73.5%
Taylor expanded in t around 0 66.2%
times-frac69.3%
Simplified69.3%
Taylor expanded in b around 0 58.0%
times-frac62.6%
Simplified62.6%
Taylor expanded in y around 0 31.6%
if 7.50000000000000033e91 < b Initial program 100.0%
Taylor expanded in y around 0 97.6%
Taylor expanded in b around inf 92.8%
neg-mul-192.8%
Simplified92.8%
associate-/l*92.8%
div-inv92.8%
div-inv92.8%
add-sqr-sqrt0.0%
sqrt-unprod8.8%
sqr-neg8.8%
sqrt-unprod8.8%
add-sqr-sqrt8.8%
exp-neg8.8%
add-sqr-sqrt0.0%
sqrt-unprod92.8%
sqr-neg92.8%
sqrt-unprod92.8%
add-sqr-sqrt92.8%
Applied egg-rr92.8%
associate-*r/92.8%
*-rgt-identity92.8%
Simplified92.8%
Taylor expanded in b around 0 36.1%
Final simplification38.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -1300.0)
(/ (- x (* x b)) y)
(if (<= b -1.95e-275)
(/ 1.0 (* a (/ y x)))
(if (<= b 2.65e-222)
(/ (- x) (/ y b))
(if (<= b 3.2e+91) (* (/ x a) (/ 1.0 y)) (/ x (+ y (* y b))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1300.0) {
tmp = (x - (x * b)) / y;
} else if (b <= -1.95e-275) {
tmp = 1.0 / (a * (y / x));
} else if (b <= 2.65e-222) {
tmp = -x / (y / b);
} else if (b <= 3.2e+91) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = x / (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 <= (-1300.0d0)) then
tmp = (x - (x * b)) / y
else if (b <= (-1.95d-275)) then
tmp = 1.0d0 / (a * (y / x))
else if (b <= 2.65d-222) then
tmp = -x / (y / b)
else if (b <= 3.2d+91) then
tmp = (x / a) * (1.0d0 / y)
else
tmp = x / (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 <= -1300.0) {
tmp = (x - (x * b)) / y;
} else if (b <= -1.95e-275) {
tmp = 1.0 / (a * (y / x));
} else if (b <= 2.65e-222) {
tmp = -x / (y / b);
} else if (b <= 3.2e+91) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = x / (y + (y * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1300.0: tmp = (x - (x * b)) / y elif b <= -1.95e-275: tmp = 1.0 / (a * (y / x)) elif b <= 2.65e-222: tmp = -x / (y / b) elif b <= 3.2e+91: tmp = (x / a) * (1.0 / y) else: tmp = x / (y + (y * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1300.0) tmp = Float64(Float64(x - Float64(x * b)) / y); elseif (b <= -1.95e-275) tmp = Float64(1.0 / Float64(a * Float64(y / x))); elseif (b <= 2.65e-222) tmp = Float64(Float64(-x) / Float64(y / b)); elseif (b <= 3.2e+91) tmp = Float64(Float64(x / a) * Float64(1.0 / y)); else tmp = Float64(x / Float64(y + Float64(y * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -1300.0) tmp = (x - (x * b)) / y; elseif (b <= -1.95e-275) tmp = 1.0 / (a * (y / x)); elseif (b <= 2.65e-222) tmp = -x / (y / b); elseif (b <= 3.2e+91) tmp = (x / a) * (1.0 / y); else tmp = x / (y + (y * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1300.0], N[(N[(x - N[(x * b), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -1.95e-275], N[(1.0 / N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.65e-222], N[((-x) / N[(y / b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.2e+91], N[(N[(x / a), $MachinePrecision] * N[(1.0 / y), $MachinePrecision]), $MachinePrecision], N[(x / N[(y + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1300:\\
\;\;\;\;\frac{x - x \cdot b}{y}\\
\mathbf{elif}\;b \leq -1.95 \cdot 10^{-275}:\\
\;\;\;\;\frac{1}{a \cdot \frac{y}{x}}\\
\mathbf{elif}\;b \leq 2.65 \cdot 10^{-222}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\mathbf{elif}\;b \leq 3.2 \cdot 10^{+91}:\\
\;\;\;\;\frac{x}{a} \cdot \frac{1}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y + y \cdot b}\\
\end{array}
\end{array}
if b < -1300Initial program 100.0%
Taylor expanded in y around 0 88.9%
Taylor expanded in b around inf 76.6%
neg-mul-176.6%
Simplified76.6%
Taylor expanded in b around 0 37.4%
Taylor expanded in y around 0 37.4%
mul-1-neg37.4%
distribute-rgt-neg-out37.4%
Simplified37.4%
if -1300 < b < -1.94999999999999986e-275Initial program 99.0%
associate-*l/94.6%
*-commutative94.6%
exp-diff94.7%
+-commutative94.7%
exp-sum74.7%
associate-/l*74.7%
*-commutative74.7%
exp-to-pow75.5%
sub-neg75.5%
metadata-eval75.5%
*-commutative75.5%
exp-to-pow75.5%
Simplified75.5%
Taylor expanded in t around 0 63.0%
times-frac67.4%
Simplified67.4%
Taylor expanded in y around 0 35.9%
associate-*r*35.9%
Simplified35.9%
clear-num35.8%
inv-pow35.8%
Applied egg-rr35.8%
unpow-135.8%
associate-*r*35.8%
*-lft-identity35.8%
times-frac42.3%
/-rgt-identity42.3%
*-commutative42.3%
associate-*r/42.3%
Simplified42.3%
Taylor expanded in b around 0 35.4%
associate-*r/41.7%
Simplified41.7%
if -1.94999999999999986e-275 < b < 2.6499999999999999e-222Initial program 95.6%
Taylor expanded in y around 0 65.8%
Taylor expanded in b around inf 19.4%
neg-mul-119.4%
Simplified19.4%
Taylor expanded in b around 0 19.4%
Taylor expanded in b around inf 38.1%
mul-1-neg38.1%
associate-*r/38.1%
*-commutative38.1%
associate-/r/54.9%
distribute-neg-frac54.9%
Simplified54.9%
if 2.6499999999999999e-222 < b < 3.19999999999999989e91Initial program 97.2%
associate-*l/89.6%
*-commutative89.6%
exp-diff80.3%
+-commutative80.3%
exp-sum72.5%
associate-/l*72.5%
*-commutative72.5%
exp-to-pow73.5%
sub-neg73.5%
metadata-eval73.5%
*-commutative73.5%
exp-to-pow73.5%
Simplified73.5%
Taylor expanded in t around 0 66.2%
times-frac69.3%
Simplified69.3%
Taylor expanded in b around 0 58.0%
times-frac62.6%
Simplified62.6%
Taylor expanded in y around 0 31.6%
if 3.19999999999999989e91 < b Initial program 100.0%
Taylor expanded in y around 0 97.6%
Taylor expanded in b around inf 92.8%
neg-mul-192.8%
Simplified92.8%
associate-/l*92.8%
div-inv92.8%
div-inv92.8%
add-sqr-sqrt0.0%
sqrt-unprod8.8%
sqr-neg8.8%
sqrt-unprod8.8%
add-sqr-sqrt8.8%
exp-neg8.8%
add-sqr-sqrt0.0%
sqrt-unprod92.8%
sqr-neg92.8%
sqrt-unprod92.8%
add-sqr-sqrt92.8%
Applied egg-rr92.8%
associate-*r/92.8%
*-rgt-identity92.8%
Simplified92.8%
Taylor expanded in b around 0 36.1%
Final simplification38.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -3100.0)
(/ (- x (* x b)) y)
(if (<= b -8.8e-277)
(/ 1.0 (* a (/ y x)))
(if (<= b 5.8e-223) (/ (- x) (/ y b)) (/ x (* a (+ y (* y b))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -3100.0) {
tmp = (x - (x * b)) / y;
} else if (b <= -8.8e-277) {
tmp = 1.0 / (a * (y / x));
} else if (b <= 5.8e-223) {
tmp = -x / (y / b);
} 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 <= (-3100.0d0)) then
tmp = (x - (x * b)) / y
else if (b <= (-8.8d-277)) then
tmp = 1.0d0 / (a * (y / x))
else if (b <= 5.8d-223) then
tmp = -x / (y / b)
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 <= -3100.0) {
tmp = (x - (x * b)) / y;
} else if (b <= -8.8e-277) {
tmp = 1.0 / (a * (y / x));
} else if (b <= 5.8e-223) {
tmp = -x / (y / b);
} else {
tmp = x / (a * (y + (y * b)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -3100.0: tmp = (x - (x * b)) / y elif b <= -8.8e-277: tmp = 1.0 / (a * (y / x)) elif b <= 5.8e-223: tmp = -x / (y / b) else: tmp = x / (a * (y + (y * b))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -3100.0) tmp = Float64(Float64(x - Float64(x * b)) / y); elseif (b <= -8.8e-277) tmp = Float64(1.0 / Float64(a * Float64(y / x))); elseif (b <= 5.8e-223) tmp = Float64(Float64(-x) / Float64(y / b)); 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 <= -3100.0) tmp = (x - (x * b)) / y; elseif (b <= -8.8e-277) tmp = 1.0 / (a * (y / x)); elseif (b <= 5.8e-223) tmp = -x / (y / b); else tmp = x / (a * (y + (y * b))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -3100.0], N[(N[(x - N[(x * b), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[b, -8.8e-277], N[(1.0 / N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 5.8e-223], N[((-x) / N[(y / b), $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 -3100:\\
\;\;\;\;\frac{x - x \cdot b}{y}\\
\mathbf{elif}\;b \leq -8.8 \cdot 10^{-277}:\\
\;\;\;\;\frac{1}{a \cdot \frac{y}{x}}\\
\mathbf{elif}\;b \leq 5.8 \cdot 10^{-223}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{a \cdot \left(y + y \cdot b\right)}\\
\end{array}
\end{array}
if b < -3100Initial program 100.0%
Taylor expanded in y around 0 88.9%
Taylor expanded in b around inf 76.6%
neg-mul-176.6%
Simplified76.6%
Taylor expanded in b around 0 37.4%
Taylor expanded in y around 0 37.4%
mul-1-neg37.4%
distribute-rgt-neg-out37.4%
Simplified37.4%
if -3100 < b < -8.79999999999999983e-277Initial program 99.0%
associate-*l/94.6%
*-commutative94.6%
exp-diff94.7%
+-commutative94.7%
exp-sum74.7%
associate-/l*74.7%
*-commutative74.7%
exp-to-pow75.5%
sub-neg75.5%
metadata-eval75.5%
*-commutative75.5%
exp-to-pow75.5%
Simplified75.5%
Taylor expanded in t around 0 63.0%
times-frac67.4%
Simplified67.4%
Taylor expanded in y around 0 35.9%
associate-*r*35.9%
Simplified35.9%
clear-num35.8%
inv-pow35.8%
Applied egg-rr35.8%
unpow-135.8%
associate-*r*35.8%
*-lft-identity35.8%
times-frac42.3%
/-rgt-identity42.3%
*-commutative42.3%
associate-*r/42.3%
Simplified42.3%
Taylor expanded in b around 0 35.4%
associate-*r/41.7%
Simplified41.7%
if -8.79999999999999983e-277 < b < 5.8000000000000001e-223Initial program 95.6%
Taylor expanded in y around 0 65.8%
Taylor expanded in b around inf 19.4%
neg-mul-119.4%
Simplified19.4%
Taylor expanded in b around 0 19.4%
Taylor expanded in b around inf 38.1%
mul-1-neg38.1%
associate-*r/38.1%
*-commutative38.1%
associate-/r/54.9%
distribute-neg-frac54.9%
Simplified54.9%
if 5.8000000000000001e-223 < b Initial program 98.3%
associate-*l/86.1%
*-commutative86.1%
exp-diff69.9%
+-commutative69.9%
exp-sum64.2%
associate-/l*64.2%
*-commutative64.2%
exp-to-pow64.8%
sub-neg64.8%
metadata-eval64.8%
*-commutative64.8%
exp-to-pow64.8%
Simplified64.8%
Taylor expanded in t around 0 73.7%
times-frac70.9%
Simplified70.9%
Taylor expanded in y around 0 64.2%
associate-*r*57.6%
Simplified57.6%
Taylor expanded in b around 0 31.0%
distribute-lft-out31.0%
*-commutative31.0%
Simplified31.0%
Final simplification37.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -1.05e-287)
(- (/ x (* y a)) (* (/ x a) (/ b y)))
(if (<= b 5e-164)
(/ (- x) (/ y b))
(if (<= b 2.9e+91) (* (/ x a) (/ 1.0 y)) (/ x (+ y (* y b)))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.05e-287) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= 5e-164) {
tmp = -x / (y / b);
} else if (b <= 2.9e+91) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = x / (y + (y * b));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= (-1.05d-287)) then
tmp = (x / (y * a)) - ((x / a) * (b / y))
else if (b <= 5d-164) then
tmp = -x / (y / b)
else if (b <= 2.9d+91) then
tmp = (x / a) * (1.0d0 / y)
else
tmp = x / (y + (y * b))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -1.05e-287) {
tmp = (x / (y * a)) - ((x / a) * (b / y));
} else if (b <= 5e-164) {
tmp = -x / (y / b);
} else if (b <= 2.9e+91) {
tmp = (x / a) * (1.0 / y);
} else {
tmp = x / (y + (y * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -1.05e-287: tmp = (x / (y * a)) - ((x / a) * (b / y)) elif b <= 5e-164: tmp = -x / (y / b) elif b <= 2.9e+91: tmp = (x / a) * (1.0 / y) else: tmp = x / (y + (y * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -1.05e-287) tmp = Float64(Float64(x / Float64(y * a)) - Float64(Float64(x / a) * Float64(b / y))); elseif (b <= 5e-164) tmp = Float64(Float64(-x) / Float64(y / b)); elseif (b <= 2.9e+91) tmp = Float64(Float64(x / a) * Float64(1.0 / y)); else tmp = Float64(x / Float64(y + Float64(y * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -1.05e-287) tmp = (x / (y * a)) - ((x / a) * (b / y)); elseif (b <= 5e-164) tmp = -x / (y / b); elseif (b <= 2.9e+91) tmp = (x / a) * (1.0 / y); else tmp = x / (y + (y * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -1.05e-287], N[(N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision] - N[(N[(x / a), $MachinePrecision] * N[(b / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 5e-164], N[((-x) / N[(y / b), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.9e+91], N[(N[(x / a), $MachinePrecision] * N[(1.0 / y), $MachinePrecision]), $MachinePrecision], N[(x / N[(y + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.05 \cdot 10^{-287}:\\
\;\;\;\;\frac{x}{y \cdot a} - \frac{x}{a} \cdot \frac{b}{y}\\
\mathbf{elif}\;b \leq 5 \cdot 10^{-164}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\mathbf{elif}\;b \leq 2.9 \cdot 10^{+91}:\\
\;\;\;\;\frac{x}{a} \cdot \frac{1}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y + y \cdot b}\\
\end{array}
\end{array}
if b < -1.0499999999999999e-287Initial program 98.9%
associate-*l/94.1%
*-commutative94.1%
exp-diff66.5%
+-commutative66.5%
exp-sum58.7%
associate-/l*58.7%
*-commutative58.7%
exp-to-pow59.0%
sub-neg59.0%
metadata-eval59.0%
*-commutative59.0%
exp-to-pow59.0%
Simplified59.0%
Taylor expanded in t around 0 61.0%
times-frac58.7%
Simplified58.7%
Taylor expanded in y around 0 62.5%
associate-*r*55.4%
Simplified55.4%
Taylor expanded in b around 0 39.1%
+-commutative39.1%
mul-1-neg39.1%
unsub-neg39.1%
*-commutative39.1%
times-frac39.7%
Simplified39.7%
if -1.0499999999999999e-287 < b < 4.99999999999999962e-164Initial program 96.6%
Taylor expanded in y around 0 67.2%
Taylor expanded in b around inf 22.0%
neg-mul-122.0%
Simplified22.0%
Taylor expanded in b around 0 22.0%
Taylor expanded in b around inf 39.1%
mul-1-neg39.1%
associate-*r/36.2%
*-commutative36.2%
associate-/r/49.6%
distribute-neg-frac49.6%
Simplified49.6%
if 4.99999999999999962e-164 < b < 2.90000000000000014e91Initial program 98.5%
associate-*l/90.0%
*-commutative90.0%
exp-diff79.1%
+-commutative79.1%
exp-sum71.8%
associate-/l*71.8%
*-commutative71.8%
exp-to-pow72.8%
sub-neg72.8%
metadata-eval72.8%
*-commutative72.8%
exp-to-pow72.8%
Simplified72.8%
Taylor expanded in t around 0 62.4%
times-frac69.6%
Simplified69.6%
Taylor expanded in b around 0 52.9%
times-frac61.8%
Simplified61.8%
Taylor expanded in y around 0 32.6%
if 2.90000000000000014e91 < b Initial program 100.0%
Taylor expanded in y around 0 97.6%
Taylor expanded in b around inf 92.8%
neg-mul-192.8%
Simplified92.8%
associate-/l*92.8%
div-inv92.8%
div-inv92.8%
add-sqr-sqrt0.0%
sqrt-unprod8.8%
sqr-neg8.8%
sqrt-unprod8.8%
add-sqr-sqrt8.8%
exp-neg8.8%
add-sqr-sqrt0.0%
sqrt-unprod92.8%
sqr-neg92.8%
sqrt-unprod92.8%
add-sqr-sqrt92.8%
Applied egg-rr92.8%
associate-*r/92.8%
*-rgt-identity92.8%
Simplified92.8%
Taylor expanded in b around 0 36.1%
Final simplification38.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -5.5e-242)
(* (/ x a) (/ 1.0 y))
(if (<= t 2.4e-163)
(/ 1.0 (* a (/ y x)))
(if (<= t 9.0) (/ x (* y a)) (/ (* x (- b)) y)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -5.5e-242) {
tmp = (x / a) * (1.0 / y);
} else if (t <= 2.4e-163) {
tmp = 1.0 / (a * (y / x));
} else if (t <= 9.0) {
tmp = x / (y * a);
} else {
tmp = (x * -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 (t <= (-5.5d-242)) then
tmp = (x / a) * (1.0d0 / y)
else if (t <= 2.4d-163) then
tmp = 1.0d0 / (a * (y / x))
else if (t <= 9.0d0) then
tmp = x / (y * a)
else
tmp = (x * -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 (t <= -5.5e-242) {
tmp = (x / a) * (1.0 / y);
} else if (t <= 2.4e-163) {
tmp = 1.0 / (a * (y / x));
} else if (t <= 9.0) {
tmp = x / (y * a);
} else {
tmp = (x * -b) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -5.5e-242: tmp = (x / a) * (1.0 / y) elif t <= 2.4e-163: tmp = 1.0 / (a * (y / x)) elif t <= 9.0: tmp = x / (y * a) else: tmp = (x * -b) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -5.5e-242) tmp = Float64(Float64(x / a) * Float64(1.0 / y)); elseif (t <= 2.4e-163) tmp = Float64(1.0 / Float64(a * Float64(y / x))); elseif (t <= 9.0) tmp = Float64(x / Float64(y * a)); else tmp = Float64(Float64(x * Float64(-b)) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -5.5e-242) tmp = (x / a) * (1.0 / y); elseif (t <= 2.4e-163) tmp = 1.0 / (a * (y / x)); elseif (t <= 9.0) tmp = x / (y * a); else tmp = (x * -b) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -5.5e-242], N[(N[(x / a), $MachinePrecision] * N[(1.0 / y), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.4e-163], N[(1.0 / N[(a * N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 9.0], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(N[(x * (-b)), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5.5 \cdot 10^{-242}:\\
\;\;\;\;\frac{x}{a} \cdot \frac{1}{y}\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-163}:\\
\;\;\;\;\frac{1}{a \cdot \frac{y}{x}}\\
\mathbf{elif}\;t \leq 9:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \left(-b\right)}{y}\\
\end{array}
\end{array}
if t < -5.4999999999999998e-242Initial program 98.7%
associate-*l/88.3%
*-commutative88.3%
exp-diff69.9%
+-commutative69.9%
exp-sum58.7%
associate-/l*58.7%
*-commutative58.7%
exp-to-pow59.2%
sub-neg59.2%
metadata-eval59.2%
*-commutative59.2%
exp-to-pow59.2%
Simplified59.2%
Taylor expanded in t around 0 70.0%
times-frac69.2%
Simplified69.2%
Taylor expanded in b around 0 58.1%
times-frac59.7%
Simplified59.7%
Taylor expanded in y around 0 36.2%
if -5.4999999999999998e-242 < t < 2.4000000000000001e-163Initial program 97.1%
associate-*l/91.4%
*-commutative91.4%
exp-diff70.3%
+-commutative70.3%
exp-sum70.3%
associate-/l*70.3%
*-commutative70.3%
exp-to-pow71.0%
sub-neg71.0%
metadata-eval71.0%
*-commutative71.0%
exp-to-pow71.0%
Simplified71.0%
Taylor expanded in t around 0 73.6%
times-frac73.8%
Simplified73.8%
Taylor expanded in y around 0 67.1%
associate-*r*61.7%
Simplified61.7%
clear-num61.7%
inv-pow61.7%
Applied egg-rr61.7%
unpow-161.7%
associate-*r*67.1%
*-lft-identity67.1%
times-frac69.7%
/-rgt-identity69.7%
*-commutative69.7%
associate-*r/61.6%
Simplified61.6%
Taylor expanded in b around 0 29.2%
associate-*r/39.1%
Simplified39.1%
if 2.4000000000000001e-163 < t < 9Initial program 98.0%
Taylor expanded in y around 0 80.4%
Taylor expanded in b around 0 34.3%
*-commutative34.3%
exp-to-pow35.6%
sub-neg35.6%
metadata-eval35.6%
+-commutative35.6%
Simplified35.6%
Taylor expanded in t around 0 41.2%
if 9 < t Initial program 100.0%
Taylor expanded in y around 0 95.1%
Taylor expanded in b around inf 42.8%
neg-mul-142.8%
Simplified42.8%
Taylor expanded in b around 0 13.7%
Taylor expanded in b around inf 26.8%
Final simplification35.1%
(FPCore (x y z t a b) :precision binary64 (if (<= t 0.0033) (/ x (* y a)) (* b (/ (- x) y))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 0.0033) {
tmp = x / (y * a);
} else {
tmp = b * (-x / 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 <= 0.0033d0) then
tmp = x / (y * a)
else
tmp = b * (-x / 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 <= 0.0033) {
tmp = x / (y * a);
} else {
tmp = b * (-x / y);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= 0.0033: tmp = x / (y * a) else: tmp = b * (-x / y) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= 0.0033) tmp = Float64(x / Float64(y * a)); else tmp = Float64(b * Float64(Float64(-x) / y)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= 0.0033) tmp = x / (y * a); else tmp = b * (-x / y); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, 0.0033], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(b * N[((-x) / y), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 0.0033:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;b \cdot \frac{-x}{y}\\
\end{array}
\end{array}
if t < 0.0033Initial program 98.3%
Taylor expanded in y around 0 79.7%
Taylor expanded in b around 0 50.4%
*-commutative50.4%
exp-to-pow50.9%
sub-neg50.9%
metadata-eval50.9%
+-commutative50.9%
Simplified50.9%
Taylor expanded in t around 0 34.7%
if 0.0033 < t Initial program 100.0%
Taylor expanded in y around 0 95.1%
Taylor expanded in b around inf 42.8%
neg-mul-142.8%
Simplified42.8%
Taylor expanded in b around 0 13.7%
Taylor expanded in b around inf 26.8%
mul-1-neg26.8%
associate-*r/20.7%
distribute-rgt-neg-in20.7%
distribute-frac-neg20.7%
Simplified20.7%
Final simplification31.5%
(FPCore (x y z t a b) :precision binary64 (if (<= t 5.3) (/ x (* y a)) (/ (- x) (/ y b))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 5.3) {
tmp = x / (y * a);
} else {
tmp = -x / (y / b);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (t <= 5.3d0) then
tmp = x / (y * a)
else
tmp = -x / (y / b)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 5.3) {
tmp = x / (y * a);
} else {
tmp = -x / (y / b);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= 5.3: tmp = x / (y * a) else: tmp = -x / (y / b) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= 5.3) tmp = Float64(x / Float64(y * a)); else tmp = Float64(Float64(-x) / Float64(y / b)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= 5.3) tmp = x / (y * a); else tmp = -x / (y / b); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, 5.3], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[((-x) / N[(y / b), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 5.3:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{-x}{\frac{y}{b}}\\
\end{array}
\end{array}
if t < 5.29999999999999982Initial program 98.3%
Taylor expanded in y around 0 79.7%
Taylor expanded in b around 0 50.4%
*-commutative50.4%
exp-to-pow50.9%
sub-neg50.9%
metadata-eval50.9%
+-commutative50.9%
Simplified50.9%
Taylor expanded in t around 0 34.7%
if 5.29999999999999982 < t Initial program 100.0%
Taylor expanded in y around 0 95.1%
Taylor expanded in b around inf 42.8%
neg-mul-142.8%
Simplified42.8%
Taylor expanded in b around 0 13.7%
Taylor expanded in b around inf 26.8%
mul-1-neg26.8%
associate-*r/20.7%
*-commutative20.7%
associate-/r/20.7%
distribute-neg-frac20.7%
Simplified20.7%
Final simplification31.5%
(FPCore (x y z t a b) :precision binary64 (if (<= t 0.03) (/ x (* y a)) (/ (* x (- b)) y)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= 0.03) {
tmp = x / (y * a);
} else {
tmp = (x * -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 (t <= 0.03d0) then
tmp = x / (y * a)
else
tmp = (x * -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 (t <= 0.03) {
tmp = x / (y * a);
} else {
tmp = (x * -b) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= 0.03: tmp = x / (y * a) else: tmp = (x * -b) / y return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= 0.03) tmp = Float64(x / Float64(y * a)); else tmp = Float64(Float64(x * Float64(-b)) / y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= 0.03) tmp = x / (y * a); else tmp = (x * -b) / y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, 0.03], N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision], N[(N[(x * (-b)), $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 0.03:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \left(-b\right)}{y}\\
\end{array}
\end{array}
if t < 0.029999999999999999Initial program 98.3%
Taylor expanded in y around 0 79.7%
Taylor expanded in b around 0 50.4%
*-commutative50.4%
exp-to-pow50.9%
sub-neg50.9%
metadata-eval50.9%
+-commutative50.9%
Simplified50.9%
Taylor expanded in t around 0 34.7%
if 0.029999999999999999 < t Initial program 100.0%
Taylor expanded in y around 0 95.1%
Taylor expanded in b around inf 42.8%
neg-mul-142.8%
Simplified42.8%
Taylor expanded in b around 0 13.7%
Taylor expanded in b around inf 26.8%
Final simplification32.9%
(FPCore (x y z t a b) :precision binary64 (if (<= b 3.8e-161) (/ 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.8e-161) {
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.8d-161) 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.8e-161) {
tmp = x / (y * a);
} else {
tmp = (x / a) / y;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 3.8e-161: 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.8e-161) 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.8e-161) 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.8e-161], 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.8 \cdot 10^{-161}:\\
\;\;\;\;\frac{x}{y \cdot a}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{a}}{y}\\
\end{array}
\end{array}
if b < 3.8000000000000001e-161Initial program 98.4%
Taylor expanded in y around 0 79.6%
Taylor expanded in b around 0 60.3%
*-commutative60.3%
exp-to-pow60.6%
sub-neg60.6%
metadata-eval60.6%
+-commutative60.6%
Simplified60.6%
Taylor expanded in t around 0 34.0%
if 3.8000000000000001e-161 < b Initial program 99.1%
Taylor expanded in y around 0 89.4%
Taylor expanded in b around 0 55.7%
*-commutative55.7%
exp-to-pow56.2%
sub-neg56.2%
metadata-eval56.2%
+-commutative56.2%
Simplified56.2%
Taylor expanded in t around 0 22.9%
Final simplification29.8%
(FPCore (x y z t a b) :precision binary64 (/ x (* y a)))
double code(double x, double y, double z, double t, double a, double b) {
return x / (y * a);
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x / (y * a)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / (y * a);
}
def code(x, y, z, t, a, b): return x / (y * a)
function code(x, y, z, t, a, b) return Float64(x / Float64(y * a)) end
function tmp = code(x, y, z, t, a, b) tmp = x / (y * a); end
code[x_, y_, z_, t_, a_, b_] := N[(x / N[(y * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y \cdot a}
\end{array}
Initial program 98.7%
Taylor expanded in y around 0 83.3%
Taylor expanded in b around 0 58.6%
*-commutative58.6%
exp-to-pow59.0%
sub-neg59.0%
metadata-eval59.0%
+-commutative59.0%
Simplified59.0%
Taylor expanded in t around 0 28.0%
Final simplification28.0%
(FPCore (x y z t a b) :precision binary64 (/ x y))
double code(double x, double y, double z, double t, double a, double b) {
return x / y;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x / y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / y;
}
def code(x, y, z, t, a, b): return x / y
function code(x, y, z, t, a, b) return Float64(x / y) end
function tmp = code(x, y, z, t, a, b) tmp = x / y; end
code[x_, y_, z_, t_, a_, b_] := N[(x / y), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{y}
\end{array}
Initial program 98.7%
Taylor expanded in y around 0 83.3%
Taylor expanded in b around inf 51.8%
neg-mul-151.8%
Simplified51.8%
Taylor expanded in b around 0 17.0%
Final simplification17.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (pow a (- t 1.0)))
(t_2 (/ (* x (/ t_1 y)) (- (+ b 1.0) (* y (log z))))))
(if (< t -0.8845848504127471)
t_2
(if (< t 852031.2288374073)
(/ (* (/ x y) t_1) (exp (- b (* (log z) y))))
t_2))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = pow(a, (t - 1.0));
double t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * log(z)));
double tmp;
if (t < -0.8845848504127471) {
tmp = t_2;
} else if (t < 852031.2288374073) {
tmp = ((x / y) * t_1) / exp((b - (log(z) * y)));
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = a ** (t - 1.0d0)
t_2 = (x * (t_1 / y)) / ((b + 1.0d0) - (y * log(z)))
if (t < (-0.8845848504127471d0)) then
tmp = t_2
else if (t < 852031.2288374073d0) then
tmp = ((x / y) * t_1) / exp((b - (log(z) * y)))
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = Math.pow(a, (t - 1.0));
double t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * Math.log(z)));
double tmp;
if (t < -0.8845848504127471) {
tmp = t_2;
} else if (t < 852031.2288374073) {
tmp = ((x / y) * t_1) / Math.exp((b - (Math.log(z) * y)));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = math.pow(a, (t - 1.0)) t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * math.log(z))) tmp = 0 if t < -0.8845848504127471: tmp = t_2 elif t < 852031.2288374073: tmp = ((x / y) * t_1) / math.exp((b - (math.log(z) * y))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b) t_1 = a ^ Float64(t - 1.0) t_2 = Float64(Float64(x * Float64(t_1 / y)) / Float64(Float64(b + 1.0) - Float64(y * log(z)))) tmp = 0.0 if (t < -0.8845848504127471) tmp = t_2; elseif (t < 852031.2288374073) tmp = Float64(Float64(Float64(x / y) * t_1) / exp(Float64(b - Float64(log(z) * y)))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = a ^ (t - 1.0); t_2 = (x * (t_1 / y)) / ((b + 1.0) - (y * log(z))); tmp = 0.0; if (t < -0.8845848504127471) tmp = t_2; elseif (t < 852031.2288374073) tmp = ((x / y) * t_1) / exp((b - (log(z) * y))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[Power[a, N[(t - 1.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * N[(t$95$1 / y), $MachinePrecision]), $MachinePrecision] / N[(N[(b + 1.0), $MachinePrecision] - N[(y * N[Log[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -0.8845848504127471], t$95$2, If[Less[t, 852031.2288374073], N[(N[(N[(x / y), $MachinePrecision] * t$95$1), $MachinePrecision] / N[Exp[N[(b - N[(N[Log[z], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {a}^{\left(t - 1\right)}\\
t_2 := \frac{x \cdot \frac{t\_1}{y}}{\left(b + 1\right) - y \cdot \log z}\\
\mathbf{if}\;t < -0.8845848504127471:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t < 852031.2288374073:\\
\;\;\;\;\frac{\frac{x}{y} \cdot t\_1}{e^{b - \log z \cdot y}}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024034
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:precision binary64
:herbie-target
(if (< t -0.8845848504127471) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z)))) (if (< t 852031.2288374073) (/ (* (/ x y) (pow a (- t 1.0))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))