
(FPCore (x y z t a b) :precision binary64 (* x (exp (+ (* y (- (log z) t)) (* a (- (log (- 1.0 z)) b))))))
double code(double x, double y, double z, double t, double a, double b) {
return x * exp(((y * (log(z) - t)) + (a * (log((1.0 - z)) - 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 * exp(((y * (log(z) - t)) + (a * (log((1.0d0 - z)) - b))))
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)) + (a * (Math.log((1.0 - z)) - b))));
}
def code(x, y, z, t, a, b): return x * math.exp(((y * (math.log(z) - t)) + (a * (math.log((1.0 - z)) - b))))
function code(x, y, z, t, a, b) return Float64(x * exp(Float64(Float64(y * Float64(log(z) - t)) + Float64(a * Float64(log(Float64(1.0 - z)) - b))))) end
function tmp = code(x, y, z, t, a, b) tmp = x * exp(((y * (log(z) - t)) + (a * (log((1.0 - z)) - b)))); end
code[x_, y_, z_, t_, a_, b_] := N[(x * N[Exp[N[(N[(y * N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision] + N[(a * N[(N[Log[N[(1.0 - z), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot e^{y \cdot \left(\log z - t\right) + a \cdot \left(\log \left(1 - z\right) - b\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 26 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (* x (exp (+ (* y (- (log z) t)) (* a (- (log (- 1.0 z)) b))))))
double code(double x, double y, double z, double t, double a, double b) {
return x * exp(((y * (log(z) - t)) + (a * (log((1.0 - z)) - 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 * exp(((y * (log(z) - t)) + (a * (log((1.0d0 - z)) - b))))
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)) + (a * (Math.log((1.0 - z)) - b))));
}
def code(x, y, z, t, a, b): return x * math.exp(((y * (math.log(z) - t)) + (a * (math.log((1.0 - z)) - b))))
function code(x, y, z, t, a, b) return Float64(x * exp(Float64(Float64(y * Float64(log(z) - t)) + Float64(a * Float64(log(Float64(1.0 - z)) - b))))) end
function tmp = code(x, y, z, t, a, b) tmp = x * exp(((y * (log(z) - t)) + (a * (log((1.0 - z)) - b)))); end
code[x_, y_, z_, t_, a_, b_] := N[(x * N[Exp[N[(N[(y * N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision] + N[(a * N[(N[Log[N[(1.0 - z), $MachinePrecision]], $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot e^{y \cdot \left(\log z - t\right) + a \cdot \left(\log \left(1 - z\right) - b\right)}
\end{array}
(FPCore (x y z t a b) :precision binary64 (/ x (exp (+ (* a (+ b z)) (* y (- t (log z)))))))
double code(double x, double y, double z, double t, double a, double b) {
return x / exp(((a * (b + z)) + (y * (t - log(z)))));
}
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(((a * (b + z)) + (y * (t - log(z)))))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x / Math.exp(((a * (b + z)) + (y * (t - Math.log(z)))));
}
def code(x, y, z, t, a, b): return x / math.exp(((a * (b + z)) + (y * (t - math.log(z)))))
function code(x, y, z, t, a, b) return Float64(x / exp(Float64(Float64(a * Float64(b + z)) + Float64(y * Float64(t - log(z)))))) end
function tmp = code(x, y, z, t, a, b) tmp = x / exp(((a * (b + z)) + (y * (t - log(z))))); end
code[x_, y_, z_, t_, a_, b_] := N[(x / N[Exp[N[(N[(a * N[(b + z), $MachinePrecision]), $MachinePrecision] + N[(y * N[(t - N[Log[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{e^{a \cdot \left(b + z\right) + y \cdot \left(t - \log z\right)}}
\end{array}
Initial program 95.8%
Simplified95.8%
Taylor expanded in z around 0
associate-+r+N/A
+-lowering-+.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6499.2%
Simplified99.2%
(FPCore (x y z t a b) :precision binary64 (if (<= a -1.3e+82) (/ x (exp (* a (- b (log (- 1.0 z)))))) (if (<= a 1.02e+93) (/ x (exp (* y (- t (log z))))) (/ x (exp (* a b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (a <= -1.3e+82) {
tmp = x / exp((a * (b - log((1.0 - z)))));
} else if (a <= 1.02e+93) {
tmp = x / exp((y * (t - log(z))));
} else {
tmp = x / exp((a * b));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (a <= (-1.3d+82)) then
tmp = x / exp((a * (b - log((1.0d0 - z)))))
else if (a <= 1.02d+93) then
tmp = x / exp((y * (t - log(z))))
else
tmp = x / exp((a * b))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (a <= -1.3e+82) {
tmp = x / Math.exp((a * (b - Math.log((1.0 - z)))));
} else if (a <= 1.02e+93) {
tmp = x / Math.exp((y * (t - Math.log(z))));
} else {
tmp = x / Math.exp((a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if a <= -1.3e+82: tmp = x / math.exp((a * (b - math.log((1.0 - z))))) elif a <= 1.02e+93: tmp = x / math.exp((y * (t - math.log(z)))) else: tmp = x / math.exp((a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (a <= -1.3e+82) tmp = Float64(x / exp(Float64(a * Float64(b - log(Float64(1.0 - z)))))); elseif (a <= 1.02e+93) tmp = Float64(x / exp(Float64(y * Float64(t - log(z))))); else tmp = Float64(x / exp(Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (a <= -1.3e+82) tmp = x / exp((a * (b - log((1.0 - z))))); elseif (a <= 1.02e+93) tmp = x / exp((y * (t - log(z)))); else tmp = x / exp((a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[a, -1.3e+82], N[(x / N[Exp[N[(a * N[(b - N[Log[N[(1.0 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.02e+93], N[(x / N[Exp[N[(y * N[(t - N[Log[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(x / N[Exp[N[(a * b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.3 \cdot 10^{+82}:\\
\;\;\;\;\frac{x}{e^{a \cdot \left(b - \log \left(1 - z\right)\right)}}\\
\mathbf{elif}\;a \leq 1.02 \cdot 10^{+93}:\\
\;\;\;\;\frac{x}{e^{y \cdot \left(t - \log z\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\end{array}
\end{array}
if a < -1.2999999999999999e82Initial program 93.3%
Simplified93.3%
Taylor expanded in y around 0
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f64N/A
--lowering--.f6480.0%
Simplified80.0%
if -1.2999999999999999e82 < a < 1.0200000000000001e93Initial program 98.3%
Simplified98.4%
Taylor expanded in a around 0
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6489.2%
Simplified89.2%
if 1.0200000000000001e93 < a Initial program 88.5%
Simplified88.5%
Taylor expanded in b around inf
*-lowering-*.f6477.4%
Simplified77.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (exp (* a b)))))
(if (<= a -8.2e+80)
t_1
(if (<= a 9.2e+93) (/ x (exp (* y (- t (log z))))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / exp((a * b));
double tmp;
if (a <= -8.2e+80) {
tmp = t_1;
} else if (a <= 9.2e+93) {
tmp = x / exp((y * (t - log(z))));
} 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 / exp((a * b))
if (a <= (-8.2d+80)) then
tmp = t_1
else if (a <= 9.2d+93) then
tmp = x / exp((y * (t - log(z))))
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.exp((a * b));
double tmp;
if (a <= -8.2e+80) {
tmp = t_1;
} else if (a <= 9.2e+93) {
tmp = x / Math.exp((y * (t - Math.log(z))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x / math.exp((a * b)) tmp = 0 if a <= -8.2e+80: tmp = t_1 elif a <= 9.2e+93: tmp = x / math.exp((y * (t - math.log(z)))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x / exp(Float64(a * b))) tmp = 0.0 if (a <= -8.2e+80) tmp = t_1; elseif (a <= 9.2e+93) tmp = Float64(x / exp(Float64(y * Float64(t - log(z))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x / exp((a * b)); tmp = 0.0; if (a <= -8.2e+80) tmp = t_1; elseif (a <= 9.2e+93) tmp = x / exp((y * (t - log(z)))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[Exp[N[(a * b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -8.2e+80], t$95$1, If[LessEqual[a, 9.2e+93], N[(x / N[Exp[N[(y * N[(t - N[Log[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{e^{a \cdot b}}\\
\mathbf{if}\;a \leq -8.2 \cdot 10^{+80}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 9.2 \cdot 10^{+93}:\\
\;\;\;\;\frac{x}{e^{y \cdot \left(t - \log z\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -8.20000000000000003e80 or 9.2000000000000006e93 < a Initial program 90.9%
Simplified90.9%
Taylor expanded in b around inf
*-lowering-*.f6477.6%
Simplified77.6%
if -8.20000000000000003e80 < a < 9.2000000000000006e93Initial program 98.3%
Simplified98.4%
Taylor expanded in a around 0
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6489.2%
Simplified89.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (exp (* y t)))))
(if (<= t -2.3e+110)
t_1
(if (<= t -7.1e-156)
(/ x (exp (* a b)))
(if (<= t 4.8e-5) (/ x (/ 1.0 (pow z y))) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / exp((y * t));
double tmp;
if (t <= -2.3e+110) {
tmp = t_1;
} else if (t <= -7.1e-156) {
tmp = x / exp((a * b));
} else if (t <= 4.8e-5) {
tmp = x / (1.0 / pow(z, 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 / exp((y * t))
if (t <= (-2.3d+110)) then
tmp = t_1
else if (t <= (-7.1d-156)) then
tmp = x / exp((a * b))
else if (t <= 4.8d-5) then
tmp = x / (1.0d0 / (z ** y))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / Math.exp((y * t));
double tmp;
if (t <= -2.3e+110) {
tmp = t_1;
} else if (t <= -7.1e-156) {
tmp = x / Math.exp((a * b));
} else if (t <= 4.8e-5) {
tmp = x / (1.0 / Math.pow(z, y));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x / math.exp((y * t)) tmp = 0 if t <= -2.3e+110: tmp = t_1 elif t <= -7.1e-156: tmp = x / math.exp((a * b)) elif t <= 4.8e-5: tmp = x / (1.0 / math.pow(z, y)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x / exp(Float64(y * t))) tmp = 0.0 if (t <= -2.3e+110) tmp = t_1; elseif (t <= -7.1e-156) tmp = Float64(x / exp(Float64(a * b))); elseif (t <= 4.8e-5) tmp = Float64(x / Float64(1.0 / (z ^ y))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x / exp((y * t)); tmp = 0.0; if (t <= -2.3e+110) tmp = t_1; elseif (t <= -7.1e-156) tmp = x / exp((a * b)); elseif (t <= 4.8e-5) tmp = x / (1.0 / (z ^ y)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[Exp[N[(y * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -2.3e+110], t$95$1, If[LessEqual[t, -7.1e-156], N[(x / N[Exp[N[(a * b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.8e-5], N[(x / N[(1.0 / N[Power[z, y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{e^{y \cdot t}}\\
\mathbf{if}\;t \leq -2.3 \cdot 10^{+110}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -7.1 \cdot 10^{-156}:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\mathbf{elif}\;t \leq 4.8 \cdot 10^{-5}:\\
\;\;\;\;\frac{x}{\frac{1}{{z}^{y}}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -2.3e110 or 4.8000000000000001e-5 < t Initial program 97.7%
Simplified97.7%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6485.5%
Simplified85.5%
if -2.3e110 < t < -7.1000000000000005e-156Initial program 91.9%
Simplified91.9%
Taylor expanded in b around inf
*-lowering-*.f6471.4%
Simplified71.4%
if -7.1000000000000005e-156 < t < 4.8000000000000001e-5Initial program 96.3%
Simplified96.4%
Taylor expanded in a around 0
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6472.7%
Simplified72.7%
Taylor expanded in t around 0
mul-1-negN/A
exp-negN/A
/-lowering-/.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f6472.7%
Simplified72.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (exp (* y t)))))
(if (<= t -3.6e+112)
t_1
(if (<= t -1.7e-155)
(/ x (exp (* a b)))
(if (<= t 4.9e-5) (* x (pow z y)) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / exp((y * t));
double tmp;
if (t <= -3.6e+112) {
tmp = t_1;
} else if (t <= -1.7e-155) {
tmp = x / exp((a * b));
} else if (t <= 4.9e-5) {
tmp = x * pow(z, 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 / exp((y * t))
if (t <= (-3.6d+112)) then
tmp = t_1
else if (t <= (-1.7d-155)) then
tmp = x / exp((a * b))
else if (t <= 4.9d-5) then
tmp = x * (z ** y)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x / Math.exp((y * t));
double tmp;
if (t <= -3.6e+112) {
tmp = t_1;
} else if (t <= -1.7e-155) {
tmp = x / Math.exp((a * b));
} else if (t <= 4.9e-5) {
tmp = x * Math.pow(z, y);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x / math.exp((y * t)) tmp = 0 if t <= -3.6e+112: tmp = t_1 elif t <= -1.7e-155: tmp = x / math.exp((a * b)) elif t <= 4.9e-5: tmp = x * math.pow(z, y) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x / exp(Float64(y * t))) tmp = 0.0 if (t <= -3.6e+112) tmp = t_1; elseif (t <= -1.7e-155) tmp = Float64(x / exp(Float64(a * b))); elseif (t <= 4.9e-5) tmp = Float64(x * (z ^ y)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x / exp((y * t)); tmp = 0.0; if (t <= -3.6e+112) tmp = t_1; elseif (t <= -1.7e-155) tmp = x / exp((a * b)); elseif (t <= 4.9e-5) tmp = x * (z ^ y); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x / N[Exp[N[(y * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3.6e+112], t$95$1, If[LessEqual[t, -1.7e-155], N[(x / N[Exp[N[(a * b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.9e-5], N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{e^{y \cdot t}}\\
\mathbf{if}\;t \leq -3.6 \cdot 10^{+112}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -1.7 \cdot 10^{-155}:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\mathbf{elif}\;t \leq 4.9 \cdot 10^{-5}:\\
\;\;\;\;x \cdot {z}^{y}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -3.6e112 or 4.9e-5 < t Initial program 97.7%
Simplified97.7%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6485.5%
Simplified85.5%
if -3.6e112 < t < -1.7e-155Initial program 91.9%
Simplified91.9%
Taylor expanded in b around inf
*-lowering-*.f6471.4%
Simplified71.4%
if -1.7e-155 < t < 4.9e-5Initial program 96.3%
Simplified96.4%
Taylor expanded in a around 0
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6472.7%
Simplified72.7%
Taylor expanded in t around 0
mul-1-negN/A
exp-negN/A
/-lowering-/.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f6472.7%
Simplified72.7%
clear-numN/A
associate-/r/N/A
remove-double-divN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6472.7%
Applied egg-rr72.7%
Final simplification77.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -1.15e+169)
(* x (+ 1.0 (* t (- (* y (* y (* t 0.5))) y))))
(if (<= t -1600000000.0)
(/
x
(+
1.0
(*
y
(+
t
(*
y
(+ (* 0.16666666666666666 (* t (* y (* t t)))) (* 0.5 (* t t))))))))
(if (<= t 6.9e+177) (* x (pow z y)) (* x (* (* t 0.5) (* t (* y y))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.15e+169) {
tmp = x * (1.0 + (t * ((y * (y * (t * 0.5))) - y)));
} else if (t <= -1600000000.0) {
tmp = x / (1.0 + (y * (t + (y * ((0.16666666666666666 * (t * (y * (t * t)))) + (0.5 * (t * t)))))));
} else if (t <= 6.9e+177) {
tmp = x * pow(z, y);
} else {
tmp = x * ((t * 0.5) * (t * (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 (t <= (-1.15d+169)) then
tmp = x * (1.0d0 + (t * ((y * (y * (t * 0.5d0))) - y)))
else if (t <= (-1600000000.0d0)) then
tmp = x / (1.0d0 + (y * (t + (y * ((0.16666666666666666d0 * (t * (y * (t * t)))) + (0.5d0 * (t * t)))))))
else if (t <= 6.9d+177) then
tmp = x * (z ** y)
else
tmp = x * ((t * 0.5d0) * (t * (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 (t <= -1.15e+169) {
tmp = x * (1.0 + (t * ((y * (y * (t * 0.5))) - y)));
} else if (t <= -1600000000.0) {
tmp = x / (1.0 + (y * (t + (y * ((0.16666666666666666 * (t * (y * (t * t)))) + (0.5 * (t * t)))))));
} else if (t <= 6.9e+177) {
tmp = x * Math.pow(z, y);
} else {
tmp = x * ((t * 0.5) * (t * (y * y)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.15e+169: tmp = x * (1.0 + (t * ((y * (y * (t * 0.5))) - y))) elif t <= -1600000000.0: tmp = x / (1.0 + (y * (t + (y * ((0.16666666666666666 * (t * (y * (t * t)))) + (0.5 * (t * t))))))) elif t <= 6.9e+177: tmp = x * math.pow(z, y) else: tmp = x * ((t * 0.5) * (t * (y * y))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.15e+169) tmp = Float64(x * Float64(1.0 + Float64(t * Float64(Float64(y * Float64(y * Float64(t * 0.5))) - y)))); elseif (t <= -1600000000.0) tmp = Float64(x / Float64(1.0 + Float64(y * Float64(t + Float64(y * Float64(Float64(0.16666666666666666 * Float64(t * Float64(y * Float64(t * t)))) + Float64(0.5 * Float64(t * t)))))))); elseif (t <= 6.9e+177) tmp = Float64(x * (z ^ y)); else tmp = Float64(x * Float64(Float64(t * 0.5) * Float64(t * Float64(y * y)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.15e+169) tmp = x * (1.0 + (t * ((y * (y * (t * 0.5))) - y))); elseif (t <= -1600000000.0) tmp = x / (1.0 + (y * (t + (y * ((0.16666666666666666 * (t * (y * (t * t)))) + (0.5 * (t * t))))))); elseif (t <= 6.9e+177) tmp = x * (z ^ y); else tmp = x * ((t * 0.5) * (t * (y * y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.15e+169], N[(x * N[(1.0 + N[(t * N[(N[(y * N[(y * N[(t * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, -1600000000.0], N[(x / N[(1.0 + N[(y * N[(t + N[(y * N[(N[(0.16666666666666666 * N[(t * N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6.9e+177], N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.15 \cdot 10^{+169}:\\
\;\;\;\;x \cdot \left(1 + t \cdot \left(y \cdot \left(y \cdot \left(t \cdot 0.5\right)\right) - y\right)\right)\\
\mathbf{elif}\;t \leq -1600000000:\\
\;\;\;\;\frac{x}{1 + y \cdot \left(t + y \cdot \left(0.16666666666666666 \cdot \left(t \cdot \left(y \cdot \left(t \cdot t\right)\right)\right) + 0.5 \cdot \left(t \cdot t\right)\right)\right)}\\
\mathbf{elif}\;t \leq 6.9 \cdot 10^{+177}:\\
\;\;\;\;x \cdot {z}^{y}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\right)\\
\end{array}
\end{array}
if t < -1.15e169Initial program 96.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6490.7%
Simplified90.7%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6462.3%
Simplified62.3%
if -1.15e169 < t < -1.6e9Initial program 93.9%
Simplified94.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6467.6%
Simplified67.6%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.4%
Simplified58.4%
if -1.6e9 < t < 6.89999999999999971e177Initial program 96.0%
Simplified96.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6472.1%
Simplified72.1%
Taylor expanded in t around 0
mul-1-negN/A
exp-negN/A
/-lowering-/.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f6468.6%
Simplified68.6%
clear-numN/A
associate-/r/N/A
remove-double-divN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6468.6%
Applied egg-rr68.6%
if 6.89999999999999971e177 < t Initial program 96.2%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6478.1%
Simplified78.1%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6433.1%
Simplified33.1%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6462.1%
Simplified62.1%
Final simplification65.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (pow z y))))
(if (<= y -2.8e+22)
t_1
(if (<= y 8800000000000.0) (/ x (exp (* a b))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * pow(z, y);
double tmp;
if (y <= -2.8e+22) {
tmp = t_1;
} else if (y <= 8800000000000.0) {
tmp = x / exp((a * b));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = x * (z ** y)
if (y <= (-2.8d+22)) then
tmp = t_1
else if (y <= 8800000000000.0d0) then
tmp = x / exp((a * b))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * Math.pow(z, y);
double tmp;
if (y <= -2.8e+22) {
tmp = t_1;
} else if (y <= 8800000000000.0) {
tmp = x / Math.exp((a * b));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * math.pow(z, y) tmp = 0 if y <= -2.8e+22: tmp = t_1 elif y <= 8800000000000.0: tmp = x / math.exp((a * b)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * (z ^ y)) tmp = 0.0 if (y <= -2.8e+22) tmp = t_1; elseif (y <= 8800000000000.0) tmp = Float64(x / exp(Float64(a * b))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * (z ^ y); tmp = 0.0; if (y <= -2.8e+22) tmp = t_1; elseif (y <= 8800000000000.0) tmp = x / exp((a * b)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.8e+22], t$95$1, If[LessEqual[y, 8800000000000.0], N[(x / N[Exp[N[(a * b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot {z}^{y}\\
\mathbf{if}\;y \leq -2.8 \cdot 10^{+22}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 8800000000000:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.8e22 or 8.8e12 < y Initial program 98.0%
Simplified98.0%
Taylor expanded in a around 0
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6494.1%
Simplified94.1%
Taylor expanded in t around 0
mul-1-negN/A
exp-negN/A
/-lowering-/.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f6474.4%
Simplified74.4%
clear-numN/A
associate-/r/N/A
remove-double-divN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6474.4%
Applied egg-rr74.4%
if -2.8e22 < y < 8.8e12Initial program 94.4%
Simplified94.5%
Taylor expanded in b around inf
*-lowering-*.f6473.0%
Simplified73.0%
Final simplification73.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (* t t))) (t_2 (* (* t 0.5) (* t (* y y)))))
(if (<= y -9.2e+219)
(* x t_2)
(if (<= y -8e+22)
(/ (/ (* x 2.0) t_1) y)
(if (<= y -1e-21)
(* (* 0.5 (* t t)) (* x (* y y)))
(if (<= y 7e-41)
(/ x (+ 1.0 (* y (* 0.5 t_1))))
(if (<= y 5.6e+143)
(* (* x (* t (* t t))) (* (* y (* y y)) -0.16666666666666666))
(/ x t_2))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * (t * t);
double t_2 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_2;
} else if (y <= -8e+22) {
tmp = ((x * 2.0) / t_1) / y;
} else if (y <= -1e-21) {
tmp = (0.5 * (t * t)) * (x * (y * y));
} else if (y <= 7e-41) {
tmp = x / (1.0 + (y * (0.5 * t_1)));
} else if (y <= 5.6e+143) {
tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666);
} else {
tmp = x / 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 = y * (t * t)
t_2 = (t * 0.5d0) * (t * (y * y))
if (y <= (-9.2d+219)) then
tmp = x * t_2
else if (y <= (-8d+22)) then
tmp = ((x * 2.0d0) / t_1) / y
else if (y <= (-1d-21)) then
tmp = (0.5d0 * (t * t)) * (x * (y * y))
else if (y <= 7d-41) then
tmp = x / (1.0d0 + (y * (0.5d0 * t_1)))
else if (y <= 5.6d+143) then
tmp = (x * (t * (t * t))) * ((y * (y * y)) * (-0.16666666666666666d0))
else
tmp = x / 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 = y * (t * t);
double t_2 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_2;
} else if (y <= -8e+22) {
tmp = ((x * 2.0) / t_1) / y;
} else if (y <= -1e-21) {
tmp = (0.5 * (t * t)) * (x * (y * y));
} else if (y <= 7e-41) {
tmp = x / (1.0 + (y * (0.5 * t_1)));
} else if (y <= 5.6e+143) {
tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666);
} else {
tmp = x / t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * (t * t) t_2 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -9.2e+219: tmp = x * t_2 elif y <= -8e+22: tmp = ((x * 2.0) / t_1) / y elif y <= -1e-21: tmp = (0.5 * (t * t)) * (x * (y * y)) elif y <= 7e-41: tmp = x / (1.0 + (y * (0.5 * t_1))) elif y <= 5.6e+143: tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666) else: tmp = x / t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y * Float64(t * t)) t_2 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -9.2e+219) tmp = Float64(x * t_2); elseif (y <= -8e+22) tmp = Float64(Float64(Float64(x * 2.0) / t_1) / y); elseif (y <= -1e-21) tmp = Float64(Float64(0.5 * Float64(t * t)) * Float64(x * Float64(y * y))); elseif (y <= 7e-41) tmp = Float64(x / Float64(1.0 + Float64(y * Float64(0.5 * t_1)))); elseif (y <= 5.6e+143) tmp = Float64(Float64(x * Float64(t * Float64(t * t))) * Float64(Float64(y * Float64(y * y)) * -0.16666666666666666)); else tmp = Float64(x / t_2); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * (t * t); t_2 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -9.2e+219) tmp = x * t_2; elseif (y <= -8e+22) tmp = ((x * 2.0) / t_1) / y; elseif (y <= -1e-21) tmp = (0.5 * (t * t)) * (x * (y * y)); elseif (y <= 7e-41) tmp = x / (1.0 + (y * (0.5 * t_1))); elseif (y <= 5.6e+143) tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666); else tmp = x / t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9.2e+219], N[(x * t$95$2), $MachinePrecision], If[LessEqual[y, -8e+22], N[(N[(N[(x * 2.0), $MachinePrecision] / t$95$1), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, -1e-21], N[(N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision] * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 7e-41], N[(x / N[(1.0 + N[(y * N[(0.5 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5.6e+143], N[(N[(x * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision], N[(x / t$95$2), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(t \cdot t\right)\\
t_2 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -9.2 \cdot 10^{+219}:\\
\;\;\;\;x \cdot t\_2\\
\mathbf{elif}\;y \leq -8 \cdot 10^{+22}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{t\_1}}{y}\\
\mathbf{elif}\;y \leq -1 \cdot 10^{-21}:\\
\;\;\;\;\left(0.5 \cdot \left(t \cdot t\right)\right) \cdot \left(x \cdot \left(y \cdot y\right)\right)\\
\mathbf{elif}\;y \leq 7 \cdot 10^{-41}:\\
\;\;\;\;\frac{x}{1 + y \cdot \left(0.5 \cdot t\_1\right)}\\
\mathbf{elif}\;y \leq 5.6 \cdot 10^{+143}:\\
\;\;\;\;\left(x \cdot \left(t \cdot \left(t \cdot t\right)\right)\right) \cdot \left(\left(y \cdot \left(y \cdot y\right)\right) \cdot -0.16666666666666666\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_2}\\
\end{array}
\end{array}
if y < -9.2000000000000004e219Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -9.2000000000000004e219 < y < -8e22Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.5%
Simplified50.5%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.6%
Simplified65.6%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6468.5%
Applied egg-rr68.5%
if -8e22 < y < -9.99999999999999908e-22Initial program 91.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6436.0%
Simplified36.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6427.8%
Simplified27.8%
Taylor expanded in t around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.2%
Simplified59.2%
if -9.99999999999999908e-22 < y < 6.9999999999999999e-41Initial program 93.8%
Simplified93.8%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6463.3%
Simplified63.3%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.1%
Simplified53.1%
Taylor expanded in t around inf
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.1%
Simplified53.1%
if 6.9999999999999999e-41 < y < 5.59999999999999996e143Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6451.4%
Simplified51.4%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified25.1%
Taylor expanded in y around inf
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.3%
Simplified52.3%
if 5.59999999999999996e143 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification60.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (* t t))) (t_2 (* (* t 0.5) (* t (* y y)))))
(if (<= y -9.2e+219)
(* x t_2)
(if (<= y -3.9e+22)
(/ (/ (* x 2.0) t_1) y)
(if (<= y -2.35e-21)
(* (* 0.5 (* t t)) (* x (* y y)))
(if (<= y 2.3e-119)
(/ x (+ 1.0 (* y (* 0.5 t_1))))
(if (<= y 4e+146)
(* x (* y (* y (* y (* t (* (* t t) -0.16666666666666666))))))
(/ x t_2))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * (t * t);
double t_2 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_2;
} else if (y <= -3.9e+22) {
tmp = ((x * 2.0) / t_1) / y;
} else if (y <= -2.35e-21) {
tmp = (0.5 * (t * t)) * (x * (y * y));
} else if (y <= 2.3e-119) {
tmp = x / (1.0 + (y * (0.5 * t_1)));
} else if (y <= 4e+146) {
tmp = x * (y * (y * (y * (t * ((t * t) * -0.16666666666666666)))));
} else {
tmp = x / 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 = y * (t * t)
t_2 = (t * 0.5d0) * (t * (y * y))
if (y <= (-9.2d+219)) then
tmp = x * t_2
else if (y <= (-3.9d+22)) then
tmp = ((x * 2.0d0) / t_1) / y
else if (y <= (-2.35d-21)) then
tmp = (0.5d0 * (t * t)) * (x * (y * y))
else if (y <= 2.3d-119) then
tmp = x / (1.0d0 + (y * (0.5d0 * t_1)))
else if (y <= 4d+146) then
tmp = x * (y * (y * (y * (t * ((t * t) * (-0.16666666666666666d0))))))
else
tmp = x / 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 = y * (t * t);
double t_2 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_2;
} else if (y <= -3.9e+22) {
tmp = ((x * 2.0) / t_1) / y;
} else if (y <= -2.35e-21) {
tmp = (0.5 * (t * t)) * (x * (y * y));
} else if (y <= 2.3e-119) {
tmp = x / (1.0 + (y * (0.5 * t_1)));
} else if (y <= 4e+146) {
tmp = x * (y * (y * (y * (t * ((t * t) * -0.16666666666666666)))));
} else {
tmp = x / t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * (t * t) t_2 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -9.2e+219: tmp = x * t_2 elif y <= -3.9e+22: tmp = ((x * 2.0) / t_1) / y elif y <= -2.35e-21: tmp = (0.5 * (t * t)) * (x * (y * y)) elif y <= 2.3e-119: tmp = x / (1.0 + (y * (0.5 * t_1))) elif y <= 4e+146: tmp = x * (y * (y * (y * (t * ((t * t) * -0.16666666666666666))))) else: tmp = x / t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y * Float64(t * t)) t_2 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -9.2e+219) tmp = Float64(x * t_2); elseif (y <= -3.9e+22) tmp = Float64(Float64(Float64(x * 2.0) / t_1) / y); elseif (y <= -2.35e-21) tmp = Float64(Float64(0.5 * Float64(t * t)) * Float64(x * Float64(y * y))); elseif (y <= 2.3e-119) tmp = Float64(x / Float64(1.0 + Float64(y * Float64(0.5 * t_1)))); elseif (y <= 4e+146) tmp = Float64(x * Float64(y * Float64(y * Float64(y * Float64(t * Float64(Float64(t * t) * -0.16666666666666666)))))); else tmp = Float64(x / t_2); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * (t * t); t_2 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -9.2e+219) tmp = x * t_2; elseif (y <= -3.9e+22) tmp = ((x * 2.0) / t_1) / y; elseif (y <= -2.35e-21) tmp = (0.5 * (t * t)) * (x * (y * y)); elseif (y <= 2.3e-119) tmp = x / (1.0 + (y * (0.5 * t_1))); elseif (y <= 4e+146) tmp = x * (y * (y * (y * (t * ((t * t) * -0.16666666666666666))))); else tmp = x / t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9.2e+219], N[(x * t$95$2), $MachinePrecision], If[LessEqual[y, -3.9e+22], N[(N[(N[(x * 2.0), $MachinePrecision] / t$95$1), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, -2.35e-21], N[(N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision] * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.3e-119], N[(x / N[(1.0 + N[(y * N[(0.5 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4e+146], N[(x * N[(y * N[(y * N[(y * N[(t * N[(N[(t * t), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / t$95$2), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(t \cdot t\right)\\
t_2 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -9.2 \cdot 10^{+219}:\\
\;\;\;\;x \cdot t\_2\\
\mathbf{elif}\;y \leq -3.9 \cdot 10^{+22}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{t\_1}}{y}\\
\mathbf{elif}\;y \leq -2.35 \cdot 10^{-21}:\\
\;\;\;\;\left(0.5 \cdot \left(t \cdot t\right)\right) \cdot \left(x \cdot \left(y \cdot y\right)\right)\\
\mathbf{elif}\;y \leq 2.3 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + y \cdot \left(0.5 \cdot t\_1\right)}\\
\mathbf{elif}\;y \leq 4 \cdot 10^{+146}:\\
\;\;\;\;x \cdot \left(y \cdot \left(y \cdot \left(y \cdot \left(t \cdot \left(\left(t \cdot t\right) \cdot -0.16666666666666666\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_2}\\
\end{array}
\end{array}
if y < -9.2000000000000004e219Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -9.2000000000000004e219 < y < -3.90000000000000021e22Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.5%
Simplified50.5%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.6%
Simplified65.6%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6468.5%
Applied egg-rr68.5%
if -3.90000000000000021e22 < y < -2.35000000000000015e-21Initial program 91.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6436.0%
Simplified36.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6427.8%
Simplified27.8%
Taylor expanded in t around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.2%
Simplified59.2%
if -2.35000000000000015e-21 < y < 2.29999999999999993e-119Initial program 93.2%
Simplified93.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6461.8%
Simplified61.8%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.8%
Simplified54.8%
Taylor expanded in t around inf
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.8%
Simplified54.8%
if 2.29999999999999993e-119 < y < 3.99999999999999973e146Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.2%
Simplified57.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified32.3%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.3%
Simplified48.3%
if 3.99999999999999973e146 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification60.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* y (* t t)))
(t_2 (* (* t 0.5) (* t (* y y))))
(t_3 (* 0.5 (* t t))))
(if (<= y -4.5e+220)
(* x t_2)
(if (<= y -5e+24)
(/ (/ (* x 2.0) t_1) y)
(if (<= y -2.35e-21)
(* t_3 (* x (* y y)))
(if (<= y 1.96e-119)
(/ x (+ 1.0 (* y (* 0.5 t_1))))
(if (<= y 1.7e+146) (* y (* y (* x t_3))) (/ x t_2))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = y * (t * t);
double t_2 = (t * 0.5) * (t * (y * y));
double t_3 = 0.5 * (t * t);
double tmp;
if (y <= -4.5e+220) {
tmp = x * t_2;
} else if (y <= -5e+24) {
tmp = ((x * 2.0) / t_1) / y;
} else if (y <= -2.35e-21) {
tmp = t_3 * (x * (y * y));
} else if (y <= 1.96e-119) {
tmp = x / (1.0 + (y * (0.5 * t_1)));
} else if (y <= 1.7e+146) {
tmp = y * (y * (x * t_3));
} else {
tmp = x / 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) :: t_3
real(8) :: tmp
t_1 = y * (t * t)
t_2 = (t * 0.5d0) * (t * (y * y))
t_3 = 0.5d0 * (t * t)
if (y <= (-4.5d+220)) then
tmp = x * t_2
else if (y <= (-5d+24)) then
tmp = ((x * 2.0d0) / t_1) / y
else if (y <= (-2.35d-21)) then
tmp = t_3 * (x * (y * y))
else if (y <= 1.96d-119) then
tmp = x / (1.0d0 + (y * (0.5d0 * t_1)))
else if (y <= 1.7d+146) then
tmp = y * (y * (x * t_3))
else
tmp = x / 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 = y * (t * t);
double t_2 = (t * 0.5) * (t * (y * y));
double t_3 = 0.5 * (t * t);
double tmp;
if (y <= -4.5e+220) {
tmp = x * t_2;
} else if (y <= -5e+24) {
tmp = ((x * 2.0) / t_1) / y;
} else if (y <= -2.35e-21) {
tmp = t_3 * (x * (y * y));
} else if (y <= 1.96e-119) {
tmp = x / (1.0 + (y * (0.5 * t_1)));
} else if (y <= 1.7e+146) {
tmp = y * (y * (x * t_3));
} else {
tmp = x / t_2;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = y * (t * t) t_2 = (t * 0.5) * (t * (y * y)) t_3 = 0.5 * (t * t) tmp = 0 if y <= -4.5e+220: tmp = x * t_2 elif y <= -5e+24: tmp = ((x * 2.0) / t_1) / y elif y <= -2.35e-21: tmp = t_3 * (x * (y * y)) elif y <= 1.96e-119: tmp = x / (1.0 + (y * (0.5 * t_1))) elif y <= 1.7e+146: tmp = y * (y * (x * t_3)) else: tmp = x / t_2 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(y * Float64(t * t)) t_2 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) t_3 = Float64(0.5 * Float64(t * t)) tmp = 0.0 if (y <= -4.5e+220) tmp = Float64(x * t_2); elseif (y <= -5e+24) tmp = Float64(Float64(Float64(x * 2.0) / t_1) / y); elseif (y <= -2.35e-21) tmp = Float64(t_3 * Float64(x * Float64(y * y))); elseif (y <= 1.96e-119) tmp = Float64(x / Float64(1.0 + Float64(y * Float64(0.5 * t_1)))); elseif (y <= 1.7e+146) tmp = Float64(y * Float64(y * Float64(x * t_3))); else tmp = Float64(x / t_2); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = y * (t * t); t_2 = (t * 0.5) * (t * (y * y)); t_3 = 0.5 * (t * t); tmp = 0.0; if (y <= -4.5e+220) tmp = x * t_2; elseif (y <= -5e+24) tmp = ((x * 2.0) / t_1) / y; elseif (y <= -2.35e-21) tmp = t_3 * (x * (y * y)); elseif (y <= 1.96e-119) tmp = x / (1.0 + (y * (0.5 * t_1))); elseif (y <= 1.7e+146) tmp = y * (y * (x * t_3)); else tmp = x / t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -4.5e+220], N[(x * t$95$2), $MachinePrecision], If[LessEqual[y, -5e+24], N[(N[(N[(x * 2.0), $MachinePrecision] / t$95$1), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, -2.35e-21], N[(t$95$3 * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.96e-119], N[(x / N[(1.0 + N[(y * N[(0.5 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.7e+146], N[(y * N[(y * N[(x * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / t$95$2), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(t \cdot t\right)\\
t_2 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
t_3 := 0.5 \cdot \left(t \cdot t\right)\\
\mathbf{if}\;y \leq -4.5 \cdot 10^{+220}:\\
\;\;\;\;x \cdot t\_2\\
\mathbf{elif}\;y \leq -5 \cdot 10^{+24}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{t\_1}}{y}\\
\mathbf{elif}\;y \leq -2.35 \cdot 10^{-21}:\\
\;\;\;\;t\_3 \cdot \left(x \cdot \left(y \cdot y\right)\right)\\
\mathbf{elif}\;y \leq 1.96 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + y \cdot \left(0.5 \cdot t\_1\right)}\\
\mathbf{elif}\;y \leq 1.7 \cdot 10^{+146}:\\
\;\;\;\;y \cdot \left(y \cdot \left(x \cdot t\_3\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_2}\\
\end{array}
\end{array}
if y < -4.50000000000000011e220Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -4.50000000000000011e220 < y < -5.00000000000000045e24Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.5%
Simplified50.5%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.6%
Simplified65.6%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6468.5%
Applied egg-rr68.5%
if -5.00000000000000045e24 < y < -2.35000000000000015e-21Initial program 91.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6436.0%
Simplified36.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6427.8%
Simplified27.8%
Taylor expanded in t around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.2%
Simplified59.2%
if -2.35000000000000015e-21 < y < 1.9600000000000001e-119Initial program 93.2%
Simplified93.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6461.8%
Simplified61.8%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.8%
Simplified54.8%
Taylor expanded in t around inf
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.8%
Simplified54.8%
if 1.9600000000000001e-119 < y < 1.69999999999999995e146Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.2%
Simplified57.2%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6432.5%
Simplified32.5%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr26.5%
Taylor expanded in t around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.9%
Simplified44.9%
if 1.69999999999999995e146 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification59.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (* t 0.5) (* t (* y y)))))
(if (<= y -1.56e+220)
(* x t_1)
(if (<= y -3.65e+22)
(/ (/ (* x 2.0) (* y (* t t))) y)
(if (<= y 7e-41)
(/ x (+ 1.0 (* b (+ a (* 0.5 (* b (* a a)))))))
(if (<= y 4.6e+144)
(* (* x (* t (* t t))) (* (* y (* y y)) -0.16666666666666666))
(/ x t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -1.56e+220) {
tmp = x * t_1;
} else if (y <= -3.65e+22) {
tmp = ((x * 2.0) / (y * (t * t))) / y;
} else if (y <= 7e-41) {
tmp = x / (1.0 + (b * (a + (0.5 * (b * (a * a))))));
} else if (y <= 4.6e+144) {
tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666);
} else {
tmp = x / 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 = (t * 0.5d0) * (t * (y * y))
if (y <= (-1.56d+220)) then
tmp = x * t_1
else if (y <= (-3.65d+22)) then
tmp = ((x * 2.0d0) / (y * (t * t))) / y
else if (y <= 7d-41) then
tmp = x / (1.0d0 + (b * (a + (0.5d0 * (b * (a * a))))))
else if (y <= 4.6d+144) then
tmp = (x * (t * (t * t))) * ((y * (y * y)) * (-0.16666666666666666d0))
else
tmp = x / 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 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -1.56e+220) {
tmp = x * t_1;
} else if (y <= -3.65e+22) {
tmp = ((x * 2.0) / (y * (t * t))) / y;
} else if (y <= 7e-41) {
tmp = x / (1.0 + (b * (a + (0.5 * (b * (a * a))))));
} else if (y <= 4.6e+144) {
tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666);
} else {
tmp = x / t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -1.56e+220: tmp = x * t_1 elif y <= -3.65e+22: tmp = ((x * 2.0) / (y * (t * t))) / y elif y <= 7e-41: tmp = x / (1.0 + (b * (a + (0.5 * (b * (a * a)))))) elif y <= 4.6e+144: tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666) else: tmp = x / t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -1.56e+220) tmp = Float64(x * t_1); elseif (y <= -3.65e+22) tmp = Float64(Float64(Float64(x * 2.0) / Float64(y * Float64(t * t))) / y); elseif (y <= 7e-41) tmp = Float64(x / Float64(1.0 + Float64(b * Float64(a + Float64(0.5 * Float64(b * Float64(a * a))))))); elseif (y <= 4.6e+144) tmp = Float64(Float64(x * Float64(t * Float64(t * t))) * Float64(Float64(y * Float64(y * y)) * -0.16666666666666666)); else tmp = Float64(x / t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -1.56e+220) tmp = x * t_1; elseif (y <= -3.65e+22) tmp = ((x * 2.0) / (y * (t * t))) / y; elseif (y <= 7e-41) tmp = x / (1.0 + (b * (a + (0.5 * (b * (a * a)))))); elseif (y <= 4.6e+144) tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666); else tmp = x / t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.56e+220], N[(x * t$95$1), $MachinePrecision], If[LessEqual[y, -3.65e+22], N[(N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, 7e-41], N[(x / N[(1.0 + N[(b * N[(a + N[(0.5 * N[(b * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.6e+144], N[(N[(x * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision], N[(x / t$95$1), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -1.56 \cdot 10^{+220}:\\
\;\;\;\;x \cdot t\_1\\
\mathbf{elif}\;y \leq -3.65 \cdot 10^{+22}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{y \cdot \left(t \cdot t\right)}}{y}\\
\mathbf{elif}\;y \leq 7 \cdot 10^{-41}:\\
\;\;\;\;\frac{x}{1 + b \cdot \left(a + 0.5 \cdot \left(b \cdot \left(a \cdot a\right)\right)\right)}\\
\mathbf{elif}\;y \leq 4.6 \cdot 10^{+144}:\\
\;\;\;\;\left(x \cdot \left(t \cdot \left(t \cdot t\right)\right)\right) \cdot \left(\left(y \cdot \left(y \cdot y\right)\right) \cdot -0.16666666666666666\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_1}\\
\end{array}
\end{array}
if y < -1.56e220Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -1.56e220 < y < -3.6499999999999999e22Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.5%
Simplified50.5%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.6%
Simplified65.6%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6468.5%
Applied egg-rr68.5%
if -3.6499999999999999e22 < y < 6.9999999999999999e-41Initial program 93.6%
Simplified93.6%
Taylor expanded in b around inf
*-lowering-*.f6475.0%
Simplified75.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.9%
Simplified56.9%
if 6.9999999999999999e-41 < y < 4.6000000000000003e144Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6451.4%
Simplified51.4%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified25.1%
Taylor expanded in y around inf
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.3%
Simplified52.3%
if 4.6000000000000003e144 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification62.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (* t 0.5) (* t (* y y)))))
(if (<= y -1.1e+220)
(* x t_1)
(if (<= y -3.65e+22)
(/ (/ (* x 2.0) (* y (* t t))) y)
(if (<= y 7e-41)
(/ x (+ 1.0 (* a (+ b (* a (* b (* b 0.5)))))))
(if (<= y 3.5e+144)
(* (* x (* t (* t t))) (* (* y (* y y)) -0.16666666666666666))
(/ x t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -1.1e+220) {
tmp = x * t_1;
} else if (y <= -3.65e+22) {
tmp = ((x * 2.0) / (y * (t * t))) / y;
} else if (y <= 7e-41) {
tmp = x / (1.0 + (a * (b + (a * (b * (b * 0.5))))));
} else if (y <= 3.5e+144) {
tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666);
} else {
tmp = x / 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 = (t * 0.5d0) * (t * (y * y))
if (y <= (-1.1d+220)) then
tmp = x * t_1
else if (y <= (-3.65d+22)) then
tmp = ((x * 2.0d0) / (y * (t * t))) / y
else if (y <= 7d-41) then
tmp = x / (1.0d0 + (a * (b + (a * (b * (b * 0.5d0))))))
else if (y <= 3.5d+144) then
tmp = (x * (t * (t * t))) * ((y * (y * y)) * (-0.16666666666666666d0))
else
tmp = x / 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 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -1.1e+220) {
tmp = x * t_1;
} else if (y <= -3.65e+22) {
tmp = ((x * 2.0) / (y * (t * t))) / y;
} else if (y <= 7e-41) {
tmp = x / (1.0 + (a * (b + (a * (b * (b * 0.5))))));
} else if (y <= 3.5e+144) {
tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666);
} else {
tmp = x / t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -1.1e+220: tmp = x * t_1 elif y <= -3.65e+22: tmp = ((x * 2.0) / (y * (t * t))) / y elif y <= 7e-41: tmp = x / (1.0 + (a * (b + (a * (b * (b * 0.5)))))) elif y <= 3.5e+144: tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666) else: tmp = x / t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -1.1e+220) tmp = Float64(x * t_1); elseif (y <= -3.65e+22) tmp = Float64(Float64(Float64(x * 2.0) / Float64(y * Float64(t * t))) / y); elseif (y <= 7e-41) tmp = Float64(x / Float64(1.0 + Float64(a * Float64(b + Float64(a * Float64(b * Float64(b * 0.5))))))); elseif (y <= 3.5e+144) tmp = Float64(Float64(x * Float64(t * Float64(t * t))) * Float64(Float64(y * Float64(y * y)) * -0.16666666666666666)); else tmp = Float64(x / t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -1.1e+220) tmp = x * t_1; elseif (y <= -3.65e+22) tmp = ((x * 2.0) / (y * (t * t))) / y; elseif (y <= 7e-41) tmp = x / (1.0 + (a * (b + (a * (b * (b * 0.5)))))); elseif (y <= 3.5e+144) tmp = (x * (t * (t * t))) * ((y * (y * y)) * -0.16666666666666666); else tmp = x / t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.1e+220], N[(x * t$95$1), $MachinePrecision], If[LessEqual[y, -3.65e+22], N[(N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, 7e-41], N[(x / N[(1.0 + N[(a * N[(b + N[(a * N[(b * N[(b * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.5e+144], N[(N[(x * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision], N[(x / t$95$1), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -1.1 \cdot 10^{+220}:\\
\;\;\;\;x \cdot t\_1\\
\mathbf{elif}\;y \leq -3.65 \cdot 10^{+22}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{y \cdot \left(t \cdot t\right)}}{y}\\
\mathbf{elif}\;y \leq 7 \cdot 10^{-41}:\\
\;\;\;\;\frac{x}{1 + a \cdot \left(b + a \cdot \left(b \cdot \left(b \cdot 0.5\right)\right)\right)}\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{+144}:\\
\;\;\;\;\left(x \cdot \left(t \cdot \left(t \cdot t\right)\right)\right) \cdot \left(\left(y \cdot \left(y \cdot y\right)\right) \cdot -0.16666666666666666\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_1}\\
\end{array}
\end{array}
if y < -1.09999999999999995e220Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -1.09999999999999995e220 < y < -3.6499999999999999e22Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.5%
Simplified50.5%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6465.6%
Simplified65.6%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6468.5%
Applied egg-rr68.5%
if -3.6499999999999999e22 < y < 6.9999999999999999e-41Initial program 93.6%
Simplified93.6%
Taylor expanded in b around inf
*-lowering-*.f6475.0%
Simplified75.0%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f6453.3%
Simplified53.3%
if 6.9999999999999999e-41 < y < 3.4999999999999998e144Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6451.4%
Simplified51.4%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified25.1%
Taylor expanded in y around inf
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.3%
Simplified52.3%
if 3.4999999999999998e144 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification60.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (* t 0.5) (* t (* y y)))))
(if (<= y -9.2e+219)
(* x t_1)
(if (<= y -5.2e-31)
(/ (/ (* x 2.0) (* y (* t t))) y)
(if (<= y 1.05e-119)
(/ x (+ 1.0 (* a b)))
(if (<= y 5.4e+144) (* y (* y (* x (* 0.5 (* t t))))) (/ x t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -5.2e-31) {
tmp = ((x * 2.0) / (y * (t * t))) / y;
} else if (y <= 1.05e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 5.4e+144) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / 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 = (t * 0.5d0) * (t * (y * y))
if (y <= (-9.2d+219)) then
tmp = x * t_1
else if (y <= (-5.2d-31)) then
tmp = ((x * 2.0d0) / (y * (t * t))) / y
else if (y <= 1.05d-119) then
tmp = x / (1.0d0 + (a * b))
else if (y <= 5.4d+144) then
tmp = y * (y * (x * (0.5d0 * (t * t))))
else
tmp = x / 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 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -5.2e-31) {
tmp = ((x * 2.0) / (y * (t * t))) / y;
} else if (y <= 1.05e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 5.4e+144) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -9.2e+219: tmp = x * t_1 elif y <= -5.2e-31: tmp = ((x * 2.0) / (y * (t * t))) / y elif y <= 1.05e-119: tmp = x / (1.0 + (a * b)) elif y <= 5.4e+144: tmp = y * (y * (x * (0.5 * (t * t)))) else: tmp = x / t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -9.2e+219) tmp = Float64(x * t_1); elseif (y <= -5.2e-31) tmp = Float64(Float64(Float64(x * 2.0) / Float64(y * Float64(t * t))) / y); elseif (y <= 1.05e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); elseif (y <= 5.4e+144) tmp = Float64(y * Float64(y * Float64(x * Float64(0.5 * Float64(t * t))))); else tmp = Float64(x / t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -9.2e+219) tmp = x * t_1; elseif (y <= -5.2e-31) tmp = ((x * 2.0) / (y * (t * t))) / y; elseif (y <= 1.05e-119) tmp = x / (1.0 + (a * b)); elseif (y <= 5.4e+144) tmp = y * (y * (x * (0.5 * (t * t)))); else tmp = x / t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9.2e+219], N[(x * t$95$1), $MachinePrecision], If[LessEqual[y, -5.2e-31], N[(N[(N[(x * 2.0), $MachinePrecision] / N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], If[LessEqual[y, 1.05e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5.4e+144], N[(y * N[(y * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / t$95$1), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -9.2 \cdot 10^{+219}:\\
\;\;\;\;x \cdot t\_1\\
\mathbf{elif}\;y \leq -5.2 \cdot 10^{-31}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{y \cdot \left(t \cdot t\right)}}{y}\\
\mathbf{elif}\;y \leq 1.05 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{elif}\;y \leq 5.4 \cdot 10^{+144}:\\
\;\;\;\;y \cdot \left(y \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_1}\\
\end{array}
\end{array}
if y < -9.2000000000000004e219Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -9.2000000000000004e219 < y < -5.19999999999999991e-31Initial program 98.2%
Simplified98.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.2%
Simplified50.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.7%
Simplified54.7%
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6456.7%
Applied egg-rr56.7%
if -5.19999999999999991e-31 < y < 1.05e-119Initial program 92.9%
Simplified92.9%
Taylor expanded in b around inf
*-lowering-*.f6480.4%
Simplified80.4%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.5%
Simplified54.5%
if 1.05e-119 < y < 5.4000000000000003e144Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.2%
Simplified57.2%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6432.5%
Simplified32.5%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr26.5%
Taylor expanded in t around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.9%
Simplified44.9%
if 5.4000000000000003e144 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification57.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (* t 0.5) (* t (* y y)))))
(if (<= y -9.2e+219)
(* x t_1)
(if (<= y -6.2e-31)
(/ (/ (* x 2.0) (* t t)) (* y y))
(if (<= y 2.4e-119)
(/ x (+ 1.0 (* a b)))
(if (<= y 2.9e+146) (* y (* y (* x (* 0.5 (* t t))))) (/ x t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -6.2e-31) {
tmp = ((x * 2.0) / (t * t)) / (y * y);
} else if (y <= 2.4e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 2.9e+146) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / 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 = (t * 0.5d0) * (t * (y * y))
if (y <= (-9.2d+219)) then
tmp = x * t_1
else if (y <= (-6.2d-31)) then
tmp = ((x * 2.0d0) / (t * t)) / (y * y)
else if (y <= 2.4d-119) then
tmp = x / (1.0d0 + (a * b))
else if (y <= 2.9d+146) then
tmp = y * (y * (x * (0.5d0 * (t * t))))
else
tmp = x / 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 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -6.2e-31) {
tmp = ((x * 2.0) / (t * t)) / (y * y);
} else if (y <= 2.4e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 2.9e+146) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -9.2e+219: tmp = x * t_1 elif y <= -6.2e-31: tmp = ((x * 2.0) / (t * t)) / (y * y) elif y <= 2.4e-119: tmp = x / (1.0 + (a * b)) elif y <= 2.9e+146: tmp = y * (y * (x * (0.5 * (t * t)))) else: tmp = x / t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -9.2e+219) tmp = Float64(x * t_1); elseif (y <= -6.2e-31) tmp = Float64(Float64(Float64(x * 2.0) / Float64(t * t)) / Float64(y * y)); elseif (y <= 2.4e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); elseif (y <= 2.9e+146) tmp = Float64(y * Float64(y * Float64(x * Float64(0.5 * Float64(t * t))))); else tmp = Float64(x / t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -9.2e+219) tmp = x * t_1; elseif (y <= -6.2e-31) tmp = ((x * 2.0) / (t * t)) / (y * y); elseif (y <= 2.4e-119) tmp = x / (1.0 + (a * b)); elseif (y <= 2.9e+146) tmp = y * (y * (x * (0.5 * (t * t)))); else tmp = x / t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9.2e+219], N[(x * t$95$1), $MachinePrecision], If[LessEqual[y, -6.2e-31], N[(N[(N[(x * 2.0), $MachinePrecision] / N[(t * t), $MachinePrecision]), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.9e+146], N[(y * N[(y * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / t$95$1), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -9.2 \cdot 10^{+219}:\\
\;\;\;\;x \cdot t\_1\\
\mathbf{elif}\;y \leq -6.2 \cdot 10^{-31}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{t \cdot t}}{y \cdot y}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{elif}\;y \leq 2.9 \cdot 10^{+146}:\\
\;\;\;\;y \cdot \left(y \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_1}\\
\end{array}
\end{array}
if y < -9.2000000000000004e219Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -9.2000000000000004e219 < y < -6.19999999999999999e-31Initial program 98.2%
Simplified98.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.2%
Simplified50.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.7%
Simplified54.7%
if -6.19999999999999999e-31 < y < 2.40000000000000009e-119Initial program 92.9%
Simplified92.9%
Taylor expanded in b around inf
*-lowering-*.f6480.4%
Simplified80.4%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.5%
Simplified54.5%
if 2.40000000000000009e-119 < y < 2.8999999999999998e146Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.2%
Simplified57.2%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6432.5%
Simplified32.5%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr26.5%
Taylor expanded in t around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.9%
Simplified44.9%
if 2.8999999999999998e146 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification57.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (* t 0.5) (* t (* y y)))))
(if (<= y -9.2e+219)
(* x t_1)
(if (<= y -1.02e-30)
(/ (* x (/ 2.0 (* t t))) (* y y))
(if (<= y 2.35e-119)
(/ x (+ 1.0 (* a b)))
(if (<= y 3.5e+146) (* y (* y (* x (* 0.5 (* t t))))) (/ x t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -1.02e-30) {
tmp = (x * (2.0 / (t * t))) / (y * y);
} else if (y <= 2.35e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 3.5e+146) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / 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 = (t * 0.5d0) * (t * (y * y))
if (y <= (-9.2d+219)) then
tmp = x * t_1
else if (y <= (-1.02d-30)) then
tmp = (x * (2.0d0 / (t * t))) / (y * y)
else if (y <= 2.35d-119) then
tmp = x / (1.0d0 + (a * b))
else if (y <= 3.5d+146) then
tmp = y * (y * (x * (0.5d0 * (t * t))))
else
tmp = x / 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 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -1.02e-30) {
tmp = (x * (2.0 / (t * t))) / (y * y);
} else if (y <= 2.35e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 3.5e+146) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -9.2e+219: tmp = x * t_1 elif y <= -1.02e-30: tmp = (x * (2.0 / (t * t))) / (y * y) elif y <= 2.35e-119: tmp = x / (1.0 + (a * b)) elif y <= 3.5e+146: tmp = y * (y * (x * (0.5 * (t * t)))) else: tmp = x / t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -9.2e+219) tmp = Float64(x * t_1); elseif (y <= -1.02e-30) tmp = Float64(Float64(x * Float64(2.0 / Float64(t * t))) / Float64(y * y)); elseif (y <= 2.35e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); elseif (y <= 3.5e+146) tmp = Float64(y * Float64(y * Float64(x * Float64(0.5 * Float64(t * t))))); else tmp = Float64(x / t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -9.2e+219) tmp = x * t_1; elseif (y <= -1.02e-30) tmp = (x * (2.0 / (t * t))) / (y * y); elseif (y <= 2.35e-119) tmp = x / (1.0 + (a * b)); elseif (y <= 3.5e+146) tmp = y * (y * (x * (0.5 * (t * t)))); else tmp = x / t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9.2e+219], N[(x * t$95$1), $MachinePrecision], If[LessEqual[y, -1.02e-30], N[(N[(x * N[(2.0 / N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(y * y), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.35e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.5e+146], N[(y * N[(y * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / t$95$1), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -9.2 \cdot 10^{+219}:\\
\;\;\;\;x \cdot t\_1\\
\mathbf{elif}\;y \leq -1.02 \cdot 10^{-30}:\\
\;\;\;\;\frac{x \cdot \frac{2}{t \cdot t}}{y \cdot y}\\
\mathbf{elif}\;y \leq 2.35 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{+146}:\\
\;\;\;\;y \cdot \left(y \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_1}\\
\end{array}
\end{array}
if y < -9.2000000000000004e219Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -9.2000000000000004e219 < y < -1.0199999999999999e-30Initial program 98.2%
Simplified98.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.2%
Simplified50.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.7%
Simplified54.7%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6454.7%
Applied egg-rr54.7%
if -1.0199999999999999e-30 < y < 2.35000000000000001e-119Initial program 92.9%
Simplified92.9%
Taylor expanded in b around inf
*-lowering-*.f6480.4%
Simplified80.4%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.5%
Simplified54.5%
if 2.35000000000000001e-119 < y < 3.5000000000000001e146Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.2%
Simplified57.2%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6432.5%
Simplified32.5%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr26.5%
Taylor expanded in t around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.9%
Simplified44.9%
if 3.5000000000000001e146 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification57.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (* t 0.5) (* t (* y y)))))
(if (<= y -9.2e+219)
(* x t_1)
(if (<= y -2.05e-30)
(* (/ 2.0 (* t t)) (/ x (* y y)))
(if (<= y 2.05e-119)
(/ x (+ 1.0 (* a b)))
(if (<= y 6e+145) (* y (* y (* x (* 0.5 (* t t))))) (/ x t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -2.05e-30) {
tmp = (2.0 / (t * t)) * (x / (y * y));
} else if (y <= 2.05e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 6e+145) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / 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 = (t * 0.5d0) * (t * (y * y))
if (y <= (-9.2d+219)) then
tmp = x * t_1
else if (y <= (-2.05d-30)) then
tmp = (2.0d0 / (t * t)) * (x / (y * y))
else if (y <= 2.05d-119) then
tmp = x / (1.0d0 + (a * b))
else if (y <= 6d+145) then
tmp = y * (y * (x * (0.5d0 * (t * t))))
else
tmp = x / 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 = (t * 0.5) * (t * (y * y));
double tmp;
if (y <= -9.2e+219) {
tmp = x * t_1;
} else if (y <= -2.05e-30) {
tmp = (2.0 / (t * t)) * (x / (y * y));
} else if (y <= 2.05e-119) {
tmp = x / (1.0 + (a * b));
} else if (y <= 6e+145) {
tmp = y * (y * (x * (0.5 * (t * t))));
} else {
tmp = x / t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (t * 0.5) * (t * (y * y)) tmp = 0 if y <= -9.2e+219: tmp = x * t_1 elif y <= -2.05e-30: tmp = (2.0 / (t * t)) * (x / (y * y)) elif y <= 2.05e-119: tmp = x / (1.0 + (a * b)) elif y <= 6e+145: tmp = y * (y * (x * (0.5 * (t * t)))) else: tmp = x / t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(t * 0.5) * Float64(t * Float64(y * y))) tmp = 0.0 if (y <= -9.2e+219) tmp = Float64(x * t_1); elseif (y <= -2.05e-30) tmp = Float64(Float64(2.0 / Float64(t * t)) * Float64(x / Float64(y * y))); elseif (y <= 2.05e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); elseif (y <= 6e+145) tmp = Float64(y * Float64(y * Float64(x * Float64(0.5 * Float64(t * t))))); else tmp = Float64(x / t_1); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (t * 0.5) * (t * (y * y)); tmp = 0.0; if (y <= -9.2e+219) tmp = x * t_1; elseif (y <= -2.05e-30) tmp = (2.0 / (t * t)) * (x / (y * y)); elseif (y <= 2.05e-119) tmp = x / (1.0 + (a * b)); elseif (y <= 6e+145) tmp = y * (y * (x * (0.5 * (t * t)))); else tmp = x / t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9.2e+219], N[(x * t$95$1), $MachinePrecision], If[LessEqual[y, -2.05e-30], N[(N[(2.0 / N[(t * t), $MachinePrecision]), $MachinePrecision] * N[(x / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.05e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6e+145], N[(y * N[(y * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / t$95$1), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\\
\mathbf{if}\;y \leq -9.2 \cdot 10^{+219}:\\
\;\;\;\;x \cdot t\_1\\
\mathbf{elif}\;y \leq -2.05 \cdot 10^{-30}:\\
\;\;\;\;\frac{2}{t \cdot t} \cdot \frac{x}{y \cdot y}\\
\mathbf{elif}\;y \leq 2.05 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{elif}\;y \leq 6 \cdot 10^{+145}:\\
\;\;\;\;y \cdot \left(y \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t\_1}\\
\end{array}
\end{array}
if y < -9.2000000000000004e219Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -9.2000000000000004e219 < y < -2.0500000000000002e-30Initial program 98.2%
Simplified98.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.2%
Simplified50.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.7%
Simplified54.7%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6449.1%
Applied egg-rr49.1%
if -2.0500000000000002e-30 < y < 2.0500000000000001e-119Initial program 92.9%
Simplified92.9%
Taylor expanded in b around inf
*-lowering-*.f6480.4%
Simplified80.4%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.5%
Simplified54.5%
if 2.0500000000000001e-119 < y < 6.0000000000000005e145Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.2%
Simplified57.2%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6432.5%
Simplified32.5%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr26.5%
Taylor expanded in t around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.9%
Simplified44.9%
if 6.0000000000000005e145 < y Initial program 96.3%
Simplified96.3%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in y around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6485.4%
Simplified85.4%
Final simplification56.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -9.2e+219)
(* x (* (* t 0.5) (* t (* y y))))
(if (<= y -3.8e-31)
(* (/ 2.0 (* t t)) (/ x (* y y)))
(if (<= y 2.3e-119)
(/ x (+ 1.0 (* a b)))
(* y (* y (* x (* 0.5 (* t t)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -9.2e+219) {
tmp = x * ((t * 0.5) * (t * (y * y)));
} else if (y <= -3.8e-31) {
tmp = (2.0 / (t * t)) * (x / (y * y));
} else if (y <= 2.3e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = y * (y * (x * (0.5 * (t * t))));
}
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 <= (-9.2d+219)) then
tmp = x * ((t * 0.5d0) * (t * (y * y)))
else if (y <= (-3.8d-31)) then
tmp = (2.0d0 / (t * t)) * (x / (y * y))
else if (y <= 2.3d-119) then
tmp = x / (1.0d0 + (a * b))
else
tmp = y * (y * (x * (0.5d0 * (t * t))))
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 <= -9.2e+219) {
tmp = x * ((t * 0.5) * (t * (y * y)));
} else if (y <= -3.8e-31) {
tmp = (2.0 / (t * t)) * (x / (y * y));
} else if (y <= 2.3e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = y * (y * (x * (0.5 * (t * t))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -9.2e+219: tmp = x * ((t * 0.5) * (t * (y * y))) elif y <= -3.8e-31: tmp = (2.0 / (t * t)) * (x / (y * y)) elif y <= 2.3e-119: tmp = x / (1.0 + (a * b)) else: tmp = y * (y * (x * (0.5 * (t * t)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -9.2e+219) tmp = Float64(x * Float64(Float64(t * 0.5) * Float64(t * Float64(y * y)))); elseif (y <= -3.8e-31) tmp = Float64(Float64(2.0 / Float64(t * t)) * Float64(x / Float64(y * y))); elseif (y <= 2.3e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(y * Float64(y * Float64(x * Float64(0.5 * Float64(t * t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -9.2e+219) tmp = x * ((t * 0.5) * (t * (y * y))); elseif (y <= -3.8e-31) tmp = (2.0 / (t * t)) * (x / (y * y)); elseif (y <= 2.3e-119) tmp = x / (1.0 + (a * b)); else tmp = y * (y * (x * (0.5 * (t * t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -9.2e+219], N[(x * N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -3.8e-31], N[(N[(2.0 / N[(t * t), $MachinePrecision]), $MachinePrecision] * N[(x / N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.3e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y * N[(y * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.2 \cdot 10^{+219}:\\
\;\;\;\;x \cdot \left(\left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\right)\\
\mathbf{elif}\;y \leq -3.8 \cdot 10^{-31}:\\
\;\;\;\;\frac{2}{t \cdot t} \cdot \frac{x}{y \cdot y}\\
\mathbf{elif}\;y \leq 2.3 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(y \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -9.2000000000000004e219Initial program 93.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.3%
Simplified57.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6463.5%
Simplified63.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
if -9.2000000000000004e219 < y < -3.8e-31Initial program 98.2%
Simplified98.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6450.2%
Simplified50.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.2%
Simplified32.2%
Taylor expanded in y around inf
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.7%
Simplified54.7%
associate-/l/N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6449.1%
Applied egg-rr49.1%
if -3.8e-31 < y < 2.29999999999999993e-119Initial program 92.9%
Simplified92.9%
Taylor expanded in b around inf
*-lowering-*.f6480.4%
Simplified80.4%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.5%
Simplified54.5%
if 2.29999999999999993e-119 < y Initial program 98.7%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6459.5%
Simplified59.5%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6427.1%
Simplified27.1%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr17.5%
Taylor expanded in t around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.8%
Simplified48.8%
Final simplification53.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -9.5e+131)
(/ x (+ 1.0 (* y t)))
(if (<= y 7.5e+14)
(- x (* a (* x b)))
(if (<= y 4.1e+191) (* x (* t (- 0.0 y))) (* y (- (/ x y) (* x t)))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -9.5e+131) {
tmp = x / (1.0 + (y * t));
} else if (y <= 7.5e+14) {
tmp = x - (a * (x * b));
} else if (y <= 4.1e+191) {
tmp = x * (t * (0.0 - y));
} else {
tmp = y * ((x / y) - (x * t));
}
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 <= (-9.5d+131)) then
tmp = x / (1.0d0 + (y * t))
else if (y <= 7.5d+14) then
tmp = x - (a * (x * b))
else if (y <= 4.1d+191) then
tmp = x * (t * (0.0d0 - y))
else
tmp = y * ((x / y) - (x * t))
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 <= -9.5e+131) {
tmp = x / (1.0 + (y * t));
} else if (y <= 7.5e+14) {
tmp = x - (a * (x * b));
} else if (y <= 4.1e+191) {
tmp = x * (t * (0.0 - y));
} else {
tmp = y * ((x / y) - (x * t));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -9.5e+131: tmp = x / (1.0 + (y * t)) elif y <= 7.5e+14: tmp = x - (a * (x * b)) elif y <= 4.1e+191: tmp = x * (t * (0.0 - y)) else: tmp = y * ((x / y) - (x * t)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -9.5e+131) tmp = Float64(x / Float64(1.0 + Float64(y * t))); elseif (y <= 7.5e+14) tmp = Float64(x - Float64(a * Float64(x * b))); elseif (y <= 4.1e+191) tmp = Float64(x * Float64(t * Float64(0.0 - y))); else tmp = Float64(y * Float64(Float64(x / y) - Float64(x * t))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -9.5e+131) tmp = x / (1.0 + (y * t)); elseif (y <= 7.5e+14) tmp = x - (a * (x * b)); elseif (y <= 4.1e+191) tmp = x * (t * (0.0 - y)); else tmp = y * ((x / y) - (x * t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -9.5e+131], N[(x / N[(1.0 + N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 7.5e+14], N[(x - N[(a * N[(x * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.1e+191], N[(x * N[(t * N[(0.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y * N[(N[(x / y), $MachinePrecision] - N[(x * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.5 \cdot 10^{+131}:\\
\;\;\;\;\frac{x}{1 + y \cdot t}\\
\mathbf{elif}\;y \leq 7.5 \cdot 10^{+14}:\\
\;\;\;\;x - a \cdot \left(x \cdot b\right)\\
\mathbf{elif}\;y \leq 4.1 \cdot 10^{+191}:\\
\;\;\;\;x \cdot \left(t \cdot \left(0 - y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(\frac{x}{y} - x \cdot t\right)\\
\end{array}
\end{array}
if y < -9.50000000000000015e131Initial program 96.6%
Simplified96.6%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6466.5%
Simplified66.5%
Taylor expanded in y around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6432.7%
Simplified32.7%
if -9.50000000000000015e131 < y < 7.5e14Initial program 95.2%
Simplified95.2%
Taylor expanded in b around inf
*-lowering-*.f6469.2%
Simplified69.2%
Taylor expanded in a around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6443.4%
Simplified43.4%
if 7.5e14 < y < 4.0999999999999999e191Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6457.1%
Simplified57.1%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6421.2%
Simplified21.2%
Taylor expanded in y around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6430.9%
Simplified30.9%
if 4.0999999999999999e191 < y Initial program 94.7%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6454.0%
Simplified54.0%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6413.5%
Simplified13.5%
Taylor expanded in y around inf
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6443.4%
Simplified43.4%
Final simplification40.8%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -1.95e-22)
(* x (* (* t 0.5) (* t (* y y))))
(if (<= y 2.1e-119)
(/ x (+ 1.0 (* a b)))
(* y (* y (* x (* 0.5 (* t t))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.95e-22) {
tmp = x * ((t * 0.5) * (t * (y * y)));
} else if (y <= 2.1e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = y * (y * (x * (0.5 * (t * t))));
}
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.95d-22)) then
tmp = x * ((t * 0.5d0) * (t * (y * y)))
else if (y <= 2.1d-119) then
tmp = x / (1.0d0 + (a * b))
else
tmp = y * (y * (x * (0.5d0 * (t * t))))
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.95e-22) {
tmp = x * ((t * 0.5) * (t * (y * y)));
} else if (y <= 2.1e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = y * (y * (x * (0.5 * (t * t))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -1.95e-22: tmp = x * ((t * 0.5) * (t * (y * y))) elif y <= 2.1e-119: tmp = x / (1.0 + (a * b)) else: tmp = y * (y * (x * (0.5 * (t * t)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -1.95e-22) tmp = Float64(x * Float64(Float64(t * 0.5) * Float64(t * Float64(y * y)))); elseif (y <= 2.1e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(y * Float64(y * Float64(x * Float64(0.5 * Float64(t * t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -1.95e-22) tmp = x * ((t * 0.5) * (t * (y * y))); elseif (y <= 2.1e-119) tmp = x / (1.0 + (a * b)); else tmp = y * (y * (x * (0.5 * (t * t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -1.95e-22], N[(x * N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.1e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y * N[(y * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.95 \cdot 10^{-22}:\\
\;\;\;\;x \cdot \left(\left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\right)\\
\mathbf{elif}\;y \leq 2.1 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(y \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -1.94999999999999999e-22Initial program 97.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6449.5%
Simplified49.5%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6431.5%
Simplified31.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6442.8%
Simplified42.8%
if -1.94999999999999999e-22 < y < 2.1e-119Initial program 93.2%
Simplified93.2%
Taylor expanded in b around inf
*-lowering-*.f6479.5%
Simplified79.5%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.0%
Simplified54.0%
if 2.1e-119 < y Initial program 98.7%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6459.5%
Simplified59.5%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6427.1%
Simplified27.1%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr17.5%
Taylor expanded in t around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.8%
Simplified48.8%
Final simplification49.6%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -1.9e-22)
(* x (* (* t 0.5) (* t (* y y))))
(if (<= y 1.55e-119)
(/ x (+ 1.0 (* a b)))
(* x (* y (* y (* 0.5 (* t t))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.9e-22) {
tmp = x * ((t * 0.5) * (t * (y * y)));
} else if (y <= 1.55e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = x * (y * (y * (0.5 * (t * t))));
}
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.9d-22)) then
tmp = x * ((t * 0.5d0) * (t * (y * y)))
else if (y <= 1.55d-119) then
tmp = x / (1.0d0 + (a * b))
else
tmp = x * (y * (y * (0.5d0 * (t * t))))
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.9e-22) {
tmp = x * ((t * 0.5) * (t * (y * y)));
} else if (y <= 1.55e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = x * (y * (y * (0.5 * (t * t))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -1.9e-22: tmp = x * ((t * 0.5) * (t * (y * y))) elif y <= 1.55e-119: tmp = x / (1.0 + (a * b)) else: tmp = x * (y * (y * (0.5 * (t * t)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -1.9e-22) tmp = Float64(x * Float64(Float64(t * 0.5) * Float64(t * Float64(y * y)))); elseif (y <= 1.55e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(x * Float64(y * Float64(y * Float64(0.5 * Float64(t * t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -1.9e-22) tmp = x * ((t * 0.5) * (t * (y * y))); elseif (y <= 1.55e-119) tmp = x / (1.0 + (a * b)); else tmp = x * (y * (y * (0.5 * (t * t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -1.9e-22], N[(x * N[(N[(t * 0.5), $MachinePrecision] * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.55e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(y * N[(y * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.9 \cdot 10^{-22}:\\
\;\;\;\;x \cdot \left(\left(t \cdot 0.5\right) \cdot \left(t \cdot \left(y \cdot y\right)\right)\right)\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(y \cdot \left(y \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -1.90000000000000012e-22Initial program 97.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6449.5%
Simplified49.5%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6431.5%
Simplified31.5%
Taylor expanded in t around inf
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6442.8%
Simplified42.8%
if -1.90000000000000012e-22 < y < 1.54999999999999989e-119Initial program 93.2%
Simplified93.2%
Taylor expanded in b around inf
*-lowering-*.f6479.5%
Simplified79.5%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.0%
Simplified54.0%
if 1.54999999999999989e-119 < y Initial program 98.7%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6459.5%
Simplified59.5%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6427.1%
Simplified27.1%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr17.5%
Taylor expanded in t around inf
unpow2N/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.6%
Simplified42.6%
Final simplification47.7%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* x (* y (* y (* 0.5 (* t t))))))) (if (<= y -1.66e-21) t_1 (if (<= y 1.08e-119) (/ x (+ 1.0 (* a b))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (y * (y * (0.5 * (t * t))));
double tmp;
if (y <= -1.66e-21) {
tmp = t_1;
} else if (y <= 1.08e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: tmp
t_1 = x * (y * (y * (0.5d0 * (t * t))))
if (y <= (-1.66d-21)) then
tmp = t_1
else if (y <= 1.08d-119) then
tmp = x / (1.0d0 + (a * b))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * (y * (y * (0.5 * (t * t))));
double tmp;
if (y <= -1.66e-21) {
tmp = t_1;
} else if (y <= 1.08e-119) {
tmp = x / (1.0 + (a * b));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * (y * (y * (0.5 * (t * t)))) tmp = 0 if y <= -1.66e-21: tmp = t_1 elif y <= 1.08e-119: tmp = x / (1.0 + (a * b)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(y * Float64(y * Float64(0.5 * Float64(t * t))))) tmp = 0.0 if (y <= -1.66e-21) tmp = t_1; elseif (y <= 1.08e-119) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * (y * (y * (0.5 * (t * t)))); tmp = 0.0; if (y <= -1.66e-21) tmp = t_1; elseif (y <= 1.08e-119) tmp = x / (1.0 + (a * b)); 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[(y * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.66e-21], t$95$1, If[LessEqual[y, 1.08e-119], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(y \cdot \left(y \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\right)\\
\mathbf{if}\;y \leq -1.66 \cdot 10^{-21}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.08 \cdot 10^{-119}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.65999999999999993e-21 or 1.0799999999999999e-119 < y Initial program 97.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6455.0%
Simplified55.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
--lowering--.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6429.1%
Simplified29.1%
*-commutativeN/A
flip--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr16.5%
Taylor expanded in t around inf
unpow2N/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.5%
Simplified38.5%
if -1.65999999999999993e-21 < y < 1.0799999999999999e-119Initial program 93.2%
Simplified93.2%
Taylor expanded in b around inf
*-lowering-*.f6479.5%
Simplified79.5%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6454.0%
Simplified54.0%
Final simplification45.3%
(FPCore (x y z t a b) :precision binary64 (if (<= y -4.2e+132) (* x (- 1.0 (* y t))) (if (<= y 8800000000000.0) (- x (* a (* x b))) (* x (* t (- 0.0 y))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -4.2e+132) {
tmp = x * (1.0 - (y * t));
} else if (y <= 8800000000000.0) {
tmp = x - (a * (x * b));
} else {
tmp = x * (t * (0.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 (y <= (-4.2d+132)) then
tmp = x * (1.0d0 - (y * t))
else if (y <= 8800000000000.0d0) then
tmp = x - (a * (x * b))
else
tmp = x * (t * (0.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 (y <= -4.2e+132) {
tmp = x * (1.0 - (y * t));
} else if (y <= 8800000000000.0) {
tmp = x - (a * (x * b));
} else {
tmp = x * (t * (0.0 - y));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -4.2e+132: tmp = x * (1.0 - (y * t)) elif y <= 8800000000000.0: tmp = x - (a * (x * b)) else: tmp = x * (t * (0.0 - y)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -4.2e+132) tmp = Float64(x * Float64(1.0 - Float64(y * t))); elseif (y <= 8800000000000.0) tmp = Float64(x - Float64(a * Float64(x * b))); else tmp = Float64(x * Float64(t * Float64(0.0 - y))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -4.2e+132) tmp = x * (1.0 - (y * t)); elseif (y <= 8800000000000.0) tmp = x - (a * (x * b)); else tmp = x * (t * (0.0 - y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -4.2e+132], N[(x * N[(1.0 - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8800000000000.0], N[(x - N[(a * N[(x * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(t * N[(0.0 - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.2 \cdot 10^{+132}:\\
\;\;\;\;x \cdot \left(1 - y \cdot t\right)\\
\mathbf{elif}\;y \leq 8800000000000:\\
\;\;\;\;x - a \cdot \left(x \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(t \cdot \left(0 - y\right)\right)\\
\end{array}
\end{array}
if y < -4.19999999999999987e132Initial program 96.4%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6465.3%
Simplified65.3%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6430.1%
Simplified30.1%
if -4.19999999999999987e132 < y < 8.8e12Initial program 95.2%
Simplified95.2%
Taylor expanded in b around inf
*-lowering-*.f6468.8%
Simplified68.8%
Taylor expanded in a around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6443.1%
Simplified43.1%
if 8.8e12 < y Initial program 97.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6455.8%
Simplified55.8%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6418.0%
Simplified18.0%
Taylor expanded in y around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6423.7%
Simplified23.7%
Final simplification38.2%
(FPCore (x y z t a b) :precision binary64 (if (<= x 5.2e+169) (/ x (+ 1.0 (* y t))) (- x (* a (* x b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= 5.2e+169) {
tmp = x / (1.0 + (y * t));
} else {
tmp = x - (a * (x * 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 (x <= 5.2d+169) then
tmp = x / (1.0d0 + (y * t))
else
tmp = x - (a * (x * 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 (x <= 5.2e+169) {
tmp = x / (1.0 + (y * t));
} else {
tmp = x - (a * (x * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if x <= 5.2e+169: tmp = x / (1.0 + (y * t)) else: tmp = x - (a * (x * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (x <= 5.2e+169) tmp = Float64(x / Float64(1.0 + Float64(y * t))); else tmp = Float64(x - Float64(a * Float64(x * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (x <= 5.2e+169) tmp = x / (1.0 + (y * t)); else tmp = x - (a * (x * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, 5.2e+169], N[(x / N[(1.0 + N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(a * N[(x * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.2 \cdot 10^{+169}:\\
\;\;\;\;\frac{x}{1 + y \cdot t}\\
\mathbf{else}:\\
\;\;\;\;x - a \cdot \left(x \cdot b\right)\\
\end{array}
\end{array}
if x < 5.19999999999999999e169Initial program 96.2%
Simplified96.2%
Taylor expanded in t around inf
*-commutativeN/A
*-lowering-*.f6456.9%
Simplified56.9%
Taylor expanded in y around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6430.8%
Simplified30.8%
if 5.19999999999999999e169 < x Initial program 93.2%
Simplified93.2%
Taylor expanded in b around inf
*-lowering-*.f6473.0%
Simplified73.0%
Taylor expanded in a around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6452.6%
Simplified52.6%
Final simplification33.3%
(FPCore (x y z t a b) :precision binary64 (if (<= x 2.7e-15) (/ x (+ 1.0 (* a b))) (- x (* a (* x b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (x <= 2.7e-15) {
tmp = x / (1.0 + (a * b));
} else {
tmp = x - (a * (x * 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 (x <= 2.7d-15) then
tmp = x / (1.0d0 + (a * b))
else
tmp = x - (a * (x * 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 (x <= 2.7e-15) {
tmp = x / (1.0 + (a * b));
} else {
tmp = x - (a * (x * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if x <= 2.7e-15: tmp = x / (1.0 + (a * b)) else: tmp = x - (a * (x * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (x <= 2.7e-15) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(x - Float64(a * Float64(x * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (x <= 2.7e-15) tmp = x / (1.0 + (a * b)); else tmp = x - (a * (x * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[x, 2.7e-15], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(a * N[(x * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 2.7 \cdot 10^{-15}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;x - a \cdot \left(x \cdot b\right)\\
\end{array}
\end{array}
if x < 2.70000000000000009e-15Initial program 97.0%
Simplified97.0%
Taylor expanded in b around inf
*-lowering-*.f6455.9%
Simplified55.9%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6430.4%
Simplified30.4%
if 2.70000000000000009e-15 < x Initial program 92.3%
Simplified92.3%
Taylor expanded in b around inf
*-lowering-*.f6463.4%
Simplified63.4%
Taylor expanded in a around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6434.5%
Simplified34.5%
Final simplification31.4%
(FPCore (x y z t a b) :precision binary64 (* x (- 1.0 (* y t))))
double code(double x, double y, double z, double t, double a, double b) {
return x * (1.0 - (y * t));
}
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 * (1.0d0 - (y * t))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x * (1.0 - (y * t));
}
def code(x, y, z, t, a, b): return x * (1.0 - (y * t))
function code(x, y, z, t, a, b) return Float64(x * Float64(1.0 - Float64(y * t))) end
function tmp = code(x, y, z, t, a, b) tmp = x * (1.0 - (y * t)); end
code[x_, y_, z_, t_, a_, b_] := N[(x * N[(1.0 - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 - y \cdot t\right)
\end{array}
Initial program 95.8%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6458.0%
Simplified58.0%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6430.7%
Simplified30.7%
(FPCore (x y z t a b) :precision binary64 x)
double code(double x, double y, double z, double t, double a, double b) {
return x;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
code = x
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x;
}
def code(x, y, z, t, a, b): return x
function code(x, y, z, t, a, b) return x end
function tmp = code(x, y, z, t, a, b) tmp = x; end
code[x_, y_, z_, t_, a_, b_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 95.8%
Simplified95.8%
Taylor expanded in b around inf
*-lowering-*.f6457.8%
Simplified57.8%
Taylor expanded in a around 0
Simplified23.9%
herbie shell --seed 2024161
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaApprox from math-functions-0.1.5.2, B"
:precision binary64
(* x (exp (+ (* y (- (log z) t)) (* a (- (log (- 1.0 z)) b))))))