
(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 22 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 (fma y i (fma (+ b -0.5) (log c) (+ z (fma x (log y) (+ t a))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return fma(y, i, fma((b + -0.5), log(c), (z + fma(x, log(y), (t + a)))));
}
function code(x, y, z, t, a, b, c, i) return fma(y, i, fma(Float64(b + -0.5), log(c), Float64(z + fma(x, log(y), Float64(t + a))))) end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(y * i + N[(N[(b + -0.5), $MachinePrecision] * N[Log[c], $MachinePrecision] + N[(z + N[(x * N[Log[y], $MachinePrecision] + N[(t + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(y, i, \mathsf{fma}\left(b + -0.5, \log c, z + \mathsf{fma}\left(x, \log y, t + a\right)\right)\right)
\end{array}
Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
+-commutative99.9%
associate-+r+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-def99.9%
+-commutative99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
+-commutative99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x y z t a b c i) :precision binary64 (+ (+ (+ t a) (fma x (log y) z)) (+ (* (+ 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 ((t + a) + fma(x, log(y), z)) + (((b + -0.5) * log(c)) + (y * i));
}
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(t + a) + fma(x, log(y), z)) + Float64(Float64(Float64(b + -0.5) * log(c)) + Float64(y * i))) end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(t + a), $MachinePrecision] + N[(x * N[Log[y], $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(b + -0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(t + a\right) + \mathsf{fma}\left(x, \log y, z\right)\right) + \left(\left(b + -0.5\right) \cdot \log c + y \cdot i\right)
\end{array}
Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x y z t a b c i) :precision binary64 (+ (* y i) (+ (+ a (+ t (+ z (* x (log y))))) (* (log c) (- b 0.5)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (y * i) + ((a + (t + (z + (x * log(y))))) + (log(c) * (b - 0.5)));
}
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 = (y * i) + ((a + (t + (z + (x * log(y))))) + (log(c) * (b - 0.5d0)))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (y * i) + ((a + (t + (z + (x * Math.log(y))))) + (Math.log(c) * (b - 0.5)));
}
def code(x, y, z, t, a, b, c, i): return (y * i) + ((a + (t + (z + (x * math.log(y))))) + (math.log(c) * (b - 0.5)))
function code(x, y, z, t, a, b, c, i) return Float64(Float64(y * i) + Float64(Float64(a + Float64(t + Float64(z + Float64(x * log(y))))) + Float64(log(c) * Float64(b - 0.5)))) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = (y * i) + ((a + (t + (z + (x * log(y))))) + (log(c) * (b - 0.5))); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(y * i), $MachinePrecision] + N[(N[(a + N[(t + N[(z + N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
y \cdot i + \left(\left(a + \left(t + \left(z + x \cdot \log y\right)\right)\right) + \log c \cdot \left(b - 0.5\right)\right)
\end{array}
Initial program 99.9%
Final simplification99.9%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= x -2.1e+200) (not (<= x 2.4e+86))) (+ (+ (+ t a) (fma x (log y) z)) (* y i)) (+ a (+ t (+ z (fma y i (* (+ 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 ((x <= -2.1e+200) || !(x <= 2.4e+86)) {
tmp = ((t + a) + fma(x, log(y), z)) + (y * i);
} else {
tmp = a + (t + (z + fma(y, i, ((b + -0.5) * log(c)))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((x <= -2.1e+200) || !(x <= 2.4e+86)) tmp = Float64(Float64(Float64(t + a) + fma(x, log(y), z)) + Float64(y * i)); else tmp = Float64(a + Float64(t + Float64(z + fma(y, i, Float64(Float64(b + -0.5) * log(c)))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[x, -2.1e+200], N[Not[LessEqual[x, 2.4e+86]], $MachinePrecision]], N[(N[(N[(t + a), $MachinePrecision] + N[(x * N[Log[y], $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision], N[(a + N[(t + N[(z + N[(y * i + N[(N[(b + -0.5), $MachinePrecision] * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.1 \cdot 10^{+200} \lor \neg \left(x \leq 2.4 \cdot 10^{+86}\right):\\
\;\;\;\;\left(\left(t + a\right) + \mathsf{fma}\left(x, \log y, z\right)\right) + y \cdot i\\
\mathbf{else}:\\
\;\;\;\;a + \left(t + \left(z + \mathsf{fma}\left(y, i, \left(b + -0.5\right) \cdot \log c\right)\right)\right)\\
\end{array}
\end{array}
if x < -2.09999999999999997e200 or 2.4e86 < x Initial program 99.7%
associate-+l+99.7%
associate-+l+99.7%
fma-def99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in b around 0 89.3%
Taylor expanded in i around inf 89.3%
*-commutative89.3%
Simplified89.3%
if -2.09999999999999997e200 < x < 2.4e86Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 98.7%
*-commutative98.7%
sub-neg98.7%
metadata-eval98.7%
+-commutative98.7%
distribute-rgt-out98.7%
+-commutative98.7%
distribute-rgt-in98.7%
fma-udef98.7%
+-commutative98.7%
Simplified98.7%
Final simplification96.4%
(FPCore (x y z t a b c i) :precision binary64 (+ (* y i) (+ (+ a (+ t (+ z (* x (log y))))) (* b (log c)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (y * i) + ((a + (t + (z + (x * log(y))))) + (b * log(c)));
}
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 = (y * i) + ((a + (t + (z + (x * log(y))))) + (b * log(c)))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (y * i) + ((a + (t + (z + (x * Math.log(y))))) + (b * Math.log(c)));
}
def code(x, y, z, t, a, b, c, i): return (y * i) + ((a + (t + (z + (x * math.log(y))))) + (b * math.log(c)))
function code(x, y, z, t, a, b, c, i) return Float64(Float64(y * i) + Float64(Float64(a + Float64(t + Float64(z + Float64(x * log(y))))) + Float64(b * log(c)))) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = (y * i) + ((a + (t + (z + (x * log(y))))) + (b * log(c))); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(y * i), $MachinePrecision] + N[(N[(a + N[(t + N[(z + N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
y \cdot i + \left(\left(a + \left(t + \left(z + x \cdot \log y\right)\right)\right) + b \cdot \log c\right)
\end{array}
Initial program 99.9%
Taylor expanded in b around inf 97.9%
*-commutative97.9%
Simplified97.9%
Final simplification97.9%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= x -2.1e+200) (not (<= x 2.25e+86))) (+ (+ (+ t a) (fma x (log y) z)) (* y i)) (+ (* y i) (+ a (+ t (+ z (* (log c) (- b 0.5))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((x <= -2.1e+200) || !(x <= 2.25e+86)) {
tmp = ((t + a) + fma(x, log(y), z)) + (y * i);
} else {
tmp = (y * i) + (a + (t + (z + (log(c) * (b - 0.5)))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((x <= -2.1e+200) || !(x <= 2.25e+86)) tmp = Float64(Float64(Float64(t + a) + fma(x, log(y), z)) + Float64(y * i)); else tmp = Float64(Float64(y * i) + Float64(a + Float64(t + Float64(z + Float64(log(c) * Float64(b - 0.5)))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[x, -2.1e+200], N[Not[LessEqual[x, 2.25e+86]], $MachinePrecision]], N[(N[(N[(t + a), $MachinePrecision] + N[(x * N[Log[y], $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + N[(t + N[(z + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.1 \cdot 10^{+200} \lor \neg \left(x \leq 2.25 \cdot 10^{+86}\right):\\
\;\;\;\;\left(\left(t + a\right) + \mathsf{fma}\left(x, \log y, z\right)\right) + y \cdot i\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + \left(t + \left(z + \log c \cdot \left(b - 0.5\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < -2.09999999999999997e200 or 2.24999999999999996e86 < x Initial program 99.7%
associate-+l+99.7%
associate-+l+99.7%
fma-def99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in b around 0 89.3%
Taylor expanded in i around inf 89.3%
*-commutative89.3%
Simplified89.3%
if -2.09999999999999997e200 < x < 2.24999999999999996e86Initial program 99.9%
Taylor expanded in x around 0 98.7%
Final simplification96.3%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= x -7.5e+233) (not (<= x 3.2e+86))) (+ (* y i) (+ a (+ t (* x (log y))))) (+ (* y i) (+ a (+ t (+ z (* (log c) (- b 0.5))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((x <= -7.5e+233) || !(x <= 3.2e+86)) {
tmp = (y * i) + (a + (t + (x * log(y))));
} else {
tmp = (y * i) + (a + (t + (z + (log(c) * (b - 0.5)))));
}
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+233)) .or. (.not. (x <= 3.2d+86))) then
tmp = (y * i) + (a + (t + (x * log(y))))
else
tmp = (y * i) + (a + (t + (z + (log(c) * (b - 0.5d0)))))
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+233) || !(x <= 3.2e+86)) {
tmp = (y * i) + (a + (t + (x * Math.log(y))));
} else {
tmp = (y * i) + (a + (t + (z + (Math.log(c) * (b - 0.5)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (x <= -7.5e+233) or not (x <= 3.2e+86): tmp = (y * i) + (a + (t + (x * math.log(y)))) else: tmp = (y * i) + (a + (t + (z + (math.log(c) * (b - 0.5))))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((x <= -7.5e+233) || !(x <= 3.2e+86)) tmp = Float64(Float64(y * i) + Float64(a + Float64(t + Float64(x * log(y))))); else tmp = Float64(Float64(y * i) + Float64(a + Float64(t + Float64(z + Float64(log(c) * Float64(b - 0.5)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((x <= -7.5e+233) || ~((x <= 3.2e+86))) tmp = (y * i) + (a + (t + (x * log(y)))); else tmp = (y * i) + (a + (t + (z + (log(c) * (b - 0.5))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[x, -7.5e+233], N[Not[LessEqual[x, 3.2e+86]], $MachinePrecision]], N[(N[(y * i), $MachinePrecision] + N[(a + N[(t + N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + N[(t + N[(z + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7.5 \cdot 10^{+233} \lor \neg \left(x \leq 3.2 \cdot 10^{+86}\right):\\
\;\;\;\;y \cdot i + \left(a + \left(t + x \cdot \log y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + \left(t + \left(z + \log c \cdot \left(b - 0.5\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < -7.4999999999999997e233 or 3.2e86 < x Initial program 99.7%
associate-+l+99.7%
associate-+l+99.7%
fma-def99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in b around 0 88.0%
Taylor expanded in i around inf 88.0%
*-commutative88.0%
Simplified88.0%
Taylor expanded in z around 0 79.2%
+-commutative79.2%
Simplified79.2%
if -7.4999999999999997e233 < x < 3.2e86Initial program 99.9%
Taylor expanded in x around 0 97.7%
Final simplification93.6%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= x -4.3e+201) (not (<= x 3.2e+86))) (+ (* y i) (+ a (+ t (* x (log y))))) (+ (* y i) (+ a (+ z (* (log c) (- b 0.5)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((x <= -4.3e+201) || !(x <= 3.2e+86)) {
tmp = (y * i) + (a + (t + (x * log(y))));
} else {
tmp = (y * i) + (a + (z + (log(c) * (b - 0.5))));
}
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 <= (-4.3d+201)) .or. (.not. (x <= 3.2d+86))) then
tmp = (y * i) + (a + (t + (x * log(y))))
else
tmp = (y * i) + (a + (z + (log(c) * (b - 0.5d0))))
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 <= -4.3e+201) || !(x <= 3.2e+86)) {
tmp = (y * i) + (a + (t + (x * Math.log(y))));
} else {
tmp = (y * i) + (a + (z + (Math.log(c) * (b - 0.5))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (x <= -4.3e+201) or not (x <= 3.2e+86): tmp = (y * i) + (a + (t + (x * math.log(y)))) else: tmp = (y * i) + (a + (z + (math.log(c) * (b - 0.5)))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((x <= -4.3e+201) || !(x <= 3.2e+86)) tmp = Float64(Float64(y * i) + Float64(a + Float64(t + Float64(x * log(y))))); else tmp = Float64(Float64(y * i) + Float64(a + Float64(z + Float64(log(c) * Float64(b - 0.5))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((x <= -4.3e+201) || ~((x <= 3.2e+86))) tmp = (y * i) + (a + (t + (x * log(y)))); else tmp = (y * i) + (a + (z + (log(c) * (b - 0.5)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[x, -4.3e+201], N[Not[LessEqual[x, 3.2e+86]], $MachinePrecision]], N[(N[(y * i), $MachinePrecision] + N[(a + N[(t + N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -4.3 \cdot 10^{+201} \lor \neg \left(x \leq 3.2 \cdot 10^{+86}\right):\\
\;\;\;\;y \cdot i + \left(a + \left(t + x \cdot \log y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + \left(z + \log c \cdot \left(b - 0.5\right)\right)\right)\\
\end{array}
\end{array}
if x < -4.2999999999999999e201 or 3.2e86 < x Initial program 99.7%
associate-+l+99.7%
associate-+l+99.7%
fma-def99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in b around 0 89.3%
Taylor expanded in i around inf 89.3%
*-commutative89.3%
Simplified89.3%
Taylor expanded in z around 0 76.9%
+-commutative76.9%
Simplified76.9%
if -4.2999999999999999e201 < x < 3.2e86Initial program 99.9%
Taylor expanded in x around 0 98.7%
Taylor expanded in t around 0 80.0%
Final simplification79.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* (log c) (- b 0.5))))
(if (<= z -5.3e+182)
(+ (* y i) (+ a (+ z t)))
(if (<= z -8.5e+136) (+ a (+ z t_1)) (+ (* y i) (+ a (+ t t_1)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = log(c) * (b - 0.5);
double tmp;
if (z <= -5.3e+182) {
tmp = (y * i) + (a + (z + t));
} else if (z <= -8.5e+136) {
tmp = a + (z + t_1);
} else {
tmp = (y * i) + (a + (t + 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 = log(c) * (b - 0.5d0)
if (z <= (-5.3d+182)) then
tmp = (y * i) + (a + (z + t))
else if (z <= (-8.5d+136)) then
tmp = a + (z + t_1)
else
tmp = (y * i) + (a + (t + 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 = Math.log(c) * (b - 0.5);
double tmp;
if (z <= -5.3e+182) {
tmp = (y * i) + (a + (z + t));
} else if (z <= -8.5e+136) {
tmp = a + (z + t_1);
} else {
tmp = (y * i) + (a + (t + t_1));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = math.log(c) * (b - 0.5) tmp = 0 if z <= -5.3e+182: tmp = (y * i) + (a + (z + t)) elif z <= -8.5e+136: tmp = a + (z + t_1) else: tmp = (y * i) + (a + (t + t_1)) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(log(c) * Float64(b - 0.5)) tmp = 0.0 if (z <= -5.3e+182) tmp = Float64(Float64(y * i) + Float64(a + Float64(z + t))); elseif (z <= -8.5e+136) tmp = Float64(a + Float64(z + t_1)); else tmp = Float64(Float64(y * i) + Float64(a + Float64(t + t_1))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = log(c) * (b - 0.5); tmp = 0.0; if (z <= -5.3e+182) tmp = (y * i) + (a + (z + t)); elseif (z <= -8.5e+136) tmp = a + (z + t_1); else tmp = (y * i) + (a + (t + t_1)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -5.3e+182], N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -8.5e+136], N[(a + N[(z + t$95$1), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + N[(t + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log c \cdot \left(b - 0.5\right)\\
\mathbf{if}\;z \leq -5.3 \cdot 10^{+182}:\\
\;\;\;\;y \cdot i + \left(a + \left(z + t\right)\right)\\
\mathbf{elif}\;z \leq -8.5 \cdot 10^{+136}:\\
\;\;\;\;a + \left(z + t_1\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + \left(t + t_1\right)\right)\\
\end{array}
\end{array}
if z < -5.3e182Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 95.7%
Taylor expanded in i around inf 95.7%
*-commutative95.7%
Simplified95.7%
Taylor expanded in x around 0 92.5%
+-commutative92.5%
Simplified92.5%
if -5.3e182 < z < -8.49999999999999966e136Initial program 99.8%
Taylor expanded in x around 0 79.6%
Taylor expanded in t around 0 69.2%
Taylor expanded in y around 0 57.0%
if -8.49999999999999966e136 < z Initial program 99.8%
Taylor expanded in x around 0 84.7%
Taylor expanded in z around 0 71.7%
Final simplification73.6%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* (log c) (- b 0.5))))
(if (<= z -5.5e+181)
(+ (* y i) (+ a (+ z t)))
(if (<= z -2.4e+135) (+ a (+ z t_1)) (+ (* y i) (+ a t_1))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = log(c) * (b - 0.5);
double tmp;
if (z <= -5.5e+181) {
tmp = (y * i) + (a + (z + t));
} else if (z <= -2.4e+135) {
tmp = a + (z + t_1);
} else {
tmp = (y * i) + (a + 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 = log(c) * (b - 0.5d0)
if (z <= (-5.5d+181)) then
tmp = (y * i) + (a + (z + t))
else if (z <= (-2.4d+135)) then
tmp = a + (z + t_1)
else
tmp = (y * i) + (a + 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 = Math.log(c) * (b - 0.5);
double tmp;
if (z <= -5.5e+181) {
tmp = (y * i) + (a + (z + t));
} else if (z <= -2.4e+135) {
tmp = a + (z + t_1);
} else {
tmp = (y * i) + (a + t_1);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = math.log(c) * (b - 0.5) tmp = 0 if z <= -5.5e+181: tmp = (y * i) + (a + (z + t)) elif z <= -2.4e+135: tmp = a + (z + t_1) else: tmp = (y * i) + (a + t_1) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(log(c) * Float64(b - 0.5)) tmp = 0.0 if (z <= -5.5e+181) tmp = Float64(Float64(y * i) + Float64(a + Float64(z + t))); elseif (z <= -2.4e+135) tmp = Float64(a + Float64(z + t_1)); else tmp = Float64(Float64(y * i) + Float64(a + t_1)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = log(c) * (b - 0.5); tmp = 0.0; if (z <= -5.5e+181) tmp = (y * i) + (a + (z + t)); elseif (z <= -2.4e+135) tmp = a + (z + t_1); else tmp = (y * i) + (a + t_1); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -5.5e+181], N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -2.4e+135], N[(a + N[(z + t$95$1), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + t$95$1), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log c \cdot \left(b - 0.5\right)\\
\mathbf{if}\;z \leq -5.5 \cdot 10^{+181}:\\
\;\;\;\;y \cdot i + \left(a + \left(z + t\right)\right)\\
\mathbf{elif}\;z \leq -2.4 \cdot 10^{+135}:\\
\;\;\;\;a + \left(z + t_1\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + t_1\right)\\
\end{array}
\end{array}
if z < -5.49999999999999991e181Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 95.7%
Taylor expanded in i around inf 95.7%
*-commutative95.7%
Simplified95.7%
Taylor expanded in x around 0 92.5%
+-commutative92.5%
Simplified92.5%
if -5.49999999999999991e181 < z < -2.39999999999999997e135Initial program 99.8%
Taylor expanded in x around 0 79.6%
Taylor expanded in t around 0 69.2%
Taylor expanded in y around 0 57.0%
if -2.39999999999999997e135 < z Initial program 99.8%
Taylor expanded in x around 0 84.7%
Taylor expanded in t around 0 68.1%
Taylor expanded in z around 0 55.2%
Final simplification59.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* y i) (+ a (+ z t)))))
(if (<= i -1.02e+26)
t_1
(if (<= i 2.9e-11)
(+ a (+ t (+ z (* x (log y)))))
(if (<= i 1.32e+26) (+ a (* (log c) (- b 0.5))) 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) + (a + (z + t));
double tmp;
if (i <= -1.02e+26) {
tmp = t_1;
} else if (i <= 2.9e-11) {
tmp = a + (t + (z + (x * log(y))));
} else if (i <= 1.32e+26) {
tmp = a + (log(c) * (b - 0.5));
} 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) + (a + (z + t))
if (i <= (-1.02d+26)) then
tmp = t_1
else if (i <= 2.9d-11) then
tmp = a + (t + (z + (x * log(y))))
else if (i <= 1.32d+26) then
tmp = a + (log(c) * (b - 0.5d0))
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) + (a + (z + t));
double tmp;
if (i <= -1.02e+26) {
tmp = t_1;
} else if (i <= 2.9e-11) {
tmp = a + (t + (z + (x * Math.log(y))));
} else if (i <= 1.32e+26) {
tmp = a + (Math.log(c) * (b - 0.5));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * i) + (a + (z + t)) tmp = 0 if i <= -1.02e+26: tmp = t_1 elif i <= 2.9e-11: tmp = a + (t + (z + (x * math.log(y)))) elif i <= 1.32e+26: tmp = a + (math.log(c) * (b - 0.5)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(y * i) + Float64(a + Float64(z + t))) tmp = 0.0 if (i <= -1.02e+26) tmp = t_1; elseif (i <= 2.9e-11) tmp = Float64(a + Float64(t + Float64(z + Float64(x * log(y))))); elseif (i <= 1.32e+26) tmp = Float64(a + Float64(log(c) * Float64(b - 0.5))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (y * i) + (a + (z + t)); tmp = 0.0; if (i <= -1.02e+26) tmp = t_1; elseif (i <= 2.9e-11) tmp = a + (t + (z + (x * log(y)))); elseif (i <= 1.32e+26) tmp = a + (log(c) * (b - 0.5)); 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[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[i, -1.02e+26], t$95$1, If[LessEqual[i, 2.9e-11], N[(a + N[(t + N[(z + N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 1.32e+26], N[(a + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot i + \left(a + \left(z + t\right)\right)\\
\mathbf{if}\;i \leq -1.02 \cdot 10^{+26}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;i \leq 2.9 \cdot 10^{-11}:\\
\;\;\;\;a + \left(t + \left(z + x \cdot \log y\right)\right)\\
\mathbf{elif}\;i \leq 1.32 \cdot 10^{+26}:\\
\;\;\;\;a + \log c \cdot \left(b - 0.5\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if i < -1.0200000000000001e26 or 1.32e26 < i Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 87.6%
Taylor expanded in i around inf 87.6%
*-commutative87.6%
Simplified87.6%
Taylor expanded in x around 0 78.3%
+-commutative78.3%
Simplified78.3%
if -1.0200000000000001e26 < i < 2.9e-11Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
fma-def99.8%
sub-neg99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in b around 0 80.8%
Taylor expanded in i around inf 77.9%
*-commutative77.9%
Simplified77.9%
Taylor expanded in y around 0 75.6%
if 2.9e-11 < i < 1.32e26Initial program 99.7%
Taylor expanded in x around 0 89.7%
Taylor expanded in t around 0 89.6%
Taylor expanded in y around 0 80.2%
Taylor expanded in z around 0 70.8%
Final simplification76.7%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= i -2.7e+26) (not (<= i 2.5e+115))) (+ (* y i) (+ a (+ z t))) (+ a (+ z (* (log c) (- b 0.5))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((i <= -2.7e+26) || !(i <= 2.5e+115)) {
tmp = (y * i) + (a + (z + t));
} else {
tmp = a + (z + (log(c) * (b - 0.5)));
}
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 ((i <= (-2.7d+26)) .or. (.not. (i <= 2.5d+115))) then
tmp = (y * i) + (a + (z + t))
else
tmp = a + (z + (log(c) * (b - 0.5d0)))
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 ((i <= -2.7e+26) || !(i <= 2.5e+115)) {
tmp = (y * i) + (a + (z + t));
} else {
tmp = a + (z + (Math.log(c) * (b - 0.5)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (i <= -2.7e+26) or not (i <= 2.5e+115): tmp = (y * i) + (a + (z + t)) else: tmp = a + (z + (math.log(c) * (b - 0.5))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((i <= -2.7e+26) || !(i <= 2.5e+115)) tmp = Float64(Float64(y * i) + Float64(a + Float64(z + t))); else tmp = Float64(a + Float64(z + Float64(log(c) * Float64(b - 0.5)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((i <= -2.7e+26) || ~((i <= 2.5e+115))) tmp = (y * i) + (a + (z + t)); else tmp = a + (z + (log(c) * (b - 0.5))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[i, -2.7e+26], N[Not[LessEqual[i, 2.5e+115]], $MachinePrecision]], N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a + N[(z + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -2.7 \cdot 10^{+26} \lor \neg \left(i \leq 2.5 \cdot 10^{+115}\right):\\
\;\;\;\;y \cdot i + \left(a + \left(z + t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;a + \left(z + \log c \cdot \left(b - 0.5\right)\right)\\
\end{array}
\end{array}
if i < -2.7e26 or 2.50000000000000004e115 < i Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 89.6%
Taylor expanded in i around inf 89.6%
*-commutative89.6%
Simplified89.6%
Taylor expanded in x around 0 79.8%
+-commutative79.8%
Simplified79.8%
if -2.7e26 < i < 2.50000000000000004e115Initial program 99.8%
Taylor expanded in x around 0 82.9%
Taylor expanded in t around 0 63.9%
Taylor expanded in y around 0 60.4%
Final simplification68.5%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= t -8.5e-253) (not (<= t 7.1e-177))) (+ (* y i) (+ a (+ z t))) (+ a (* (log c) (- b 0.5)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((t <= -8.5e-253) || !(t <= 7.1e-177)) {
tmp = (y * i) + (a + (z + t));
} else {
tmp = a + (log(c) * (b - 0.5));
}
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 ((t <= (-8.5d-253)) .or. (.not. (t <= 7.1d-177))) then
tmp = (y * i) + (a + (z + t))
else
tmp = a + (log(c) * (b - 0.5d0))
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 ((t <= -8.5e-253) || !(t <= 7.1e-177)) {
tmp = (y * i) + (a + (z + t));
} else {
tmp = a + (Math.log(c) * (b - 0.5));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (t <= -8.5e-253) or not (t <= 7.1e-177): tmp = (y * i) + (a + (z + t)) else: tmp = a + (math.log(c) * (b - 0.5)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((t <= -8.5e-253) || !(t <= 7.1e-177)) tmp = Float64(Float64(y * i) + Float64(a + Float64(z + t))); else tmp = Float64(a + Float64(log(c) * Float64(b - 0.5))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((t <= -8.5e-253) || ~((t <= 7.1e-177))) tmp = (y * i) + (a + (z + t)); else tmp = a + (log(c) * (b - 0.5)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[t, -8.5e-253], N[Not[LessEqual[t, 7.1e-177]], $MachinePrecision]], N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -8.5 \cdot 10^{-253} \lor \neg \left(t \leq 7.1 \cdot 10^{-177}\right):\\
\;\;\;\;y \cdot i + \left(a + \left(z + t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;a + \log c \cdot \left(b - 0.5\right)\\
\end{array}
\end{array}
if t < -8.4999999999999999e-253 or 7.10000000000000043e-177 < t Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 84.4%
Taylor expanded in i around inf 83.5%
*-commutative83.5%
Simplified83.5%
Taylor expanded in x around 0 70.2%
+-commutative70.2%
Simplified70.2%
if -8.4999999999999999e-253 < t < 7.10000000000000043e-177Initial program 99.8%
Taylor expanded in x around 0 83.2%
Taylor expanded in t around 0 83.2%
Taylor expanded in y around 0 75.0%
Taylor expanded in z around 0 58.2%
Final simplification68.3%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= x -8.8e+234) (not (<= x 2.2e+86))) (+ (* y i) (* x (log y))) (+ (* y i) (+ a (+ z t)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((x <= -8.8e+234) || !(x <= 2.2e+86)) {
tmp = (y * i) + (x * log(y));
} else {
tmp = (y * i) + (a + (z + t));
}
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 <= (-8.8d+234)) .or. (.not. (x <= 2.2d+86))) then
tmp = (y * i) + (x * log(y))
else
tmp = (y * i) + (a + (z + t))
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 <= -8.8e+234) || !(x <= 2.2e+86)) {
tmp = (y * i) + (x * Math.log(y));
} else {
tmp = (y * i) + (a + (z + t));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (x <= -8.8e+234) or not (x <= 2.2e+86): tmp = (y * i) + (x * math.log(y)) else: tmp = (y * i) + (a + (z + t)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((x <= -8.8e+234) || !(x <= 2.2e+86)) tmp = Float64(Float64(y * i) + Float64(x * log(y))); else tmp = Float64(Float64(y * i) + Float64(a + Float64(z + t))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((x <= -8.8e+234) || ~((x <= 2.2e+86))) tmp = (y * i) + (x * log(y)); else tmp = (y * i) + (a + (z + t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[x, -8.8e+234], N[Not[LessEqual[x, 2.2e+86]], $MachinePrecision]], N[(N[(y * i), $MachinePrecision] + N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -8.8 \cdot 10^{+234} \lor \neg \left(x \leq 2.2 \cdot 10^{+86}\right):\\
\;\;\;\;y \cdot i + x \cdot \log y\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + \left(z + t\right)\right)\\
\end{array}
\end{array}
if x < -8.8000000000000003e234 or 2.20000000000000003e86 < x Initial program 99.7%
Taylor expanded in x around inf 69.2%
if -8.8000000000000003e234 < x < 2.20000000000000003e86Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 81.8%
Taylor expanded in i around inf 79.5%
*-commutative79.5%
Simplified79.5%
Taylor expanded in x around 0 77.5%
+-commutative77.5%
Simplified77.5%
Final simplification75.7%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= b -3.5e+155) (not (<= b 2.55e+178))) (+ a (* b (log c))) (+ (* y i) (+ a (+ z t)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((b <= -3.5e+155) || !(b <= 2.55e+178)) {
tmp = a + (b * log(c));
} else {
tmp = (y * i) + (a + (z + t));
}
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 ((b <= (-3.5d+155)) .or. (.not. (b <= 2.55d+178))) then
tmp = a + (b * log(c))
else
tmp = (y * i) + (a + (z + t))
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 ((b <= -3.5e+155) || !(b <= 2.55e+178)) {
tmp = a + (b * Math.log(c));
} else {
tmp = (y * i) + (a + (z + t));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (b <= -3.5e+155) or not (b <= 2.55e+178): tmp = a + (b * math.log(c)) else: tmp = (y * i) + (a + (z + t)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((b <= -3.5e+155) || !(b <= 2.55e+178)) tmp = Float64(a + Float64(b * log(c))); else tmp = Float64(Float64(y * i) + Float64(a + Float64(z + t))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((b <= -3.5e+155) || ~((b <= 2.55e+178))) tmp = a + (b * log(c)); else tmp = (y * i) + (a + (z + t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[b, -3.5e+155], N[Not[LessEqual[b, 2.55e+178]], $MachinePrecision]], N[(a + N[(b * N[Log[c], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.5 \cdot 10^{+155} \lor \neg \left(b \leq 2.55 \cdot 10^{+178}\right):\\
\;\;\;\;a + b \cdot \log c\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + \left(a + \left(z + t\right)\right)\\
\end{array}
\end{array}
if b < -3.49999999999999985e155 or 2.5499999999999998e178 < b Initial program 99.8%
Taylor expanded in x around 0 96.5%
Taylor expanded in t around 0 85.8%
Taylor expanded in y around 0 71.0%
Taylor expanded in b around inf 59.6%
*-commutative59.6%
Simplified59.6%
if -3.49999999999999985e155 < b < 2.5499999999999998e178Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 96.0%
Taylor expanded in i around inf 93.7%
*-commutative93.7%
Simplified93.7%
Taylor expanded in x around 0 76.2%
+-commutative76.2%
Simplified76.2%
Final simplification72.0%
(FPCore (x y z t a b c i)
:precision binary64
(if (or (<= i -9e+192)
(not
(or (<= i -2.05e+161)
(and (not (<= i -2.45e+114)) (<= i 4.7e+112)))))
(* y i)
(+ z a)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((i <= -9e+192) || !((i <= -2.05e+161) || (!(i <= -2.45e+114) && (i <= 4.7e+112)))) {
tmp = y * i;
} else {
tmp = 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) :: tmp
if ((i <= (-9d+192)) .or. (.not. (i <= (-2.05d+161)) .or. (.not. (i <= (-2.45d+114))) .and. (i <= 4.7d+112))) then
tmp = y * i
else
tmp = 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 tmp;
if ((i <= -9e+192) || !((i <= -2.05e+161) || (!(i <= -2.45e+114) && (i <= 4.7e+112)))) {
tmp = y * i;
} else {
tmp = z + a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (i <= -9e+192) or not ((i <= -2.05e+161) or (not (i <= -2.45e+114) and (i <= 4.7e+112))): tmp = y * i else: tmp = z + a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((i <= -9e+192) || !((i <= -2.05e+161) || (!(i <= -2.45e+114) && (i <= 4.7e+112)))) tmp = Float64(y * i); else tmp = Float64(z + a); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((i <= -9e+192) || ~(((i <= -2.05e+161) || (~((i <= -2.45e+114)) && (i <= 4.7e+112))))) tmp = y * i; else tmp = z + a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[i, -9e+192], N[Not[Or[LessEqual[i, -2.05e+161], And[N[Not[LessEqual[i, -2.45e+114]], $MachinePrecision], LessEqual[i, 4.7e+112]]]], $MachinePrecision]], N[(y * i), $MachinePrecision], N[(z + a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -9 \cdot 10^{+192} \lor \neg \left(i \leq -2.05 \cdot 10^{+161} \lor \neg \left(i \leq -2.45 \cdot 10^{+114}\right) \land i \leq 4.7 \cdot 10^{+112}\right):\\
\;\;\;\;y \cdot i\\
\mathbf{else}:\\
\;\;\;\;z + a\\
\end{array}
\end{array}
if i < -9e192 or -2.0500000000000001e161 < i < -2.45e114 or 4.69999999999999997e112 < i Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 59.7%
*-commutative59.7%
Simplified59.7%
if -9e192 < i < -2.0500000000000001e161 or -2.45e114 < i < 4.69999999999999997e112Initial program 99.8%
Taylor expanded in x around 0 83.8%
Taylor expanded in t around 0 65.6%
Taylor expanded in y around 0 59.8%
Taylor expanded in z around inf 38.9%
Final simplification45.0%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= i -8.6e+29) (not (<= i 4.2e+107))) (+ a (* y i)) (+ z a)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((i <= -8.6e+29) || !(i <= 4.2e+107)) {
tmp = a + (y * i);
} else {
tmp = 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) :: tmp
if ((i <= (-8.6d+29)) .or. (.not. (i <= 4.2d+107))) then
tmp = a + (y * i)
else
tmp = 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 tmp;
if ((i <= -8.6e+29) || !(i <= 4.2e+107)) {
tmp = a + (y * i);
} else {
tmp = z + a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (i <= -8.6e+29) or not (i <= 4.2e+107): tmp = a + (y * i) else: tmp = z + a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((i <= -8.6e+29) || !(i <= 4.2e+107)) tmp = Float64(a + Float64(y * i)); else tmp = Float64(z + a); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((i <= -8.6e+29) || ~((i <= 4.2e+107))) tmp = a + (y * i); else tmp = z + a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[i, -8.6e+29], N[Not[LessEqual[i, 4.2e+107]], $MachinePrecision]], N[(a + N[(y * i), $MachinePrecision]), $MachinePrecision], N[(z + a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -8.6 \cdot 10^{+29} \lor \neg \left(i \leq 4.2 \cdot 10^{+107}\right):\\
\;\;\;\;a + y \cdot i\\
\mathbf{else}:\\
\;\;\;\;z + a\\
\end{array}
\end{array}
if i < -8.6000000000000006e29 or 4.1999999999999999e107 < i Initial program 99.9%
Taylor expanded in a around inf 57.0%
if -8.6000000000000006e29 < i < 4.1999999999999999e107Initial program 99.8%
Taylor expanded in x around 0 82.8%
Taylor expanded in t around 0 63.7%
Taylor expanded in y around 0 60.2%
Taylor expanded in z around inf 37.5%
Final simplification45.7%
(FPCore (x y z t a b c i) :precision binary64 (+ (* y i) (+ a (+ z t))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (y * i) + (a + (z + t));
}
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 = (y * i) + (a + (z + t))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (y * i) + (a + (z + t));
}
def code(x, y, z, t, a, b, c, i): return (y * i) + (a + (z + t))
function code(x, y, z, t, a, b, c, i) return Float64(Float64(y * i) + Float64(a + Float64(z + t))) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = (y * i) + (a + (z + t)); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(y * i), $MachinePrecision] + N[(a + N[(z + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
y \cdot i + \left(a + \left(z + t\right)\right)
\end{array}
Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
fma-def99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in b around 0 82.9%
Taylor expanded in i around inf 81.2%
*-commutative81.2%
Simplified81.2%
Taylor expanded in x around 0 67.4%
+-commutative67.4%
Simplified67.4%
Final simplification67.4%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 2.85e-123) z (if (<= a 8.5e+156) (* y i) a)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 2.85e-123) {
tmp = z;
} else if (a <= 8.5e+156) {
tmp = y * i;
} 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 (a <= 2.85d-123) then
tmp = z
else if (a <= 8.5d+156) then
tmp = y * i
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 (a <= 2.85e-123) {
tmp = z;
} else if (a <= 8.5e+156) {
tmp = y * i;
} else {
tmp = a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 2.85e-123: tmp = z elif a <= 8.5e+156: tmp = y * i else: tmp = a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 2.85e-123) tmp = z; elseif (a <= 8.5e+156) tmp = Float64(y * i); else tmp = a; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (a <= 2.85e-123) tmp = z; elseif (a <= 8.5e+156) tmp = y * i; else tmp = a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 2.85e-123], z, If[LessEqual[a, 8.5e+156], N[(y * i), $MachinePrecision], a]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 2.85 \cdot 10^{-123}:\\
\;\;\;\;z\\
\mathbf{elif}\;a \leq 8.5 \cdot 10^{+156}:\\
\;\;\;\;y \cdot i\\
\mathbf{else}:\\
\;\;\;\;a\\
\end{array}
\end{array}
if a < 2.85000000000000014e-123Initial program 99.8%
Taylor expanded in x around 0 85.6%
Taylor expanded in t around 0 69.2%
Taylor expanded in y around 0 47.9%
Taylor expanded in z around inf 18.0%
if 2.85000000000000014e-123 < a < 8.49999999999999948e156Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
fma-def99.8%
sub-neg99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in y around inf 28.1%
*-commutative28.1%
Simplified28.1%
if 8.49999999999999948e156 < a Initial program 100.0%
Taylor expanded in x around 0 88.2%
Taylor expanded in t around 0 84.6%
Taylor expanded in y around 0 70.9%
Taylor expanded in a around inf 61.1%
Final simplification24.8%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 2.25e+151) (+ z (* y i)) (+ a (* y i))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 2.25e+151) {
tmp = z + (y * i);
} 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 (a <= 2.25d+151) then
tmp = z + (y * i)
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 (a <= 2.25e+151) {
tmp = z + (y * i);
} else {
tmp = a + (y * i);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 2.25e+151: tmp = z + (y * i) else: tmp = a + (y * i) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 2.25e+151) tmp = Float64(z + Float64(y * i)); 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 (a <= 2.25e+151) tmp = z + (y * i); else tmp = a + (y * i); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 2.25e+151], N[(z + N[(y * i), $MachinePrecision]), $MachinePrecision], N[(a + N[(y * i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 2.25 \cdot 10^{+151}:\\
\;\;\;\;z + y \cdot i\\
\mathbf{else}:\\
\;\;\;\;a + y \cdot i\\
\end{array}
\end{array}
if a < 2.2499999999999999e151Initial program 99.8%
Taylor expanded in z around inf 40.6%
if 2.2499999999999999e151 < a Initial program 100.0%
Taylor expanded in a around inf 76.3%
Final simplification44.7%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 4.6e+151) z a))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 4.6e+151) {
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 (a <= 4.6d+151) 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 (a <= 4.6e+151) {
tmp = z;
} else {
tmp = a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 4.6e+151: tmp = z else: tmp = a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 4.6e+151) 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 (a <= 4.6e+151) tmp = z; else tmp = a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 4.6e+151], z, a]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 4.6 \cdot 10^{+151}:\\
\;\;\;\;z\\
\mathbf{else}:\\
\;\;\;\;a\\
\end{array}
\end{array}
if a < 4.6000000000000002e151Initial program 99.8%
Taylor expanded in x around 0 85.6%
Taylor expanded in t around 0 68.8%
Taylor expanded in y around 0 46.6%
Taylor expanded in z around inf 18.6%
if 4.6000000000000002e151 < a Initial program 100.0%
Taylor expanded in x around 0 89.0%
Taylor expanded in t around 0 85.7%
Taylor expanded in y around 0 69.5%
Taylor expanded in a around inf 60.3%
Final simplification23.5%
(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.9%
Taylor expanded in x around 0 86.0%
Taylor expanded in t around 0 70.8%
Taylor expanded in y around 0 49.3%
Taylor expanded in a around inf 16.2%
Final simplification16.2%
herbie shell --seed 2024011
(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)))