
(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 19 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 97.7%
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-+.f6499.6%
Simplified99.6%
Final simplification99.6%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (exp (* y (- (log z) t))))))
(if (<= y -9.5e-19)
t_1
(if (<= y 1.15e-64) (* x (exp (* a (- (log (- 1.0 z)) 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 <= -9.5e-19) {
tmp = t_1;
} else if (y <= 1.15e-64) {
tmp = x * exp((a * (log((1.0 - z)) - 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 <= (-9.5d-19)) then
tmp = t_1
else if (y <= 1.15d-64) then
tmp = x * exp((a * (log((1.0d0 - z)) - 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 <= -9.5e-19) {
tmp = t_1;
} else if (y <= 1.15e-64) {
tmp = x * Math.exp((a * (Math.log((1.0 - z)) - 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 <= -9.5e-19: tmp = t_1 elif y <= 1.15e-64: tmp = x * math.exp((a * (math.log((1.0 - z)) - 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 <= -9.5e-19) tmp = t_1; elseif (y <= 1.15e-64) tmp = Float64(x * exp(Float64(a * Float64(log(Float64(1.0 - z)) - 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 <= -9.5e-19) tmp = t_1; elseif (y <= 1.15e-64) tmp = x * exp((a * (log((1.0 - z)) - 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, -9.5e-19], t$95$1, If[LessEqual[y, 1.15e-64], N[(x * N[Exp[N[(a * N[(N[Log[N[(1.0 - z), $MachinePrecision]], $MachinePrecision] - 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 -9.5 \cdot 10^{-19}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.15 \cdot 10^{-64}:\\
\;\;\;\;x \cdot e^{a \cdot \left(\log \left(1 - z\right) - b\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -9.4999999999999995e-19 or 1.1500000000000001e-64 < y Initial program 98.6%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6488.0%
Simplified88.0%
if -9.4999999999999995e-19 < y < 1.1500000000000001e-64Initial program 96.6%
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--.f6486.5%
Simplified86.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (exp (* y (- (log z) t))))))
(if (<= y -2.3e-18)
t_1
(if (<= y 3.35e-59) (* 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-18) {
tmp = t_1;
} else if (y <= 3.35e-59) {
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-18)) then
tmp = t_1
else if (y <= 3.35d-59) 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-18) {
tmp = t_1;
} else if (y <= 3.35e-59) {
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-18: tmp = t_1 elif y <= 3.35e-59: 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-18) tmp = t_1; elseif (y <= 3.35e-59) 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-18) tmp = t_1; elseif (y <= 3.35e-59) 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-18], t$95$1, If[LessEqual[y, 3.35e-59], 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^{-18}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.35 \cdot 10^{-59}:\\
\;\;\;\;x \cdot e^{0 - a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.3000000000000001e-18 or 3.35e-59 < y Initial program 98.6%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6488.0%
Simplified88.0%
if -2.3000000000000001e-18 < y < 3.35e-59Initial program 96.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6485.7%
Simplified85.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (pow z y))) (t_2 (/ x (exp (* y t)))))
(if (<= y -3.4e+168)
t_2
(if (<= y -1.02e-14)
t_1
(if (<= y 1.22e-59)
(* x (exp (- 0.0 (* a b))))
(if (<= y 7.0) t_2 t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * pow(z, y);
double t_2 = x / exp((y * t));
double tmp;
if (y <= -3.4e+168) {
tmp = t_2;
} else if (y <= -1.02e-14) {
tmp = t_1;
} else if (y <= 1.22e-59) {
tmp = x * exp((0.0 - (a * b)));
} else if (y <= 7.0) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = x * (z ** y)
t_2 = x / exp((y * t))
if (y <= (-3.4d+168)) then
tmp = t_2
else if (y <= (-1.02d-14)) then
tmp = t_1
else if (y <= 1.22d-59) then
tmp = x * exp((0.0d0 - (a * b)))
else if (y <= 7.0d0) then
tmp = t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * Math.pow(z, y);
double t_2 = x / Math.exp((y * t));
double tmp;
if (y <= -3.4e+168) {
tmp = t_2;
} else if (y <= -1.02e-14) {
tmp = t_1;
} else if (y <= 1.22e-59) {
tmp = x * Math.exp((0.0 - (a * b)));
} else if (y <= 7.0) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * math.pow(z, y) t_2 = x / math.exp((y * t)) tmp = 0 if y <= -3.4e+168: tmp = t_2 elif y <= -1.02e-14: tmp = t_1 elif y <= 1.22e-59: tmp = x * math.exp((0.0 - (a * b))) elif y <= 7.0: tmp = t_2 else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * (z ^ y)) t_2 = Float64(x / exp(Float64(y * t))) tmp = 0.0 if (y <= -3.4e+168) tmp = t_2; elseif (y <= -1.02e-14) tmp = t_1; elseif (y <= 1.22e-59) tmp = Float64(x * exp(Float64(0.0 - Float64(a * b)))); elseif (y <= 7.0) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * (z ^ y); t_2 = x / exp((y * t)); tmp = 0.0; if (y <= -3.4e+168) tmp = t_2; elseif (y <= -1.02e-14) tmp = t_1; elseif (y <= 1.22e-59) tmp = x * exp((0.0 - (a * b))); elseif (y <= 7.0) tmp = t_2; 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]}, Block[{t$95$2 = N[(x / N[Exp[N[(y * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.4e+168], t$95$2, If[LessEqual[y, -1.02e-14], t$95$1, If[LessEqual[y, 1.22e-59], N[(x * N[Exp[N[(0.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 7.0], t$95$2, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot {z}^{y}\\
t_2 := \frac{x}{e^{y \cdot t}}\\
\mathbf{if}\;y \leq -3.4 \cdot 10^{+168}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.02 \cdot 10^{-14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.22 \cdot 10^{-59}:\\
\;\;\;\;x \cdot e^{0 - a \cdot b}\\
\mathbf{elif}\;y \leq 7:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.40000000000000003e168 or 1.22e-59 < y < 7Initial program 99.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6483.3%
Simplified83.3%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6483.3%
Applied egg-rr83.3%
if -3.40000000000000003e168 < y < -1.02e-14 or 7 < y Initial program 98.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.9%
Simplified87.9%
Taylor expanded in t around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6471.9%
Simplified71.9%
if -1.02e-14 < y < 1.22e-59Initial program 96.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6485.5%
Simplified85.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (pow z y))) (t_2 (/ x (exp (* y t)))))
(if (<= y -1.85e+167)
t_2
(if (<= y -1.02e-14)
t_1
(if (<= y 5.6e-59) (/ x (exp (* a b))) (if (<= y 3.8) t_2 t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * pow(z, y);
double t_2 = x / exp((y * t));
double tmp;
if (y <= -1.85e+167) {
tmp = t_2;
} else if (y <= -1.02e-14) {
tmp = t_1;
} else if (y <= 5.6e-59) {
tmp = x / exp((a * b));
} else if (y <= 3.8) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = x * (z ** y)
t_2 = x / exp((y * t))
if (y <= (-1.85d+167)) then
tmp = t_2
else if (y <= (-1.02d-14)) then
tmp = t_1
else if (y <= 5.6d-59) then
tmp = x / exp((a * b))
else if (y <= 3.8d0) then
tmp = t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x * Math.pow(z, y);
double t_2 = x / Math.exp((y * t));
double tmp;
if (y <= -1.85e+167) {
tmp = t_2;
} else if (y <= -1.02e-14) {
tmp = t_1;
} else if (y <= 5.6e-59) {
tmp = x / Math.exp((a * b));
} else if (y <= 3.8) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x * math.pow(z, y) t_2 = x / math.exp((y * t)) tmp = 0 if y <= -1.85e+167: tmp = t_2 elif y <= -1.02e-14: tmp = t_1 elif y <= 5.6e-59: tmp = x / math.exp((a * b)) elif y <= 3.8: tmp = t_2 else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x * (z ^ y)) t_2 = Float64(x / exp(Float64(y * t))) tmp = 0.0 if (y <= -1.85e+167) tmp = t_2; elseif (y <= -1.02e-14) tmp = t_1; elseif (y <= 5.6e-59) tmp = Float64(x / exp(Float64(a * b))); elseif (y <= 3.8) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x * (z ^ y); t_2 = x / exp((y * t)); tmp = 0.0; if (y <= -1.85e+167) tmp = t_2; elseif (y <= -1.02e-14) tmp = t_1; elseif (y <= 5.6e-59) tmp = x / exp((a * b)); elseif (y <= 3.8) tmp = t_2; 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]}, Block[{t$95$2 = N[(x / N[Exp[N[(y * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.85e+167], t$95$2, If[LessEqual[y, -1.02e-14], t$95$1, If[LessEqual[y, 5.6e-59], N[(x / N[Exp[N[(a * b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.8], t$95$2, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot {z}^{y}\\
t_2 := \frac{x}{e^{y \cdot t}}\\
\mathbf{if}\;y \leq -1.85 \cdot 10^{+167}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.02 \cdot 10^{-14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 5.6 \cdot 10^{-59}:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\mathbf{elif}\;y \leq 3.8:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.85e167 or 5.59999999999999961e-59 < y < 3.7999999999999998Initial program 99.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6483.3%
Simplified83.3%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6483.3%
Applied egg-rr83.3%
if -1.85e167 < y < -1.02e-14 or 3.7999999999999998 < y Initial program 98.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.9%
Simplified87.9%
Taylor expanded in t around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6471.9%
Simplified71.9%
if -1.02e-14 < y < 5.59999999999999961e-59Initial program 96.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6485.5%
Simplified85.5%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6485.5%
Applied egg-rr85.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* x (pow z y))))
(if (<= y -1.02e-14)
t_1
(if (<= y -5.8e-145)
(/
x
(+
1.0
(*
a
(+
b
(*
a
(+
(* 0.16666666666666666 (* a (* b (* b b))))
(* (* b b) 0.5)))))))
(if (<= y 0.17)
(+
x
(*
x
(*
a
(- (* a (* (* b b) (+ 0.5 (* (* a b) -0.16666666666666666)))) 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 <= -1.02e-14) {
tmp = t_1;
} else if (y <= -5.8e-145) {
tmp = x / (1.0 + (a * (b + (a * ((0.16666666666666666 * (a * (b * (b * b)))) + ((b * b) * 0.5))))));
} else if (y <= 0.17) {
tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - 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 <= (-1.02d-14)) then
tmp = t_1
else if (y <= (-5.8d-145)) then
tmp = x / (1.0d0 + (a * (b + (a * ((0.16666666666666666d0 * (a * (b * (b * b)))) + ((b * b) * 0.5d0))))))
else if (y <= 0.17d0) then
tmp = x + (x * (a * ((a * ((b * b) * (0.5d0 + ((a * b) * (-0.16666666666666666d0))))) - 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 <= -1.02e-14) {
tmp = t_1;
} else if (y <= -5.8e-145) {
tmp = x / (1.0 + (a * (b + (a * ((0.16666666666666666 * (a * (b * (b * b)))) + ((b * b) * 0.5))))));
} else if (y <= 0.17) {
tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - 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 <= -1.02e-14: tmp = t_1 elif y <= -5.8e-145: tmp = x / (1.0 + (a * (b + (a * ((0.16666666666666666 * (a * (b * (b * b)))) + ((b * b) * 0.5)))))) elif y <= 0.17: tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - 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 <= -1.02e-14) tmp = t_1; elseif (y <= -5.8e-145) tmp = Float64(x / Float64(1.0 + Float64(a * Float64(b + Float64(a * Float64(Float64(0.16666666666666666 * Float64(a * Float64(b * Float64(b * b)))) + Float64(Float64(b * b) * 0.5))))))); elseif (y <= 0.17) tmp = Float64(x + Float64(x * Float64(a * Float64(Float64(a * Float64(Float64(b * b) * Float64(0.5 + Float64(Float64(a * b) * -0.16666666666666666)))) - 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 <= -1.02e-14) tmp = t_1; elseif (y <= -5.8e-145) tmp = x / (1.0 + (a * (b + (a * ((0.16666666666666666 * (a * (b * (b * b)))) + ((b * b) * 0.5)))))); elseif (y <= 0.17) tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - 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, -1.02e-14], t$95$1, If[LessEqual[y, -5.8e-145], N[(x / N[(1.0 + N[(a * N[(b + N[(a * N[(N[(0.16666666666666666 * N[(a * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(b * b), $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 0.17], N[(x + N[(x * N[(a * N[(N[(a * N[(N[(b * b), $MachinePrecision] * N[(0.5 + N[(N[(a * b), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot {z}^{y}\\
\mathbf{if}\;y \leq -1.02 \cdot 10^{-14}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -5.8 \cdot 10^{-145}:\\
\;\;\;\;\frac{x}{1 + a \cdot \left(b + a \cdot \left(0.16666666666666666 \cdot \left(a \cdot \left(b \cdot \left(b \cdot b\right)\right)\right) + \left(b \cdot b\right) \cdot 0.5\right)\right)}\\
\mathbf{elif}\;y \leq 0.17:\\
\;\;\;\;x + x \cdot \left(a \cdot \left(a \cdot \left(\left(b \cdot b\right) \cdot \left(0.5 + \left(a \cdot b\right) \cdot -0.16666666666666666\right)\right) - b\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.02e-14 or 0.170000000000000012 < y Initial program 98.4%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6488.4%
Simplified88.4%
Taylor expanded in t around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6469.8%
Simplified69.8%
if -1.02e-14 < y < -5.79999999999999968e-145Initial program 96.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6479.7%
Simplified79.7%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6479.7%
Applied egg-rr79.7%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.8%
Simplified69.8%
if -5.79999999999999968e-145 < y < 0.170000000000000012Initial program 97.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6481.4%
Simplified81.4%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
Simplified54.6%
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Applied egg-rr58.2%
Final simplification65.3%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* x (pow z y)))) (if (<= y -9.5e-6) t_1 (if (<= y 3.2e+25) (/ 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 <= -9.5e-6) {
tmp = t_1;
} else if (y <= 3.2e+25) {
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 <= (-9.5d-6)) then
tmp = t_1
else if (y <= 3.2d+25) 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 <= -9.5e-6) {
tmp = t_1;
} else if (y <= 3.2e+25) {
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 <= -9.5e-6: tmp = t_1 elif y <= 3.2e+25: 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 <= -9.5e-6) tmp = t_1; elseif (y <= 3.2e+25) 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 <= -9.5e-6) tmp = t_1; elseif (y <= 3.2e+25) 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, -9.5e-6], t$95$1, If[LessEqual[y, 3.2e+25], 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 -9.5 \cdot 10^{-6}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.2 \cdot 10^{+25}:\\
\;\;\;\;\frac{x}{e^{a \cdot b}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -9.5000000000000005e-6 or 3.1999999999999999e25 < y Initial program 98.3%
Taylor expanded in a around 0
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
log-lowering-log.f6490.1%
Simplified90.1%
Taylor expanded in t around 0
*-lowering-*.f64N/A
pow-lowering-pow.f6470.2%
Simplified70.2%
if -9.5000000000000005e-6 < y < 3.1999999999999999e25Initial program 97.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.1%
Simplified80.1%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6480.1%
Applied egg-rr80.1%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -1.6e+45)
(* 0.5 (* t (* t (* x (* y y)))))
(if (<= y 0.00245)
(+
x
(*
x
(* a (- (* a (* (* b b) (+ 0.5 (* (* a b) -0.16666666666666666)))) b))))
(* x (* b (* b (* 0.5 (* a a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.6e+45) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 0.00245) {
tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - b)));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-1.6d+45)) then
tmp = 0.5d0 * (t * (t * (x * (y * y))))
else if (y <= 0.00245d0) then
tmp = x + (x * (a * ((a * ((b * b) * (0.5d0 + ((a * b) * (-0.16666666666666666d0))))) - b)))
else
tmp = x * (b * (b * (0.5d0 * (a * a))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.6e+45) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 0.00245) {
tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - b)));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -1.6e+45: tmp = 0.5 * (t * (t * (x * (y * y)))) elif y <= 0.00245: tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - b))) else: tmp = x * (b * (b * (0.5 * (a * a)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -1.6e+45) tmp = Float64(0.5 * Float64(t * Float64(t * Float64(x * Float64(y * y))))); elseif (y <= 0.00245) tmp = Float64(x + Float64(x * Float64(a * Float64(Float64(a * Float64(Float64(b * b) * Float64(0.5 + Float64(Float64(a * b) * -0.16666666666666666)))) - b)))); else tmp = Float64(x * Float64(b * Float64(b * Float64(0.5 * Float64(a * a))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -1.6e+45) tmp = 0.5 * (t * (t * (x * (y * y)))); elseif (y <= 0.00245) tmp = x + (x * (a * ((a * ((b * b) * (0.5 + ((a * b) * -0.16666666666666666)))) - b))); else tmp = x * (b * (b * (0.5 * (a * a)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -1.6e+45], N[(0.5 * N[(t * N[(t * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 0.00245], N[(x + N[(x * N[(a * N[(N[(a * N[(N[(b * b), $MachinePrecision] * N[(0.5 + N[(N[(a * b), $MachinePrecision] * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(0.5 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.6 \cdot 10^{+45}:\\
\;\;\;\;0.5 \cdot \left(t \cdot \left(t \cdot \left(x \cdot \left(y \cdot y\right)\right)\right)\right)\\
\mathbf{elif}\;y \leq 0.00245:\\
\;\;\;\;x + x \cdot \left(a \cdot \left(a \cdot \left(\left(b \cdot b\right) \cdot \left(0.5 + \left(a \cdot b\right) \cdot -0.16666666666666666\right)\right) - b\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(0.5 \cdot \left(a \cdot a\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -1.6000000000000001e45Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6463.3%
Simplified63.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
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.0%
Simplified38.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.7%
Simplified58.7%
if -1.6000000000000001e45 < y < 0.0024499999999999999Initial program 97.2%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6477.8%
Simplified77.8%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
Simplified51.0%
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
+-lowering-+.f64N/A
Applied egg-rr55.0%
if 0.0024499999999999999 < y Initial program 97.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6435.3%
Simplified35.3%
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-*.f6418.9%
Simplified18.9%
*-commutativeN/A
flip3--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr0.8%
Taylor expanded in b around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.2%
Simplified51.2%
Final simplification54.7%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -3.2e+45)
(* 0.5 (* t (* t (* x (* y y)))))
(if (<= y 4.4e-105)
(* x (+ 1.0 (* a (- (* 0.5 (* a (* b b))) b))))
(* x (* b (* b (* 0.5 (* a a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -3.2e+45) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 4.4e-105) {
tmp = x * (1.0 + (a * ((0.5 * (a * (b * b))) - b)));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-3.2d+45)) then
tmp = 0.5d0 * (t * (t * (x * (y * y))))
else if (y <= 4.4d-105) then
tmp = x * (1.0d0 + (a * ((0.5d0 * (a * (b * b))) - b)))
else
tmp = x * (b * (b * (0.5d0 * (a * a))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -3.2e+45) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 4.4e-105) {
tmp = x * (1.0 + (a * ((0.5 * (a * (b * b))) - b)));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -3.2e+45: tmp = 0.5 * (t * (t * (x * (y * y)))) elif y <= 4.4e-105: tmp = x * (1.0 + (a * ((0.5 * (a * (b * b))) - b))) else: tmp = x * (b * (b * (0.5 * (a * a)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -3.2e+45) tmp = Float64(0.5 * Float64(t * Float64(t * Float64(x * Float64(y * y))))); elseif (y <= 4.4e-105) tmp = Float64(x * Float64(1.0 + Float64(a * Float64(Float64(0.5 * Float64(a * Float64(b * b))) - b)))); else tmp = Float64(x * Float64(b * Float64(b * Float64(0.5 * Float64(a * a))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -3.2e+45) tmp = 0.5 * (t * (t * (x * (y * y)))); elseif (y <= 4.4e-105) tmp = x * (1.0 + (a * ((0.5 * (a * (b * b))) - b))); else tmp = x * (b * (b * (0.5 * (a * a)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -3.2e+45], N[(0.5 * N[(t * N[(t * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.4e-105], N[(x * N[(1.0 + N[(a * N[(N[(0.5 * N[(a * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(0.5 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.2 \cdot 10^{+45}:\\
\;\;\;\;0.5 \cdot \left(t \cdot \left(t \cdot \left(x \cdot \left(y \cdot y\right)\right)\right)\right)\\
\mathbf{elif}\;y \leq 4.4 \cdot 10^{-105}:\\
\;\;\;\;x \cdot \left(1 + a \cdot \left(0.5 \cdot \left(a \cdot \left(b \cdot b\right)\right) - b\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(0.5 \cdot \left(a \cdot a\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -3.2000000000000003e45Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6463.3%
Simplified63.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
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.0%
Simplified38.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.7%
Simplified58.7%
if -3.2000000000000003e45 < y < 4.40000000000000008e-105Initial program 96.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.7%
Simplified80.7%
Taylor expanded in a around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.8%
Simplified56.8%
if 4.40000000000000008e-105 < y Initial program 97.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6440.7%
Simplified40.7%
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-*.f6421.1%
Simplified21.1%
*-commutativeN/A
flip3--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr3.2%
Taylor expanded in b around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.2%
Simplified47.2%
Final simplification54.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -1.35e+43)
(* 0.5 (* t (* t (* x (* y y)))))
(if (<= y 4.3e-105)
(* x (+ 1.0 (* b (* 0.5 (* b (* a a))))))
(* x (* b (* b (* 0.5 (* a a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.35e+43) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 4.3e-105) {
tmp = x * (1.0 + (b * (0.5 * (b * (a * a)))));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-1.35d+43)) then
tmp = 0.5d0 * (t * (t * (x * (y * y))))
else if (y <= 4.3d-105) then
tmp = x * (1.0d0 + (b * (0.5d0 * (b * (a * a)))))
else
tmp = x * (b * (b * (0.5d0 * (a * a))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.35e+43) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 4.3e-105) {
tmp = x * (1.0 + (b * (0.5 * (b * (a * a)))));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -1.35e+43: tmp = 0.5 * (t * (t * (x * (y * y)))) elif y <= 4.3e-105: tmp = x * (1.0 + (b * (0.5 * (b * (a * a))))) else: tmp = x * (b * (b * (0.5 * (a * a)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -1.35e+43) tmp = Float64(0.5 * Float64(t * Float64(t * Float64(x * Float64(y * y))))); elseif (y <= 4.3e-105) tmp = Float64(x * Float64(1.0 + Float64(b * Float64(0.5 * Float64(b * Float64(a * a)))))); else tmp = Float64(x * Float64(b * Float64(b * Float64(0.5 * Float64(a * a))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -1.35e+43) tmp = 0.5 * (t * (t * (x * (y * y)))); elseif (y <= 4.3e-105) tmp = x * (1.0 + (b * (0.5 * (b * (a * a))))); else tmp = x * (b * (b * (0.5 * (a * a)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -1.35e+43], N[(0.5 * N[(t * N[(t * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.3e-105], N[(x * N[(1.0 + N[(b * N[(0.5 * N[(b * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(0.5 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.35 \cdot 10^{+43}:\\
\;\;\;\;0.5 \cdot \left(t \cdot \left(t \cdot \left(x \cdot \left(y \cdot y\right)\right)\right)\right)\\
\mathbf{elif}\;y \leq 4.3 \cdot 10^{-105}:\\
\;\;\;\;x \cdot \left(1 + b \cdot \left(0.5 \cdot \left(b \cdot \left(a \cdot a\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(0.5 \cdot \left(a \cdot a\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -1.3500000000000001e43Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6463.3%
Simplified63.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
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.0%
Simplified38.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6458.7%
Simplified58.7%
if -1.3500000000000001e43 < y < 4.29999999999999964e-105Initial program 96.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.7%
Simplified80.7%
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-*.f6454.5%
Simplified54.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.5%
Simplified53.5%
if 4.29999999999999964e-105 < y Initial program 97.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6440.7%
Simplified40.7%
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-*.f6421.1%
Simplified21.1%
*-commutativeN/A
flip3--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr3.2%
Taylor expanded in b around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.2%
Simplified47.2%
Final simplification52.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -9.6e+86)
(* x (- 1.0 (* y t)))
(if (<= y -4.9e-17)
(* b (- (/ x b) (* x a)))
(if (<= y 8e+15) (/ x (+ 1.0 (* a b))) (- 0.0 (* x (* a b)))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -9.6e+86) {
tmp = x * (1.0 - (y * t));
} else if (y <= -4.9e-17) {
tmp = b * ((x / b) - (x * a));
} else if (y <= 8e+15) {
tmp = x / (1.0 + (a * b));
} else {
tmp = 0.0 - (x * (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 <= (-9.6d+86)) then
tmp = x * (1.0d0 - (y * t))
else if (y <= (-4.9d-17)) then
tmp = b * ((x / b) - (x * a))
else if (y <= 8d+15) then
tmp = x / (1.0d0 + (a * b))
else
tmp = 0.0d0 - (x * (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 <= -9.6e+86) {
tmp = x * (1.0 - (y * t));
} else if (y <= -4.9e-17) {
tmp = b * ((x / b) - (x * a));
} else if (y <= 8e+15) {
tmp = x / (1.0 + (a * b));
} else {
tmp = 0.0 - (x * (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -9.6e+86: tmp = x * (1.0 - (y * t)) elif y <= -4.9e-17: tmp = b * ((x / b) - (x * a)) elif y <= 8e+15: tmp = x / (1.0 + (a * b)) else: tmp = 0.0 - (x * (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -9.6e+86) tmp = Float64(x * Float64(1.0 - Float64(y * t))); elseif (y <= -4.9e-17) tmp = Float64(b * Float64(Float64(x / b) - Float64(x * a))); elseif (y <= 8e+15) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(0.0 - Float64(x * Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -9.6e+86) tmp = x * (1.0 - (y * t)); elseif (y <= -4.9e-17) tmp = b * ((x / b) - (x * a)); elseif (y <= 8e+15) tmp = x / (1.0 + (a * b)); else tmp = 0.0 - (x * (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -9.6e+86], N[(x * N[(1.0 - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.9e-17], N[(b * N[(N[(x / b), $MachinePrecision] - N[(x * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8e+15], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(x * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.6 \cdot 10^{+86}:\\
\;\;\;\;x \cdot \left(1 - y \cdot t\right)\\
\mathbf{elif}\;y \leq -4.9 \cdot 10^{-17}:\\
\;\;\;\;b \cdot \left(\frac{x}{b} - x \cdot a\right)\\
\mathbf{elif}\;y \leq 8 \cdot 10^{+15}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;0 - x \cdot \left(a \cdot b\right)\\
\end{array}
\end{array}
if y < -9.6000000000000001e86Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6469.3%
Simplified69.3%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6428.9%
Simplified28.9%
if -9.6000000000000001e86 < y < -4.90000000000000012e-17Initial program 100.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6436.5%
Simplified36.5%
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-*.f6417.9%
Simplified17.9%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6413.2%
Simplified13.2%
Taylor expanded in b around inf
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6435.3%
Simplified35.3%
if -4.90000000000000012e-17 < y < 8e15Initial program 97.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6481.5%
Simplified81.5%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6481.5%
Applied egg-rr81.5%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6449.2%
Simplified49.2%
if 8e15 < y Initial program 96.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6432.7%
Simplified32.7%
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-*.f6418.4%
Simplified18.4%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6415.4%
Simplified15.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6434.0%
Simplified34.0%
Final simplification41.1%
(FPCore (x y z t a b)
:precision binary64
(if (<= y -7.5e-15)
(* 0.5 (* t (* t (* x (* y y)))))
(if (<= y 1.05e-104)
(/ x (+ 1.0 (* a b)))
(* x (* b (* b (* 0.5 (* a a))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -7.5e-15) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 1.05e-104) {
tmp = x / (1.0 + (a * b));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8) :: tmp
if (y <= (-7.5d-15)) then
tmp = 0.5d0 * (t * (t * (x * (y * y))))
else if (y <= 1.05d-104) then
tmp = x / (1.0d0 + (a * b))
else
tmp = x * (b * (b * (0.5d0 * (a * a))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -7.5e-15) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 1.05e-104) {
tmp = x / (1.0 + (a * b));
} else {
tmp = x * (b * (b * (0.5 * (a * a))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -7.5e-15: tmp = 0.5 * (t * (t * (x * (y * y)))) elif y <= 1.05e-104: tmp = x / (1.0 + (a * b)) else: tmp = x * (b * (b * (0.5 * (a * a)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -7.5e-15) tmp = Float64(0.5 * Float64(t * Float64(t * Float64(x * Float64(y * y))))); elseif (y <= 1.05e-104) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(x * Float64(b * Float64(b * Float64(0.5 * Float64(a * a))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -7.5e-15) tmp = 0.5 * (t * (t * (x * (y * y)))); elseif (y <= 1.05e-104) tmp = x / (1.0 + (a * b)); else tmp = x * (b * (b * (0.5 * (a * a)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -7.5e-15], N[(0.5 * N[(t * N[(t * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.05e-104], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(b * N[(b * N[(0.5 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -7.5 \cdot 10^{-15}:\\
\;\;\;\;0.5 \cdot \left(t \cdot \left(t \cdot \left(x \cdot \left(y \cdot y\right)\right)\right)\right)\\
\mathbf{elif}\;y \leq 1.05 \cdot 10^{-104}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(b \cdot \left(b \cdot \left(0.5 \cdot \left(a \cdot a\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -7.4999999999999996e-15Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6459.6%
Simplified59.6%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.0%
Simplified36.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.1%
Simplified51.1%
if -7.4999999999999996e-15 < y < 1.04999999999999999e-104Initial program 96.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6485.6%
Simplified85.6%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6485.6%
Applied egg-rr85.6%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6453.7%
Simplified53.7%
if 1.04999999999999999e-104 < y Initial program 97.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6440.0%
Simplified40.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
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6421.4%
Simplified21.4%
*-commutativeN/A
flip3--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr3.3%
Taylor expanded in b around inf
associate-*r*N/A
unpow2N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.7%
Simplified47.7%
Final simplification51.1%
(FPCore (x y z t a b) :precision binary64 (if (<= y -9.2e-18) (* 0.5 (* t (* t (* x (* y y))))) (if (<= y 7.5e-7) (/ x (+ 1.0 (* a b))) (* 0.5 (* a (* a (* x (* b b))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -9.2e-18) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 7.5e-7) {
tmp = x / (1.0 + (a * b));
} else {
tmp = 0.5 * (a * (a * (x * (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 <= (-9.2d-18)) then
tmp = 0.5d0 * (t * (t * (x * (y * y))))
else if (y <= 7.5d-7) then
tmp = x / (1.0d0 + (a * b))
else
tmp = 0.5d0 * (a * (a * (x * (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 <= -9.2e-18) {
tmp = 0.5 * (t * (t * (x * (y * y))));
} else if (y <= 7.5e-7) {
tmp = x / (1.0 + (a * b));
} else {
tmp = 0.5 * (a * (a * (x * (b * b))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -9.2e-18: tmp = 0.5 * (t * (t * (x * (y * y)))) elif y <= 7.5e-7: tmp = x / (1.0 + (a * b)) else: tmp = 0.5 * (a * (a * (x * (b * b)))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -9.2e-18) tmp = Float64(0.5 * Float64(t * Float64(t * Float64(x * Float64(y * y))))); elseif (y <= 7.5e-7) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(0.5 * Float64(a * Float64(a * Float64(x * Float64(b * b))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -9.2e-18) tmp = 0.5 * (t * (t * (x * (y * y)))); elseif (y <= 7.5e-7) tmp = x / (1.0 + (a * b)); else tmp = 0.5 * (a * (a * (x * (b * b)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -9.2e-18], N[(0.5 * N[(t * N[(t * N[(x * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 7.5e-7], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.5 * N[(a * N[(a * N[(x * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.2 \cdot 10^{-18}:\\
\;\;\;\;0.5 \cdot \left(t \cdot \left(t \cdot \left(x \cdot \left(y \cdot y\right)\right)\right)\right)\\
\mathbf{elif}\;y \leq 7.5 \cdot 10^{-7}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(a \cdot \left(a \cdot \left(x \cdot \left(b \cdot b\right)\right)\right)\right)\\
\end{array}
\end{array}
if y < -9.2000000000000004e-18Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6459.6%
Simplified59.6%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
mul-1-negN/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6436.0%
Simplified36.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.1%
Simplified51.1%
if -9.2000000000000004e-18 < y < 7.5000000000000002e-7Initial program 96.9%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6482.3%
Simplified82.3%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6482.3%
Applied egg-rr82.3%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6451.0%
Simplified51.0%
if 7.5000000000000002e-7 < y Initial program 97.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6434.8%
Simplified34.8%
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-*.f6418.7%
Simplified18.7%
*-commutativeN/A
flip3--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr0.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.4%
Simplified50.4%
Final simplification50.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* 0.5 (* a (* a (* x (* b b)))))))
(if (<= y -22000000000000.0)
t_1
(if (<= y 3.5e-6) (/ 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.5 * (a * (a * (x * (b * b))));
double tmp;
if (y <= -22000000000000.0) {
tmp = t_1;
} else if (y <= 3.5e-6) {
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.5d0 * (a * (a * (x * (b * b))))
if (y <= (-22000000000000.0d0)) then
tmp = t_1
else if (y <= 3.5d-6) 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.5 * (a * (a * (x * (b * b))));
double tmp;
if (y <= -22000000000000.0) {
tmp = t_1;
} else if (y <= 3.5e-6) {
tmp = x / (1.0 + (a * b));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = 0.5 * (a * (a * (x * (b * b)))) tmp = 0 if y <= -22000000000000.0: tmp = t_1 elif y <= 3.5e-6: tmp = x / (1.0 + (a * b)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(0.5 * Float64(a * Float64(a * Float64(x * Float64(b * b))))) tmp = 0.0 if (y <= -22000000000000.0) tmp = t_1; elseif (y <= 3.5e-6) 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.5 * (a * (a * (x * (b * b)))); tmp = 0.0; if (y <= -22000000000000.0) tmp = t_1; elseif (y <= 3.5e-6) 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.5 * N[(a * N[(a * N[(x * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -22000000000000.0], t$95$1, If[LessEqual[y, 3.5e-6], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.5 \cdot \left(a \cdot \left(a \cdot \left(x \cdot \left(b \cdot b\right)\right)\right)\right)\\
\mathbf{if}\;y \leq -22000000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{-6}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.2e13 or 3.49999999999999995e-6 < y Initial program 98.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6436.5%
Simplified36.5%
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-*.f6420.9%
Simplified20.9%
*-commutativeN/A
flip3--N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr5.7%
Taylor expanded in b around inf
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6443.0%
Simplified43.0%
if -2.2e13 < y < 3.49999999999999995e-6Initial program 97.0%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6480.7%
Simplified80.7%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6480.6%
Applied egg-rr80.6%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6449.6%
Simplified49.6%
Final simplification46.4%
(FPCore (x y z t a b) :precision binary64 (if (<= y -1.7e+206) (* x (- 1.0 (* y t))) (if (<= y 9e+15) (/ x (+ 1.0 (* a b))) (- 0.0 (* x (* a b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -1.7e+206) {
tmp = x * (1.0 - (y * t));
} else if (y <= 9e+15) {
tmp = x / (1.0 + (a * b));
} else {
tmp = 0.0 - (x * (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.7d+206)) then
tmp = x * (1.0d0 - (y * t))
else if (y <= 9d+15) then
tmp = x / (1.0d0 + (a * b))
else
tmp = 0.0d0 - (x * (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.7e+206) {
tmp = x * (1.0 - (y * t));
} else if (y <= 9e+15) {
tmp = x / (1.0 + (a * b));
} else {
tmp = 0.0 - (x * (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -1.7e+206: tmp = x * (1.0 - (y * t)) elif y <= 9e+15: tmp = x / (1.0 + (a * b)) else: tmp = 0.0 - (x * (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -1.7e+206) tmp = Float64(x * Float64(1.0 - Float64(y * t))); elseif (y <= 9e+15) tmp = Float64(x / Float64(1.0 + Float64(a * b))); else tmp = Float64(0.0 - Float64(x * Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -1.7e+206) tmp = x * (1.0 - (y * t)); elseif (y <= 9e+15) tmp = x / (1.0 + (a * b)); else tmp = 0.0 - (x * (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -1.7e+206], N[(x * N[(1.0 - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 9e+15], N[(x / N[(1.0 + N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(x * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.7 \cdot 10^{+206}:\\
\;\;\;\;x \cdot \left(1 - y \cdot t\right)\\
\mathbf{elif}\;y \leq 9 \cdot 10^{+15}:\\
\;\;\;\;\frac{x}{1 + a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;0 - x \cdot \left(a \cdot b\right)\\
\end{array}
\end{array}
if y < -1.69999999999999999e206Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6486.2%
Simplified86.2%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6439.6%
Simplified39.6%
if -1.69999999999999999e206 < y < 9e15Initial program 97.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6472.5%
Simplified72.5%
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6472.5%
Applied egg-rr72.5%
Taylor expanded in a around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6441.2%
Simplified41.2%
if 9e15 < y Initial program 96.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6432.7%
Simplified32.7%
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-*.f6418.4%
Simplified18.4%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6415.4%
Simplified15.4%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6434.0%
Simplified34.0%
Final simplification39.3%
(FPCore (x y z t a b) :precision binary64 (if (<= y -8e+68) (* x (- 1.0 (* y t))) (if (<= y 1.8e-27) (* x (- 1.0 (* a b))) (- 0.0 (* x (* a b))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -8e+68) {
tmp = x * (1.0 - (y * t));
} else if (y <= 1.8e-27) {
tmp = x * (1.0 - (a * b));
} else {
tmp = 0.0 - (x * (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 <= (-8d+68)) then
tmp = x * (1.0d0 - (y * t))
else if (y <= 1.8d-27) then
tmp = x * (1.0d0 - (a * b))
else
tmp = 0.0d0 - (x * (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 <= -8e+68) {
tmp = x * (1.0 - (y * t));
} else if (y <= 1.8e-27) {
tmp = x * (1.0 - (a * b));
} else {
tmp = 0.0 - (x * (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -8e+68: tmp = x * (1.0 - (y * t)) elif y <= 1.8e-27: tmp = x * (1.0 - (a * b)) else: tmp = 0.0 - (x * (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -8e+68) tmp = Float64(x * Float64(1.0 - Float64(y * t))); elseif (y <= 1.8e-27) tmp = Float64(x * Float64(1.0 - Float64(a * b))); else tmp = Float64(0.0 - Float64(x * Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -8e+68) tmp = x * (1.0 - (y * t)); elseif (y <= 1.8e-27) tmp = x * (1.0 - (a * b)); else tmp = 0.0 - (x * (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -8e+68], N[(x * N[(1.0 - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.8e-27], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(x * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8 \cdot 10^{+68}:\\
\;\;\;\;x \cdot \left(1 - y \cdot t\right)\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{-27}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;0 - x \cdot \left(a \cdot b\right)\\
\end{array}
\end{array}
if y < -7.99999999999999962e68Initial program 100.0%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6467.1%
Simplified67.1%
Taylor expanded in y around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6427.5%
Simplified27.5%
if -7.99999999999999962e68 < y < 1.7999999999999999e-27Initial program 97.2%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6478.7%
Simplified78.7%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6444.5%
Simplified44.5%
if 1.7999999999999999e-27 < y Initial program 97.2%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6434.4%
Simplified34.4%
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-*.f6417.8%
Simplified17.8%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6413.8%
Simplified13.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6431.3%
Simplified31.3%
Final simplification37.9%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (- 0.0 (* x (* a b))))) (if (<= y -22000000000000.0) t_1 (if (<= y 1.08e-100) x t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = 0.0 - (x * (a * b));
double tmp;
if (y <= -22000000000000.0) {
tmp = t_1;
} else if (y <= 1.08e-100) {
tmp = x;
} 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 - (x * (a * b))
if (y <= (-22000000000000.0d0)) then
tmp = t_1
else if (y <= 1.08d-100) then
tmp = x
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 - (x * (a * b));
double tmp;
if (y <= -22000000000000.0) {
tmp = t_1;
} else if (y <= 1.08e-100) {
tmp = x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = 0.0 - (x * (a * b)) tmp = 0 if y <= -22000000000000.0: tmp = t_1 elif y <= 1.08e-100: tmp = x else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(0.0 - Float64(x * Float64(a * b))) tmp = 0.0 if (y <= -22000000000000.0) tmp = t_1; elseif (y <= 1.08e-100) tmp = x; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = 0.0 - (x * (a * b)); tmp = 0.0; if (y <= -22000000000000.0) tmp = t_1; elseif (y <= 1.08e-100) tmp = x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(0.0 - N[(x * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -22000000000000.0], t$95$1, If[LessEqual[y, 1.08e-100], x, t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0 - x \cdot \left(a \cdot b\right)\\
\mathbf{if}\;y \leq -22000000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.08 \cdot 10^{-100}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.2e13 or 1.0800000000000001e-100 < y Initial program 98.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6439.6%
Simplified39.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-*.f6422.4%
Simplified22.4%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6416.1%
Simplified16.1%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6425.8%
Simplified25.8%
if -2.2e13 < y < 1.0800000000000001e-100Initial program 96.6%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6483.0%
Simplified83.0%
Taylor expanded in a around 0
Simplified38.7%
Final simplification31.6%
(FPCore (x y z t a b) :precision binary64 (if (<= y 6.6e-29) (* x (- 1.0 (* a b))) (- 0.0 (* x (* a b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= 6.6e-29) {
tmp = x * (1.0 - (a * b));
} else {
tmp = 0.0 - (x * (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.6d-29) then
tmp = x * (1.0d0 - (a * b))
else
tmp = 0.0d0 - (x * (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.6e-29) {
tmp = x * (1.0 - (a * b));
} else {
tmp = 0.0 - (x * (a * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= 6.6e-29: tmp = x * (1.0 - (a * b)) else: tmp = 0.0 - (x * (a * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= 6.6e-29) tmp = Float64(x * Float64(1.0 - Float64(a * b))); else tmp = Float64(0.0 - Float64(x * Float64(a * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= 6.6e-29) tmp = x * (1.0 - (a * b)); else tmp = 0.0 - (x * (a * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, 6.6e-29], N[(x * N[(1.0 - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 - N[(x * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 6.6 \cdot 10^{-29}:\\
\;\;\;\;x \cdot \left(1 - a \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;0 - x \cdot \left(a \cdot b\right)\\
\end{array}
\end{array}
if y < 6.60000000000000055e-29Initial program 97.9%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6468.6%
Simplified68.6%
Taylor expanded in a around 0
neg-mul-1N/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f6437.2%
Simplified37.2%
if 6.60000000000000055e-29 < y Initial program 97.2%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6434.4%
Simplified34.4%
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-*.f6417.8%
Simplified17.8%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6413.8%
Simplified13.8%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6431.3%
Simplified31.3%
Final simplification35.5%
(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 97.7%
Taylor expanded in b around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6459.1%
Simplified59.1%
Taylor expanded in a around 0
Simplified20.1%
herbie shell --seed 2024160
(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))))))