
(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 15 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%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* (log c) (- b 0.5))) (t_2 (+ z t_1)))
(if (<= t_1 -1e+49)
(+ a t_2)
(if (<= t_1 2e+215)
(+ a (+ t (+ z (+ (* x (log y)) (* y i)))))
(+ (* y i) t_2)))))
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 t_2 = z + t_1;
double tmp;
if (t_1 <= -1e+49) {
tmp = a + t_2;
} else if (t_1 <= 2e+215) {
tmp = a + (t + (z + ((x * log(y)) + (y * i))));
} else {
tmp = (y * i) + 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 = log(c) * (b - 0.5d0)
t_2 = z + t_1
if (t_1 <= (-1d+49)) then
tmp = a + t_2
else if (t_1 <= 2d+215) then
tmp = a + (t + (z + ((x * log(y)) + (y * i))))
else
tmp = (y * i) + 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 = Math.log(c) * (b - 0.5);
double t_2 = z + t_1;
double tmp;
if (t_1 <= -1e+49) {
tmp = a + t_2;
} else if (t_1 <= 2e+215) {
tmp = a + (t + (z + ((x * Math.log(y)) + (y * i))));
} else {
tmp = (y * i) + t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = math.log(c) * (b - 0.5) t_2 = z + t_1 tmp = 0 if t_1 <= -1e+49: tmp = a + t_2 elif t_1 <= 2e+215: tmp = a + (t + (z + ((x * math.log(y)) + (y * i)))) else: tmp = (y * i) + t_2 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(log(c) * Float64(b - 0.5)) t_2 = Float64(z + t_1) tmp = 0.0 if (t_1 <= -1e+49) tmp = Float64(a + t_2); elseif (t_1 <= 2e+215) tmp = Float64(a + Float64(t + Float64(z + Float64(Float64(x * log(y)) + Float64(y * i))))); else tmp = Float64(Float64(y * i) + t_2); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = log(c) * (b - 0.5); t_2 = z + t_1; tmp = 0.0; if (t_1 <= -1e+49) tmp = a + t_2; elseif (t_1 <= 2e+215) tmp = a + (t + (z + ((x * log(y)) + (y * i)))); else tmp = (y * i) + t_2; 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]}, Block[{t$95$2 = N[(z + t$95$1), $MachinePrecision]}, If[LessEqual[t$95$1, -1e+49], N[(a + t$95$2), $MachinePrecision], If[LessEqual[t$95$1, 2e+215], N[(a + N[(t + N[(z + N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * i), $MachinePrecision] + t$95$2), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log c \cdot \left(b - 0.5\right)\\
t_2 := z + t\_1\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{+49}:\\
\;\;\;\;a + t\_2\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{+215}:\\
\;\;\;\;a + \left(t + \left(z + \left(x \cdot \log y + y \cdot i\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot i + t\_2\\
\end{array}
\end{array}
if (*.f64 (-.f64 b #s(literal 1/2 binary64)) (log.f64 c)) < -9.99999999999999946e48Initial program 99.8%
associate-+l+99.8%
associate-+l+99.8%
+-commutative99.8%
associate-+l+99.8%
+-commutative99.8%
associate-+l+99.8%
+-commutative99.8%
associate-+l+99.8%
+-commutative99.8%
fma-define99.8%
+-commutative99.8%
fma-define99.8%
Simplified99.8%
Taylor expanded in i around inf 60.5%
associate-+r+60.5%
associate-+r+60.5%
associate-/l*60.5%
sub-neg60.5%
metadata-eval60.5%
associate-/l*60.4%
+-commutative60.4%
Simplified60.4%
Taylor expanded in t around 0 54.3%
Taylor expanded in i around 0 75.7%
Taylor expanded in x around 0 65.9%
if -9.99999999999999946e48 < (*.f64 (-.f64 b #s(literal 1/2 binary64)) (log.f64 c)) < 1.99999999999999981e215Initial program 99.9%
Taylor expanded in b around inf 97.3%
*-commutative97.3%
Simplified97.3%
Taylor expanded in b around 0 94.1%
if 1.99999999999999981e215 < (*.f64 (-.f64 b #s(literal 1/2 binary64)) (log.f64 c)) Initial program 99.6%
add-cube-cbrt99.5%
pow399.5%
Applied egg-rr99.5%
Taylor expanded in z around inf 74.1%
Final simplification86.5%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= x -3.6e+95) (not (<= x 3.6e+174))) (+ a (+ t (+ z (+ (* x (log y)) (* y i))))) (fma 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 <= -3.6e+95) || !(x <= 3.6e+174)) {
tmp = a + (t + (z + ((x * log(y)) + (y * i))));
} else {
tmp = fma(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 <= -3.6e+95) || !(x <= 3.6e+174)) tmp = Float64(a + Float64(t + Float64(z + Float64(Float64(x * log(y)) + Float64(y * i))))); else tmp = fma(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, -3.6e+95], N[Not[LessEqual[x, 3.6e+174]], $MachinePrecision]], N[(a + N[(t + N[(z + N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y * i + 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 -3.6 \cdot 10^{+95} \lor \neg \left(x \leq 3.6 \cdot 10^{+174}\right):\\
\;\;\;\;a + \left(t + \left(z + \left(x \cdot \log y + y \cdot i\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y, i, a + \left(t + \left(z + \log c \cdot \left(b - 0.5\right)\right)\right)\right)\\
\end{array}
\end{array}
if x < -3.59999999999999978e95 or 3.6000000000000002e174 < x Initial program 99.9%
Taylor expanded in b around inf 99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in b around 0 94.8%
if -3.59999999999999978e95 < x < 3.6000000000000002e174Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around 0 99.4%
Final simplification98.3%
(FPCore (x y z t a b c i) :precision binary64 (+ (+ (+ a (+ t (+ z (* x (log y))))) (* (log c) (- b 0.5))) (* y i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((a + (t + (z + (x * log(y))))) + (log(c) * (b - 0.5))) + (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 = ((a + (t + (z + (x * log(y))))) + (log(c) * (b - 0.5d0))) + (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 ((a + (t + (z + (x * Math.log(y))))) + (Math.log(c) * (b - 0.5))) + (y * i);
}
def code(x, y, z, t, a, b, c, i): return ((a + (t + (z + (x * math.log(y))))) + (math.log(c) * (b - 0.5))) + (y * i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(a + Float64(t + Float64(z + Float64(x * log(y))))) + Float64(log(c) * Float64(b - 0.5))) + Float64(y * i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = ((a + (t + (z + (x * log(y))))) + (log(c) * (b - 0.5))) + (y * i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(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] + N[(y * i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(a + \left(t + \left(z + x \cdot \log y\right)\right)\right) + \log c \cdot \left(b - 0.5\right)\right) + y \cdot i
\end{array}
Initial program 99.9%
Final simplification99.9%
(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 98.0%
*-commutative98.0%
Simplified98.0%
Final simplification98.0%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= x -2.3e+96) (not (<= x 1.45e+174))) (+ a (+ t (+ z (+ (* x (log y)) (* y i))))) (+ a (+ t (+ z (+ (* (log c) (- b 0.5)) (* y i)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((x <= -2.3e+96) || !(x <= 1.45e+174)) {
tmp = a + (t + (z + ((x * log(y)) + (y * i))));
} else {
tmp = a + (t + (z + ((log(c) * (b - 0.5)) + (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 ((x <= (-2.3d+96)) .or. (.not. (x <= 1.45d+174))) then
tmp = a + (t + (z + ((x * log(y)) + (y * i))))
else
tmp = a + (t + (z + ((log(c) * (b - 0.5d0)) + (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 ((x <= -2.3e+96) || !(x <= 1.45e+174)) {
tmp = a + (t + (z + ((x * Math.log(y)) + (y * i))));
} else {
tmp = a + (t + (z + ((Math.log(c) * (b - 0.5)) + (y * i))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (x <= -2.3e+96) or not (x <= 1.45e+174): tmp = a + (t + (z + ((x * math.log(y)) + (y * i)))) else: tmp = a + (t + (z + ((math.log(c) * (b - 0.5)) + (y * i)))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((x <= -2.3e+96) || !(x <= 1.45e+174)) tmp = Float64(a + Float64(t + Float64(z + Float64(Float64(x * log(y)) + Float64(y * i))))); else tmp = Float64(a + Float64(t + Float64(z + Float64(Float64(log(c) * Float64(b - 0.5)) + Float64(y * i))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((x <= -2.3e+96) || ~((x <= 1.45e+174))) tmp = a + (t + (z + ((x * log(y)) + (y * i)))); else tmp = a + (t + (z + ((log(c) * (b - 0.5)) + (y * i)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[x, -2.3e+96], N[Not[LessEqual[x, 1.45e+174]], $MachinePrecision]], N[(a + N[(t + N[(z + N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a + N[(t + N[(z + N[(N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision] + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.3 \cdot 10^{+96} \lor \neg \left(x \leq 1.45 \cdot 10^{+174}\right):\\
\;\;\;\;a + \left(t + \left(z + \left(x \cdot \log y + y \cdot i\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;a + \left(t + \left(z + \left(\log c \cdot \left(b - 0.5\right) + y \cdot i\right)\right)\right)\\
\end{array}
\end{array}
if x < -2.30000000000000015e96 or 1.45e174 < x Initial program 99.9%
Taylor expanded in b around inf 99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in b around 0 94.8%
if -2.30000000000000015e96 < x < 1.45e174Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around 0 99.4%
Final simplification98.2%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= i -7.2e+167) (not (<= i 4.4e+166))) (* a (+ (* i (/ y a)) 1.0)) (+ 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 ((i <= -7.2e+167) || !(i <= 4.4e+166)) {
tmp = a * ((i * (y / a)) + 1.0);
} else {
tmp = 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 ((i <= (-7.2d+167)) .or. (.not. (i <= 4.4d+166))) then
tmp = a * ((i * (y / a)) + 1.0d0)
else
tmp = 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 ((i <= -7.2e+167) || !(i <= 4.4e+166)) {
tmp = a * ((i * (y / a)) + 1.0);
} else {
tmp = a + (t + (z + (Math.log(c) * (b - 0.5))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (i <= -7.2e+167) or not (i <= 4.4e+166): tmp = a * ((i * (y / a)) + 1.0) else: tmp = 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 ((i <= -7.2e+167) || !(i <= 4.4e+166)) tmp = Float64(a * Float64(Float64(i * Float64(y / a)) + 1.0)); else tmp = 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 ((i <= -7.2e+167) || ~((i <= 4.4e+166))) tmp = a * ((i * (y / a)) + 1.0); else tmp = 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[i, -7.2e+167], N[Not[LessEqual[i, 4.4e+166]], $MachinePrecision]], N[(a * N[(N[(i * N[(y / a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(a + N[(t + N[(z + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -7.2 \cdot 10^{+167} \lor \neg \left(i \leq 4.4 \cdot 10^{+166}\right):\\
\;\;\;\;a \cdot \left(i \cdot \frac{y}{a} + 1\right)\\
\mathbf{else}:\\
\;\;\;\;a + \left(t + \left(z + \log c \cdot \left(b - 0.5\right)\right)\right)\\
\end{array}
\end{array}
if i < -7.20000000000000049e167 or 4.3999999999999998e166 < i Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define100.0%
+-commutative100.0%
fma-define100.0%
Simplified100.0%
Taylor expanded in a around -inf 81.4%
Taylor expanded in i around inf 66.4%
associate-/l*66.3%
Simplified66.3%
if -7.20000000000000049e167 < i < 4.3999999999999998e166Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in x around 0 86.6%
Taylor expanded in y around 0 75.7%
Final simplification73.5%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= i -1.05e+168) (not (<= i 6e+166))) (* a (+ (* i (/ y a)) 1.0)) (+ 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 <= -1.05e+168) || !(i <= 6e+166)) {
tmp = a * ((i * (y / a)) + 1.0);
} 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 <= (-1.05d+168)) .or. (.not. (i <= 6d+166))) then
tmp = a * ((i * (y / a)) + 1.0d0)
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 <= -1.05e+168) || !(i <= 6e+166)) {
tmp = a * ((i * (y / a)) + 1.0);
} 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 <= -1.05e+168) or not (i <= 6e+166): tmp = a * ((i * (y / a)) + 1.0) 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 <= -1.05e+168) || !(i <= 6e+166)) tmp = Float64(a * Float64(Float64(i * Float64(y / a)) + 1.0)); 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 <= -1.05e+168) || ~((i <= 6e+166))) tmp = a * ((i * (y / a)) + 1.0); 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, -1.05e+168], N[Not[LessEqual[i, 6e+166]], $MachinePrecision]], N[(a * N[(N[(i * N[(y / a), $MachinePrecision]), $MachinePrecision] + 1.0), $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 -1.05 \cdot 10^{+168} \lor \neg \left(i \leq 6 \cdot 10^{+166}\right):\\
\;\;\;\;a \cdot \left(i \cdot \frac{y}{a} + 1\right)\\
\mathbf{else}:\\
\;\;\;\;a + \left(z + \log c \cdot \left(b - 0.5\right)\right)\\
\end{array}
\end{array}
if i < -1.05000000000000001e168 or 5.99999999999999997e166 < i Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define100.0%
+-commutative100.0%
fma-define100.0%
Simplified100.0%
Taylor expanded in a around -inf 81.4%
Taylor expanded in i around inf 66.4%
associate-/l*66.3%
Simplified66.3%
if -1.05000000000000001e168 < i < 5.99999999999999997e166Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in i around inf 55.7%
associate-+r+55.7%
associate-+r+55.7%
associate-/l*55.7%
sub-neg55.7%
metadata-eval55.7%
associate-/l*55.7%
+-commutative55.7%
Simplified55.7%
Taylor expanded in t around 0 50.4%
Taylor expanded in i around 0 75.5%
Taylor expanded in x around 0 62.7%
Final simplification63.5%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 2.9e+149) (+ (* y i) (+ z (* (log c) (- b 0.5)))) (* a (+ (* i (/ y a)) 1.0))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 2.9e+149) {
tmp = (y * i) + (z + (log(c) * (b - 0.5)));
} else {
tmp = a * ((i * (y / a)) + 1.0);
}
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.9d+149) then
tmp = (y * i) + (z + (log(c) * (b - 0.5d0)))
else
tmp = a * ((i * (y / a)) + 1.0d0)
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.9e+149) {
tmp = (y * i) + (z + (Math.log(c) * (b - 0.5)));
} else {
tmp = a * ((i * (y / a)) + 1.0);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 2.9e+149: tmp = (y * i) + (z + (math.log(c) * (b - 0.5))) else: tmp = a * ((i * (y / a)) + 1.0) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 2.9e+149) tmp = Float64(Float64(y * i) + Float64(z + Float64(log(c) * Float64(b - 0.5)))); else tmp = Float64(a * Float64(Float64(i * Float64(y / a)) + 1.0)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (a <= 2.9e+149) tmp = (y * i) + (z + (log(c) * (b - 0.5))); else tmp = a * ((i * (y / a)) + 1.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 2.9e+149], N[(N[(y * i), $MachinePrecision] + N[(z + N[(N[Log[c], $MachinePrecision] * N[(b - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a * N[(N[(i * N[(y / a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 2.9 \cdot 10^{+149}:\\
\;\;\;\;y \cdot i + \left(z + \log c \cdot \left(b - 0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(i \cdot \frac{y}{a} + 1\right)\\
\end{array}
\end{array}
if a < 2.9000000000000002e149Initial program 99.9%
add-cube-cbrt99.6%
pow399.7%
Applied egg-rr99.7%
Taylor expanded in z around inf 64.9%
if 2.9000000000000002e149 < a Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in a around -inf 99.9%
Taylor expanded in i around inf 83.3%
associate-/l*83.3%
Simplified83.3%
Final simplification66.8%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 34000000.0) z (if (<= a 2.25e+145) (* y i) (* a (+ 1.0 (/ z a))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 34000000.0) {
tmp = z;
} else if (a <= 2.25e+145) {
tmp = y * i;
} else {
tmp = a * (1.0 + (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 (a <= 34000000.0d0) then
tmp = z
else if (a <= 2.25d+145) then
tmp = y * i
else
tmp = a * (1.0d0 + (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 (a <= 34000000.0) {
tmp = z;
} else if (a <= 2.25e+145) {
tmp = y * i;
} else {
tmp = a * (1.0 + (z / a));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 34000000.0: tmp = z elif a <= 2.25e+145: tmp = y * i else: tmp = a * (1.0 + (z / a)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 34000000.0) tmp = z; elseif (a <= 2.25e+145) tmp = Float64(y * i); else tmp = Float64(a * Float64(1.0 + Float64(z / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (a <= 34000000.0) tmp = z; elseif (a <= 2.25e+145) tmp = y * i; else tmp = a * (1.0 + (z / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 34000000.0], z, If[LessEqual[a, 2.25e+145], N[(y * i), $MachinePrecision], N[(a * N[(1.0 + N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 34000000:\\
\;\;\;\;z\\
\mathbf{elif}\;a \leq 2.25 \cdot 10^{+145}:\\
\;\;\;\;y \cdot i\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(1 + \frac{z}{a}\right)\\
\end{array}
\end{array}
if a < 3.4e7Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in z around inf 21.7%
if 3.4e7 < a < 2.2499999999999999e145Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in y around inf 39.0%
*-commutative39.0%
Simplified39.0%
if 2.2499999999999999e145 < a Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in a around -inf 99.9%
Taylor expanded in z around inf 64.3%
Final simplification28.3%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 25500000.0) z (if (<= a 9.4e+148) (* y i) (* a (+ 1.0 (/ t a))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 25500000.0) {
tmp = z;
} else if (a <= 9.4e+148) {
tmp = y * i;
} else {
tmp = a * (1.0 + (t / 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 <= 25500000.0d0) then
tmp = z
else if (a <= 9.4d+148) then
tmp = y * i
else
tmp = a * (1.0d0 + (t / 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 <= 25500000.0) {
tmp = z;
} else if (a <= 9.4e+148) {
tmp = y * i;
} else {
tmp = a * (1.0 + (t / a));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 25500000.0: tmp = z elif a <= 9.4e+148: tmp = y * i else: tmp = a * (1.0 + (t / a)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 25500000.0) tmp = z; elseif (a <= 9.4e+148) tmp = Float64(y * i); else tmp = Float64(a * Float64(1.0 + Float64(t / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (a <= 25500000.0) tmp = z; elseif (a <= 9.4e+148) tmp = y * i; else tmp = a * (1.0 + (t / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 25500000.0], z, If[LessEqual[a, 9.4e+148], N[(y * i), $MachinePrecision], N[(a * N[(1.0 + N[(t / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 25500000:\\
\;\;\;\;z\\
\mathbf{elif}\;a \leq 9.4 \cdot 10^{+148}:\\
\;\;\;\;y \cdot i\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(1 + \frac{t}{a}\right)\\
\end{array}
\end{array}
if a < 2.55e7Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in z around inf 21.7%
if 2.55e7 < a < 9.3999999999999994e148Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in y around inf 37.8%
*-commutative37.8%
Simplified37.8%
if 9.3999999999999994e148 < a Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in a around -inf 99.9%
Taylor expanded in t around inf 61.2%
Final simplification27.7%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 27500000.0) z (* a (+ (* i (/ y a)) 1.0))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 27500000.0) {
tmp = z;
} else {
tmp = a * ((i * (y / a)) + 1.0);
}
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 <= 27500000.0d0) then
tmp = z
else
tmp = a * ((i * (y / a)) + 1.0d0)
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 <= 27500000.0) {
tmp = z;
} else {
tmp = a * ((i * (y / a)) + 1.0);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 27500000.0: tmp = z else: tmp = a * ((i * (y / a)) + 1.0) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 27500000.0) tmp = z; else tmp = Float64(a * Float64(Float64(i * Float64(y / a)) + 1.0)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (a <= 27500000.0) tmp = z; else tmp = a * ((i * (y / a)) + 1.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 27500000.0], z, N[(a * N[(N[(i * N[(y / a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 27500000:\\
\;\;\;\;z\\
\mathbf{else}:\\
\;\;\;\;a \cdot \left(i \cdot \frac{y}{a} + 1\right)\\
\end{array}
\end{array}
if a < 2.75e7Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in z around inf 21.7%
if 2.75e7 < a Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in a around -inf 99.9%
Taylor expanded in i around inf 64.0%
associate-/l*64.0%
Simplified64.0%
Final simplification31.2%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 22000000.0) z (if (<= a 2.4e+148) (* 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 <= 22000000.0) {
tmp = z;
} else if (a <= 2.4e+148) {
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 <= 22000000.0d0) then
tmp = z
else if (a <= 2.4d+148) 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 <= 22000000.0) {
tmp = z;
} else if (a <= 2.4e+148) {
tmp = y * i;
} else {
tmp = a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 22000000.0: tmp = z elif a <= 2.4e+148: tmp = y * i else: tmp = a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 22000000.0) tmp = z; elseif (a <= 2.4e+148) 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 <= 22000000.0) tmp = z; elseif (a <= 2.4e+148) tmp = y * i; else tmp = a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 22000000.0], z, If[LessEqual[a, 2.4e+148], N[(y * i), $MachinePrecision], a]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 22000000:\\
\;\;\;\;z\\
\mathbf{elif}\;a \leq 2.4 \cdot 10^{+148}:\\
\;\;\;\;y \cdot i\\
\mathbf{else}:\\
\;\;\;\;a\\
\end{array}
\end{array}
if a < 2.2e7Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in z around inf 21.7%
if 2.2e7 < a < 2.39999999999999995e148Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in y around inf 37.8%
*-commutative37.8%
Simplified37.8%
if 2.39999999999999995e148 < a Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in a around inf 61.2%
(FPCore (x y z t a b c i) :precision binary64 (if (<= a 3.3e+102) z a))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (a <= 3.3e+102) {
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 <= 3.3d+102) 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 <= 3.3e+102) {
tmp = z;
} else {
tmp = a;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if a <= 3.3e+102: tmp = z else: tmp = a return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (a <= 3.3e+102) 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 <= 3.3e+102) tmp = z; else tmp = a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[a, 3.3e+102], z, a]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 3.3 \cdot 10^{+102}:\\
\;\;\;\;z\\
\mathbf{else}:\\
\;\;\;\;a\\
\end{array}
\end{array}
if a < 3.29999999999999999e102Initial program 99.9%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in z around inf 20.8%
if 3.29999999999999999e102 < a Initial program 100.0%
associate-+l+100.0%
associate-+l+100.0%
+-commutative100.0%
associate-+l+100.0%
+-commutative100.0%
associate-+l+100.0%
+-commutative100.0%
associate-+l+100.0%
+-commutative100.0%
fma-define100.0%
+-commutative100.0%
fma-define100.0%
Simplified100.0%
Taylor expanded in a around inf 50.6%
(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%
associate-+l+99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
+-commutative99.9%
fma-define99.9%
Simplified99.9%
Taylor expanded in a around inf 16.5%
herbie shell --seed 2024186
(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)))