
(FPCore (x y z t a) :precision binary64 (+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))
double code(double x, double y, double z, double t, double a) {
return ((log((x + y)) + log(z)) - t) + ((a - 0.5) * log(t));
}
real(8) function code(x, y, z, t, a)
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
code = ((log((x + y)) + log(z)) - t) + ((a - 0.5d0) * log(t))
end function
public static double code(double x, double y, double z, double t, double a) {
return ((Math.log((x + y)) + Math.log(z)) - t) + ((a - 0.5) * Math.log(t));
}
def code(x, y, z, t, a): return ((math.log((x + y)) + math.log(z)) - t) + ((a - 0.5) * math.log(t))
function code(x, y, z, t, a) return Float64(Float64(Float64(log(Float64(x + y)) + log(z)) - t) + Float64(Float64(a - 0.5) * log(t))) end
function tmp = code(x, y, z, t, a) tmp = ((log((x + y)) + log(z)) - t) + ((a - 0.5) * log(t)); end
code[x_, y_, z_, t_, a_] := N[(N[(N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[Log[z], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision] + N[(N[(a - 0.5), $MachinePrecision] * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))
double code(double x, double y, double z, double t, double a) {
return ((log((x + y)) + log(z)) - t) + ((a - 0.5) * log(t));
}
real(8) function code(x, y, z, t, a)
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
code = ((log((x + y)) + log(z)) - t) + ((a - 0.5d0) * log(t))
end function
public static double code(double x, double y, double z, double t, double a) {
return ((Math.log((x + y)) + Math.log(z)) - t) + ((a - 0.5) * Math.log(t));
}
def code(x, y, z, t, a): return ((math.log((x + y)) + math.log(z)) - t) + ((a - 0.5) * math.log(t))
function code(x, y, z, t, a) return Float64(Float64(Float64(log(Float64(x + y)) + log(z)) - t) + Float64(Float64(a - 0.5) * log(t))) end
function tmp = code(x, y, z, t, a) tmp = ((log((x + y)) + log(z)) - t) + ((a - 0.5) * log(t)); end
code[x_, y_, z_, t_, a_] := N[(N[(N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[Log[z], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision] + N[(N[(a - 0.5), $MachinePrecision] * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t
\end{array}
(FPCore (x y z t a) :precision binary64 (fma (+ a -0.5) (log t) (+ (log (+ x y)) (- (log z) t))))
double code(double x, double y, double z, double t, double a) {
return fma((a + -0.5), log(t), (log((x + y)) + (log(z) - t)));
}
function code(x, y, z, t, a) return fma(Float64(a + -0.5), log(t), Float64(log(Float64(x + y)) + Float64(log(z) - t))) end
code[x_, y_, z_, t_, a_] := N[(N[(a + -0.5), $MachinePrecision] * N[Log[t], $MachinePrecision] + N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(a + -0.5, \log t, \log \left(x + y\right) + \left(\log z - t\right)\right)
\end{array}
Initial program 99.7%
+-commutative99.7%
fma-def99.7%
sub-neg99.7%
metadata-eval99.7%
associate--l+99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x y z t a) :precision binary64 (if (<= t 210.0) (+ (log z) (+ (log (+ x y)) (* (log t) (- a 0.5)))) (- (* a (log t)) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 210.0) {
tmp = log(z) + (log((x + y)) + (log(t) * (a - 0.5)));
} else {
tmp = (a * log(t)) - t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if (t <= 210.0d0) then
tmp = log(z) + (log((x + y)) + (log(t) * (a - 0.5d0)))
else
tmp = (a * log(t)) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 210.0) {
tmp = Math.log(z) + (Math.log((x + y)) + (Math.log(t) * (a - 0.5)));
} else {
tmp = (a * Math.log(t)) - t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 210.0: tmp = math.log(z) + (math.log((x + y)) + (math.log(t) * (a - 0.5))) else: tmp = (a * math.log(t)) - t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 210.0) tmp = Float64(log(z) + Float64(log(Float64(x + y)) + Float64(log(t) * Float64(a - 0.5)))); else tmp = Float64(Float64(a * log(t)) - t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 210.0) tmp = log(z) + (log((x + y)) + (log(t) * (a - 0.5))); else tmp = (a * log(t)) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 210.0], N[(N[Log[z], $MachinePrecision] + N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[t], $MachinePrecision] * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 210:\\
\;\;\;\;\log z + \left(\log \left(x + y\right) + \log t \cdot \left(a - 0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;a \cdot \log t - t\\
\end{array}
\end{array}
if t < 210Initial program 99.5%
associate--l+99.5%
sub-neg99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in t around 0 98.3%
if 210 < t Initial program 99.9%
associate-+l-99.9%
associate--l+99.9%
sub-neg99.9%
+-commutative99.9%
*-commutative99.9%
distribute-rgt-neg-in99.9%
fma-udef99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
associate-+r-99.9%
fma-udef99.9%
associate--r+99.9%
+-commutative99.9%
sum-log81.4%
Applied egg-rr81.4%
Taylor expanded in x around inf 71.8%
mul-1-neg71.8%
log-rec71.8%
remove-double-neg71.8%
log-prod55.3%
*-commutative55.3%
Simplified55.3%
Taylor expanded in a around inf 98.8%
*-commutative98.8%
Simplified98.8%
Final simplification98.6%
(FPCore (x y z t a) :precision binary64 (if (<= t 210.0) (+ (+ (log z) (log y)) (* (log t) (- a 0.5))) (- (* a (log t)) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 210.0) {
tmp = (log(z) + log(y)) + (log(t) * (a - 0.5));
} else {
tmp = (a * log(t)) - t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if (t <= 210.0d0) then
tmp = (log(z) + log(y)) + (log(t) * (a - 0.5d0))
else
tmp = (a * log(t)) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 210.0) {
tmp = (Math.log(z) + Math.log(y)) + (Math.log(t) * (a - 0.5));
} else {
tmp = (a * Math.log(t)) - t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 210.0: tmp = (math.log(z) + math.log(y)) + (math.log(t) * (a - 0.5)) else: tmp = (a * math.log(t)) - t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 210.0) tmp = Float64(Float64(log(z) + log(y)) + Float64(log(t) * Float64(a - 0.5))); else tmp = Float64(Float64(a * log(t)) - t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 210.0) tmp = (log(z) + log(y)) + (log(t) * (a - 0.5)); else tmp = (a * log(t)) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 210.0], N[(N[(N[Log[z], $MachinePrecision] + N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(N[Log[t], $MachinePrecision] * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 210:\\
\;\;\;\;\left(\log z + \log y\right) + \log t \cdot \left(a - 0.5\right)\\
\mathbf{else}:\\
\;\;\;\;a \cdot \log t - t\\
\end{array}
\end{array}
if t < 210Initial program 99.5%
associate-+l-99.5%
associate--l+99.5%
sub-neg99.5%
+-commutative99.5%
*-commutative99.5%
distribute-rgt-neg-in99.5%
fma-udef99.5%
sub-neg99.5%
+-commutative99.5%
distribute-neg-in99.5%
metadata-eval99.5%
metadata-eval99.5%
unsub-neg99.5%
Simplified99.5%
Taylor expanded in x around 0 70.6%
Taylor expanded in t around 0 70.1%
if 210 < t Initial program 99.9%
associate-+l-99.9%
associate--l+99.9%
sub-neg99.9%
+-commutative99.9%
*-commutative99.9%
distribute-rgt-neg-in99.9%
fma-udef99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
associate-+r-99.9%
fma-udef99.9%
associate--r+99.9%
+-commutative99.9%
sum-log81.4%
Applied egg-rr81.4%
Taylor expanded in x around inf 71.8%
mul-1-neg71.8%
log-rec71.8%
remove-double-neg71.8%
log-prod55.3%
*-commutative55.3%
Simplified55.3%
Taylor expanded in a around inf 98.8%
*-commutative98.8%
Simplified98.8%
Final simplification83.7%
(FPCore (x y z t a) :precision binary64 (+ (+ (log (+ x y)) (- (log z) t)) (* (+ a -0.5) (log t))))
double code(double x, double y, double z, double t, double a) {
return (log((x + y)) + (log(z) - t)) + ((a + -0.5) * log(t));
}
real(8) function code(x, y, z, t, a)
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
code = (log((x + y)) + (log(z) - t)) + ((a + (-0.5d0)) * log(t))
end function
public static double code(double x, double y, double z, double t, double a) {
return (Math.log((x + y)) + (Math.log(z) - t)) + ((a + -0.5) * Math.log(t));
}
def code(x, y, z, t, a): return (math.log((x + y)) + (math.log(z) - t)) + ((a + -0.5) * math.log(t))
function code(x, y, z, t, a) return Float64(Float64(log(Float64(x + y)) + Float64(log(z) - t)) + Float64(Float64(a + -0.5) * log(t))) end
function tmp = code(x, y, z, t, a) tmp = (log((x + y)) + (log(z) - t)) + ((a + -0.5) * log(t)); end
code[x_, y_, z_, t_, a_] := N[(N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision] + N[(N[(a + -0.5), $MachinePrecision] * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\log \left(x + y\right) + \left(\log z - t\right)\right) + \left(a + -0.5\right) \cdot \log t
\end{array}
Initial program 99.7%
associate--l+99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x y z t a) :precision binary64 (+ (+ (log z) (log y)) (- (* (log t) (- a 0.5)) t)))
double code(double x, double y, double z, double t, double a) {
return (log(z) + log(y)) + ((log(t) * (a - 0.5)) - t);
}
real(8) function code(x, y, z, t, a)
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
code = (log(z) + log(y)) + ((log(t) * (a - 0.5d0)) - t)
end function
public static double code(double x, double y, double z, double t, double a) {
return (Math.log(z) + Math.log(y)) + ((Math.log(t) * (a - 0.5)) - t);
}
def code(x, y, z, t, a): return (math.log(z) + math.log(y)) + ((math.log(t) * (a - 0.5)) - t)
function code(x, y, z, t, a) return Float64(Float64(log(z) + log(y)) + Float64(Float64(log(t) * Float64(a - 0.5)) - t)) end
function tmp = code(x, y, z, t, a) tmp = (log(z) + log(y)) + ((log(t) * (a - 0.5)) - t); end
code[x_, y_, z_, t_, a_] := N[(N[(N[Log[z], $MachinePrecision] + N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(N[(N[Log[t], $MachinePrecision] * N[(a - 0.5), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\log z + \log y\right) + \left(\log t \cdot \left(a - 0.5\right) - t\right)
\end{array}
Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 72.4%
Final simplification72.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (* a (log t)) t)))
(if (<= a -1.4)
t_1
(if (<= a 7.2e-72)
(+ (log (+ x y)) (- (log z) t))
(if (<= a 6.2e-33) (log (/ z (/ (pow t (- 0.5 a)) y))) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (a * log(t)) - t;
double tmp;
if (a <= -1.4) {
tmp = t_1;
} else if (a <= 7.2e-72) {
tmp = log((x + y)) + (log(z) - t);
} else if (a <= 6.2e-33) {
tmp = log((z / (pow(t, (0.5 - a)) / y)));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: t_1
real(8) :: tmp
t_1 = (a * log(t)) - t
if (a <= (-1.4d0)) then
tmp = t_1
else if (a <= 7.2d-72) then
tmp = log((x + y)) + (log(z) - t)
else if (a <= 6.2d-33) then
tmp = log((z / ((t ** (0.5d0 - a)) / y)))
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 t_1 = (a * Math.log(t)) - t;
double tmp;
if (a <= -1.4) {
tmp = t_1;
} else if (a <= 7.2e-72) {
tmp = Math.log((x + y)) + (Math.log(z) - t);
} else if (a <= 6.2e-33) {
tmp = Math.log((z / (Math.pow(t, (0.5 - a)) / y)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (a * math.log(t)) - t tmp = 0 if a <= -1.4: tmp = t_1 elif a <= 7.2e-72: tmp = math.log((x + y)) + (math.log(z) - t) elif a <= 6.2e-33: tmp = math.log((z / (math.pow(t, (0.5 - a)) / y))) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(a * log(t)) - t) tmp = 0.0 if (a <= -1.4) tmp = t_1; elseif (a <= 7.2e-72) tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); elseif (a <= 6.2e-33) tmp = log(Float64(z / Float64((t ^ Float64(0.5 - a)) / y))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (a * log(t)) - t; tmp = 0.0; if (a <= -1.4) tmp = t_1; elseif (a <= 7.2e-72) tmp = log((x + y)) + (log(z) - t); elseif (a <= 6.2e-33) tmp = log((z / ((t ^ (0.5 - a)) / y))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]}, If[LessEqual[a, -1.4], t$95$1, If[LessEqual[a, 7.2e-72], N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 6.2e-33], N[Log[N[(z / N[(N[Power[t, N[(0.5 - a), $MachinePrecision]], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \log t - t\\
\mathbf{if}\;a \leq -1.4:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 7.2 \cdot 10^{-72}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\mathbf{elif}\;a \leq 6.2 \cdot 10^{-33}:\\
\;\;\;\;\log \left(\frac{z}{\frac{{t}^{\left(0.5 - a\right)}}{y}}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.3999999999999999 or 6.19999999999999994e-33 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.6%
sub-neg99.6%
+-commutative99.6%
*-commutative99.6%
distribute-rgt-neg-in99.6%
fma-udef99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
metadata-eval99.6%
metadata-eval99.6%
unsub-neg99.6%
Simplified99.6%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.9%
Applied egg-rr80.9%
Taylor expanded in x around inf 72.0%
mul-1-neg72.0%
log-rec72.0%
remove-double-neg72.0%
log-prod56.2%
*-commutative56.2%
Simplified56.2%
Taylor expanded in a around inf 98.2%
*-commutative98.2%
Simplified98.2%
if -1.3999999999999999 < a < 7.2e-72Initial program 99.8%
associate-+l-99.8%
associate--l+99.8%
sub-neg99.8%
+-commutative99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
fma-udef99.8%
sub-neg99.8%
+-commutative99.8%
distribute-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
unsub-neg99.8%
Simplified99.8%
Taylor expanded in t around inf 60.0%
if 7.2e-72 < a < 6.19999999999999994e-33Initial program 98.9%
associate-+l-98.9%
associate--l+99.2%
sub-neg99.2%
+-commutative99.2%
*-commutative99.2%
distribute-rgt-neg-in99.2%
fma-udef99.2%
sub-neg99.2%
+-commutative99.2%
distribute-neg-in99.2%
metadata-eval99.2%
metadata-eval99.2%
unsub-neg99.2%
Simplified99.2%
Taylor expanded in x around 0 90.1%
Taylor expanded in t around 0 80.2%
Taylor expanded in y around inf 80.2%
*-commutative80.2%
mul-1-neg80.2%
log-rec80.2%
remove-double-neg80.2%
log-prod61.1%
log-pow61.1%
log-div61.1%
associate-/l*61.1%
Simplified61.1%
Final simplification82.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (* a (log t)) t)))
(if (<= a -1.85)
t_1
(if (<= a 6.5e-72)
(+ (log (+ x y)) (- (log z) t))
(if (<= a 5.2e-33) (log (/ (* y z) (pow t (- 0.5 a)))) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (a * log(t)) - t;
double tmp;
if (a <= -1.85) {
tmp = t_1;
} else if (a <= 6.5e-72) {
tmp = log((x + y)) + (log(z) - t);
} else if (a <= 5.2e-33) {
tmp = log(((y * z) / pow(t, (0.5 - a))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: t_1
real(8) :: tmp
t_1 = (a * log(t)) - t
if (a <= (-1.85d0)) then
tmp = t_1
else if (a <= 6.5d-72) then
tmp = log((x + y)) + (log(z) - t)
else if (a <= 5.2d-33) then
tmp = log(((y * z) / (t ** (0.5d0 - a))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (a * Math.log(t)) - t;
double tmp;
if (a <= -1.85) {
tmp = t_1;
} else if (a <= 6.5e-72) {
tmp = Math.log((x + y)) + (Math.log(z) - t);
} else if (a <= 5.2e-33) {
tmp = Math.log(((y * z) / Math.pow(t, (0.5 - a))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (a * math.log(t)) - t tmp = 0 if a <= -1.85: tmp = t_1 elif a <= 6.5e-72: tmp = math.log((x + y)) + (math.log(z) - t) elif a <= 5.2e-33: tmp = math.log(((y * z) / math.pow(t, (0.5 - a)))) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(a * log(t)) - t) tmp = 0.0 if (a <= -1.85) tmp = t_1; elseif (a <= 6.5e-72) tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); elseif (a <= 5.2e-33) tmp = log(Float64(Float64(y * z) / (t ^ Float64(0.5 - a)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (a * log(t)) - t; tmp = 0.0; if (a <= -1.85) tmp = t_1; elseif (a <= 6.5e-72) tmp = log((x + y)) + (log(z) - t); elseif (a <= 5.2e-33) tmp = log(((y * z) / (t ^ (0.5 - a)))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]}, If[LessEqual[a, -1.85], t$95$1, If[LessEqual[a, 6.5e-72], N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 5.2e-33], N[Log[N[(N[(y * z), $MachinePrecision] / N[Power[t, N[(0.5 - a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \log t - t\\
\mathbf{if}\;a \leq -1.85:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 6.5 \cdot 10^{-72}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\mathbf{elif}\;a \leq 5.2 \cdot 10^{-33}:\\
\;\;\;\;\log \left(\frac{y \cdot z}{{t}^{\left(0.5 - a\right)}}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.8500000000000001 or 5.19999999999999988e-33 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.6%
sub-neg99.6%
+-commutative99.6%
*-commutative99.6%
distribute-rgt-neg-in99.6%
fma-udef99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
metadata-eval99.6%
metadata-eval99.6%
unsub-neg99.6%
Simplified99.6%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.9%
Applied egg-rr80.9%
Taylor expanded in x around inf 72.0%
mul-1-neg72.0%
log-rec72.0%
remove-double-neg72.0%
log-prod56.2%
*-commutative56.2%
Simplified56.2%
Taylor expanded in a around inf 98.2%
*-commutative98.2%
Simplified98.2%
if -1.8500000000000001 < a < 6.4999999999999997e-72Initial program 99.8%
associate-+l-99.8%
associate--l+99.8%
sub-neg99.8%
+-commutative99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
fma-udef99.8%
sub-neg99.8%
+-commutative99.8%
distribute-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
unsub-neg99.8%
Simplified99.8%
Taylor expanded in t around inf 60.0%
if 6.4999999999999997e-72 < a < 5.19999999999999988e-33Initial program 98.9%
associate-+l-98.9%
associate--l+99.2%
sub-neg99.2%
+-commutative99.2%
*-commutative99.2%
distribute-rgt-neg-in99.2%
fma-udef99.2%
sub-neg99.2%
+-commutative99.2%
distribute-neg-in99.2%
metadata-eval99.2%
metadata-eval99.2%
unsub-neg99.2%
Simplified99.2%
Taylor expanded in x around 0 90.1%
Taylor expanded in t around 0 80.2%
+-commutative80.2%
log-prod61.1%
add-log-exp61.1%
diff-log61.1%
exp-to-pow61.1%
Applied egg-rr61.1%
Final simplification82.3%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -880000000000.0) (not (<= a 2.1e-32))) (- (* a (log t)) t) (- (+ (log (* (+ x y) z)) (* (log t) (- a 0.5))) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -880000000000.0) || !(a <= 2.1e-32)) {
tmp = (a * log(t)) - t;
} else {
tmp = (log(((x + y) * z)) + (log(t) * (a - 0.5))) - t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if ((a <= (-880000000000.0d0)) .or. (.not. (a <= 2.1d-32))) then
tmp = (a * log(t)) - t
else
tmp = (log(((x + y) * z)) + (log(t) * (a - 0.5d0))) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -880000000000.0) || !(a <= 2.1e-32)) {
tmp = (a * Math.log(t)) - t;
} else {
tmp = (Math.log(((x + y) * z)) + (Math.log(t) * (a - 0.5))) - t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -880000000000.0) or not (a <= 2.1e-32): tmp = (a * math.log(t)) - t else: tmp = (math.log(((x + y) * z)) + (math.log(t) * (a - 0.5))) - t return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -880000000000.0) || !(a <= 2.1e-32)) tmp = Float64(Float64(a * log(t)) - t); else tmp = Float64(Float64(log(Float64(Float64(x + y) * z)) + Float64(log(t) * Float64(a - 0.5))) - t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -880000000000.0) || ~((a <= 2.1e-32))) tmp = (a * log(t)) - t; else tmp = (log(((x + y) * z)) + (log(t) * (a - 0.5))) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -880000000000.0], N[Not[LessEqual[a, 2.1e-32]], $MachinePrecision]], N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[(N[Log[N[(N[(x + y), $MachinePrecision] * z), $MachinePrecision]], $MachinePrecision] + N[(N[Log[t], $MachinePrecision] * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -880000000000 \lor \neg \left(a \leq 2.1 \cdot 10^{-32}\right):\\
\;\;\;\;a \cdot \log t - t\\
\mathbf{else}:\\
\;\;\;\;\left(\log \left(\left(x + y\right) \cdot z\right) + \log t \cdot \left(a - 0.5\right)\right) - t\\
\end{array}
\end{array}
if a < -8.8e11 or 2.0999999999999999e-32 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.7%
Applied egg-rr80.7%
Taylor expanded in x around inf 72.3%
mul-1-neg72.3%
log-rec72.3%
remove-double-neg72.3%
log-prod56.2%
*-commutative56.2%
Simplified56.2%
Taylor expanded in a around inf 98.9%
*-commutative98.9%
Simplified98.9%
if -8.8e11 < a < 2.0999999999999999e-32Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.1%
Applied egg-rr80.1%
Final simplification90.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (* a (log t)) t)))
(if (<= a -1.4)
t_1
(if (<= a 1.45e-78)
(- (+ (log z) (log y)) t)
(if (<= a 6.2e-33) (+ (log (* y z)) (* -0.5 (log t))) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (a * log(t)) - t;
double tmp;
if (a <= -1.4) {
tmp = t_1;
} else if (a <= 1.45e-78) {
tmp = (log(z) + log(y)) - t;
} else if (a <= 6.2e-33) {
tmp = log((y * z)) + (-0.5 * log(t));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: t_1
real(8) :: tmp
t_1 = (a * log(t)) - t
if (a <= (-1.4d0)) then
tmp = t_1
else if (a <= 1.45d-78) then
tmp = (log(z) + log(y)) - t
else if (a <= 6.2d-33) then
tmp = log((y * z)) + ((-0.5d0) * log(t))
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 t_1 = (a * Math.log(t)) - t;
double tmp;
if (a <= -1.4) {
tmp = t_1;
} else if (a <= 1.45e-78) {
tmp = (Math.log(z) + Math.log(y)) - t;
} else if (a <= 6.2e-33) {
tmp = Math.log((y * z)) + (-0.5 * Math.log(t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (a * math.log(t)) - t tmp = 0 if a <= -1.4: tmp = t_1 elif a <= 1.45e-78: tmp = (math.log(z) + math.log(y)) - t elif a <= 6.2e-33: tmp = math.log((y * z)) + (-0.5 * math.log(t)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(a * log(t)) - t) tmp = 0.0 if (a <= -1.4) tmp = t_1; elseif (a <= 1.45e-78) tmp = Float64(Float64(log(z) + log(y)) - t); elseif (a <= 6.2e-33) tmp = Float64(log(Float64(y * z)) + Float64(-0.5 * log(t))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (a * log(t)) - t; tmp = 0.0; if (a <= -1.4) tmp = t_1; elseif (a <= 1.45e-78) tmp = (log(z) + log(y)) - t; elseif (a <= 6.2e-33) tmp = log((y * z)) + (-0.5 * log(t)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]}, If[LessEqual[a, -1.4], t$95$1, If[LessEqual[a, 1.45e-78], N[(N[(N[Log[z], $MachinePrecision] + N[Log[y], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], If[LessEqual[a, 6.2e-33], N[(N[Log[N[(y * z), $MachinePrecision]], $MachinePrecision] + N[(-0.5 * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \log t - t\\
\mathbf{if}\;a \leq -1.4:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 1.45 \cdot 10^{-78}:\\
\;\;\;\;\left(\log z + \log y\right) - t\\
\mathbf{elif}\;a \leq 6.2 \cdot 10^{-33}:\\
\;\;\;\;\log \left(y \cdot z\right) + -0.5 \cdot \log t\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.3999999999999999 or 6.19999999999999994e-33 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.6%
sub-neg99.6%
+-commutative99.6%
*-commutative99.6%
distribute-rgt-neg-in99.6%
fma-udef99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
metadata-eval99.6%
metadata-eval99.6%
unsub-neg99.6%
Simplified99.6%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.9%
Applied egg-rr80.9%
Taylor expanded in x around inf 72.0%
mul-1-neg72.0%
log-rec72.0%
remove-double-neg72.0%
log-prod56.2%
*-commutative56.2%
Simplified56.2%
Taylor expanded in a around inf 98.2%
*-commutative98.2%
Simplified98.2%
if -1.3999999999999999 < a < 1.45e-78Initial program 99.8%
associate-+l-99.8%
associate--l+99.8%
sub-neg99.8%
+-commutative99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
fma-udef99.8%
sub-neg99.8%
+-commutative99.8%
distribute-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
unsub-neg99.8%
Simplified99.8%
Taylor expanded in x around 0 64.4%
Taylor expanded in t around inf 43.9%
if 1.45e-78 < a < 6.19999999999999994e-33Initial program 98.9%
associate-+l-98.9%
associate--l+99.2%
sub-neg99.2%
+-commutative99.2%
*-commutative99.2%
distribute-rgt-neg-in99.2%
fma-udef99.2%
sub-neg99.2%
+-commutative99.2%
distribute-neg-in99.2%
metadata-eval99.2%
metadata-eval99.2%
unsub-neg99.2%
Simplified99.2%
Taylor expanded in x around 0 90.1%
Taylor expanded in t around 0 80.2%
Taylor expanded in a around 0 80.2%
cancel-sign-sub-inv80.2%
log-prod61.1%
metadata-eval61.1%
*-commutative61.1%
Simplified61.1%
Final simplification76.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (* a (log t)) t)))
(if (<= a -1.4)
t_1
(if (<= a 1.05e-73)
(+ (log (+ x y)) (- (log z) t))
(if (<= a 5.3e-33) (+ (log (* y z)) (* -0.5 (log t))) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (a * log(t)) - t;
double tmp;
if (a <= -1.4) {
tmp = t_1;
} else if (a <= 1.05e-73) {
tmp = log((x + y)) + (log(z) - t);
} else if (a <= 5.3e-33) {
tmp = log((y * z)) + (-0.5 * log(t));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: t_1
real(8) :: tmp
t_1 = (a * log(t)) - t
if (a <= (-1.4d0)) then
tmp = t_1
else if (a <= 1.05d-73) then
tmp = log((x + y)) + (log(z) - t)
else if (a <= 5.3d-33) then
tmp = log((y * z)) + ((-0.5d0) * log(t))
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 t_1 = (a * Math.log(t)) - t;
double tmp;
if (a <= -1.4) {
tmp = t_1;
} else if (a <= 1.05e-73) {
tmp = Math.log((x + y)) + (Math.log(z) - t);
} else if (a <= 5.3e-33) {
tmp = Math.log((y * z)) + (-0.5 * Math.log(t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (a * math.log(t)) - t tmp = 0 if a <= -1.4: tmp = t_1 elif a <= 1.05e-73: tmp = math.log((x + y)) + (math.log(z) - t) elif a <= 5.3e-33: tmp = math.log((y * z)) + (-0.5 * math.log(t)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(a * log(t)) - t) tmp = 0.0 if (a <= -1.4) tmp = t_1; elseif (a <= 1.05e-73) tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); elseif (a <= 5.3e-33) tmp = Float64(log(Float64(y * z)) + Float64(-0.5 * log(t))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (a * log(t)) - t; tmp = 0.0; if (a <= -1.4) tmp = t_1; elseif (a <= 1.05e-73) tmp = log((x + y)) + (log(z) - t); elseif (a <= 5.3e-33) tmp = log((y * z)) + (-0.5 * log(t)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]}, If[LessEqual[a, -1.4], t$95$1, If[LessEqual[a, 1.05e-73], N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 5.3e-33], N[(N[Log[N[(y * z), $MachinePrecision]], $MachinePrecision] + N[(-0.5 * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \log t - t\\
\mathbf{if}\;a \leq -1.4:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 1.05 \cdot 10^{-73}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\mathbf{elif}\;a \leq 5.3 \cdot 10^{-33}:\\
\;\;\;\;\log \left(y \cdot z\right) + -0.5 \cdot \log t\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.3999999999999999 or 5.29999999999999968e-33 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.6%
sub-neg99.6%
+-commutative99.6%
*-commutative99.6%
distribute-rgt-neg-in99.6%
fma-udef99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
metadata-eval99.6%
metadata-eval99.6%
unsub-neg99.6%
Simplified99.6%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.9%
Applied egg-rr80.9%
Taylor expanded in x around inf 72.0%
mul-1-neg72.0%
log-rec72.0%
remove-double-neg72.0%
log-prod56.2%
*-commutative56.2%
Simplified56.2%
Taylor expanded in a around inf 98.2%
*-commutative98.2%
Simplified98.2%
if -1.3999999999999999 < a < 1.0499999999999999e-73Initial program 99.8%
associate-+l-99.8%
associate--l+99.8%
sub-neg99.8%
+-commutative99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
fma-udef99.8%
sub-neg99.8%
+-commutative99.8%
distribute-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
unsub-neg99.8%
Simplified99.8%
Taylor expanded in t around inf 60.0%
if 1.0499999999999999e-73 < a < 5.29999999999999968e-33Initial program 98.9%
associate-+l-98.9%
associate--l+99.2%
sub-neg99.2%
+-commutative99.2%
*-commutative99.2%
distribute-rgt-neg-in99.2%
fma-udef99.2%
sub-neg99.2%
+-commutative99.2%
distribute-neg-in99.2%
metadata-eval99.2%
metadata-eval99.2%
unsub-neg99.2%
Simplified99.2%
Taylor expanded in x around 0 90.1%
Taylor expanded in t around 0 80.2%
Taylor expanded in a around 0 80.2%
cancel-sign-sub-inv80.2%
log-prod61.1%
metadata-eval61.1%
*-commutative61.1%
Simplified61.1%
Final simplification82.3%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -220000000000.0) (not (<= a 2.1e-32))) (- (* a (log t)) t) (+ (log (* y z)) (- (* (log t) (- a 0.5)) t))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -220000000000.0) || !(a <= 2.1e-32)) {
tmp = (a * log(t)) - t;
} else {
tmp = log((y * z)) + ((log(t) * (a - 0.5)) - t);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if ((a <= (-220000000000.0d0)) .or. (.not. (a <= 2.1d-32))) then
tmp = (a * log(t)) - t
else
tmp = log((y * z)) + ((log(t) * (a - 0.5d0)) - t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -220000000000.0) || !(a <= 2.1e-32)) {
tmp = (a * Math.log(t)) - t;
} else {
tmp = Math.log((y * z)) + ((Math.log(t) * (a - 0.5)) - t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -220000000000.0) or not (a <= 2.1e-32): tmp = (a * math.log(t)) - t else: tmp = math.log((y * z)) + ((math.log(t) * (a - 0.5)) - t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -220000000000.0) || !(a <= 2.1e-32)) tmp = Float64(Float64(a * log(t)) - t); else tmp = Float64(log(Float64(y * z)) + Float64(Float64(log(t) * Float64(a - 0.5)) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -220000000000.0) || ~((a <= 2.1e-32))) tmp = (a * log(t)) - t; else tmp = log((y * z)) + ((log(t) * (a - 0.5)) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -220000000000.0], N[Not[LessEqual[a, 2.1e-32]], $MachinePrecision]], N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[Log[N[(y * z), $MachinePrecision]], $MachinePrecision] + N[(N[(N[Log[t], $MachinePrecision] * N[(a - 0.5), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -220000000000 \lor \neg \left(a \leq 2.1 \cdot 10^{-32}\right):\\
\;\;\;\;a \cdot \log t - t\\
\mathbf{else}:\\
\;\;\;\;\log \left(y \cdot z\right) + \left(\log t \cdot \left(a - 0.5\right) - t\right)\\
\end{array}
\end{array}
if a < -2.2e11 or 2.0999999999999999e-32 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.7%
Applied egg-rr80.7%
Taylor expanded in x around inf 72.3%
mul-1-neg72.3%
log-rec72.3%
remove-double-neg72.3%
log-prod56.2%
*-commutative56.2%
Simplified56.2%
Taylor expanded in a around inf 98.9%
*-commutative98.9%
Simplified98.9%
if -2.2e11 < a < 2.0999999999999999e-32Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 66.7%
Taylor expanded in y around inf 66.7%
mul-1-neg66.7%
log-rec66.7%
remove-double-neg66.7%
log-prod52.7%
Simplified52.7%
Final simplification79.2%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -2e-9) (not (<= a 1.22e-32))) (- (* a (log t)) t) (- (log (* y z)) (+ t (* (log t) 0.5)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -2e-9) || !(a <= 1.22e-32)) {
tmp = (a * log(t)) - t;
} else {
tmp = log((y * z)) - (t + (log(t) * 0.5));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if ((a <= (-2d-9)) .or. (.not. (a <= 1.22d-32))) then
tmp = (a * log(t)) - t
else
tmp = log((y * z)) - (t + (log(t) * 0.5d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -2e-9) || !(a <= 1.22e-32)) {
tmp = (a * Math.log(t)) - t;
} else {
tmp = Math.log((y * z)) - (t + (Math.log(t) * 0.5));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -2e-9) or not (a <= 1.22e-32): tmp = (a * math.log(t)) - t else: tmp = math.log((y * z)) - (t + (math.log(t) * 0.5)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -2e-9) || !(a <= 1.22e-32)) tmp = Float64(Float64(a * log(t)) - t); else tmp = Float64(log(Float64(y * z)) - Float64(t + Float64(log(t) * 0.5))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -2e-9) || ~((a <= 1.22e-32))) tmp = (a * log(t)) - t; else tmp = log((y * z)) - (t + (log(t) * 0.5)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -2e-9], N[Not[LessEqual[a, 1.22e-32]], $MachinePrecision]], N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[Log[N[(y * z), $MachinePrecision]], $MachinePrecision] - N[(t + N[(N[Log[t], $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2 \cdot 10^{-9} \lor \neg \left(a \leq 1.22 \cdot 10^{-32}\right):\\
\;\;\;\;a \cdot \log t - t\\
\mathbf{else}:\\
\;\;\;\;\log \left(y \cdot z\right) - \left(t + \log t \cdot 0.5\right)\\
\end{array}
\end{array}
if a < -2.00000000000000012e-9 or 1.22e-32 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.5%
Applied egg-rr80.5%
Taylor expanded in x around inf 71.2%
mul-1-neg71.2%
log-rec71.2%
remove-double-neg71.2%
log-prod55.5%
*-commutative55.5%
Simplified55.5%
Taylor expanded in a around inf 98.1%
*-commutative98.1%
Simplified98.1%
if -2.00000000000000012e-9 < a < 1.22e-32Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 65.5%
add-cube-cbrt65.3%
pow365.3%
Applied egg-rr65.3%
rem-cube-cbrt65.5%
*-commutative65.5%
flip--65.5%
associate-*l/65.5%
metadata-eval65.5%
pow265.5%
Applied egg-rr65.5%
Taylor expanded in a around 0 65.5%
log-prod51.9%
*-commutative51.9%
Simplified51.9%
Final simplification79.1%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -1.6) (not (<= a 2.05))) (- (* a (log t)) t) (- (+ (log z) (log y)) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -1.6) || !(a <= 2.05)) {
tmp = (a * log(t)) - t;
} else {
tmp = (log(z) + log(y)) - t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if ((a <= (-1.6d0)) .or. (.not. (a <= 2.05d0))) then
tmp = (a * log(t)) - t
else
tmp = (log(z) + log(y)) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -1.6) || !(a <= 2.05)) {
tmp = (a * Math.log(t)) - t;
} else {
tmp = (Math.log(z) + Math.log(y)) - t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -1.6) or not (a <= 2.05): tmp = (a * math.log(t)) - t else: tmp = (math.log(z) + math.log(y)) - t return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -1.6) || !(a <= 2.05)) tmp = Float64(Float64(a * log(t)) - t); else tmp = Float64(Float64(log(z) + log(y)) - t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -1.6) || ~((a <= 2.05))) tmp = (a * log(t)) - t; else tmp = (log(z) + log(y)) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -1.6], N[Not[LessEqual[a, 2.05]], $MachinePrecision]], N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[(N[Log[z], $MachinePrecision] + N[Log[y], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.6 \lor \neg \left(a \leq 2.05\right):\\
\;\;\;\;a \cdot \log t - t\\
\mathbf{else}:\\
\;\;\;\;\left(\log z + \log y\right) - t\\
\end{array}
\end{array}
if a < -1.6000000000000001 or 2.0499999999999998 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log81.2%
Applied egg-rr81.2%
Taylor expanded in x around inf 71.7%
mul-1-neg71.7%
log-rec71.7%
remove-double-neg71.7%
log-prod56.2%
*-commutative56.2%
Simplified56.2%
Taylor expanded in a around inf 99.3%
*-commutative99.3%
Simplified99.3%
if -1.6000000000000001 < a < 2.0499999999999998Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 66.4%
Taylor expanded in t around inf 42.4%
Final simplification74.9%
(FPCore (x y z t a) :precision binary64 (if (or (<= a -4.1e+69) (not (<= a 3.4e+16))) (* a (log t)) (- (log y) t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -4.1e+69) || !(a <= 3.4e+16)) {
tmp = a * log(t);
} else {
tmp = log(y) - t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if ((a <= (-4.1d+69)) .or. (.not. (a <= 3.4d+16))) then
tmp = a * log(t)
else
tmp = log(y) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a <= -4.1e+69) || !(a <= 3.4e+16)) {
tmp = a * Math.log(t);
} else {
tmp = Math.log(y) - t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a <= -4.1e+69) or not (a <= 3.4e+16): tmp = a * math.log(t) else: tmp = math.log(y) - t return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((a <= -4.1e+69) || !(a <= 3.4e+16)) tmp = Float64(a * log(t)); else tmp = Float64(log(y) - t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a <= -4.1e+69) || ~((a <= 3.4e+16))) tmp = a * log(t); else tmp = log(y) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[a, -4.1e+69], N[Not[LessEqual[a, 3.4e+16]], $MachinePrecision]], N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision], N[(N[Log[y], $MachinePrecision] - t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -4.1 \cdot 10^{+69} \lor \neg \left(a \leq 3.4 \cdot 10^{+16}\right):\\
\;\;\;\;a \cdot \log t\\
\mathbf{else}:\\
\;\;\;\;\log y - t\\
\end{array}
\end{array}
if a < -4.0999999999999999e69 or 3.4e16 < a Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 76.3%
Taylor expanded in a around inf 81.5%
if -4.0999999999999999e69 < a < 3.4e16Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in t around inf 55.4%
neg-mul-155.4%
Simplified55.4%
Taylor expanded in x around 0 41.9%
Final simplification62.3%
(FPCore (x y z t a) :precision binary64 (if (<= t 246.0) (log (+ x y)) (- t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 246.0) {
tmp = log((x + y));
} else {
tmp = -t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if (t <= 246.0d0) then
tmp = log((x + y))
else
tmp = -t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 246.0) {
tmp = Math.log((x + y));
} else {
tmp = -t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 246.0: tmp = math.log((x + y)) else: tmp = -t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 246.0) tmp = log(Float64(x + y)); else tmp = Float64(-t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 246.0) tmp = log((x + y)); else tmp = -t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 246.0], N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision], (-t)]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 246:\\
\;\;\;\;\log \left(x + y\right)\\
\mathbf{else}:\\
\;\;\;\;-t\\
\end{array}
\end{array}
if t < 246Initial program 99.5%
associate-+l-99.5%
associate--l+99.5%
sub-neg99.5%
+-commutative99.5%
*-commutative99.5%
distribute-rgt-neg-in99.5%
fma-udef99.5%
sub-neg99.5%
+-commutative99.5%
distribute-neg-in99.5%
metadata-eval99.5%
metadata-eval99.5%
unsub-neg99.5%
Simplified99.5%
Taylor expanded in t around inf 7.9%
neg-mul-17.9%
Simplified7.9%
Taylor expanded in t around 0 7.9%
+-commutative7.9%
Simplified7.9%
if 246 < t Initial program 99.9%
associate-+l-99.9%
associate--l+99.9%
sub-neg99.9%
+-commutative99.9%
*-commutative99.9%
distribute-rgt-neg-in99.9%
fma-udef99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in x around 0 74.3%
Taylor expanded in t around inf 68.5%
mul-1-neg68.5%
Simplified68.5%
Final simplification36.3%
(FPCore (x y z t a) :precision binary64 (if (<= t 2.6e+58) (* a (log t)) (- t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 2.6e+58) {
tmp = a * log(t);
} else {
tmp = -t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
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) :: tmp
if (t <= 2.6d+58) then
tmp = a * log(t)
else
tmp = -t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 2.6e+58) {
tmp = a * Math.log(t);
} else {
tmp = -t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 2.6e+58: tmp = a * math.log(t) else: tmp = -t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 2.6e+58) tmp = Float64(a * log(t)); else tmp = Float64(-t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 2.6e+58) tmp = a * log(t); else tmp = -t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 2.6e+58], N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision], (-t)]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 2.6 \cdot 10^{+58}:\\
\;\;\;\;a \cdot \log t\\
\mathbf{else}:\\
\;\;\;\;-t\\
\end{array}
\end{array}
if t < 2.59999999999999988e58Initial program 99.5%
associate-+l-99.5%
associate--l+99.5%
sub-neg99.5%
+-commutative99.5%
*-commutative99.5%
distribute-rgt-neg-in99.5%
fma-udef99.5%
sub-neg99.5%
+-commutative99.5%
distribute-neg-in99.5%
metadata-eval99.5%
metadata-eval99.5%
unsub-neg99.5%
Simplified99.5%
Taylor expanded in x around 0 71.5%
Taylor expanded in a around inf 60.4%
if 2.59999999999999988e58 < t Initial program 99.9%
associate-+l-99.9%
associate--l+99.9%
sub-neg99.9%
+-commutative99.9%
*-commutative99.9%
distribute-rgt-neg-in99.9%
fma-udef99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in x around 0 73.9%
Taylor expanded in t around inf 77.1%
mul-1-neg77.1%
Simplified77.1%
Final simplification66.9%
(FPCore (x y z t a) :precision binary64 (- (* a (log t)) t))
double code(double x, double y, double z, double t, double a) {
return (a * log(t)) - t;
}
real(8) function code(x, y, z, t, a)
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
code = (a * log(t)) - t
end function
public static double code(double x, double y, double z, double t, double a) {
return (a * Math.log(t)) - t;
}
def code(x, y, z, t, a): return (a * math.log(t)) - t
function code(x, y, z, t, a) return Float64(Float64(a * log(t)) - t) end
function tmp = code(x, y, z, t, a) tmp = (a * log(t)) - t; end
code[x_, y_, z_, t_, a_] := N[(N[(a * N[Log[t], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
a \cdot \log t - t
\end{array}
Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
associate-+r-99.7%
fma-udef99.7%
associate--r+99.7%
+-commutative99.7%
sum-log80.4%
Applied egg-rr80.4%
Taylor expanded in x around inf 67.3%
mul-1-neg67.3%
log-rec67.3%
remove-double-neg67.3%
log-prod53.0%
*-commutative53.0%
Simplified53.0%
Taylor expanded in a around inf 77.6%
*-commutative77.6%
Simplified77.6%
Final simplification77.6%
(FPCore (x y z t a) :precision binary64 (- t))
double code(double x, double y, double z, double t, double a) {
return -t;
}
real(8) function code(x, y, z, t, a)
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
code = -t
end function
public static double code(double x, double y, double z, double t, double a) {
return -t;
}
def code(x, y, z, t, a): return -t
function code(x, y, z, t, a) return Float64(-t) end
function tmp = code(x, y, z, t, a) tmp = -t; end
code[x_, y_, z_, t_, a_] := (-t)
\begin{array}{l}
\\
-t
\end{array}
Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-udef99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around 0 72.4%
Taylor expanded in t around inf 33.5%
mul-1-neg33.5%
Simplified33.5%
Final simplification33.5%
(FPCore (x y z t a) :precision binary64 (+ (log (+ x y)) (+ (- (log z) t) (* (- a 0.5) (log t)))))
double code(double x, double y, double z, double t, double a) {
return log((x + y)) + ((log(z) - t) + ((a - 0.5) * log(t)));
}
real(8) function code(x, y, z, t, a)
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
code = log((x + y)) + ((log(z) - t) + ((a - 0.5d0) * log(t)))
end function
public static double code(double x, double y, double z, double t, double a) {
return Math.log((x + y)) + ((Math.log(z) - t) + ((a - 0.5) * Math.log(t)));
}
def code(x, y, z, t, a): return math.log((x + y)) + ((math.log(z) - t) + ((a - 0.5) * math.log(t)))
function code(x, y, z, t, a) return Float64(log(Float64(x + y)) + Float64(Float64(log(z) - t) + Float64(Float64(a - 0.5) * log(t)))) end
function tmp = code(x, y, z, t, a) tmp = log((x + y)) + ((log(z) - t) + ((a - 0.5) * log(t))); end
code[x_, y_, z_, t_, a_] := N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision] + N[(N[(a - 0.5), $MachinePrecision] * N[Log[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\log \left(x + y\right) + \left(\left(\log z - t\right) + \left(a - 0.5\right) \cdot \log t\right)
\end{array}
herbie shell --seed 2024031
(FPCore (x y z t a)
:name "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(+ (log (+ x y)) (+ (- (log z) t) (* (- a 0.5) (log t))))
(+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))