
(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 22 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 (- (* y (- (log z) t)) (* a (+ z b))))))
double code(double x, double y, double z, double t, double a, double b) {
return x * exp(((y * (log(z) - t)) - (a * (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 * (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 * (z + b))));
}
def code(x, y, z, t, a, b): return x * math.exp(((y * (math.log(z) - t)) - (a * (z + b))))
function code(x, y, z, t, a, b) return Float64(x * exp(Float64(Float64(y * Float64(log(z) - t)) - Float64(a * Float64(z + b))))) end
function tmp = code(x, y, z, t, a, b) tmp = x * exp(((y * (log(z) - t)) - (a * (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[(z + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot e^{y \cdot \left(\log z - t\right) - a \cdot \left(z + b\right)}
\end{array}
Initial program 96.0%
Taylor expanded in z around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f64N/A
associate-*r*N/A
associate-*r*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
+-lowering-+.f6498.8%
Simplified98.8%
Final simplification98.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (exp (* y (- (log z) t))))))
(if (<= y -2.3e-80)
t_1
(if (<= y 0.00045) (* x (exp (- 0.0 (* a b)))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * exp((y * (log(z) - t)));
double tmp;
if (y <= -2.3e-80) {
tmp = t_1;
} else if (y <= 0.00045) {
tmp = x * exp((0.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 * exp((y * (log(z) - t)))
if (y <= (-2.3d-80)) then
tmp = t_1
else if (y <= 0.00045d0) then
tmp = x * exp((0.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 * Math.exp((y * (Math.log(z) - t)));
double tmp;
if (y <= -2.3e-80) {
tmp = t_1;
} else if (y <= 0.00045) {
tmp = x * Math.exp((0.0 - (a * b)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * math.exp((y * (math.log(z) - t))) tmp = 0 if y <= -2.3e-80: tmp = t_1 elif y <= 0.00045: tmp = x * math.exp((0.0 - (a * b))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * exp(Float64(y * Float64(log(z) - t)))) tmp = 0.0 if (y <= -2.3e-80) tmp = t_1; elseif (y <= 0.00045) tmp = Float64(x * exp(Float64(0.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 * exp((y * (log(z) - t))); tmp = 0.0; if (y <= -2.3e-80) tmp = t_1; elseif (y <= 0.00045) tmp = x * exp((0.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[Exp[N[(y * N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.3e-80], t$95$1, If[LessEqual[y, 0.00045], N[(x * N[Exp[N[(0.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot e^{y \cdot \left(\log z - t\right)}\\
\mathbf{if}\;y \leq -2.3 \cdot 10^{-80}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 0.00045:\\
\;\;\;\;x \cdot e^{0 - a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.2999999999999998e-80 or 4.4999999999999999e-4 < y Initial program 96.6%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6485.8%
Simplified85.8%
if -2.2999999999999998e-80 < y < 4.4999999999999999e-4Initial program 95.3%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6481.8%
Simplified81.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (/ x (exp (* y t)))))
(if (<= t -440.0)
t_1
(if (<= t 8e-216)
(* x (pow z y))
(if (<= t 2e+44) (/ x (exp (* a b))) 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 <= -440.0) {
tmp = t_1;
} else if (t <= 8e-216) {
tmp = x * pow(z, y);
} else if (t <= 2e+44) {
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 / exp((y * t))
if (t <= (-440.0d0)) then
tmp = t_1
else if (t <= 8d-216) then
tmp = x * (z ** y)
else if (t <= 2d+44) 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.exp((y * t));
double tmp;
if (t <= -440.0) {
tmp = t_1;
} else if (t <= 8e-216) {
tmp = x * Math.pow(z, y);
} else if (t <= 2e+44) {
tmp = x / Math.exp((a * b));
} 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 <= -440.0: tmp = t_1 elif t <= 8e-216: tmp = x * math.pow(z, y) elif t <= 2e+44: tmp = x / math.exp((a * b)) 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 <= -440.0) tmp = t_1; elseif (t <= 8e-216) tmp = Float64(x * (z ^ y)); elseif (t <= 2e+44) 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 / exp((y * t)); tmp = 0.0; if (t <= -440.0) tmp = t_1; elseif (t <= 8e-216) tmp = x * (z ^ y); elseif (t <= 2e+44) 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[Exp[N[(y * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -440.0], t$95$1, If[LessEqual[t, 8e-216], N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2e+44], N[(x / N[Exp[N[(a * b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{e^{y \cdot t}}\\
\mathbf{if}\;t \leq -440:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 8 \cdot 10^{-216}:\\
\;\;\;\;x \cdot {z}^{y}\\
\mathbf{elif}\;t \leq 2 \cdot 10^{+44}:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -440 or 2.0000000000000002e44 < t Initial program 94.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6479.9%
Simplified79.9%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6479.9%
Applied egg-rr79.9%
if -440 < t < 8.0000000000000003e-216Initial program 97.6%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6467.9%
Simplified67.9%
Taylor expanded in t around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6467.9%
Simplified67.9%
if 8.0000000000000003e-216 < t < 2.0000000000000002e44Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6476.3%
Simplified76.3%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6476.3%
Applied egg-rr76.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* 0.5 (* a a))))
(if (<= a -2.2e+211)
(* x (+ 1.0 (* b (- (* 0.5 (* b (* a a))) a))))
(if (<= a -1.26e+33)
(/
x
(+
1.0
(*
b
(+ a (* b (+ (* 0.16666666666666666 (* b (* a (* a a)))) t_1))))))
(if (<= a 3.6e+61)
(* x (pow z y))
(/ x (+ 1.0 (* b (+ a (* b t_1))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = 0.5 * (a * a);
double tmp;
if (a <= -2.2e+211) {
tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a)));
} else if (a <= -1.26e+33) {
tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * (a * (a * a)))) + t_1)))));
} else if (a <= 3.6e+61) {
tmp = x * pow(z, y);
} else {
tmp = x / (1.0 + (b * (a + (b * 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 = 0.5d0 * (a * a)
if (a <= (-2.2d+211)) then
tmp = x * (1.0d0 + (b * ((0.5d0 * (b * (a * a))) - a)))
else if (a <= (-1.26d+33)) then
tmp = x / (1.0d0 + (b * (a + (b * ((0.16666666666666666d0 * (b * (a * (a * a)))) + t_1)))))
else if (a <= 3.6d+61) then
tmp = x * (z ** y)
else
tmp = x / (1.0d0 + (b * (a + (b * 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 = 0.5 * (a * a);
double tmp;
if (a <= -2.2e+211) {
tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a)));
} else if (a <= -1.26e+33) {
tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * (a * (a * a)))) + t_1)))));
} else if (a <= 3.6e+61) {
tmp = x * Math.pow(z, y);
} else {
tmp = x / (1.0 + (b * (a + (b * t_1))));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = 0.5 * (a * a) tmp = 0 if a <= -2.2e+211: tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a))) elif a <= -1.26e+33: tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * (a * (a * a)))) + t_1))))) elif a <= 3.6e+61: tmp = x * math.pow(z, y) else: tmp = x / (1.0 + (b * (a + (b * t_1)))) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(0.5 * Float64(a * a)) tmp = 0.0 if (a <= -2.2e+211) tmp = Float64(x * Float64(1.0 + Float64(b * Float64(Float64(0.5 * Float64(b * Float64(a * a))) - a)))); elseif (a <= -1.26e+33) tmp = Float64(x / Float64(1.0 + Float64(b * Float64(a + Float64(b * Float64(Float64(0.16666666666666666 * Float64(b * Float64(a * Float64(a * a)))) + t_1)))))); elseif (a <= 3.6e+61) tmp = Float64(x * (z ^ y)); else tmp = Float64(x / Float64(1.0 + Float64(b * Float64(a + Float64(b * t_1))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = 0.5 * (a * a); tmp = 0.0; if (a <= -2.2e+211) tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a))); elseif (a <= -1.26e+33) tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * (a * (a * a)))) + t_1))))); elseif (a <= 3.6e+61) tmp = x * (z ^ y); else tmp = x / (1.0 + (b * (a + (b * t_1)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(0.5 * N[(a * a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.2e+211], N[(x * N[(1.0 + N[(b * N[(N[(0.5 * N[(b * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, -1.26e+33], N[(x / N[(1.0 + N[(b * N[(a + N[(b * N[(N[(0.16666666666666666 * N[(b * N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 3.6e+61], N[(x * N[Power[z, y], $MachinePrecision]), $MachinePrecision], N[(x / N[(1.0 + N[(b * N[(a + N[(b * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.5 \cdot \left(a \cdot a\right)\\
\mathbf{if}\;a \leq -2.2 \cdot 10^{+211}:\\
\;\;\;\;x \cdot \left(1 + b \cdot \left(0.5 \cdot \left(b \cdot \left(a \cdot a\right)\right) - a\right)\right)\\
\mathbf{elif}\;a \leq -1.26 \cdot 10^{+33}:\\
\;\;\;\;\frac{x}{1 + b \cdot \left(a + b \cdot \left(0.16666666666666666 \cdot \left(b \cdot \left(a \cdot \left(a \cdot a\right)\right)\right) + t\_1\right)\right)}\\
\mathbf{elif}\;a \leq 3.6 \cdot 10^{+61}:\\
\;\;\;\;x \cdot {z}^{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{1 + b \cdot \left(a + b \cdot t\_1\right)}\\
\end{array}
\end{array}
if a < -2.20000000000000004e211Initial program 90.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6476.6%
Simplified76.6%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.2%
Simplified67.2%
if -2.20000000000000004e211 < a < -1.26e33Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6481.6%
Simplified81.6%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6481.6%
Applied egg-rr81.6%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.9%
Simplified53.9%
if -1.26e33 < a < 3.6000000000000001e61Initial program 100.0%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6487.1%
Simplified87.1%
Taylor expanded in t around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6465.0%
Simplified65.0%
if 3.6000000000000001e61 < a Initial program 82.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6473.8%
Simplified73.8%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6473.8%
Applied egg-rr73.8%
Taylor expanded in a around 0
distribute-rgt-inN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
Simplified61.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* x (pow z y)))) (if (<= y -0.000225) t_1 (if (<= y 8.5e+24) (/ 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 <= -0.000225) {
tmp = t_1;
} else if (y <= 8.5e+24) {
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 <= (-0.000225d0)) then
tmp = t_1
else if (y <= 8.5d+24) 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 <= -0.000225) {
tmp = t_1;
} else if (y <= 8.5e+24) {
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 <= -0.000225: tmp = t_1 elif y <= 8.5e+24: 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 <= -0.000225) tmp = t_1; elseif (y <= 8.5e+24) 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 <= -0.000225) tmp = t_1; elseif (y <= 8.5e+24) 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, -0.000225], t$95$1, If[LessEqual[y, 8.5e+24], 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 -0.000225:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 8.5 \cdot 10^{+24}:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.2499999999999999e-4 or 8.49999999999999959e24 < y Initial program 96.1%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6485.9%
Simplified85.9%
Taylor expanded in t around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6466.4%
Simplified66.4%
if -2.2499999999999999e-4 < y < 8.49999999999999959e24Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6478.4%
Simplified78.4%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6478.4%
Applied egg-rr78.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* a (* a a))))
(if (<= y -3.2e+76)
(* x (+ 1.0 (* t (- (* 0.5 (* t (* y y))) y))))
(if (<= y 1.55e-202)
(/
x
(+
1.0
(*
b
(+ a (* b (+ (* 0.16666666666666666 (* b t_1)) (* 0.5 (* a a))))))))
(* x (* b (* b (* b (* t_1 -0.16666666666666666)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = a * (a * a);
double tmp;
if (y <= -3.2e+76) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 1.55e-202) {
tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * t_1)) + (0.5 * (a * a)))))));
} else {
tmp = x * (b * (b * (b * (t_1 * -0.16666666666666666))));
}
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 = a * (a * a)
if (y <= (-3.2d+76)) then
tmp = x * (1.0d0 + (t * ((0.5d0 * (t * (y * y))) - y)))
else if (y <= 1.55d-202) then
tmp = x / (1.0d0 + (b * (a + (b * ((0.16666666666666666d0 * (b * t_1)) + (0.5d0 * (a * a)))))))
else
tmp = x * (b * (b * (b * (t_1 * (-0.16666666666666666d0)))))
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 = a * (a * a);
double tmp;
if (y <= -3.2e+76) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 1.55e-202) {
tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * t_1)) + (0.5 * (a * a)))))));
} else {
tmp = x * (b * (b * (b * (t_1 * -0.16666666666666666))));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = a * (a * a) tmp = 0 if y <= -3.2e+76: tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))) elif y <= 1.55e-202: tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * t_1)) + (0.5 * (a * a))))))) else: tmp = x * (b * (b * (b * (t_1 * -0.16666666666666666)))) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(a * Float64(a * a)) tmp = 0.0 if (y <= -3.2e+76) tmp = Float64(x * Float64(1.0 + Float64(t * Float64(Float64(0.5 * Float64(t * Float64(y * y))) - y)))); elseif (y <= 1.55e-202) tmp = Float64(x / Float64(1.0 + Float64(b * Float64(a + Float64(b * Float64(Float64(0.16666666666666666 * Float64(b * t_1)) + Float64(0.5 * Float64(a * a)))))))); else tmp = Float64(x * Float64(b * Float64(b * Float64(b * Float64(t_1 * -0.16666666666666666))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = a * (a * a); tmp = 0.0; if (y <= -3.2e+76) tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))); elseif (y <= 1.55e-202) tmp = x / (1.0 + (b * (a + (b * ((0.16666666666666666 * (b * t_1)) + (0.5 * (a * a))))))); else tmp = x * (b * (b * (b * (t_1 * -0.16666666666666666)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.2e+76], N[(x * N[(1.0 + N[(t * N[(N[(0.5 * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.55e-202], N[(x / N[(1.0 + N[(b * N[(a + N[(b * N[(N[(0.16666666666666666 * N[(b * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(0.5 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(b * N[(t$95$1 * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \left(a \cdot a\right)\\
\mathbf{if}\;y \leq -3.2 \cdot 10^{+76}:\\
\;\;\;\;x \cdot \left(1 + t \cdot \left(0.5 \cdot \left(t \cdot \left(y \cdot y\right)\right) - y\right)\right)\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-202}:\\
\;\;\;\;\frac{x}{1 + b \cdot \left(a + b \cdot \left(0.16666666666666666 \cdot \left(b \cdot t\_1\right) + 0.5 \cdot \left(a \cdot a\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(b \cdot \left(t\_1 \cdot -0.16666666666666666\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -3.19999999999999976e76Initial program 98.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6462.3%
Simplified62.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.1%
Simplified60.1%
if -3.19999999999999976e76 < y < 1.55e-202Initial program 95.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6475.4%
Applied egg-rr75.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.7%
Simplified55.7%
if 1.55e-202 < y Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified26.5%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.9%
Simplified51.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (* (* y y) (* 0.5 (* t t))))))
(if (<= y -8.2e-30)
t_1
(if (<= y 5e-18)
(* x (- 1.0 (* a b)))
(if (<= y 9e+157) t_1 (* x (- 0.0 (* a b))))))))
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 <= -8.2e-30) {
tmp = t_1;
} else if (y <= 5e-18) {
tmp = x * (1.0 - (a * b));
} else if (y <= 9e+157) {
tmp = t_1;
} else {
tmp = x * (0.0 - (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) :: t_1
real(8) :: tmp
t_1 = x * ((y * y) * (0.5d0 * (t * t)))
if (y <= (-8.2d-30)) then
tmp = t_1
else if (y <= 5d-18) then
tmp = x * (1.0d0 - (a * b))
else if (y <= 9d+157) then
tmp = t_1
else
tmp = x * (0.0d0 - (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 t_1 = x * ((y * y) * (0.5 * (t * t)));
double tmp;
if (y <= -8.2e-30) {
tmp = t_1;
} else if (y <= 5e-18) {
tmp = x * (1.0 - (a * b));
} else if (y <= 9e+157) {
tmp = t_1;
} else {
tmp = x * (0.0 - (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * ((y * y) * (0.5 * (t * t))) tmp = 0 if y <= -8.2e-30: tmp = t_1 elif y <= 5e-18: tmp = x * (1.0 - (a * b)) elif y <= 9e+157: tmp = t_1 else: tmp = x * (0.0 - (a * b)) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * Float64(Float64(y * y) * Float64(0.5 * Float64(t * t)))) tmp = 0.0 if (y <= -8.2e-30) tmp = t_1; elseif (y <= 5e-18) tmp = Float64(x * Float64(1.0 - Float64(a * b))); elseif (y <= 9e+157) tmp = t_1; else tmp = Float64(x * Float64(0.0 - Float64(a * b))); 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 <= -8.2e-30) tmp = t_1; elseif (y <= 5e-18) tmp = x * (1.0 - (a * b)); elseif (y <= 9e+157) tmp = t_1; else tmp = x * (0.0 - (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x * N[(N[(y * y), $MachinePrecision] * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -8.2e-30], t$95$1, If[LessEqual[y, 5e-18], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 9e+157], t$95$1, N[(x * N[(0.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(\left(y \cdot y\right) \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\\
\mathbf{if}\;y \leq -8.2 \cdot 10^{-30}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 5 \cdot 10^{-18}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\mathbf{elif}\;y \leq 9 \cdot 10^{+157}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0 - a \cdot b\right)\\
\end{array}
\end{array}
if y < -8.2000000000000007e-30 or 5.00000000000000036e-18 < y < 8.9999999999999997e157Initial program 96.3%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6455.7%
Simplified55.7%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6433.0%
Simplified33.0%
Taylor expanded in y around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6438.2%
Simplified38.2%
Taylor expanded in t around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.2%
Simplified44.2%
if -8.2000000000000007e-30 < y < 5.00000000000000036e-18Initial program 96.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.5%
Simplified80.5%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6443.9%
Simplified43.9%
if 8.9999999999999997e157 < y Initial program 93.3%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6425.9%
Simplified25.9%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6413.0%
Simplified13.0%
Taylor expanded in a around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6443.2%
Simplified43.2%
Final simplification44.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -8e+77)
(* x (+ 1.0 (* t (- (* 0.5 (* t (* y y))) y))))
(if (<= y 1.55e-202)
(/ x (+ 1.0 (* b (+ a (* b (* 0.5 (* a a)))))))
(* x (* b (* b (* b (* (* a (* a a)) -0.16666666666666666))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -8e+77) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 1.55e-202) {
tmp = x / (1.0 + (b * (a + (b * (0.5 * (a * a))))));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-8d+77)) then
tmp = x * (1.0d0 + (t * ((0.5d0 * (t * (y * y))) - y)))
else if (y <= 1.55d-202) then
tmp = x / (1.0d0 + (b * (a + (b * (0.5d0 * (a * a))))))
else
tmp = x * (b * (b * (b * ((a * (a * a)) * (-0.16666666666666666d0)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -8e+77) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 1.55e-202) {
tmp = x / (1.0 + (b * (a + (b * (0.5 * (a * a))))));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -8e+77: tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))) elif y <= 1.55e-202: tmp = x / (1.0 + (b * (a + (b * (0.5 * (a * a)))))) else: tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -8e+77) tmp = Float64(x * Float64(1.0 + Float64(t * Float64(Float64(0.5 * Float64(t * Float64(y * y))) - y)))); elseif (y <= 1.55e-202) tmp = Float64(x / Float64(1.0 + Float64(b * Float64(a + Float64(b * Float64(0.5 * Float64(a * a))))))); else tmp = Float64(x * Float64(b * Float64(b * Float64(b * Float64(Float64(a * Float64(a * a)) * -0.16666666666666666))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -8e+77) tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))); elseif (y <= 1.55e-202) tmp = x / (1.0 + (b * (a + (b * (0.5 * (a * a)))))); else tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -8e+77], N[(x * N[(1.0 + N[(t * N[(N[(0.5 * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.55e-202], N[(x / N[(1.0 + N[(b * N[(a + N[(b * N[(0.5 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(b * N[(N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8 \cdot 10^{+77}:\\
\;\;\;\;x \cdot \left(1 + t \cdot \left(0.5 \cdot \left(t \cdot \left(y \cdot y\right)\right) - y\right)\right)\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-202}:\\
\;\;\;\;\frac{x}{1 + b \cdot \left(a + b \cdot \left(0.5 \cdot \left(a \cdot a\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(b \cdot \left(\left(a \cdot \left(a \cdot a\right)\right) \cdot -0.16666666666666666\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -7.99999999999999986e77Initial program 98.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6462.3%
Simplified62.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.1%
Simplified60.1%
if -7.99999999999999986e77 < y < 1.55e-202Initial program 95.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6475.4%
Applied egg-rr75.4%
Taylor expanded in a around 0
distribute-rgt-inN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-inN/A
+-lowering-+.f64N/A
Simplified53.0%
if 1.55e-202 < y Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified26.5%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.9%
Simplified51.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -5.6e+79)
(* x (+ 1.0 (* t (- (* 0.5 (* t (* y y))) y))))
(if (<= y 1.55e-202)
(/ x (+ 1.0 (* a (+ b (* 0.5 (* a (* b b)))))))
(* x (* b (* b (* b (* (* a (* a a)) -0.16666666666666666))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -5.6e+79) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 1.55e-202) {
tmp = x / (1.0 + (a * (b + (0.5 * (a * (b * b))))));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-5.6d+79)) then
tmp = x * (1.0d0 + (t * ((0.5d0 * (t * (y * y))) - y)))
else if (y <= 1.55d-202) then
tmp = x / (1.0d0 + (a * (b + (0.5d0 * (a * (b * b))))))
else
tmp = x * (b * (b * (b * ((a * (a * a)) * (-0.16666666666666666d0)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -5.6e+79) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 1.55e-202) {
tmp = x / (1.0 + (a * (b + (0.5 * (a * (b * b))))));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -5.6e+79: tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))) elif y <= 1.55e-202: tmp = x / (1.0 + (a * (b + (0.5 * (a * (b * b)))))) else: tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -5.6e+79) tmp = Float64(x * Float64(1.0 + Float64(t * Float64(Float64(0.5 * Float64(t * Float64(y * y))) - y)))); elseif (y <= 1.55e-202) tmp = Float64(x / Float64(1.0 + Float64(a * Float64(b + Float64(0.5 * Float64(a * Float64(b * b))))))); else tmp = Float64(x * Float64(b * Float64(b * Float64(b * Float64(Float64(a * Float64(a * a)) * -0.16666666666666666))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -5.6e+79) tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))); elseif (y <= 1.55e-202) tmp = x / (1.0 + (a * (b + (0.5 * (a * (b * b)))))); else tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -5.6e+79], N[(x * N[(1.0 + N[(t * N[(N[(0.5 * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.55e-202], N[(x / N[(1.0 + N[(a * N[(b + N[(0.5 * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(b * N[(N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.6 \cdot 10^{+79}:\\
\;\;\;\;x \cdot \left(1 + t \cdot \left(0.5 \cdot \left(t \cdot \left(y \cdot y\right)\right) - y\right)\right)\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-202}:\\
\;\;\;\;\frac{x}{1 + a \cdot \left(b + 0.5 \cdot \left(a \cdot \left(b \cdot b\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(b \cdot \left(\left(a \cdot \left(a \cdot a\right)\right) \cdot -0.16666666666666666\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -5.6000000000000002e79Initial program 98.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6462.3%
Simplified62.3%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6460.1%
Simplified60.1%
if -5.6000000000000002e79 < y < 1.55e-202Initial program 95.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6475.4%
Simplified75.4%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6475.4%
Applied egg-rr75.4%
exp-prodN/A
sqr-powN/A
pow2N/A
pow-lowering-pow.f64N/A
pow-expN/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6475.4%
Applied egg-rr75.4%
Taylor expanded in a around 0
distribute-rgt-inN/A
*-commutativeN/A
associate-+r+N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
associate-*l*N/A
associate-+r+N/A
distribute-rgt-inN/A
+-lowering-+.f64N/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
Simplified50.9%
if 1.55e-202 < y Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified26.5%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.9%
Simplified51.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -8.2e+71)
(* x (+ 1.0 (* t (- (* 0.5 (* t (* y y))) y))))
(if (<= y 8.5e-192)
(* x (+ 1.0 (* b (- (* 0.5 (* b (* a a))) a))))
(* x (* b (* b (* b (* (* a (* a a)) -0.16666666666666666))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -8.2e+71) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 8.5e-192) {
tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a)));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
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 <= (-8.2d+71)) then
tmp = x * (1.0d0 + (t * ((0.5d0 * (t * (y * y))) - y)))
else if (y <= 8.5d-192) then
tmp = x * (1.0d0 + (b * ((0.5d0 * (b * (a * a))) - a)))
else
tmp = x * (b * (b * (b * ((a * (a * a)) * (-0.16666666666666666d0)))))
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 <= -8.2e+71) {
tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y)));
} else if (y <= 8.5e-192) {
tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a)));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -8.2e+71: tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))) elif y <= 8.5e-192: tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a))) else: tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -8.2e+71) tmp = Float64(x * Float64(1.0 + Float64(t * Float64(Float64(0.5 * Float64(t * Float64(y * y))) - y)))); elseif (y <= 8.5e-192) tmp = Float64(x * Float64(1.0 + Float64(b * Float64(Float64(0.5 * Float64(b * Float64(a * a))) - a)))); else tmp = Float64(x * Float64(b * Float64(b * Float64(b * Float64(Float64(a * Float64(a * a)) * -0.16666666666666666))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -8.2e+71) tmp = x * (1.0 + (t * ((0.5 * (t * (y * y))) - y))); elseif (y <= 8.5e-192) tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a))); else tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -8.2e+71], N[(x * N[(1.0 + N[(t * N[(N[(0.5 * N[(t * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.5e-192], N[(x * N[(1.0 + N[(b * N[(N[(0.5 * N[(b * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(b * N[(N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8.2 \cdot 10^{+71}:\\
\;\;\;\;x \cdot \left(1 + t \cdot \left(0.5 \cdot \left(t \cdot \left(y \cdot y\right)\right) - y\right)\right)\\
\mathbf{elif}\;y \leq 8.5 \cdot 10^{-192}:\\
\;\;\;\;x \cdot \left(1 + b \cdot \left(0.5 \cdot \left(b \cdot \left(a \cdot a\right)\right) - a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(b \cdot \left(\left(a \cdot \left(a \cdot a\right)\right) \cdot -0.16666666666666666\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -8.2000000000000004e71Initial program 98.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6461.8%
Simplified61.8%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.8%
Simplified57.8%
if -8.2000000000000004e71 < y < 8.49999999999999985e-192Initial program 95.2%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6475.9%
Simplified75.9%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.1%
Simplified48.1%
if 8.49999999999999985e-192 < y Initial program 95.9%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6448.9%
Simplified48.9%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified25.3%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.0%
Simplified52.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -5.7e+25)
(* x (* (* y y) (* 0.5 (* t t))))
(if (<= y 8.5e-192)
(* x (+ 1.0 (* b (- (* 0.5 (* b (* a a))) a))))
(* x (* b (* b (* b (* (* a (* a a)) -0.16666666666666666))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -5.7e+25) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 8.5e-192) {
tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a)));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-5.7d+25)) then
tmp = x * ((y * y) * (0.5d0 * (t * t)))
else if (y <= 8.5d-192) then
tmp = x * (1.0d0 + (b * ((0.5d0 * (b * (a * a))) - a)))
else
tmp = x * (b * (b * (b * ((a * (a * a)) * (-0.16666666666666666d0)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -5.7e+25) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 8.5e-192) {
tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a)));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -5.7e+25: tmp = x * ((y * y) * (0.5 * (t * t))) elif y <= 8.5e-192: tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a))) else: tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -5.7e+25) tmp = Float64(x * Float64(Float64(y * y) * Float64(0.5 * Float64(t * t)))); elseif (y <= 8.5e-192) tmp = Float64(x * Float64(1.0 + Float64(b * Float64(Float64(0.5 * Float64(b * Float64(a * a))) - a)))); else tmp = Float64(x * Float64(b * Float64(b * Float64(b * Float64(Float64(a * Float64(a * a)) * -0.16666666666666666))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -5.7e+25) tmp = x * ((y * y) * (0.5 * (t * t))); elseif (y <= 8.5e-192) tmp = x * (1.0 + (b * ((0.5 * (b * (a * a))) - a))); else tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -5.7e+25], N[(x * N[(N[(y * y), $MachinePrecision] * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.5e-192], N[(x * N[(1.0 + N[(b * N[(N[(0.5 * N[(b * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(b * N[(N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.7 \cdot 10^{+25}:\\
\;\;\;\;x \cdot \left(\left(y \cdot y\right) \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\\
\mathbf{elif}\;y \leq 8.5 \cdot 10^{-192}:\\
\;\;\;\;x \cdot \left(1 + b \cdot \left(0.5 \cdot \left(b \cdot \left(a \cdot a\right)\right) - a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(b \cdot \left(\left(a \cdot \left(a \cdot a\right)\right) \cdot -0.16666666666666666\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -5.6999999999999996e25Initial program 96.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6461.5%
Simplified61.5%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6435.3%
Simplified35.3%
Taylor expanded in y around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6444.0%
Simplified44.0%
Taylor expanded in t around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.9%
Simplified46.9%
if -5.6999999999999996e25 < y < 8.49999999999999985e-192Initial program 95.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6479.9%
Simplified79.9%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.8%
Simplified51.8%
if 8.49999999999999985e-192 < y Initial program 95.9%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6448.9%
Simplified48.9%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified25.3%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.0%
Simplified52.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -6.6e-33)
(* x (* (* y y) (* 0.5 (* t t))))
(if (<= y 1.55e-202)
(* x (+ 1.0 (* a (- (* (* b b) (* a 0.5)) b))))
(* x (* b (* b (* b (* (* a (* a a)) -0.16666666666666666))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -6.6e-33) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 1.55e-202) {
tmp = x * (1.0 + (a * (((b * b) * (a * 0.5)) - b)));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
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 <= (-6.6d-33)) then
tmp = x * ((y * y) * (0.5d0 * (t * t)))
else if (y <= 1.55d-202) then
tmp = x * (1.0d0 + (a * (((b * b) * (a * 0.5d0)) - b)))
else
tmp = x * (b * (b * (b * ((a * (a * a)) * (-0.16666666666666666d0)))))
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 <= -6.6e-33) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 1.55e-202) {
tmp = x * (1.0 + (a * (((b * b) * (a * 0.5)) - b)));
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -6.6e-33: tmp = x * ((y * y) * (0.5 * (t * t))) elif y <= 1.55e-202: tmp = x * (1.0 + (a * (((b * b) * (a * 0.5)) - b))) else: tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -6.6e-33) tmp = Float64(x * Float64(Float64(y * y) * Float64(0.5 * Float64(t * t)))); elseif (y <= 1.55e-202) tmp = Float64(x * Float64(1.0 + Float64(a * Float64(Float64(Float64(b * b) * Float64(a * 0.5)) - b)))); else tmp = Float64(x * Float64(b * Float64(b * Float64(b * Float64(Float64(a * Float64(a * a)) * -0.16666666666666666))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -6.6e-33) tmp = x * ((y * y) * (0.5 * (t * t))); elseif (y <= 1.55e-202) tmp = x * (1.0 + (a * (((b * b) * (a * 0.5)) - b))); else tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -6.6e-33], N[(x * N[(N[(y * y), $MachinePrecision] * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.55e-202], N[(x * N[(1.0 + N[(a * N[(N[(N[(b * b), $MachinePrecision] * N[(a * 0.5), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(b * N[(N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.6 \cdot 10^{-33}:\\
\;\;\;\;x \cdot \left(\left(y \cdot y\right) \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-202}:\\
\;\;\;\;x \cdot \left(1 + a \cdot \left(\left(b \cdot b\right) \cdot \left(a \cdot 0.5\right) - b\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(b \cdot \left(\left(a \cdot \left(a \cdot a\right)\right) \cdot -0.16666666666666666\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -6.6000000000000005e-33Initial program 97.4%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6458.8%
Simplified58.8%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6433.0%
Simplified33.0%
Taylor expanded in y around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6440.3%
Simplified40.3%
Taylor expanded in t around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6443.7%
Simplified43.7%
if -6.6000000000000005e-33 < y < 1.55e-202Initial program 94.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6484.1%
Simplified84.1%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.6%
Simplified51.6%
if 1.55e-202 < y Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified26.5%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.9%
Simplified51.9%
Final simplification49.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -2.9e-30)
(* x (* (* y y) (* 0.5 (* t t))))
(if (<= y 1.55e-202)
(/ x (+ (* a b) 1.0))
(* x (* b (* b (* b (* (* a (* a a)) -0.16666666666666666))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -2.9e-30) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 1.55e-202) {
tmp = x / ((a * b) + 1.0);
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
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 <= (-2.9d-30)) then
tmp = x * ((y * y) * (0.5d0 * (t * t)))
else if (y <= 1.55d-202) then
tmp = x / ((a * b) + 1.0d0)
else
tmp = x * (b * (b * (b * ((a * (a * a)) * (-0.16666666666666666d0)))))
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 <= -2.9e-30) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 1.55e-202) {
tmp = x / ((a * b) + 1.0);
} else {
tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -2.9e-30: tmp = x * ((y * y) * (0.5 * (t * t))) elif y <= 1.55e-202: tmp = x / ((a * b) + 1.0) else: tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -2.9e-30) tmp = Float64(x * Float64(Float64(y * y) * Float64(0.5 * Float64(t * t)))); elseif (y <= 1.55e-202) tmp = Float64(x / Float64(Float64(a * b) + 1.0)); else tmp = Float64(x * Float64(b * Float64(b * Float64(b * Float64(Float64(a * Float64(a * a)) * -0.16666666666666666))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -2.9e-30) tmp = x * ((y * y) * (0.5 * (t * t))); elseif (y <= 1.55e-202) tmp = x / ((a * b) + 1.0); else tmp = x * (b * (b * (b * ((a * (a * a)) * -0.16666666666666666)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -2.9e-30], N[(x * N[(N[(y * y), $MachinePrecision] * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.55e-202], N[(x / N[(N[(a * b), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(b * N[(N[(a * N[(a * a), $MachinePrecision]), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.9 \cdot 10^{-30}:\\
\;\;\;\;x \cdot \left(\left(y \cdot y\right) \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\\
\mathbf{elif}\;y \leq 1.55 \cdot 10^{-202}:\\
\;\;\;\;\frac{x}{a \cdot b + 1}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(b \cdot \left(\left(a \cdot \left(a \cdot a\right)\right) \cdot -0.16666666666666666\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -2.89999999999999989e-30Initial program 97.3%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6458.3%
Simplified58.3%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.0%
Simplified32.0%
Taylor expanded in y around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6439.5%
Simplified39.5%
Taylor expanded in t around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.2%
Simplified44.2%
if -2.89999999999999989e-30 < y < 1.55e-202Initial program 94.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6484.3%
Simplified84.3%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6484.2%
Applied egg-rr84.2%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6450.0%
Simplified50.0%
if 1.55e-202 < y Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6450.0%
Simplified50.0%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified26.5%
Taylor expanded in b around inf
associate-*r*N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/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
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.9%
Simplified51.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -4e-29)
(* x (* (* y y) (* 0.5 (* t t))))
(if (<= y 2.4e-215)
(/ x (+ (* a b) 1.0))
(* x (* a (* a (* a (* -0.16666666666666666 (* b (* b b))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -4e-29) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 2.4e-215) {
tmp = x / ((a * b) + 1.0);
} else {
tmp = x * (a * (a * (a * (-0.16666666666666666 * (b * (b * b))))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-4d-29)) then
tmp = x * ((y * y) * (0.5d0 * (t * t)))
else if (y <= 2.4d-215) then
tmp = x / ((a * b) + 1.0d0)
else
tmp = x * (a * (a * (a * ((-0.16666666666666666d0) * (b * (b * b))))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -4e-29) {
tmp = x * ((y * y) * (0.5 * (t * t)));
} else if (y <= 2.4e-215) {
tmp = x / ((a * b) + 1.0);
} else {
tmp = x * (a * (a * (a * (-0.16666666666666666 * (b * (b * b))))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -4e-29: tmp = x * ((y * y) * (0.5 * (t * t))) elif y <= 2.4e-215: tmp = x / ((a * b) + 1.0) else: tmp = x * (a * (a * (a * (-0.16666666666666666 * (b * (b * b)))))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -4e-29) tmp = Float64(x * Float64(Float64(y * y) * Float64(0.5 * Float64(t * t)))); elseif (y <= 2.4e-215) tmp = Float64(x / Float64(Float64(a * b) + 1.0)); else tmp = Float64(x * Float64(a * Float64(a * Float64(a * Float64(-0.16666666666666666 * Float64(b * Float64(b * b))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -4e-29) tmp = x * ((y * y) * (0.5 * (t * t))); elseif (y <= 2.4e-215) tmp = x / ((a * b) + 1.0); else tmp = x * (a * (a * (a * (-0.16666666666666666 * (b * (b * b)))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -4e-29], N[(x * N[(N[(y * y), $MachinePrecision] * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-215], N[(x / N[(N[(a * b), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(x * N[(a * N[(a * N[(a * N[(-0.16666666666666666 * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4 \cdot 10^{-29}:\\
\;\;\;\;x \cdot \left(\left(y \cdot y\right) \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-215}:\\
\;\;\;\;\frac{x}{a \cdot b + 1}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(a \cdot \left(a \cdot \left(a \cdot \left(-0.16666666666666666 \cdot \left(b \cdot \left(b \cdot b\right)\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -3.99999999999999977e-29Initial program 97.3%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6458.3%
Simplified58.3%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.0%
Simplified32.0%
Taylor expanded in y around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6439.5%
Simplified39.5%
Taylor expanded in t around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.2%
Simplified44.2%
if -3.99999999999999977e-29 < y < 2.4000000000000001e-215Initial program 95.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6485.5%
Simplified85.5%
sub0-negN/A
exp-negN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6485.5%
Applied egg-rr85.5%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6452.4%
Simplified52.4%
if 2.4000000000000001e-215 < y Initial program 95.3%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6451.1%
Simplified51.1%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified26.0%
Taylor expanded in a around inf
cube-multN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
Simplified20.6%
Taylor expanded in a around inf
*-commutativeN/A
unpow2N/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.2%
Simplified47.2%
Final simplification47.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* 0.5 (* t t))))
(if (<= y -6.2e-30)
(* x (* (* y y) t_1))
(if (<= y 1.5e-17) (* x (- 1.0 (* a b))) (* x (* y (* y t_1)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = 0.5 * (t * t);
double tmp;
if (y <= -6.2e-30) {
tmp = x * ((y * y) * t_1);
} else if (y <= 1.5e-17) {
tmp = x * (1.0 - (a * b));
} else {
tmp = x * (y * (y * 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 = 0.5d0 * (t * t)
if (y <= (-6.2d-30)) then
tmp = x * ((y * y) * t_1)
else if (y <= 1.5d-17) then
tmp = x * (1.0d0 - (a * b))
else
tmp = x * (y * (y * 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 = 0.5 * (t * t);
double tmp;
if (y <= -6.2e-30) {
tmp = x * ((y * y) * t_1);
} else if (y <= 1.5e-17) {
tmp = x * (1.0 - (a * b));
} else {
tmp = x * (y * (y * t_1));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = 0.5 * (t * t) tmp = 0 if y <= -6.2e-30: tmp = x * ((y * y) * t_1) elif y <= 1.5e-17: tmp = x * (1.0 - (a * b)) else: tmp = x * (y * (y * t_1)) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(0.5 * Float64(t * t)) tmp = 0.0 if (y <= -6.2e-30) tmp = Float64(x * Float64(Float64(y * y) * t_1)); elseif (y <= 1.5e-17) tmp = Float64(x * Float64(1.0 - Float64(a * b))); else tmp = Float64(x * Float64(y * Float64(y * t_1))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = 0.5 * (t * t); tmp = 0.0; if (y <= -6.2e-30) tmp = x * ((y * y) * t_1); elseif (y <= 1.5e-17) tmp = x * (1.0 - (a * b)); else tmp = x * (y * (y * t_1)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6.2e-30], N[(x * N[(N[(y * y), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.5e-17], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(y * N[(y * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.5 \cdot \left(t \cdot t\right)\\
\mathbf{if}\;y \leq -6.2 \cdot 10^{-30}:\\
\;\;\;\;x \cdot \left(\left(y \cdot y\right) \cdot t\_1\right)\\
\mathbf{elif}\;y \leq 1.5 \cdot 10^{-17}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(y \cdot \left(y \cdot t\_1\right)\right)\\
\end{array}
\end{array}
if y < -6.19999999999999982e-30Initial program 97.3%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6458.3%
Simplified58.3%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.0%
Simplified32.0%
Taylor expanded in y around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6439.5%
Simplified39.5%
Taylor expanded in t around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6444.2%
Simplified44.2%
if -6.19999999999999982e-30 < y < 1.50000000000000003e-17Initial program 96.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.5%
Simplified80.5%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6443.9%
Simplified43.9%
if 1.50000000000000003e-17 < y Initial program 93.7%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6455.2%
Simplified55.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6430.5%
Simplified30.5%
Taylor expanded in y around inf
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
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-*.f6445.9%
Simplified45.9%
Final simplification44.5%
(FPCore (x y z t a b) :precision binary64 (if (<= y -6.2e+22) (- x (* t (* x y))) (if (<= y 1.1e-17) (* x (- 1.0 (* a b))) (* x (- 0.0 (* a b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -6.2e+22) {
tmp = x - (t * (x * y));
} else if (y <= 1.1e-17) {
tmp = x * (1.0 - (a * b));
} else {
tmp = x * (0.0 - (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 (y <= (-6.2d+22)) then
tmp = x - (t * (x * y))
else if (y <= 1.1d-17) then
tmp = x * (1.0d0 - (a * b))
else
tmp = x * (0.0d0 - (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 (y <= -6.2e+22) {
tmp = x - (t * (x * y));
} else if (y <= 1.1e-17) {
tmp = x * (1.0 - (a * b));
} else {
tmp = x * (0.0 - (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -6.2e+22: tmp = x - (t * (x * y)) elif y <= 1.1e-17: tmp = x * (1.0 - (a * b)) else: tmp = x * (0.0 - (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -6.2e+22) tmp = Float64(x - Float64(t * Float64(x * y))); elseif (y <= 1.1e-17) tmp = Float64(x * Float64(1.0 - Float64(a * b))); else tmp = Float64(x * Float64(0.0 - Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -6.2e+22) tmp = x - (t * (x * y)); elseif (y <= 1.1e-17) tmp = x * (1.0 - (a * b)); else tmp = x * (0.0 - (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -6.2e+22], N[(x - N[(t * N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.1e-17], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(0.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.2 \cdot 10^{+22}:\\
\;\;\;\;x - t \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{-17}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0 - a \cdot b\right)\\
\end{array}
\end{array}
if y < -6.2000000000000004e22Initial program 96.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6460.5%
Simplified60.5%
Taylor expanded in y around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6427.0%
Simplified27.0%
if -6.2000000000000004e22 < y < 1.1e-17Initial program 96.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6478.4%
Simplified78.4%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6440.9%
Simplified40.9%
if 1.1e-17 < y Initial program 93.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6436.9%
Simplified36.9%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f649.2%
Simplified9.2%
Taylor expanded in a around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6429.8%
Simplified29.8%
Final simplification34.7%
(FPCore (x y z t a b) :precision binary64 (if (<= y -1.45e-28) (- 0.0 (* y (* x t))) (if (<= y 1.8e-16) (* x (- 1.0 (* a b))) (* x (- 0.0 (* a b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.45e-28) {
tmp = 0.0 - (y * (x * t));
} else if (y <= 1.8e-16) {
tmp = x * (1.0 - (a * b));
} else {
tmp = x * (0.0 - (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 (y <= (-1.45d-28)) then
tmp = 0.0d0 - (y * (x * t))
else if (y <= 1.8d-16) then
tmp = x * (1.0d0 - (a * b))
else
tmp = x * (0.0d0 - (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 (y <= -1.45e-28) {
tmp = 0.0 - (y * (x * t));
} else if (y <= 1.8e-16) {
tmp = x * (1.0 - (a * b));
} else {
tmp = x * (0.0 - (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -1.45e-28: tmp = 0.0 - (y * (x * t)) elif y <= 1.8e-16: tmp = x * (1.0 - (a * b)) else: tmp = x * (0.0 - (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -1.45e-28) tmp = Float64(0.0 - Float64(y * Float64(x * t))); elseif (y <= 1.8e-16) tmp = Float64(x * Float64(1.0 - Float64(a * b))); else tmp = Float64(x * Float64(0.0 - Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -1.45e-28) tmp = 0.0 - (y * (x * t)); elseif (y <= 1.8e-16) tmp = x * (1.0 - (a * b)); else tmp = x * (0.0 - (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -1.45e-28], N[(0.0 - N[(y * N[(x * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.8e-16], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(0.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.45 \cdot 10^{-28}:\\
\;\;\;\;0 - y \cdot \left(x \cdot t\right)\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{-16}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0 - a \cdot b\right)\\
\end{array}
\end{array}
if y < -1.45000000000000006e-28Initial program 97.3%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6458.3%
Simplified58.3%
Taylor expanded in y around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6423.9%
Simplified23.9%
Taylor expanded in t around inf
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6419.7%
Simplified19.7%
if -1.45000000000000006e-28 < y < 1.79999999999999991e-16Initial program 96.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.5%
Simplified80.5%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6443.9%
Simplified43.9%
if 1.79999999999999991e-16 < y Initial program 93.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6436.9%
Simplified36.9%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f649.2%
Simplified9.2%
Taylor expanded in a around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6429.8%
Simplified29.8%
Final simplification33.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (- 0.0 (* y (* x t))))) (if (<= y -1.25e-28) t_1 (if (<= y 1.75e-16) (* x (- 1.0 (* a b))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = 0.0 - (y * (x * t));
double tmp;
if (y <= -1.25e-28) {
tmp = t_1;
} else if (y <= 1.75e-16) {
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 = 0.0d0 - (y * (x * t))
if (y <= (-1.25d-28)) then
tmp = t_1
else if (y <= 1.75d-16) 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 = 0.0 - (y * (x * t));
double tmp;
if (y <= -1.25e-28) {
tmp = t_1;
} else if (y <= 1.75e-16) {
tmp = x * (1.0 - (a * b));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = 0.0 - (y * (x * t)) tmp = 0 if y <= -1.25e-28: tmp = t_1 elif y <= 1.75e-16: tmp = x * (1.0 - (a * b)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(0.0 - Float64(y * Float64(x * t))) tmp = 0.0 if (y <= -1.25e-28) tmp = t_1; elseif (y <= 1.75e-16) 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 = 0.0 - (y * (x * t)); tmp = 0.0; if (y <= -1.25e-28) tmp = t_1; elseif (y <= 1.75e-16) 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[(0.0 - N[(y * N[(x * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.25e-28], t$95$1, If[LessEqual[y, 1.75e-16], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0 - y \cdot \left(x \cdot t\right)\\
\mathbf{if}\;y \leq -1.25 \cdot 10^{-28}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.75 \cdot 10^{-16}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.25e-28 or 1.75000000000000009e-16 < y Initial program 95.7%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6456.9%
Simplified56.9%
Taylor expanded in y around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6417.2%
Simplified17.2%
Taylor expanded in t around inf
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6419.7%
Simplified19.7%
if -1.25e-28 < y < 1.75000000000000009e-16Initial program 96.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.5%
Simplified80.5%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6443.9%
Simplified43.9%
Final simplification31.0%
(FPCore (x y z t a b) :precision binary64 (if (<= y 8.5e-192) (* t (- (/ x t) (* x y))) (* x (- 0.0 (* a b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 8.5e-192) {
tmp = t * ((x / t) - (x * y));
} else {
tmp = x * (0.0 - (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 (y <= 8.5d-192) then
tmp = t * ((x / t) - (x * y))
else
tmp = x * (0.0d0 - (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 (y <= 8.5e-192) {
tmp = t * ((x / t) - (x * y));
} else {
tmp = x * (0.0 - (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= 8.5e-192: tmp = t * ((x / t) - (x * y)) else: tmp = x * (0.0 - (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= 8.5e-192) tmp = Float64(t * Float64(Float64(x / t) - Float64(x * y))); else tmp = Float64(x * Float64(0.0 - Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= 8.5e-192) tmp = t * ((x / t) - (x * y)); else tmp = x * (0.0 - (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, 8.5e-192], N[(t * N[(N[(x / t), $MachinePrecision] - N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(0.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 8.5 \cdot 10^{-192}:\\
\;\;\;\;t \cdot \left(\frac{x}{t} - x \cdot y\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(0 - a \cdot b\right)\\
\end{array}
\end{array}
if y < 8.49999999999999985e-192Initial program 96.1%
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
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6434.0%
Simplified34.0%
Taylor expanded in t around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6436.7%
Simplified36.7%
if 8.49999999999999985e-192 < y Initial program 95.9%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6448.9%
Simplified48.9%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6417.5%
Simplified17.5%
Taylor expanded in a around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6428.9%
Simplified28.9%
Final simplification33.8%
(FPCore (x y z t a b) :precision binary64 (if (<= b 6400000.0) (* x (- 1.0 (* y t))) (* x (- 1.0 (* a b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 6400000.0) {
tmp = x * (1.0 - (y * t));
} else {
tmp = x * (1.0 - (a * b));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (b <= 6400000.0d0) then
tmp = x * (1.0d0 - (y * t))
else
tmp = x * (1.0d0 - (a * b))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= 6400000.0) {
tmp = x * (1.0 - (y * t));
} else {
tmp = x * (1.0 - (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= 6400000.0: tmp = x * (1.0 - (y * t)) else: tmp = x * (1.0 - (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= 6400000.0) tmp = Float64(x * Float64(1.0 - Float64(y * t))); else tmp = Float64(x * Float64(1.0 - Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= 6400000.0) tmp = x * (1.0 - (y * t)); else tmp = x * (1.0 - (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, 6400000.0], N[(x * N[(1.0 - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 6400000:\\
\;\;\;\;x \cdot \left(1 - y \cdot t\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\end{array}
\end{array}
if b < 6.4e6Initial program 95.3%
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 y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6429.3%
Simplified29.3%
if 6.4e6 < b Initial program 98.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6467.5%
Simplified67.5%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6426.7%
Simplified26.7%
(FPCore (x y z t a b) :precision binary64 (* x (- 1.0 (* a b))))
double code(double x, double y, double z, double t, double a, double b) {
return x * (1.0 - (a * 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 * (1.0d0 - (a * b))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x * (1.0 - (a * b));
}
def code(x, y, z, t, a, b): return x * (1.0 - (a * b))
function code(x, y, z, t, a, b) return Float64(x * Float64(1.0 - Float64(a * b))) end
function tmp = code(x, y, z, t, a, b) tmp = x * (1.0 - (a * b)); end
code[x_, y_, z_, t_, a_, b_] := N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 - a \cdot b\right)
\end{array}
Initial program 96.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6458.8%
Simplified58.8%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6425.1%
Simplified25.1%
(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 96.0%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6469.6%
Simplified69.6%
Taylor expanded in y around 0
Simplified18.0%
herbie shell --seed 2024152
(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))))))