
(FPCore (x y z t a b c i) :precision binary64 (+ (+ (+ (+ (+ (* x (log y)) z) t) a) (* (- b 0.5) (log c))) (* y i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((((x * log(y)) + z) + t) + a) + ((b - 0.5) * log(c))) + (y * i);
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
code = (((((x * log(y)) + z) + t) + a) + ((b - 0.5d0) * log(c))) + (y * i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((((x * Math.log(y)) + z) + t) + a) + ((b - 0.5) * Math.log(c))) + (y * i);
}
def code(x, y, z, t, a, b, c, i): return (((((x * math.log(y)) + z) + t) + a) + ((b - 0.5) * math.log(c))) + (y * i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(x * log(y)) + z) + t) + a) + Float64(Float64(b - 0.5) * log(c))) + Float64(y * i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = (((((x * log(y)) + z) + t) + a) + ((b - 0.5) * log(c))) + (y * i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + t), $MachinePrecision] + a), $MachinePrecision] + N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(x \cdot \log y + z\right) + t\right) + a\right) + \left(b - 0.5\right) \cdot \log c\right) + y \cdot i
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c i) :precision binary64 (+ (+ (+ (+ (+ (* x (log y)) z) t) a) (* (- b 0.5) (log c))) (* y i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((((x * log(y)) + z) + t) + a) + ((b - 0.5) * log(c))) + (y * i);
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
code = (((((x * log(y)) + z) + t) + a) + ((b - 0.5d0) * log(c))) + (y * i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((((x * Math.log(y)) + z) + t) + a) + ((b - 0.5) * Math.log(c))) + (y * i);
}
def code(x, y, z, t, a, b, c, i): return (((((x * math.log(y)) + z) + t) + a) + ((b - 0.5) * math.log(c))) + (y * i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(x * log(y)) + z) + t) + a) + Float64(Float64(b - 0.5) * log(c))) + Float64(y * i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = (((((x * log(y)) + z) + t) + a) + ((b - 0.5) * log(c))) + (y * i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + t), $MachinePrecision] + a), $MachinePrecision] + N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(x \cdot \log y + z\right) + t\right) + a\right) + \left(b - 0.5\right) \cdot \log c\right) + y \cdot i
\end{array}
(FPCore (x y z t a b c i) :precision binary64 (+ (+ (+ (+ (+ (* x (log y)) z) t) a) (* (- b 0.5) (log c))) (* y i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((((x * log(y)) + z) + t) + a) + ((b - 0.5) * log(c))) + (y * i);
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
code = (((((x * log(y)) + z) + t) + a) + ((b - 0.5d0) * log(c))) + (y * i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((((x * Math.log(y)) + z) + t) + a) + ((b - 0.5) * Math.log(c))) + (y * i);
}
def code(x, y, z, t, a, b, c, i): return (((((x * math.log(y)) + z) + t) + a) + ((b - 0.5) * math.log(c))) + (y * i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(x * log(y)) + z) + t) + a) + Float64(Float64(b - 0.5) * log(c))) + Float64(y * i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = (((((x * log(y)) + z) + t) + a) + ((b - 0.5) * log(c))) + (y * i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + t), $MachinePrecision] + a), $MachinePrecision] + N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(x \cdot \log y + z\right) + t\right) + a\right) + \left(b - 0.5\right) \cdot \log c\right) + y \cdot i
\end{array}
Initial program 99.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* (- b 0.5) (log c))) (t_2 (+ (* y i) (+ z t_1))))
(if (<= t_1 -2e+95)
t_2
(if (<= t_1 260.0)
(+ (+ (* x (log y)) z) (+ t a))
(if (<= t_1 5e+90) (+ (* y i) (+ (* (log c) -0.5) (+ z a))) t_2)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (b - 0.5) * log(c);
double t_2 = (y * i) + (z + t_1);
double tmp;
if (t_1 <= -2e+95) {
tmp = t_2;
} else if (t_1 <= 260.0) {
tmp = ((x * log(y)) + z) + (t + a);
} else if (t_1 <= 5e+90) {
tmp = (y * i) + ((log(c) * -0.5) + (z + a));
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (b - 0.5d0) * log(c)
t_2 = (y * i) + (z + t_1)
if (t_1 <= (-2d+95)) then
tmp = t_2
else if (t_1 <= 260.0d0) then
tmp = ((x * log(y)) + z) + (t + a)
else if (t_1 <= 5d+90) then
tmp = (y * i) + ((log(c) * (-0.5d0)) + (z + a))
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (b - 0.5) * Math.log(c);
double t_2 = (y * i) + (z + t_1);
double tmp;
if (t_1 <= -2e+95) {
tmp = t_2;
} else if (t_1 <= 260.0) {
tmp = ((x * Math.log(y)) + z) + (t + a);
} else if (t_1 <= 5e+90) {
tmp = (y * i) + ((Math.log(c) * -0.5) + (z + a));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (b - 0.5) * math.log(c) t_2 = (y * i) + (z + t_1) tmp = 0 if t_1 <= -2e+95: tmp = t_2 elif t_1 <= 260.0: tmp = ((x * math.log(y)) + z) + (t + a) elif t_1 <= 5e+90: tmp = (y * i) + ((math.log(c) * -0.5) + (z + a)) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(b - 0.5) * log(c)) t_2 = Float64(Float64(y * i) + Float64(z + t_1)) tmp = 0.0 if (t_1 <= -2e+95) tmp = t_2; elseif (t_1 <= 260.0) tmp = Float64(Float64(Float64(x * log(y)) + z) + Float64(t + a)); elseif (t_1 <= 5e+90) tmp = Float64(Float64(y * i) + Float64(Float64(log(c) * -0.5) + Float64(z + a))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (b - 0.5) * log(c); t_2 = (y * i) + (z + t_1); tmp = 0.0; if (t_1 <= -2e+95) tmp = t_2; elseif (t_1 <= 260.0) tmp = ((x * log(y)) + z) + (t + a); elseif (t_1 <= 5e+90) tmp = (y * i) + ((log(c) * -0.5) + (z + a)); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(y * i), $MachinePrecision] + N[(z + t$95$1), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e+95], t$95$2, If[LessEqual[t$95$1, 260.0], N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + N[(t + a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 5e+90], N[(N[(y * i), $MachinePrecision] + N[(N[(N[Log[c], $MachinePrecision] * -0.5), $MachinePrecision] + N[(z + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - 0.5\right) \cdot \log c\\
t_2 := y \cdot i + \left(z + t\_1\right)\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+95}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 260:\\
\;\;\;\;\left(x \cdot \log y + z\right) + \left(t + a\right)\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+90}:\\
\;\;\;\;y \cdot i + \left(\log c \cdot -0.5 + \left(z + a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (*.f64 (-.f64 b #s(literal 1/2 binary64)) (log.f64 c)) < -2.00000000000000004e95 or 5.0000000000000004e90 < (*.f64 (-.f64 b #s(literal 1/2 binary64)) (log.f64 c)) Initial program 99.8%
Taylor expanded in z around inf
Simplified73.7%
if -2.00000000000000004e95 < (*.f64 (-.f64 b #s(literal 1/2 binary64)) (log.f64 c)) < 260Initial program 99.9%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified83.8%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6479.9%
Simplified79.9%
Taylor expanded in b around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f6479.0%
Simplified79.0%
if 260 < (*.f64 (-.f64 b #s(literal 1/2 binary64)) (log.f64 c)) < 5.0000000000000004e90Initial program 99.9%
Taylor expanded in z around inf
Simplified75.0%
Taylor expanded in b around 0
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6467.2%
Simplified67.2%
Final simplification75.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* (- b 0.5) (log c))))
(if (<= i -126.0)
(+ (* y i) (+ t_1 (+ a (+ z t))))
(if (<= i 2.2e-23)
(+ (* x (log y)) (+ (+ z (+ t a)) (* (log c) (+ b -0.5))))
(+ (* y i) (+ t_1 (+ z a)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (b - 0.5) * log(c);
double tmp;
if (i <= -126.0) {
tmp = (y * i) + (t_1 + (a + (z + t)));
} else if (i <= 2.2e-23) {
tmp = (x * log(y)) + ((z + (t + a)) + (log(c) * (b + -0.5)));
} else {
tmp = (y * i) + (t_1 + (z + a));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = (b - 0.5d0) * log(c)
if (i <= (-126.0d0)) then
tmp = (y * i) + (t_1 + (a + (z + t)))
else if (i <= 2.2d-23) then
tmp = (x * log(y)) + ((z + (t + a)) + (log(c) * (b + (-0.5d0))))
else
tmp = (y * i) + (t_1 + (z + a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (b - 0.5) * Math.log(c);
double tmp;
if (i <= -126.0) {
tmp = (y * i) + (t_1 + (a + (z + t)));
} else if (i <= 2.2e-23) {
tmp = (x * Math.log(y)) + ((z + (t + a)) + (Math.log(c) * (b + -0.5)));
} else {
tmp = (y * i) + (t_1 + (z + a));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (b - 0.5) * math.log(c) tmp = 0 if i <= -126.0: tmp = (y * i) + (t_1 + (a + (z + t))) elif i <= 2.2e-23: tmp = (x * math.log(y)) + ((z + (t + a)) + (math.log(c) * (b + -0.5))) else: tmp = (y * i) + (t_1 + (z + a)) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(b - 0.5) * log(c)) tmp = 0.0 if (i <= -126.0) tmp = Float64(Float64(y * i) + Float64(t_1 + Float64(a + Float64(z + t)))); elseif (i <= 2.2e-23) tmp = Float64(Float64(x * log(y)) + Float64(Float64(z + Float64(t + a)) + Float64(log(c) * Float64(b + -0.5)))); else tmp = Float64(Float64(y * i) + Float64(t_1 + Float64(z + a))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (b - 0.5) * log(c); tmp = 0.0; if (i <= -126.0) tmp = (y * i) + (t_1 + (a + (z + t))); elseif (i <= 2.2e-23) tmp = (x * log(y)) + ((z + (t + a)) + (log(c) * (b + -0.5))); else tmp = (y * i) + (t_1 + (z + a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[i, -126.0], N[(N[(y * i), $MachinePrecision] + N[(t$95$1 + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 2.2e-23], N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(N[(z + N[(t + a), $MachinePrecision]), $MachinePrecision] + N[(N[Log[c], $MachinePrecision] * N[(b + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(t$95$1 + N[(z + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - 0.5\right) \cdot \log c\\
\mathbf{if}\;i \leq -126:\\
\;\;\;\;y \cdot i + \left(t\_1 + \left(a + \left(z + t\right)\right)\right)\\
\mathbf{elif}\;i \leq 2.2 \cdot 10^{-23}:\\
\;\;\;\;x \cdot \log y + \left(\left(z + \left(t + a\right)\right) + \log c \cdot \left(b + -0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(t\_1 + \left(z + a\right)\right)\\
\end{array}
\end{array}
if i < -126Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f6490.4%
Simplified90.4%
if -126 < i < 2.1999999999999999e-23Initial program 99.8%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified98.9%
if 2.1999999999999999e-23 < i Initial program 100.0%
Taylor expanded in z around inf
Simplified75.3%
Final simplification90.0%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= x -1.75e+152)
(+ (* x (log y)) (+ (+ z (+ t a)) (* b (log c))))
(if (<= x 7.2e+166)
(+ (* y i) (+ (* (- b 0.5) (log c)) (+ a (+ z t))))
(* x (+ (log y) (/ (+ z (* y i)) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -1.75e+152) {
tmp = (x * log(y)) + ((z + (t + a)) + (b * log(c)));
} else if (x <= 7.2e+166) {
tmp = (y * i) + (((b - 0.5) * log(c)) + (a + (z + t)));
} else {
tmp = x * (log(y) + ((z + (y * i)) / x));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (x <= (-1.75d+152)) then
tmp = (x * log(y)) + ((z + (t + a)) + (b * log(c)))
else if (x <= 7.2d+166) then
tmp = (y * i) + (((b - 0.5d0) * log(c)) + (a + (z + t)))
else
tmp = x * (log(y) + ((z + (y * i)) / x))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -1.75e+152) {
tmp = (x * Math.log(y)) + ((z + (t + a)) + (b * Math.log(c)));
} else if (x <= 7.2e+166) {
tmp = (y * i) + (((b - 0.5) * Math.log(c)) + (a + (z + t)));
} else {
tmp = x * (Math.log(y) + ((z + (y * i)) / x));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if x <= -1.75e+152: tmp = (x * math.log(y)) + ((z + (t + a)) + (b * math.log(c))) elif x <= 7.2e+166: tmp = (y * i) + (((b - 0.5) * math.log(c)) + (a + (z + t))) else: tmp = x * (math.log(y) + ((z + (y * i)) / x)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (x <= -1.75e+152) tmp = Float64(Float64(x * log(y)) + Float64(Float64(z + Float64(t + a)) + Float64(b * log(c)))); elseif (x <= 7.2e+166) tmp = Float64(Float64(y * i) + Float64(Float64(Float64(b - 0.5) * log(c)) + Float64(a + Float64(z + t)))); else tmp = Float64(x * Float64(log(y) + Float64(Float64(z + Float64(y * i)) / x))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (x <= -1.75e+152) tmp = (x * log(y)) + ((z + (t + a)) + (b * log(c))); elseif (x <= 7.2e+166) tmp = (y * i) + (((b - 0.5) * log(c)) + (a + (z + t))); else tmp = x * (log(y) + ((z + (y * i)) / x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[x, -1.75e+152], N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(N[(z + N[(t + a), $MachinePrecision]), $MachinePrecision] + N[(b * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 7.2e+166], N[(N[(y * i), $MachinePrecision] + N[(N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[y], $MachinePrecision] + N[(N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.75 \cdot 10^{+152}:\\
\;\;\;\;x \cdot \log y + \left(\left(z + \left(t + a\right)\right) + b \cdot \log c\right)\\
\mathbf{elif}\;x \leq 7.2 \cdot 10^{+166}:\\
\;\;\;\;y \cdot i + \left(\left(b - 0.5\right) \cdot \log c + \left(a + \left(z + t\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log y + \frac{z + y \cdot i}{x}\right)\\
\end{array}
\end{array}
if x < -1.74999999999999991e152Initial program 99.7%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified91.4%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6491.4%
Simplified91.4%
if -1.74999999999999991e152 < x < 7.1999999999999994e166Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f6495.9%
Simplified95.9%
if 7.1999999999999994e166 < x Initial program 99.8%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified67.4%
Taylor expanded in x around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f6465.8%
Simplified65.8%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
distribute-lft-outN/A
mul-1-negN/A
remove-double-negN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6488.7%
Simplified88.7%
Final simplification94.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* (- b 0.5) (log c))))
(if (<= x -2e+162)
(+ (* y i) (+ (* x (log y)) t_1))
(if (<= x 3.2e+167)
(+ (* y i) (+ t_1 (+ a (+ z t))))
(* x (+ (log y) (/ (+ z (* y i)) x)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (b - 0.5) * log(c);
double tmp;
if (x <= -2e+162) {
tmp = (y * i) + ((x * log(y)) + t_1);
} else if (x <= 3.2e+167) {
tmp = (y * i) + (t_1 + (a + (z + t)));
} else {
tmp = x * (log(y) + ((z + (y * i)) / x));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = (b - 0.5d0) * log(c)
if (x <= (-2d+162)) then
tmp = (y * i) + ((x * log(y)) + t_1)
else if (x <= 3.2d+167) then
tmp = (y * i) + (t_1 + (a + (z + t)))
else
tmp = x * (log(y) + ((z + (y * i)) / x))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (b - 0.5) * Math.log(c);
double tmp;
if (x <= -2e+162) {
tmp = (y * i) + ((x * Math.log(y)) + t_1);
} else if (x <= 3.2e+167) {
tmp = (y * i) + (t_1 + (a + (z + t)));
} else {
tmp = x * (Math.log(y) + ((z + (y * i)) / x));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (b - 0.5) * math.log(c) tmp = 0 if x <= -2e+162: tmp = (y * i) + ((x * math.log(y)) + t_1) elif x <= 3.2e+167: tmp = (y * i) + (t_1 + (a + (z + t))) else: tmp = x * (math.log(y) + ((z + (y * i)) / x)) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(b - 0.5) * log(c)) tmp = 0.0 if (x <= -2e+162) tmp = Float64(Float64(y * i) + Float64(Float64(x * log(y)) + t_1)); elseif (x <= 3.2e+167) tmp = Float64(Float64(y * i) + Float64(t_1 + Float64(a + Float64(z + t)))); else tmp = Float64(x * Float64(log(y) + Float64(Float64(z + Float64(y * i)) / x))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (b - 0.5) * log(c); tmp = 0.0; if (x <= -2e+162) tmp = (y * i) + ((x * log(y)) + t_1); elseif (x <= 3.2e+167) tmp = (y * i) + (t_1 + (a + (z + t))); else tmp = x * (log(y) + ((z + (y * i)) / x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -2e+162], N[(N[(y * i), $MachinePrecision] + N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.2e+167], N[(N[(y * i), $MachinePrecision] + N[(t$95$1 + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[y], $MachinePrecision] + N[(N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - 0.5\right) \cdot \log c\\
\mathbf{if}\;x \leq -2 \cdot 10^{+162}:\\
\;\;\;\;y \cdot i + \left(x \cdot \log y + t\_1\right)\\
\mathbf{elif}\;x \leq 3.2 \cdot 10^{+167}:\\
\;\;\;\;y \cdot i + \left(t\_1 + \left(a + \left(z + t\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log y + \frac{z + y \cdot i}{x}\right)\\
\end{array}
\end{array}
if x < -1.9999999999999999e162Initial program 99.7%
Taylor expanded in x around inf
*-lowering-*.f64N/A
log-lowering-log.f6482.6%
Simplified82.6%
if -1.9999999999999999e162 < x < 3.19999999999999981e167Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f6495.5%
Simplified95.5%
if 3.19999999999999981e167 < x Initial program 99.8%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified67.4%
Taylor expanded in x around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f6465.8%
Simplified65.8%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
distribute-lft-outN/A
mul-1-negN/A
remove-double-negN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6488.7%
Simplified88.7%
Final simplification93.1%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= x -7.5e+155)
(+ (* x (log y)) (+ z (* b (log c))))
(if (<= x 7.2e+166)
(+ (* y i) (+ (* (- b 0.5) (log c)) (+ a (+ z t))))
(* x (+ (log y) (/ (+ z (* y i)) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -7.5e+155) {
tmp = (x * log(y)) + (z + (b * log(c)));
} else if (x <= 7.2e+166) {
tmp = (y * i) + (((b - 0.5) * log(c)) + (a + (z + t)));
} else {
tmp = x * (log(y) + ((z + (y * i)) / x));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (x <= (-7.5d+155)) then
tmp = (x * log(y)) + (z + (b * log(c)))
else if (x <= 7.2d+166) then
tmp = (y * i) + (((b - 0.5d0) * log(c)) + (a + (z + t)))
else
tmp = x * (log(y) + ((z + (y * i)) / x))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -7.5e+155) {
tmp = (x * Math.log(y)) + (z + (b * Math.log(c)));
} else if (x <= 7.2e+166) {
tmp = (y * i) + (((b - 0.5) * Math.log(c)) + (a + (z + t)));
} else {
tmp = x * (Math.log(y) + ((z + (y * i)) / x));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if x <= -7.5e+155: tmp = (x * math.log(y)) + (z + (b * math.log(c))) elif x <= 7.2e+166: tmp = (y * i) + (((b - 0.5) * math.log(c)) + (a + (z + t))) else: tmp = x * (math.log(y) + ((z + (y * i)) / x)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (x <= -7.5e+155) tmp = Float64(Float64(x * log(y)) + Float64(z + Float64(b * log(c)))); elseif (x <= 7.2e+166) tmp = Float64(Float64(y * i) + Float64(Float64(Float64(b - 0.5) * log(c)) + Float64(a + Float64(z + t)))); else tmp = Float64(x * Float64(log(y) + Float64(Float64(z + Float64(y * i)) / x))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (x <= -7.5e+155) tmp = (x * log(y)) + (z + (b * log(c))); elseif (x <= 7.2e+166) tmp = (y * i) + (((b - 0.5) * log(c)) + (a + (z + t))); else tmp = x * (log(y) + ((z + (y * i)) / x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[x, -7.5e+155], N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(z + N[(b * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 7.2e+166], N[(N[(y * i), $MachinePrecision] + N[(N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[y], $MachinePrecision] + N[(N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7.5 \cdot 10^{+155}:\\
\;\;\;\;x \cdot \log y + \left(z + b \cdot \log c\right)\\
\mathbf{elif}\;x \leq 7.2 \cdot 10^{+166}:\\
\;\;\;\;y \cdot i + \left(\left(b - 0.5\right) \cdot \log c + \left(a + \left(z + t\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log y + \frac{z + y \cdot i}{x}\right)\\
\end{array}
\end{array}
if x < -7.4999999999999999e155Initial program 99.7%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified91.4%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6491.4%
Simplified91.4%
Taylor expanded in z around inf
Simplified80.3%
if -7.4999999999999999e155 < x < 7.1999999999999994e166Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f6495.9%
Simplified95.9%
if 7.1999999999999994e166 < x Initial program 99.8%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified67.4%
Taylor expanded in x around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f6465.8%
Simplified65.8%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
distribute-lft-outN/A
mul-1-negN/A
remove-double-negN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6488.7%
Simplified88.7%
Final simplification93.0%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= x -5.6e+159)
(+ (+ (* x (log y)) z) (+ t a))
(if (<= x 5e+167)
(+ (* y i) (+ (* (- b 0.5) (log c)) (+ a (+ z t))))
(* x (+ (log y) (/ (+ z (* y i)) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -5.6e+159) {
tmp = ((x * log(y)) + z) + (t + a);
} else if (x <= 5e+167) {
tmp = (y * i) + (((b - 0.5) * log(c)) + (a + (z + t)));
} else {
tmp = x * (log(y) + ((z + (y * i)) / x));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (x <= (-5.6d+159)) then
tmp = ((x * log(y)) + z) + (t + a)
else if (x <= 5d+167) then
tmp = (y * i) + (((b - 0.5d0) * log(c)) + (a + (z + t)))
else
tmp = x * (log(y) + ((z + (y * i)) / x))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -5.6e+159) {
tmp = ((x * Math.log(y)) + z) + (t + a);
} else if (x <= 5e+167) {
tmp = (y * i) + (((b - 0.5) * Math.log(c)) + (a + (z + t)));
} else {
tmp = x * (Math.log(y) + ((z + (y * i)) / x));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if x <= -5.6e+159: tmp = ((x * math.log(y)) + z) + (t + a) elif x <= 5e+167: tmp = (y * i) + (((b - 0.5) * math.log(c)) + (a + (z + t))) else: tmp = x * (math.log(y) + ((z + (y * i)) / x)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (x <= -5.6e+159) tmp = Float64(Float64(Float64(x * log(y)) + z) + Float64(t + a)); elseif (x <= 5e+167) tmp = Float64(Float64(y * i) + Float64(Float64(Float64(b - 0.5) * log(c)) + Float64(a + Float64(z + t)))); else tmp = Float64(x * Float64(log(y) + Float64(Float64(z + Float64(y * i)) / x))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (x <= -5.6e+159) tmp = ((x * log(y)) + z) + (t + a); elseif (x <= 5e+167) tmp = (y * i) + (((b - 0.5) * log(c)) + (a + (z + t))); else tmp = x * (log(y) + ((z + (y * i)) / x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[x, -5.6e+159], N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + N[(t + a), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 5e+167], N[(N[(y * i), $MachinePrecision] + N[(N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[y], $MachinePrecision] + N[(N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.6 \cdot 10^{+159}:\\
\;\;\;\;\left(x \cdot \log y + z\right) + \left(t + a\right)\\
\mathbf{elif}\;x \leq 5 \cdot 10^{+167}:\\
\;\;\;\;y \cdot i + \left(\left(b - 0.5\right) \cdot \log c + \left(a + \left(z + t\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log y + \frac{z + y \cdot i}{x}\right)\\
\end{array}
\end{array}
if x < -5.6000000000000002e159Initial program 99.7%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified91.4%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6491.4%
Simplified91.4%
Taylor expanded in b around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f6476.8%
Simplified76.8%
if -5.6000000000000002e159 < x < 4.9999999999999997e167Initial program 99.9%
Taylor expanded in x around 0
+-lowering-+.f6495.9%
Simplified95.9%
if 4.9999999999999997e167 < x Initial program 99.8%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified67.4%
Taylor expanded in x around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f6465.8%
Simplified65.8%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
distribute-lft-outN/A
mul-1-negN/A
remove-double-negN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6488.7%
Simplified88.7%
Final simplification92.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* y i) (+ z (* (- b 0.5) (log c))))))
(if (<= (- b 0.5) -2.5e+72)
t_1
(if (<= (- b 0.5) 1e+89) (+ (+ (* x (log y)) z) (+ t a)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (y * i) + (z + ((b - 0.5) * log(c)));
double tmp;
if ((b - 0.5) <= -2.5e+72) {
tmp = t_1;
} else if ((b - 0.5) <= 1e+89) {
tmp = ((x * log(y)) + z) + (t + a);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = (y * i) + (z + ((b - 0.5d0) * log(c)))
if ((b - 0.5d0) <= (-2.5d+72)) then
tmp = t_1
else if ((b - 0.5d0) <= 1d+89) then
tmp = ((x * log(y)) + z) + (t + a)
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 c, double i) {
double t_1 = (y * i) + (z + ((b - 0.5) * Math.log(c)));
double tmp;
if ((b - 0.5) <= -2.5e+72) {
tmp = t_1;
} else if ((b - 0.5) <= 1e+89) {
tmp = ((x * Math.log(y)) + z) + (t + a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * i) + (z + ((b - 0.5) * math.log(c))) tmp = 0 if (b - 0.5) <= -2.5e+72: tmp = t_1 elif (b - 0.5) <= 1e+89: tmp = ((x * math.log(y)) + z) + (t + a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(y * i) + Float64(z + Float64(Float64(b - 0.5) * log(c)))) tmp = 0.0 if (Float64(b - 0.5) <= -2.5e+72) tmp = t_1; elseif (Float64(b - 0.5) <= 1e+89) tmp = Float64(Float64(Float64(x * log(y)) + z) + Float64(t + a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (y * i) + (z + ((b - 0.5) * log(c))); tmp = 0.0; if ((b - 0.5) <= -2.5e+72) tmp = t_1; elseif ((b - 0.5) <= 1e+89) tmp = ((x * log(y)) + z) + (t + a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(y * i), $MachinePrecision] + N[(z + N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(b - 0.5), $MachinePrecision], -2.5e+72], t$95$1, If[LessEqual[N[(b - 0.5), $MachinePrecision], 1e+89], N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + N[(t + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot i + \left(z + \left(b - 0.5\right) \cdot \log c\right)\\
\mathbf{if}\;b - 0.5 \leq -2.5 \cdot 10^{+72}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b - 0.5 \leq 10^{+89}:\\
\;\;\;\;\left(x \cdot \log y + z\right) + \left(t + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 b #s(literal 1/2 binary64)) < -2.49999999999999996e72 or 9.99999999999999995e88 < (-.f64 b #s(literal 1/2 binary64)) Initial program 99.8%
Taylor expanded in z around inf
Simplified72.7%
if -2.49999999999999996e72 < (-.f64 b #s(literal 1/2 binary64)) < 9.99999999999999995e88Initial program 99.9%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified83.4%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6479.7%
Simplified79.7%
Taylor expanded in b around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f6477.9%
Simplified77.9%
Final simplification76.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* x (log y))))
(if (<= x -6.4e+152)
(+ (+ t_1 z) (+ t a))
(if (<= x -1.3e+98)
(+ (* y i) (* b (log c)))
(if (<= x 3.8e+173)
(+ (* (log c) (+ b -0.5)) (+ z a))
(+ t_1 (* y i)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x * log(y);
double tmp;
if (x <= -6.4e+152) {
tmp = (t_1 + z) + (t + a);
} else if (x <= -1.3e+98) {
tmp = (y * i) + (b * log(c));
} else if (x <= 3.8e+173) {
tmp = (log(c) * (b + -0.5)) + (z + a);
} else {
tmp = t_1 + (y * i);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x * log(y)
if (x <= (-6.4d+152)) then
tmp = (t_1 + z) + (t + a)
else if (x <= (-1.3d+98)) then
tmp = (y * i) + (b * log(c))
else if (x <= 3.8d+173) then
tmp = (log(c) * (b + (-0.5d0))) + (z + a)
else
tmp = t_1 + (y * i)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x * Math.log(y);
double tmp;
if (x <= -6.4e+152) {
tmp = (t_1 + z) + (t + a);
} else if (x <= -1.3e+98) {
tmp = (y * i) + (b * Math.log(c));
} else if (x <= 3.8e+173) {
tmp = (Math.log(c) * (b + -0.5)) + (z + a);
} else {
tmp = t_1 + (y * i);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x * math.log(y) tmp = 0 if x <= -6.4e+152: tmp = (t_1 + z) + (t + a) elif x <= -1.3e+98: tmp = (y * i) + (b * math.log(c)) elif x <= 3.8e+173: tmp = (math.log(c) * (b + -0.5)) + (z + a) else: tmp = t_1 + (y * i) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x * log(y)) tmp = 0.0 if (x <= -6.4e+152) tmp = Float64(Float64(t_1 + z) + Float64(t + a)); elseif (x <= -1.3e+98) tmp = Float64(Float64(y * i) + Float64(b * log(c))); elseif (x <= 3.8e+173) tmp = Float64(Float64(log(c) * Float64(b + -0.5)) + Float64(z + a)); else tmp = Float64(t_1 + Float64(y * i)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x * log(y); tmp = 0.0; if (x <= -6.4e+152) tmp = (t_1 + z) + (t + a); elseif (x <= -1.3e+98) tmp = (y * i) + (b * log(c)); elseif (x <= 3.8e+173) tmp = (log(c) * (b + -0.5)) + (z + a); else tmp = t_1 + (y * i); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -6.4e+152], N[(N[(t$95$1 + z), $MachinePrecision] + N[(t + a), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -1.3e+98], N[(N[(y * i), $MachinePrecision] + N[(b * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.8e+173], N[(N[(N[Log[c], $MachinePrecision] * N[(b + -0.5), $MachinePrecision]), $MachinePrecision] + N[(z + a), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + N[(y * i), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \log y\\
\mathbf{if}\;x \leq -6.4 \cdot 10^{+152}:\\
\;\;\;\;\left(t\_1 + z\right) + \left(t + a\right)\\
\mathbf{elif}\;x \leq -1.3 \cdot 10^{+98}:\\
\;\;\;\;y \cdot i + b \cdot \log c\\
\mathbf{elif}\;x \leq 3.8 \cdot 10^{+173}:\\
\;\;\;\;\log c \cdot \left(b + -0.5\right) + \left(z + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 + y \cdot i\\
\end{array}
\end{array}
if x < -6.40000000000000011e152Initial program 99.7%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified91.4%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6491.4%
Simplified91.4%
Taylor expanded in b around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f6476.8%
Simplified76.8%
if -6.40000000000000011e152 < x < -1.3e98Initial program 99.8%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6485.6%
Simplified85.6%
if -1.3e98 < x < 3.80000000000000011e173Initial program 99.9%
Taylor expanded in z around inf
Simplified75.1%
Taylor expanded in y around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6457.9%
Simplified57.9%
if 3.80000000000000011e173 < x Initial program 99.8%
Taylor expanded in x around inf
*-lowering-*.f64N/A
log-lowering-log.f6481.6%
Simplified81.6%
Final simplification64.6%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= x -2.5e+159)
(+ (+ (* x (log y)) z) (+ t a))
(if (<= x 1.9e+167)
(+ (* y i) (+ (* (- b 0.5) (log c)) (+ z a)))
(* x (+ (log y) (/ (+ z (* y i)) x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -2.5e+159) {
tmp = ((x * log(y)) + z) + (t + a);
} else if (x <= 1.9e+167) {
tmp = (y * i) + (((b - 0.5) * log(c)) + (z + a));
} else {
tmp = x * (log(y) + ((z + (y * i)) / x));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (x <= (-2.5d+159)) then
tmp = ((x * log(y)) + z) + (t + a)
else if (x <= 1.9d+167) then
tmp = (y * i) + (((b - 0.5d0) * log(c)) + (z + a))
else
tmp = x * (log(y) + ((z + (y * i)) / x))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (x <= -2.5e+159) {
tmp = ((x * Math.log(y)) + z) + (t + a);
} else if (x <= 1.9e+167) {
tmp = (y * i) + (((b - 0.5) * Math.log(c)) + (z + a));
} else {
tmp = x * (Math.log(y) + ((z + (y * i)) / x));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if x <= -2.5e+159: tmp = ((x * math.log(y)) + z) + (t + a) elif x <= 1.9e+167: tmp = (y * i) + (((b - 0.5) * math.log(c)) + (z + a)) else: tmp = x * (math.log(y) + ((z + (y * i)) / x)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (x <= -2.5e+159) tmp = Float64(Float64(Float64(x * log(y)) + z) + Float64(t + a)); elseif (x <= 1.9e+167) tmp = Float64(Float64(y * i) + Float64(Float64(Float64(b - 0.5) * log(c)) + Float64(z + a))); else tmp = Float64(x * Float64(log(y) + Float64(Float64(z + Float64(y * i)) / x))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (x <= -2.5e+159) tmp = ((x * log(y)) + z) + (t + a); elseif (x <= 1.9e+167) tmp = (y * i) + (((b - 0.5) * log(c)) + (z + a)); else tmp = x * (log(y) + ((z + (y * i)) / x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[x, -2.5e+159], N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + N[(t + a), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.9e+167], N[(N[(y * i), $MachinePrecision] + N[(N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision] + N[(z + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Log[y], $MachinePrecision] + N[(N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.5 \cdot 10^{+159}:\\
\;\;\;\;\left(x \cdot \log y + z\right) + \left(t + a\right)\\
\mathbf{elif}\;x \leq 1.9 \cdot 10^{+167}:\\
\;\;\;\;y \cdot i + \left(\left(b - 0.5\right) \cdot \log c + \left(z + a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\log y + \frac{z + y \cdot i}{x}\right)\\
\end{array}
\end{array}
if x < -2.50000000000000002e159Initial program 99.7%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified91.4%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6491.4%
Simplified91.4%
Taylor expanded in b around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f6476.8%
Simplified76.8%
if -2.50000000000000002e159 < x < 1.89999999999999997e167Initial program 99.9%
Taylor expanded in z around inf
Simplified75.8%
if 1.89999999999999997e167 < x Initial program 99.8%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified67.4%
Taylor expanded in x around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f6465.8%
Simplified65.8%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
distribute-lft-outN/A
mul-1-negN/A
remove-double-negN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6488.7%
Simplified88.7%
Final simplification77.4%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= z -9.5e+179)
(+ (* y i) (* z (+ (/ a z) 1.0)))
(if (<= z -7.8e+114)
(+ (+ (* x (log y)) z) (+ t a))
(+ (* y i) (+ a (* (- b 0.5) (log c)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -9.5e+179) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else if (z <= -7.8e+114) {
tmp = ((x * log(y)) + z) + (t + a);
} else {
tmp = (y * i) + (a + ((b - 0.5) * log(c)));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (z <= (-9.5d+179)) then
tmp = (y * i) + (z * ((a / z) + 1.0d0))
else if (z <= (-7.8d+114)) then
tmp = ((x * log(y)) + z) + (t + a)
else
tmp = (y * i) + (a + ((b - 0.5d0) * log(c)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -9.5e+179) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else if (z <= -7.8e+114) {
tmp = ((x * Math.log(y)) + z) + (t + a);
} else {
tmp = (y * i) + (a + ((b - 0.5) * Math.log(c)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if z <= -9.5e+179: tmp = (y * i) + (z * ((a / z) + 1.0)) elif z <= -7.8e+114: tmp = ((x * math.log(y)) + z) + (t + a) else: tmp = (y * i) + (a + ((b - 0.5) * math.log(c))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (z <= -9.5e+179) tmp = Float64(Float64(y * i) + Float64(z * Float64(Float64(a / z) + 1.0))); elseif (z <= -7.8e+114) tmp = Float64(Float64(Float64(x * log(y)) + z) + Float64(t + a)); else tmp = Float64(Float64(y * i) + Float64(a + Float64(Float64(b - 0.5) * log(c)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (z <= -9.5e+179) tmp = (y * i) + (z * ((a / z) + 1.0)); elseif (z <= -7.8e+114) tmp = ((x * log(y)) + z) + (t + a); else tmp = (y * i) + (a + ((b - 0.5) * log(c))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[z, -9.5e+179], N[(N[(y * i), $MachinePrecision] + N[(z * N[(N[(a / z), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -7.8e+114], N[(N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision] + N[(t + a), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + N[(N[(b - 0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9.5 \cdot 10^{+179}:\\
\;\;\;\;y \cdot i + z \cdot \left(\frac{a}{z} + 1\right)\\
\mathbf{elif}\;z \leq -7.8 \cdot 10^{+114}:\\
\;\;\;\;\left(x \cdot \log y + z\right) + \left(t + a\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + \left(b - 0.5\right) \cdot \log c\right)\\
\end{array}
\end{array}
if z < -9.5e179Initial program 99.9%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified99.9%
Taylor expanded in a around inf
/-lowering-/.f6470.1%
Simplified70.1%
if -9.5e179 < z < -7.8000000000000001e114Initial program 99.9%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified99.9%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6499.9%
Simplified99.9%
Taylor expanded in b around 0
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f6488.2%
Simplified88.2%
if -7.8000000000000001e114 < z Initial program 99.8%
Taylor expanded in a around inf
Simplified58.5%
Final simplification61.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* x (+ (log y) (/ (+ z (* y i)) x)))))
(if (<= x -4.5e+97)
t_1
(if (<= x 6.5e+167) (+ (* (log c) (+ b -0.5)) (+ z a)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x * (log(y) + ((z + (y * i)) / x));
double tmp;
if (x <= -4.5e+97) {
tmp = t_1;
} else if (x <= 6.5e+167) {
tmp = (log(c) * (b + -0.5)) + (z + a);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x * (log(y) + ((z + (y * i)) / x))
if (x <= (-4.5d+97)) then
tmp = t_1
else if (x <= 6.5d+167) then
tmp = (log(c) * (b + (-0.5d0))) + (z + a)
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 c, double i) {
double t_1 = x * (Math.log(y) + ((z + (y * i)) / x));
double tmp;
if (x <= -4.5e+97) {
tmp = t_1;
} else if (x <= 6.5e+167) {
tmp = (Math.log(c) * (b + -0.5)) + (z + a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x * (math.log(y) + ((z + (y * i)) / x)) tmp = 0 if x <= -4.5e+97: tmp = t_1 elif x <= 6.5e+167: tmp = (math.log(c) * (b + -0.5)) + (z + a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x * Float64(log(y) + Float64(Float64(z + Float64(y * i)) / x))) tmp = 0.0 if (x <= -4.5e+97) tmp = t_1; elseif (x <= 6.5e+167) tmp = Float64(Float64(log(c) * Float64(b + -0.5)) + Float64(z + a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x * (log(y) + ((z + (y * i)) / x)); tmp = 0.0; if (x <= -4.5e+97) tmp = t_1; elseif (x <= 6.5e+167) tmp = (log(c) * (b + -0.5)) + (z + a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x * N[(N[Log[y], $MachinePrecision] + N[(N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -4.5e+97], t$95$1, If[LessEqual[x, 6.5e+167], N[(N[(N[Log[c], $MachinePrecision] * N[(b + -0.5), $MachinePrecision]), $MachinePrecision] + N[(z + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(\log y + \frac{z + y \cdot i}{x}\right)\\
\mathbf{if}\;x \leq -4.5 \cdot 10^{+97}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 6.5 \cdot 10^{+167}:\\
\;\;\;\;\log c \cdot \left(b + -0.5\right) + \left(z + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -4.49999999999999976e97 or 6.5e167 < x Initial program 99.8%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified67.5%
Taylor expanded in x around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f6454.8%
Simplified54.8%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
distribute-lft-outN/A
mul-1-negN/A
remove-double-negN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6476.4%
Simplified76.4%
if -4.49999999999999976e97 < x < 6.5e167Initial program 99.9%
Taylor expanded in z around inf
Simplified74.9%
Taylor expanded in y around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6457.7%
Simplified57.7%
Final simplification63.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* x (log y)) (* y i))))
(if (<= x -2.8e+147)
t_1
(if (<= x 1.55e+173) (+ (* (log c) (+ b -0.5)) (+ z a)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (x * log(y)) + (y * i);
double tmp;
if (x <= -2.8e+147) {
tmp = t_1;
} else if (x <= 1.55e+173) {
tmp = (log(c) * (b + -0.5)) + (z + a);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = (x * log(y)) + (y * i)
if (x <= (-2.8d+147)) then
tmp = t_1
else if (x <= 1.55d+173) then
tmp = (log(c) * (b + (-0.5d0))) + (z + a)
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 c, double i) {
double t_1 = (x * Math.log(y)) + (y * i);
double tmp;
if (x <= -2.8e+147) {
tmp = t_1;
} else if (x <= 1.55e+173) {
tmp = (Math.log(c) * (b + -0.5)) + (z + a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (x * math.log(y)) + (y * i) tmp = 0 if x <= -2.8e+147: tmp = t_1 elif x <= 1.55e+173: tmp = (math.log(c) * (b + -0.5)) + (z + a) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(x * log(y)) + Float64(y * i)) tmp = 0.0 if (x <= -2.8e+147) tmp = t_1; elseif (x <= 1.55e+173) tmp = Float64(Float64(log(c) * Float64(b + -0.5)) + Float64(z + a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (x * log(y)) + (y * i); tmp = 0.0; if (x <= -2.8e+147) tmp = t_1; elseif (x <= 1.55e+173) tmp = (log(c) * (b + -0.5)) + (z + a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -2.8e+147], t$95$1, If[LessEqual[x, 1.55e+173], N[(N[(N[Log[c], $MachinePrecision] * N[(b + -0.5), $MachinePrecision]), $MachinePrecision] + N[(z + a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \log y + y \cdot i\\
\mathbf{if}\;x \leq -2.8 \cdot 10^{+147}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 1.55 \cdot 10^{+173}:\\
\;\;\;\;\log c \cdot \left(b + -0.5\right) + \left(z + a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -2.8000000000000001e147 or 1.55e173 < x Initial program 99.8%
Taylor expanded in x around inf
*-lowering-*.f64N/A
log-lowering-log.f6474.2%
Simplified74.2%
if -2.8000000000000001e147 < x < 1.55e173Initial program 99.9%
Taylor expanded in z around inf
Simplified75.6%
Taylor expanded in y around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f6457.0%
Simplified57.0%
Final simplification61.4%
(FPCore (x y z t a b c i) :precision binary64 (if (<= z -1.06e+40) (+ (* y i) (* z (+ (/ a z) 1.0))) (if (<= z -1.12e-172) (+ (* x (log y)) a) (+ a (* y i)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -1.06e+40) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else if (z <= -1.12e-172) {
tmp = (x * log(y)) + a;
} else {
tmp = a + (y * i);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (z <= (-1.06d+40)) then
tmp = (y * i) + (z * ((a / z) + 1.0d0))
else if (z <= (-1.12d-172)) then
tmp = (x * log(y)) + a
else
tmp = a + (y * i)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -1.06e+40) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else if (z <= -1.12e-172) {
tmp = (x * Math.log(y)) + a;
} else {
tmp = a + (y * i);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if z <= -1.06e+40: tmp = (y * i) + (z * ((a / z) + 1.0)) elif z <= -1.12e-172: tmp = (x * math.log(y)) + a else: tmp = a + (y * i) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (z <= -1.06e+40) tmp = Float64(Float64(y * i) + Float64(z * Float64(Float64(a / z) + 1.0))); elseif (z <= -1.12e-172) tmp = Float64(Float64(x * log(y)) + a); else tmp = Float64(a + Float64(y * i)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (z <= -1.06e+40) tmp = (y * i) + (z * ((a / z) + 1.0)); elseif (z <= -1.12e-172) tmp = (x * log(y)) + a; else tmp = a + (y * i); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[z, -1.06e+40], N[(N[(y * i), $MachinePrecision] + N[(z * N[(N[(a / z), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -1.12e-172], N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision], N[(a + N[(y * i), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.06 \cdot 10^{+40}:\\
\;\;\;\;y \cdot i + z \cdot \left(\frac{a}{z} + 1\right)\\
\mathbf{elif}\;z \leq -1.12 \cdot 10^{-172}:\\
\;\;\;\;x \cdot \log y + a\\
\mathbf{else}:\\
\;\;\;\;a + y \cdot i\\
\end{array}
\end{array}
if z < -1.05999999999999996e40Initial program 99.9%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified99.8%
Taylor expanded in a around inf
/-lowering-/.f6457.4%
Simplified57.4%
if -1.05999999999999996e40 < z < -1.11999999999999996e-172Initial program 99.9%
Taylor expanded in y around 0
associate-+r+N/A
associate-+r+N/A
cancel-sign-subN/A
log-recN/A
unsub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-+r+N/A
associate-+r+N/A
associate-+r+N/A
Simplified84.4%
Taylor expanded in a around inf
Simplified41.5%
if -1.11999999999999996e-172 < z Initial program 99.8%
Taylor expanded in a around inf
Simplified41.4%
Final simplification44.8%
(FPCore (x y z t a b c i) :precision binary64 (if (<= z -1.95e+38) (+ (* y i) (* z (+ (/ a z) 1.0))) (if (<= z -8.5e+23) (* x (log y)) (+ a (* y i)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -1.95e+38) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else if (z <= -8.5e+23) {
tmp = x * log(y);
} else {
tmp = a + (y * i);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (z <= (-1.95d+38)) then
tmp = (y * i) + (z * ((a / z) + 1.0d0))
else if (z <= (-8.5d+23)) then
tmp = x * log(y)
else
tmp = a + (y * i)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -1.95e+38) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else if (z <= -8.5e+23) {
tmp = x * Math.log(y);
} else {
tmp = a + (y * i);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if z <= -1.95e+38: tmp = (y * i) + (z * ((a / z) + 1.0)) elif z <= -8.5e+23: tmp = x * math.log(y) else: tmp = a + (y * i) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (z <= -1.95e+38) tmp = Float64(Float64(y * i) + Float64(z * Float64(Float64(a / z) + 1.0))); elseif (z <= -8.5e+23) tmp = Float64(x * log(y)); else tmp = Float64(a + Float64(y * i)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (z <= -1.95e+38) tmp = (y * i) + (z * ((a / z) + 1.0)); elseif (z <= -8.5e+23) tmp = x * log(y); else tmp = a + (y * i); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[z, -1.95e+38], N[(N[(y * i), $MachinePrecision] + N[(z * N[(N[(a / z), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -8.5e+23], N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision], N[(a + N[(y * i), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.95 \cdot 10^{+38}:\\
\;\;\;\;y \cdot i + z \cdot \left(\frac{a}{z} + 1\right)\\
\mathbf{elif}\;z \leq -8.5 \cdot 10^{+23}:\\
\;\;\;\;x \cdot \log y\\
\mathbf{else}:\\
\;\;\;\;a + y \cdot i\\
\end{array}
\end{array}
if z < -1.95000000000000012e38Initial program 99.9%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified99.8%
Taylor expanded in a around inf
/-lowering-/.f6457.4%
Simplified57.4%
if -1.95000000000000012e38 < z < -8.5000000000000001e23Initial program 100.0%
Taylor expanded in x around inf
*-lowering-*.f64N/A
log-lowering-log.f6439.8%
Simplified39.8%
if -8.5000000000000001e23 < z Initial program 99.8%
Taylor expanded in a around inf
Simplified39.3%
Final simplification43.1%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y 6.2e-140) (+ t a) (if (<= y 5.2e-96) z (if (<= y 6.5e+48) (+ t a) (* y i)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 6.2e-140) {
tmp = t + a;
} else if (y <= 5.2e-96) {
tmp = z;
} else if (y <= 6.5e+48) {
tmp = t + a;
} else {
tmp = y * i;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (y <= 6.2d-140) then
tmp = t + a
else if (y <= 5.2d-96) then
tmp = z
else if (y <= 6.5d+48) then
tmp = t + a
else
tmp = y * i
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 6.2e-140) {
tmp = t + a;
} else if (y <= 5.2e-96) {
tmp = z;
} else if (y <= 6.5e+48) {
tmp = t + a;
} else {
tmp = y * i;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= 6.2e-140: tmp = t + a elif y <= 5.2e-96: tmp = z elif y <= 6.5e+48: tmp = t + a else: tmp = y * i return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= 6.2e-140) tmp = Float64(t + a); elseif (y <= 5.2e-96) tmp = z; elseif (y <= 6.5e+48) tmp = Float64(t + a); else tmp = Float64(y * i); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= 6.2e-140) tmp = t + a; elseif (y <= 5.2e-96) tmp = z; elseif (y <= 6.5e+48) tmp = t + a; else tmp = y * i; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, 6.2e-140], N[(t + a), $MachinePrecision], If[LessEqual[y, 5.2e-96], z, If[LessEqual[y, 6.5e+48], N[(t + a), $MachinePrecision], N[(y * i), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 6.2 \cdot 10^{-140}:\\
\;\;\;\;t + a\\
\mathbf{elif}\;y \leq 5.2 \cdot 10^{-96}:\\
\;\;\;\;z\\
\mathbf{elif}\;y \leq 6.5 \cdot 10^{+48}:\\
\;\;\;\;t + a\\
\mathbf{else}:\\
\;\;\;\;y \cdot i\\
\end{array}
\end{array}
if y < 6.1999999999999998e-140 or 5.2000000000000003e-96 < y < 6.49999999999999972e48Initial program 99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6478.0%
Simplified78.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6475.0%
Simplified75.0%
Taylor expanded in t around inf
Simplified44.3%
if 6.1999999999999998e-140 < y < 5.2000000000000003e-96Initial program 100.0%
Taylor expanded in z around inf
Simplified28.2%
if 6.49999999999999972e48 < y Initial program 99.8%
Taylor expanded in y around inf
*-lowering-*.f6446.9%
Simplified46.9%
Final simplification44.1%
(FPCore (x y z t a b c i) :precision binary64 (if (<= z -2.35e+14) (+ (* y i) (* z (+ (/ a z) 1.0))) (+ a (* y i))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -2.35e+14) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else {
tmp = a + (y * i);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (z <= (-2.35d+14)) then
tmp = (y * i) + (z * ((a / z) + 1.0d0))
else
tmp = a + (y * i)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -2.35e+14) {
tmp = (y * i) + (z * ((a / z) + 1.0));
} else {
tmp = a + (y * i);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if z <= -2.35e+14: tmp = (y * i) + (z * ((a / z) + 1.0)) else: tmp = a + (y * i) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (z <= -2.35e+14) tmp = Float64(Float64(y * i) + Float64(z * Float64(Float64(a / z) + 1.0))); else tmp = Float64(a + Float64(y * i)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (z <= -2.35e+14) tmp = (y * i) + (z * ((a / z) + 1.0)); else tmp = a + (y * i); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[z, -2.35e+14], N[(N[(y * i), $MachinePrecision] + N[(z * N[(N[(a / z), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a + N[(y * i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.35 \cdot 10^{+14}:\\
\;\;\;\;y \cdot i + z \cdot \left(\frac{a}{z} + 1\right)\\
\mathbf{else}:\\
\;\;\;\;a + y \cdot i\\
\end{array}
\end{array}
if z < -2.35e14Initial program 99.9%
Taylor expanded in z around -inf
associate-*r*N/A
*-commutativeN/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l*N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
*-lowering-*.f64N/A
Simplified99.8%
Taylor expanded in a around inf
/-lowering-/.f6455.0%
Simplified55.0%
if -2.35e14 < z Initial program 99.8%
Taylor expanded in a around inf
Simplified39.2%
Final simplification42.8%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y 7.6e+47) (+ z a) (if (<= y 3.6e+98) (+ a (* y i)) (+ z (* y i)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 7.6e+47) {
tmp = z + a;
} else if (y <= 3.6e+98) {
tmp = a + (y * i);
} else {
tmp = z + (y * i);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (y <= 7.6d+47) then
tmp = z + a
else if (y <= 3.6d+98) then
tmp = a + (y * i)
else
tmp = z + (y * i)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 7.6e+47) {
tmp = z + a;
} else if (y <= 3.6e+98) {
tmp = a + (y * i);
} else {
tmp = z + (y * i);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= 7.6e+47: tmp = z + a elif y <= 3.6e+98: tmp = a + (y * i) else: tmp = z + (y * i) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= 7.6e+47) tmp = Float64(z + a); elseif (y <= 3.6e+98) tmp = Float64(a + Float64(y * i)); else tmp = Float64(z + Float64(y * i)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= 7.6e+47) tmp = z + a; elseif (y <= 3.6e+98) tmp = a + (y * i); else tmp = z + (y * i); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, 7.6e+47], N[(z + a), $MachinePrecision], If[LessEqual[y, 3.6e+98], N[(a + N[(y * i), $MachinePrecision]), $MachinePrecision], N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 7.6 \cdot 10^{+47}:\\
\;\;\;\;z + a\\
\mathbf{elif}\;y \leq 3.6 \cdot 10^{+98}:\\
\;\;\;\;a + y \cdot i\\
\mathbf{else}:\\
\;\;\;\;z + y \cdot i\\
\end{array}
\end{array}
if y < 7.6000000000000007e47Initial program 99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6478.2%
Simplified78.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6475.6%
Simplified75.6%
Taylor expanded in z around inf
Simplified38.9%
if 7.6000000000000007e47 < y < 3.59999999999999981e98Initial program 99.9%
Taylor expanded in a around inf
Simplified44.3%
if 3.59999999999999981e98 < y Initial program 99.8%
Taylor expanded in z around inf
Simplified54.7%
Final simplification44.5%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y 2.4e-140) a (if (<= y 4.7e+48) z (* y i))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 2.4e-140) {
tmp = a;
} else if (y <= 4.7e+48) {
tmp = z;
} else {
tmp = y * i;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (y <= 2.4d-140) then
tmp = a
else if (y <= 4.7d+48) then
tmp = z
else
tmp = y * i
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 2.4e-140) {
tmp = a;
} else if (y <= 4.7e+48) {
tmp = z;
} else {
tmp = y * i;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= 2.4e-140: tmp = a elif y <= 4.7e+48: tmp = z else: tmp = y * i return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= 2.4e-140) tmp = a; elseif (y <= 4.7e+48) tmp = z; else tmp = Float64(y * i); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= 2.4e-140) tmp = a; elseif (y <= 4.7e+48) tmp = z; else tmp = y * i; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, 2.4e-140], a, If[LessEqual[y, 4.7e+48], z, N[(y * i), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2.4 \cdot 10^{-140}:\\
\;\;\;\;a\\
\mathbf{elif}\;y \leq 4.7 \cdot 10^{+48}:\\
\;\;\;\;z\\
\mathbf{else}:\\
\;\;\;\;y \cdot i\\
\end{array}
\end{array}
if y < 2.39999999999999987e-140Initial program 99.9%
Taylor expanded in a around inf
Simplified25.9%
if 2.39999999999999987e-140 < y < 4.70000000000000012e48Initial program 99.8%
Taylor expanded in z around inf
Simplified19.5%
if 4.70000000000000012e48 < y Initial program 99.8%
Taylor expanded in y around inf
*-lowering-*.f6446.9%
Simplified46.9%
Final simplification32.3%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y 1.05e+48) (+ z a) (+ a (* y i))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 1.05e+48) {
tmp = z + a;
} else {
tmp = a + (y * i);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (y <= 1.05d+48) then
tmp = z + a
else
tmp = a + (y * i)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 1.05e+48) {
tmp = z + a;
} else {
tmp = a + (y * i);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= 1.05e+48: tmp = z + a else: tmp = a + (y * i) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= 1.05e+48) tmp = Float64(z + a); else tmp = Float64(a + Float64(y * i)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= 1.05e+48) tmp = z + a; else tmp = a + (y * i); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, 1.05e+48], N[(z + a), $MachinePrecision], N[(a + N[(y * i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.05 \cdot 10^{+48}:\\
\;\;\;\;z + a\\
\mathbf{else}:\\
\;\;\;\;a + y \cdot i\\
\end{array}
\end{array}
if y < 1.0499999999999999e48Initial program 99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6478.2%
Simplified78.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6475.6%
Simplified75.6%
Taylor expanded in z around inf
Simplified38.9%
if 1.0499999999999999e48 < y Initial program 99.8%
Taylor expanded in a around inf
Simplified52.7%
Final simplification44.5%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y 8.8e+49) (+ z a) (* y i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 8.8e+49) {
tmp = z + a;
} else {
tmp = y * i;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (y <= 8.8d+49) then
tmp = z + a
else
tmp = y * i
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= 8.8e+49) {
tmp = z + a;
} else {
tmp = y * i;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= 8.8e+49: tmp = z + a else: tmp = y * i return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= 8.8e+49) tmp = Float64(z + a); else tmp = Float64(y * i); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= 8.8e+49) tmp = z + a; else tmp = y * i; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, 8.8e+49], N[(z + a), $MachinePrecision], N[(y * i), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 8.8 \cdot 10^{+49}:\\
\;\;\;\;z + a\\
\mathbf{else}:\\
\;\;\;\;y \cdot i\\
\end{array}
\end{array}
if y < 8.8000000000000003e49Initial program 99.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6478.3%
Simplified78.3%
Taylor expanded in y around 0
+-lowering-+.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
log-lowering-log.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6475.7%
Simplified75.7%
Taylor expanded in z around inf
Simplified39.2%
if 8.8000000000000003e49 < y Initial program 99.8%
Taylor expanded in y around inf
*-lowering-*.f6446.9%
Simplified46.9%
Final simplification42.3%
(FPCore (x y z t a b c i) :precision binary64 (if (<= z -1.75e+48) z a))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -1.75e+48) {
tmp = z;
} else {
tmp = a;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
real(8) :: tmp
if (z <= (-1.75d+48)) then
tmp = z
else
tmp = a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (z <= -1.75e+48) {
tmp = z;
} else {
tmp = a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if z <= -1.75e+48: tmp = z else: tmp = a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (z <= -1.75e+48) tmp = z; else tmp = a; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (z <= -1.75e+48) tmp = z; else tmp = a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[z, -1.75e+48], z, a]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.75 \cdot 10^{+48}:\\
\;\;\;\;z\\
\mathbf{else}:\\
\;\;\;\;a\\
\end{array}
\end{array}
if z < -1.7499999999999999e48Initial program 99.9%
Taylor expanded in z around inf
Simplified32.1%
if -1.7499999999999999e48 < z Initial program 99.8%
Taylor expanded in a around inf
Simplified18.9%
(FPCore (x y z t a b c i) :precision binary64 a)
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return a;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: c
real(8), intent (in) :: i
code = a
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return a;
}
def code(x, y, z, t, a, b, c, i): return a
function code(x, y, z, t, a, b, c, i) return a end
function tmp = code(x, y, z, t, a, b, c, i) tmp = a; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := a
\begin{array}{l}
\\
a
\end{array}
Initial program 99.8%
Taylor expanded in a around inf
Simplified16.9%
herbie shell --seed 2024160
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logBeta from math-functions-0.1.5.2, B"
:precision binary64
(+ (+ (+ (+ (+ (* x (log y)) z) t) a) (* (- b 0.5) (log c))) (* y i)))