
(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 17 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 (+ (log (+ x y)) (- (log z) (fma (log t) (- 0.5 a) t))))
double code(double x, double y, double z, double t, double a) {
return log((x + y)) + (log(z) - fma(log(t), (0.5 - a), t));
}
function code(x, y, z, t, a) return Float64(log(Float64(x + y)) + Float64(log(z) - fma(log(t), Float64(0.5 - a), t))) end
code[x_, y_, z_, t_, a_] := N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - N[(N[Log[t], $MachinePrecision] * N[(0.5 - a), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\log \left(x + y\right) + \left(\log z - \mathsf{fma}\left(\log t, 0.5 - a, t\right)\right)
\end{array}
Initial program 99.6%
associate-+l-99.6%
associate--l+99.6%
sub-neg99.6%
+-commutative99.6%
*-commutative99.6%
distribute-rgt-neg-in99.6%
fma-undefine99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
metadata-eval99.6%
metadata-eval99.6%
unsub-neg99.6%
Simplified99.6%
(FPCore (x y z t a) :precision binary64 (if (<= t 0.43) (+ (log y) (- (log z) (* (log t) (- 0.5 a)))) (+ (log y) (- (+ (log z) (* (log t) a)) t))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 0.43) {
tmp = log(y) + (log(z) - (log(t) * (0.5 - a)));
} else {
tmp = log(y) + ((log(z) + (log(t) * a)) - 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 <= 0.43d0) then
tmp = log(y) + (log(z) - (log(t) * (0.5d0 - a)))
else
tmp = log(y) + ((log(z) + (log(t) * a)) - 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 <= 0.43) {
tmp = Math.log(y) + (Math.log(z) - (Math.log(t) * (0.5 - a)));
} else {
tmp = Math.log(y) + ((Math.log(z) + (Math.log(t) * a)) - t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 0.43: tmp = math.log(y) + (math.log(z) - (math.log(t) * (0.5 - a))) else: tmp = math.log(y) + ((math.log(z) + (math.log(t) * a)) - t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 0.43) tmp = Float64(log(y) + Float64(log(z) - Float64(log(t) * Float64(0.5 - a)))); else tmp = Float64(log(y) + Float64(Float64(log(z) + Float64(log(t) * a)) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 0.43) tmp = log(y) + (log(z) - (log(t) * (0.5 - a))); else tmp = log(y) + ((log(z) + (log(t) * a)) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 0.43], N[(N[Log[y], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - N[(N[Log[t], $MachinePrecision] * N[(0.5 - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Log[y], $MachinePrecision] + N[(N[(N[Log[z], $MachinePrecision] + N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 0.43:\\
\;\;\;\;\log y + \left(\log z - \log t \cdot \left(0.5 - a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\log y + \left(\left(\log z + \log t \cdot a\right) - t\right)\\
\end{array}
\end{array}
if t < 0.429999999999999993Initial program 99.4%
associate--l+99.4%
+-commutative99.4%
associate-+l+99.3%
+-commutative99.3%
fma-define99.3%
sub-neg99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around 0 69.4%
associate--l+69.4%
remove-double-neg69.4%
log-rec69.4%
mul-1-neg69.4%
mul-1-neg69.4%
log-rec69.4%
remove-double-neg69.4%
sub-neg69.4%
metadata-eval69.4%
Simplified69.4%
Taylor expanded in t around 0 68.5%
if 0.429999999999999993 < t Initial program 99.9%
associate--l+99.9%
+-commutative99.9%
associate-+l+99.9%
+-commutative99.9%
fma-define99.9%
sub-neg99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 71.5%
associate--l+71.5%
remove-double-neg71.5%
log-rec71.5%
mul-1-neg71.5%
mul-1-neg71.5%
log-rec71.5%
remove-double-neg71.5%
sub-neg71.5%
metadata-eval71.5%
Simplified71.5%
add-cube-cbrt71.3%
pow371.3%
Applied egg-rr71.3%
Taylor expanded in a around inf 71.5%
*-commutative71.5%
Simplified71.5%
Final simplification70.1%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (log t) (- 0.5 a))))
(if (<= t 1.9e-6)
(+ (log y) (- (log z) t_1))
(if (<= t 1.4e+99)
(- (- (log (* y z)) t_1) t)
(+ (log (+ x y)) (- (log z) t))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = log(t) * (0.5 - a);
double tmp;
if (t <= 1.9e-6) {
tmp = log(y) + (log(z) - t_1);
} else if (t <= 1.4e+99) {
tmp = (log((y * z)) - t_1) - t;
} else {
tmp = log((x + y)) + (log(z) - 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) :: t_1
real(8) :: tmp
t_1 = log(t) * (0.5d0 - a)
if (t <= 1.9d-6) then
tmp = log(y) + (log(z) - t_1)
else if (t <= 1.4d+99) then
tmp = (log((y * z)) - t_1) - t
else
tmp = log((x + y)) + (log(z) - t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.log(t) * (0.5 - a);
double tmp;
if (t <= 1.9e-6) {
tmp = Math.log(y) + (Math.log(z) - t_1);
} else if (t <= 1.4e+99) {
tmp = (Math.log((y * z)) - t_1) - t;
} else {
tmp = Math.log((x + y)) + (Math.log(z) - t);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.log(t) * (0.5 - a) tmp = 0 if t <= 1.9e-6: tmp = math.log(y) + (math.log(z) - t_1) elif t <= 1.4e+99: tmp = (math.log((y * z)) - t_1) - t else: tmp = math.log((x + y)) + (math.log(z) - t) return tmp
function code(x, y, z, t, a) t_1 = Float64(log(t) * Float64(0.5 - a)) tmp = 0.0 if (t <= 1.9e-6) tmp = Float64(log(y) + Float64(log(z) - t_1)); elseif (t <= 1.4e+99) tmp = Float64(Float64(log(Float64(y * z)) - t_1) - t); else tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = log(t) * (0.5 - a); tmp = 0.0; if (t <= 1.9e-6) tmp = log(y) + (log(z) - t_1); elseif (t <= 1.4e+99) tmp = (log((y * z)) - t_1) - t; else tmp = log((x + y)) + (log(z) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[Log[t], $MachinePrecision] * N[(0.5 - a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, 1.9e-6], N[(N[Log[y], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.4e+99], N[(N[(N[Log[N[(y * z), $MachinePrecision]], $MachinePrecision] - t$95$1), $MachinePrecision] - t), $MachinePrecision], N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log t \cdot \left(0.5 - a\right)\\
\mathbf{if}\;t \leq 1.9 \cdot 10^{-6}:\\
\;\;\;\;\log y + \left(\log z - t\_1\right)\\
\mathbf{elif}\;t \leq 1.4 \cdot 10^{+99}:\\
\;\;\;\;\left(\log \left(y \cdot z\right) - t\_1\right) - t\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\end{array}
\end{array}
if t < 1.9e-6Initial program 99.4%
associate--l+99.4%
+-commutative99.4%
associate-+l+99.3%
+-commutative99.3%
fma-define99.3%
sub-neg99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around 0 68.9%
associate--l+68.9%
remove-double-neg68.9%
log-rec68.9%
mul-1-neg68.9%
mul-1-neg68.9%
log-rec68.9%
remove-double-neg68.9%
sub-neg68.9%
metadata-eval68.9%
Simplified68.9%
Taylor expanded in t around 0 68.9%
if 1.9e-6 < t < 1.4e99Initial program 99.7%
associate--l+99.6%
+-commutative99.6%
associate-+l+99.7%
+-commutative99.7%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
add-sqr-sqrt30.0%
pow230.0%
Applied egg-rr24.3%
Taylor expanded in x around 0 61.2%
if 1.4e99 < 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-undefine100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
Taylor expanded in t around inf 83.6%
Final simplification73.8%
(FPCore (x y z t a) :precision binary64 (+ (- (+ (log (+ x y)) (log z)) t) (* (log t) (- a 0.5))))
double code(double x, double y, double z, double t, double a) {
return ((log((x + y)) + log(z)) - t) + (log(t) * (a - 0.5));
}
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) + (log(t) * (a - 0.5d0))
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) + (Math.log(t) * (a - 0.5));
}
def code(x, y, z, t, a): return ((math.log((x + y)) + math.log(z)) - t) + (math.log(t) * (a - 0.5))
function code(x, y, z, t, a) return Float64(Float64(Float64(log(Float64(x + y)) + log(z)) - t) + Float64(log(t) * Float64(a - 0.5))) end
function tmp = code(x, y, z, t, a) tmp = ((log((x + y)) + log(z)) - t) + (log(t) * (a - 0.5)); 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[Log[t], $MachinePrecision] * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \log t \cdot \left(a - 0.5\right)
\end{array}
Initial program 99.6%
Final simplification99.6%
(FPCore (x y z t a) :precision binary64 (+ (log y) (- (+ (log z) (* (log t) (+ a -0.5))) t)))
double code(double x, double y, double z, double t, double a) {
return log(y) + ((log(z) + (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(y) + ((log(z) + (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(y) + ((Math.log(z) + (Math.log(t) * (a + -0.5))) - t);
}
def code(x, y, z, t, a): return math.log(y) + ((math.log(z) + (math.log(t) * (a + -0.5))) - t)
function code(x, y, z, t, a) return Float64(log(y) + Float64(Float64(log(z) + Float64(log(t) * Float64(a + -0.5))) - t)) end
function tmp = code(x, y, z, t, a) tmp = log(y) + ((log(z) + (log(t) * (a + -0.5))) - t); end
code[x_, y_, z_, t_, a_] := N[(N[Log[y], $MachinePrecision] + N[(N[(N[Log[z], $MachinePrecision] + N[(N[Log[t], $MachinePrecision] * N[(a + -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\log y + \left(\left(\log z + \log t \cdot \left(a + -0.5\right)\right) - t\right)
\end{array}
Initial program 99.6%
associate--l+99.6%
+-commutative99.6%
associate-+l+99.6%
+-commutative99.6%
fma-define99.6%
sub-neg99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around 0 70.5%
associate--l+70.5%
remove-double-neg70.5%
log-rec70.5%
mul-1-neg70.5%
mul-1-neg70.5%
log-rec70.5%
remove-double-neg70.5%
sub-neg70.5%
metadata-eval70.5%
Simplified70.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (log t) a))
(t_2 (- (log (* y (* z (pow t (+ a -0.5))))) t))
(t_3 (+ (log (+ x y)) (- (log z) t))))
(if (<= a -3.8e+183)
t_1
(if (<= a -3.9e+100)
t_3
(if (<= a -1.45)
t_1
(if (<= a 9.8e-293)
t_2
(if (<= a 2.7e-179) t_3 (if (<= a 6.4) t_2 t_1))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = log(t) * a;
double t_2 = log((y * (z * pow(t, (a + -0.5))))) - t;
double t_3 = log((x + y)) + (log(z) - t);
double tmp;
if (a <= -3.8e+183) {
tmp = t_1;
} else if (a <= -3.9e+100) {
tmp = t_3;
} else if (a <= -1.45) {
tmp = t_1;
} else if (a <= 9.8e-293) {
tmp = t_2;
} else if (a <= 2.7e-179) {
tmp = t_3;
} else if (a <= 6.4) {
tmp = t_2;
} 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) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = log(t) * a
t_2 = log((y * (z * (t ** (a + (-0.5d0)))))) - t
t_3 = log((x + y)) + (log(z) - t)
if (a <= (-3.8d+183)) then
tmp = t_1
else if (a <= (-3.9d+100)) then
tmp = t_3
else if (a <= (-1.45d0)) then
tmp = t_1
else if (a <= 9.8d-293) then
tmp = t_2
else if (a <= 2.7d-179) then
tmp = t_3
else if (a <= 6.4d0) then
tmp = t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.log(t) * a;
double t_2 = Math.log((y * (z * Math.pow(t, (a + -0.5))))) - t;
double t_3 = Math.log((x + y)) + (Math.log(z) - t);
double tmp;
if (a <= -3.8e+183) {
tmp = t_1;
} else if (a <= -3.9e+100) {
tmp = t_3;
} else if (a <= -1.45) {
tmp = t_1;
} else if (a <= 9.8e-293) {
tmp = t_2;
} else if (a <= 2.7e-179) {
tmp = t_3;
} else if (a <= 6.4) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.log(t) * a t_2 = math.log((y * (z * math.pow(t, (a + -0.5))))) - t t_3 = math.log((x + y)) + (math.log(z) - t) tmp = 0 if a <= -3.8e+183: tmp = t_1 elif a <= -3.9e+100: tmp = t_3 elif a <= -1.45: tmp = t_1 elif a <= 9.8e-293: tmp = t_2 elif a <= 2.7e-179: tmp = t_3 elif a <= 6.4: tmp = t_2 else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(log(t) * a) t_2 = Float64(log(Float64(y * Float64(z * (t ^ Float64(a + -0.5))))) - t) t_3 = Float64(log(Float64(x + y)) + Float64(log(z) - t)) tmp = 0.0 if (a <= -3.8e+183) tmp = t_1; elseif (a <= -3.9e+100) tmp = t_3; elseif (a <= -1.45) tmp = t_1; elseif (a <= 9.8e-293) tmp = t_2; elseif (a <= 2.7e-179) tmp = t_3; elseif (a <= 6.4) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = log(t) * a; t_2 = log((y * (z * (t ^ (a + -0.5))))) - t; t_3 = log((x + y)) + (log(z) - t); tmp = 0.0; if (a <= -3.8e+183) tmp = t_1; elseif (a <= -3.9e+100) tmp = t_3; elseif (a <= -1.45) tmp = t_1; elseif (a <= 9.8e-293) tmp = t_2; elseif (a <= 2.7e-179) tmp = t_3; elseif (a <= 6.4) tmp = t_2; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]}, Block[{t$95$2 = N[(N[Log[N[(y * N[(z * N[Power[t, N[(a + -0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - t), $MachinePrecision]}, Block[{t$95$3 = N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -3.8e+183], t$95$1, If[LessEqual[a, -3.9e+100], t$95$3, If[LessEqual[a, -1.45], t$95$1, If[LessEqual[a, 9.8e-293], t$95$2, If[LessEqual[a, 2.7e-179], t$95$3, If[LessEqual[a, 6.4], t$95$2, t$95$1]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log t \cdot a\\
t_2 := \log \left(y \cdot \left(z \cdot {t}^{\left(a + -0.5\right)}\right)\right) - t\\
t_3 := \log \left(x + y\right) + \left(\log z - t\right)\\
\mathbf{if}\;a \leq -3.8 \cdot 10^{+183}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq -3.9 \cdot 10^{+100}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;a \leq -1.45:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 9.8 \cdot 10^{-293}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;a \leq 2.7 \cdot 10^{-179}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;a \leq 6.4:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -3.80000000000000001e183 or -3.9e100 < a < -1.44999999999999996 or 6.4000000000000004 < 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-undefine99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in a around inf 74.8%
*-commutative74.8%
Simplified74.8%
if -3.80000000000000001e183 < a < -3.9e100 or 9.8000000000000001e-293 < a < 2.69999999999999988e-179Initial program 99.7%
associate-+l-99.7%
associate--l+99.8%
sub-neg99.8%
+-commutative99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
fma-undefine99.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 69.1%
if -1.44999999999999996 < a < 9.8000000000000001e-293 or 2.69999999999999988e-179 < a < 6.4000000000000004Initial program 99.5%
associate--l+99.5%
+-commutative99.5%
associate-+l+99.5%
+-commutative99.5%
fma-define99.5%
sub-neg99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around 0 58.6%
associate--l+58.6%
remove-double-neg58.6%
log-rec58.6%
mul-1-neg58.6%
mul-1-neg58.6%
log-rec58.6%
remove-double-neg58.6%
sub-neg58.6%
metadata-eval58.6%
Simplified58.6%
associate-+r-58.6%
add-log-exp51.8%
sum-log36.3%
exp-sum36.2%
add-exp-log36.3%
exp-to-pow36.3%
Applied egg-rr36.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (log t) a)) (t_2 (+ (log (+ x y)) (- (log z) t))))
(if (<= a -3.8e+183)
t_1
(if (<= a -5.5e+100)
t_2
(if (<= a -5.5e+47)
t_1
(if (<= a 6.8e-72)
t_2
(if (<= a 2.5) (log (* z (* (+ x y) (pow t (+ a -0.5))))) t_1)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = log(t) * a;
double t_2 = log((x + y)) + (log(z) - t);
double tmp;
if (a <= -3.8e+183) {
tmp = t_1;
} else if (a <= -5.5e+100) {
tmp = t_2;
} else if (a <= -5.5e+47) {
tmp = t_1;
} else if (a <= 6.8e-72) {
tmp = t_2;
} else if (a <= 2.5) {
tmp = log((z * ((x + y) * pow(t, (a + -0.5)))));
} 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) :: t_2
real(8) :: tmp
t_1 = log(t) * a
t_2 = log((x + y)) + (log(z) - t)
if (a <= (-3.8d+183)) then
tmp = t_1
else if (a <= (-5.5d+100)) then
tmp = t_2
else if (a <= (-5.5d+47)) then
tmp = t_1
else if (a <= 6.8d-72) then
tmp = t_2
else if (a <= 2.5d0) then
tmp = log((z * ((x + y) * (t ** (a + (-0.5d0))))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.log(t) * a;
double t_2 = Math.log((x + y)) + (Math.log(z) - t);
double tmp;
if (a <= -3.8e+183) {
tmp = t_1;
} else if (a <= -5.5e+100) {
tmp = t_2;
} else if (a <= -5.5e+47) {
tmp = t_1;
} else if (a <= 6.8e-72) {
tmp = t_2;
} else if (a <= 2.5) {
tmp = Math.log((z * ((x + y) * Math.pow(t, (a + -0.5)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.log(t) * a t_2 = math.log((x + y)) + (math.log(z) - t) tmp = 0 if a <= -3.8e+183: tmp = t_1 elif a <= -5.5e+100: tmp = t_2 elif a <= -5.5e+47: tmp = t_1 elif a <= 6.8e-72: tmp = t_2 elif a <= 2.5: tmp = math.log((z * ((x + y) * math.pow(t, (a + -0.5))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(log(t) * a) t_2 = Float64(log(Float64(x + y)) + Float64(log(z) - t)) tmp = 0.0 if (a <= -3.8e+183) tmp = t_1; elseif (a <= -5.5e+100) tmp = t_2; elseif (a <= -5.5e+47) tmp = t_1; elseif (a <= 6.8e-72) tmp = t_2; elseif (a <= 2.5) tmp = log(Float64(z * Float64(Float64(x + y) * (t ^ Float64(a + -0.5))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = log(t) * a; t_2 = log((x + y)) + (log(z) - t); tmp = 0.0; if (a <= -3.8e+183) tmp = t_1; elseif (a <= -5.5e+100) tmp = t_2; elseif (a <= -5.5e+47) tmp = t_1; elseif (a <= 6.8e-72) tmp = t_2; elseif (a <= 2.5) tmp = log((z * ((x + y) * (t ^ (a + -0.5))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]}, Block[{t$95$2 = N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -3.8e+183], t$95$1, If[LessEqual[a, -5.5e+100], t$95$2, If[LessEqual[a, -5.5e+47], t$95$1, If[LessEqual[a, 6.8e-72], t$95$2, If[LessEqual[a, 2.5], N[Log[N[(z * N[(N[(x + y), $MachinePrecision] * N[Power[t, N[(a + -0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log t \cdot a\\
t_2 := \log \left(x + y\right) + \left(\log z - t\right)\\
\mathbf{if}\;a \leq -3.8 \cdot 10^{+183}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq -5.5 \cdot 10^{+100}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;a \leq -5.5 \cdot 10^{+47}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 6.8 \cdot 10^{-72}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;a \leq 2.5:\\
\;\;\;\;\log \left(z \cdot \left(\left(x + y\right) \cdot {t}^{\left(a + -0.5\right)}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -3.80000000000000001e183 or -5.5000000000000002e100 < a < -5.4999999999999998e47 or 2.5 < a Initial program 99.6%
associate-+l-99.6%
associate--l+99.6%
sub-neg99.6%
+-commutative99.6%
*-commutative99.6%
distribute-rgt-neg-in99.6%
fma-undefine99.7%
sub-neg99.7%
+-commutative99.7%
distribute-neg-in99.7%
metadata-eval99.7%
metadata-eval99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in a around inf 81.1%
*-commutative81.1%
Simplified81.1%
if -3.80000000000000001e183 < a < -5.5000000000000002e100 or -5.4999999999999998e47 < a < 6.7999999999999997e-72Initial program 99.7%
associate-+l-99.7%
associate--l+99.7%
sub-neg99.7%
+-commutative99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
fma-undefine99.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 61.7%
if 6.7999999999999997e-72 < a < 2.5Initial program 98.7%
associate--l+98.7%
+-commutative98.7%
associate-+l+98.8%
+-commutative98.8%
fma-define98.8%
sub-neg98.8%
metadata-eval98.8%
Simplified98.8%
+-commutative98.8%
fma-undefine98.8%
metadata-eval98.8%
sub-neg98.8%
associate-+r+98.7%
associate--l+98.7%
add-log-exp73.3%
+-commutative73.3%
exp-sum64.5%
Applied egg-rr65.8%
associate-*l/65.8%
*-commutative65.8%
+-commutative65.8%
Simplified65.8%
Taylor expanded in t around 0 73.9%
*-commutative73.9%
+-commutative73.9%
exp-to-pow74.7%
sub-neg74.7%
metadata-eval74.7%
Simplified74.7%
Final simplification69.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (pow t (+ a -0.5))) (t_2 (* (log t) a)))
(if (<= t 7.5e-208)
t_2
(if (<= t 4.2e-173)
(log (* (* y z) t_1))
(if (<= t 6.5e-104)
t_2
(if (<= t 2.85e-43)
(log (* y (* z t_1)))
(if (<= t 1.15e+27) t_2 (+ (log (+ x y)) (- (log z) t)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = pow(t, (a + -0.5));
double t_2 = log(t) * a;
double tmp;
if (t <= 7.5e-208) {
tmp = t_2;
} else if (t <= 4.2e-173) {
tmp = log(((y * z) * t_1));
} else if (t <= 6.5e-104) {
tmp = t_2;
} else if (t <= 2.85e-43) {
tmp = log((y * (z * t_1)));
} else if (t <= 1.15e+27) {
tmp = t_2;
} else {
tmp = log((x + y)) + (log(z) - 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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = t ** (a + (-0.5d0))
t_2 = log(t) * a
if (t <= 7.5d-208) then
tmp = t_2
else if (t <= 4.2d-173) then
tmp = log(((y * z) * t_1))
else if (t <= 6.5d-104) then
tmp = t_2
else if (t <= 2.85d-43) then
tmp = log((y * (z * t_1)))
else if (t <= 1.15d+27) then
tmp = t_2
else
tmp = log((x + y)) + (log(z) - t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.pow(t, (a + -0.5));
double t_2 = Math.log(t) * a;
double tmp;
if (t <= 7.5e-208) {
tmp = t_2;
} else if (t <= 4.2e-173) {
tmp = Math.log(((y * z) * t_1));
} else if (t <= 6.5e-104) {
tmp = t_2;
} else if (t <= 2.85e-43) {
tmp = Math.log((y * (z * t_1)));
} else if (t <= 1.15e+27) {
tmp = t_2;
} else {
tmp = Math.log((x + y)) + (Math.log(z) - t);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.pow(t, (a + -0.5)) t_2 = math.log(t) * a tmp = 0 if t <= 7.5e-208: tmp = t_2 elif t <= 4.2e-173: tmp = math.log(((y * z) * t_1)) elif t <= 6.5e-104: tmp = t_2 elif t <= 2.85e-43: tmp = math.log((y * (z * t_1))) elif t <= 1.15e+27: tmp = t_2 else: tmp = math.log((x + y)) + (math.log(z) - t) return tmp
function code(x, y, z, t, a) t_1 = t ^ Float64(a + -0.5) t_2 = Float64(log(t) * a) tmp = 0.0 if (t <= 7.5e-208) tmp = t_2; elseif (t <= 4.2e-173) tmp = log(Float64(Float64(y * z) * t_1)); elseif (t <= 6.5e-104) tmp = t_2; elseif (t <= 2.85e-43) tmp = log(Float64(y * Float64(z * t_1))); elseif (t <= 1.15e+27) tmp = t_2; else tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = t ^ (a + -0.5); t_2 = log(t) * a; tmp = 0.0; if (t <= 7.5e-208) tmp = t_2; elseif (t <= 4.2e-173) tmp = log(((y * z) * t_1)); elseif (t <= 6.5e-104) tmp = t_2; elseif (t <= 2.85e-43) tmp = log((y * (z * t_1))); elseif (t <= 1.15e+27) tmp = t_2; else tmp = log((x + y)) + (log(z) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[Power[t, N[(a + -0.5), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]}, If[LessEqual[t, 7.5e-208], t$95$2, If[LessEqual[t, 4.2e-173], N[Log[N[(N[(y * z), $MachinePrecision] * t$95$1), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 6.5e-104], t$95$2, If[LessEqual[t, 2.85e-43], N[Log[N[(y * N[(z * t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 1.15e+27], t$95$2, N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {t}^{\left(a + -0.5\right)}\\
t_2 := \log t \cdot a\\
\mathbf{if}\;t \leq 7.5 \cdot 10^{-208}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq 4.2 \cdot 10^{-173}:\\
\;\;\;\;\log \left(\left(y \cdot z\right) \cdot t\_1\right)\\
\mathbf{elif}\;t \leq 6.5 \cdot 10^{-104}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq 2.85 \cdot 10^{-43}:\\
\;\;\;\;\log \left(y \cdot \left(z \cdot t\_1\right)\right)\\
\mathbf{elif}\;t \leq 1.15 \cdot 10^{+27}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\end{array}
\end{array}
if t < 7.4999999999999999e-208 or 4.20000000000000003e-173 < t < 6.49999999999999991e-104 or 2.85e-43 < t < 1.15e27Initial program 99.4%
associate-+l-99.4%
associate--l+99.4%
sub-neg99.4%
+-commutative99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
fma-undefine99.4%
sub-neg99.4%
+-commutative99.4%
distribute-neg-in99.4%
metadata-eval99.4%
metadata-eval99.4%
unsub-neg99.4%
Simplified99.4%
Taylor expanded in a around inf 55.6%
*-commutative55.6%
Simplified55.6%
if 7.4999999999999999e-208 < t < 4.20000000000000003e-173Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.7%
+-commutative99.7%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around 0 70.7%
associate--l+70.7%
remove-double-neg70.7%
log-rec70.7%
mul-1-neg70.7%
mul-1-neg70.7%
log-rec70.7%
remove-double-neg70.7%
sub-neg70.7%
metadata-eval70.7%
Simplified70.7%
Taylor expanded in t around 0 70.7%
*-un-lft-identity70.7%
add-log-exp41.9%
sum-log42.3%
sub-neg42.3%
metadata-eval42.3%
exp-sum42.3%
add-exp-log42.3%
pow-to-exp42.4%
Applied egg-rr42.4%
*-lft-identity42.4%
associate-*r*51.7%
Simplified51.7%
if 6.49999999999999991e-104 < t < 2.85e-43Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.1%
+-commutative99.1%
fma-define99.1%
sub-neg99.1%
metadata-eval99.1%
Simplified99.1%
Taylor expanded in x around 0 59.6%
associate--l+59.6%
remove-double-neg59.6%
log-rec59.6%
mul-1-neg59.6%
mul-1-neg59.6%
log-rec59.6%
remove-double-neg59.6%
sub-neg59.6%
metadata-eval59.6%
Simplified59.6%
Taylor expanded in t around 0 59.6%
*-un-lft-identity59.6%
add-log-exp40.7%
sum-log33.7%
sub-neg33.7%
metadata-eval33.7%
exp-sum33.8%
add-exp-log33.8%
pow-to-exp33.8%
Applied egg-rr33.8%
*-lft-identity33.8%
Simplified33.8%
if 1.15e27 < 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-undefine99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in t around inf 78.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (log t) a)))
(if (<= t 8e-208)
t_1
(if (<= t 5.2e-175)
(log (* z (* y (pow t -0.5))))
(if (<= t 5.6e-105)
t_1
(if (<= t 2e-43)
(log (* y (* z (pow t (+ a -0.5)))))
(if (<= t 1.6e+21) t_1 (+ (log (+ x y)) (- (log z) t)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = log(t) * a;
double tmp;
if (t <= 8e-208) {
tmp = t_1;
} else if (t <= 5.2e-175) {
tmp = log((z * (y * pow(t, -0.5))));
} else if (t <= 5.6e-105) {
tmp = t_1;
} else if (t <= 2e-43) {
tmp = log((y * (z * pow(t, (a + -0.5)))));
} else if (t <= 1.6e+21) {
tmp = t_1;
} else {
tmp = log((x + y)) + (log(z) - 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) :: t_1
real(8) :: tmp
t_1 = log(t) * a
if (t <= 8d-208) then
tmp = t_1
else if (t <= 5.2d-175) then
tmp = log((z * (y * (t ** (-0.5d0)))))
else if (t <= 5.6d-105) then
tmp = t_1
else if (t <= 2d-43) then
tmp = log((y * (z * (t ** (a + (-0.5d0))))))
else if (t <= 1.6d+21) then
tmp = t_1
else
tmp = log((x + y)) + (log(z) - t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.log(t) * a;
double tmp;
if (t <= 8e-208) {
tmp = t_1;
} else if (t <= 5.2e-175) {
tmp = Math.log((z * (y * Math.pow(t, -0.5))));
} else if (t <= 5.6e-105) {
tmp = t_1;
} else if (t <= 2e-43) {
tmp = Math.log((y * (z * Math.pow(t, (a + -0.5)))));
} else if (t <= 1.6e+21) {
tmp = t_1;
} else {
tmp = Math.log((x + y)) + (Math.log(z) - t);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.log(t) * a tmp = 0 if t <= 8e-208: tmp = t_1 elif t <= 5.2e-175: tmp = math.log((z * (y * math.pow(t, -0.5)))) elif t <= 5.6e-105: tmp = t_1 elif t <= 2e-43: tmp = math.log((y * (z * math.pow(t, (a + -0.5))))) elif t <= 1.6e+21: tmp = t_1 else: tmp = math.log((x + y)) + (math.log(z) - t) return tmp
function code(x, y, z, t, a) t_1 = Float64(log(t) * a) tmp = 0.0 if (t <= 8e-208) tmp = t_1; elseif (t <= 5.2e-175) tmp = log(Float64(z * Float64(y * (t ^ -0.5)))); elseif (t <= 5.6e-105) tmp = t_1; elseif (t <= 2e-43) tmp = log(Float64(y * Float64(z * (t ^ Float64(a + -0.5))))); elseif (t <= 1.6e+21) tmp = t_1; else tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = log(t) * a; tmp = 0.0; if (t <= 8e-208) tmp = t_1; elseif (t <= 5.2e-175) tmp = log((z * (y * (t ^ -0.5)))); elseif (t <= 5.6e-105) tmp = t_1; elseif (t <= 2e-43) tmp = log((y * (z * (t ^ (a + -0.5))))); elseif (t <= 1.6e+21) tmp = t_1; else tmp = log((x + y)) + (log(z) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]}, If[LessEqual[t, 8e-208], t$95$1, If[LessEqual[t, 5.2e-175], N[Log[N[(z * N[(y * N[Power[t, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 5.6e-105], t$95$1, If[LessEqual[t, 2e-43], N[Log[N[(y * N[(z * N[Power[t, N[(a + -0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 1.6e+21], t$95$1, N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log t \cdot a\\
\mathbf{if}\;t \leq 8 \cdot 10^{-208}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 5.2 \cdot 10^{-175}:\\
\;\;\;\;\log \left(z \cdot \left(y \cdot {t}^{-0.5}\right)\right)\\
\mathbf{elif}\;t \leq 5.6 \cdot 10^{-105}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2 \cdot 10^{-43}:\\
\;\;\;\;\log \left(y \cdot \left(z \cdot {t}^{\left(a + -0.5\right)}\right)\right)\\
\mathbf{elif}\;t \leq 1.6 \cdot 10^{+21}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\end{array}
\end{array}
if t < 8.0000000000000008e-208 or 5.2e-175 < t < 5.6e-105 or 2.00000000000000015e-43 < t < 1.6e21Initial program 99.4%
associate-+l-99.4%
associate--l+99.4%
sub-neg99.4%
+-commutative99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
fma-undefine99.4%
sub-neg99.4%
+-commutative99.4%
distribute-neg-in99.4%
metadata-eval99.4%
metadata-eval99.4%
unsub-neg99.4%
Simplified99.4%
Taylor expanded in a around inf 55.6%
*-commutative55.6%
Simplified55.6%
if 8.0000000000000008e-208 < t < 5.2e-175Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.7%
+-commutative99.7%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around 0 70.7%
associate--l+70.7%
remove-double-neg70.7%
log-rec70.7%
mul-1-neg70.7%
mul-1-neg70.7%
log-rec70.7%
remove-double-neg70.7%
sub-neg70.7%
metadata-eval70.7%
Simplified70.7%
Taylor expanded in t around 0 70.7%
Taylor expanded in a around 0 56.3%
+-commutative56.3%
*-commutative56.3%
Simplified56.3%
*-un-lft-identity56.3%
add-log-exp34.2%
sum-log34.6%
+-commutative34.6%
exp-sum34.6%
add-exp-log34.5%
pow-to-exp34.6%
Applied egg-rr34.6%
*-lft-identity34.6%
*-commutative34.6%
associate-*r*44.0%
Simplified44.0%
if 5.6e-105 < t < 2.00000000000000015e-43Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.1%
+-commutative99.1%
fma-define99.1%
sub-neg99.1%
metadata-eval99.1%
Simplified99.1%
Taylor expanded in x around 0 59.6%
associate--l+59.6%
remove-double-neg59.6%
log-rec59.6%
mul-1-neg59.6%
mul-1-neg59.6%
log-rec59.6%
remove-double-neg59.6%
sub-neg59.6%
metadata-eval59.6%
Simplified59.6%
Taylor expanded in t around 0 59.6%
*-un-lft-identity59.6%
add-log-exp40.7%
sum-log33.7%
sub-neg33.7%
metadata-eval33.7%
exp-sum33.8%
add-exp-log33.8%
pow-to-exp33.8%
Applied egg-rr33.8%
*-lft-identity33.8%
Simplified33.8%
if 1.6e21 < 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-undefine99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in t around inf 78.9%
Final simplification64.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (log t) a)))
(if (<= t 7.3e-208)
t_1
(if (<= t 1.8e-174)
(log (* z (* y (pow t -0.5))))
(if (<= t 6.4e-103)
t_1
(if (<= t 2e-43)
(log (* y (* z (pow t -0.5))))
(if (<= t 4.4e+25) t_1 (+ (log (+ x y)) (- (log z) t)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = log(t) * a;
double tmp;
if (t <= 7.3e-208) {
tmp = t_1;
} else if (t <= 1.8e-174) {
tmp = log((z * (y * pow(t, -0.5))));
} else if (t <= 6.4e-103) {
tmp = t_1;
} else if (t <= 2e-43) {
tmp = log((y * (z * pow(t, -0.5))));
} else if (t <= 4.4e+25) {
tmp = t_1;
} else {
tmp = log((x + y)) + (log(z) - 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) :: t_1
real(8) :: tmp
t_1 = log(t) * a
if (t <= 7.3d-208) then
tmp = t_1
else if (t <= 1.8d-174) then
tmp = log((z * (y * (t ** (-0.5d0)))))
else if (t <= 6.4d-103) then
tmp = t_1
else if (t <= 2d-43) then
tmp = log((y * (z * (t ** (-0.5d0)))))
else if (t <= 4.4d+25) then
tmp = t_1
else
tmp = log((x + y)) + (log(z) - t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.log(t) * a;
double tmp;
if (t <= 7.3e-208) {
tmp = t_1;
} else if (t <= 1.8e-174) {
tmp = Math.log((z * (y * Math.pow(t, -0.5))));
} else if (t <= 6.4e-103) {
tmp = t_1;
} else if (t <= 2e-43) {
tmp = Math.log((y * (z * Math.pow(t, -0.5))));
} else if (t <= 4.4e+25) {
tmp = t_1;
} else {
tmp = Math.log((x + y)) + (Math.log(z) - t);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.log(t) * a tmp = 0 if t <= 7.3e-208: tmp = t_1 elif t <= 1.8e-174: tmp = math.log((z * (y * math.pow(t, -0.5)))) elif t <= 6.4e-103: tmp = t_1 elif t <= 2e-43: tmp = math.log((y * (z * math.pow(t, -0.5)))) elif t <= 4.4e+25: tmp = t_1 else: tmp = math.log((x + y)) + (math.log(z) - t) return tmp
function code(x, y, z, t, a) t_1 = Float64(log(t) * a) tmp = 0.0 if (t <= 7.3e-208) tmp = t_1; elseif (t <= 1.8e-174) tmp = log(Float64(z * Float64(y * (t ^ -0.5)))); elseif (t <= 6.4e-103) tmp = t_1; elseif (t <= 2e-43) tmp = log(Float64(y * Float64(z * (t ^ -0.5)))); elseif (t <= 4.4e+25) tmp = t_1; else tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = log(t) * a; tmp = 0.0; if (t <= 7.3e-208) tmp = t_1; elseif (t <= 1.8e-174) tmp = log((z * (y * (t ^ -0.5)))); elseif (t <= 6.4e-103) tmp = t_1; elseif (t <= 2e-43) tmp = log((y * (z * (t ^ -0.5)))); elseif (t <= 4.4e+25) tmp = t_1; else tmp = log((x + y)) + (log(z) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]}, If[LessEqual[t, 7.3e-208], t$95$1, If[LessEqual[t, 1.8e-174], N[Log[N[(z * N[(y * N[Power[t, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 6.4e-103], t$95$1, If[LessEqual[t, 2e-43], N[Log[N[(y * N[(z * N[Power[t, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 4.4e+25], t$95$1, N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log t \cdot a\\
\mathbf{if}\;t \leq 7.3 \cdot 10^{-208}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 1.8 \cdot 10^{-174}:\\
\;\;\;\;\log \left(z \cdot \left(y \cdot {t}^{-0.5}\right)\right)\\
\mathbf{elif}\;t \leq 6.4 \cdot 10^{-103}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2 \cdot 10^{-43}:\\
\;\;\;\;\log \left(y \cdot \left(z \cdot {t}^{-0.5}\right)\right)\\
\mathbf{elif}\;t \leq 4.4 \cdot 10^{+25}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\end{array}
\end{array}
if t < 7.30000000000000002e-208 or 1.79999999999999999e-174 < t < 6.39999999999999953e-103 or 2.00000000000000015e-43 < t < 4.4000000000000001e25Initial program 99.4%
associate-+l-99.4%
associate--l+99.4%
sub-neg99.4%
+-commutative99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
fma-undefine99.4%
sub-neg99.4%
+-commutative99.4%
distribute-neg-in99.4%
metadata-eval99.4%
metadata-eval99.4%
unsub-neg99.4%
Simplified99.4%
Taylor expanded in a around inf 55.6%
*-commutative55.6%
Simplified55.6%
if 7.30000000000000002e-208 < t < 1.79999999999999999e-174Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.7%
+-commutative99.7%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around 0 70.7%
associate--l+70.7%
remove-double-neg70.7%
log-rec70.7%
mul-1-neg70.7%
mul-1-neg70.7%
log-rec70.7%
remove-double-neg70.7%
sub-neg70.7%
metadata-eval70.7%
Simplified70.7%
Taylor expanded in t around 0 70.7%
Taylor expanded in a around 0 56.3%
+-commutative56.3%
*-commutative56.3%
Simplified56.3%
*-un-lft-identity56.3%
add-log-exp34.2%
sum-log34.6%
+-commutative34.6%
exp-sum34.6%
add-exp-log34.5%
pow-to-exp34.6%
Applied egg-rr34.6%
*-lft-identity34.6%
*-commutative34.6%
associate-*r*44.0%
Simplified44.0%
if 6.39999999999999953e-103 < t < 2.00000000000000015e-43Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.1%
+-commutative99.1%
fma-define99.1%
sub-neg99.1%
metadata-eval99.1%
Simplified99.1%
Taylor expanded in x around 0 59.6%
associate--l+59.6%
remove-double-neg59.6%
log-rec59.6%
mul-1-neg59.6%
mul-1-neg59.6%
log-rec59.6%
remove-double-neg59.6%
sub-neg59.6%
metadata-eval59.6%
Simplified59.6%
Taylor expanded in t around 0 59.6%
Taylor expanded in a around 0 40.1%
+-commutative40.1%
*-commutative40.1%
Simplified40.1%
*-un-lft-identity40.1%
add-log-exp40.1%
sum-log33.1%
+-commutative33.1%
exp-sum33.2%
add-exp-log33.2%
pow-to-exp33.2%
Applied egg-rr33.2%
*-lft-identity33.2%
Simplified33.2%
if 4.4000000000000001e25 < 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-undefine99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in t around inf 78.9%
Final simplification64.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (log t) a)))
(if (<= t 7.5e-208)
t_1
(if (<= t 6.2e-175)
(log (* z (* y (pow t -0.5))))
(if (<= t 7.2e-106)
t_1
(if (<= t 2.5e-43)
(log (* y (* z (pow t -0.5))))
(if (<= t 4.8e+22) t_1 (- t))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = log(t) * a;
double tmp;
if (t <= 7.5e-208) {
tmp = t_1;
} else if (t <= 6.2e-175) {
tmp = log((z * (y * pow(t, -0.5))));
} else if (t <= 7.2e-106) {
tmp = t_1;
} else if (t <= 2.5e-43) {
tmp = log((y * (z * pow(t, -0.5))));
} else if (t <= 4.8e+22) {
tmp = t_1;
} 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) :: t_1
real(8) :: tmp
t_1 = log(t) * a
if (t <= 7.5d-208) then
tmp = t_1
else if (t <= 6.2d-175) then
tmp = log((z * (y * (t ** (-0.5d0)))))
else if (t <= 7.2d-106) then
tmp = t_1
else if (t <= 2.5d-43) then
tmp = log((y * (z * (t ** (-0.5d0)))))
else if (t <= 4.8d+22) then
tmp = t_1
else
tmp = -t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.log(t) * a;
double tmp;
if (t <= 7.5e-208) {
tmp = t_1;
} else if (t <= 6.2e-175) {
tmp = Math.log((z * (y * Math.pow(t, -0.5))));
} else if (t <= 7.2e-106) {
tmp = t_1;
} else if (t <= 2.5e-43) {
tmp = Math.log((y * (z * Math.pow(t, -0.5))));
} else if (t <= 4.8e+22) {
tmp = t_1;
} else {
tmp = -t;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.log(t) * a tmp = 0 if t <= 7.5e-208: tmp = t_1 elif t <= 6.2e-175: tmp = math.log((z * (y * math.pow(t, -0.5)))) elif t <= 7.2e-106: tmp = t_1 elif t <= 2.5e-43: tmp = math.log((y * (z * math.pow(t, -0.5)))) elif t <= 4.8e+22: tmp = t_1 else: tmp = -t return tmp
function code(x, y, z, t, a) t_1 = Float64(log(t) * a) tmp = 0.0 if (t <= 7.5e-208) tmp = t_1; elseif (t <= 6.2e-175) tmp = log(Float64(z * Float64(y * (t ^ -0.5)))); elseif (t <= 7.2e-106) tmp = t_1; elseif (t <= 2.5e-43) tmp = log(Float64(y * Float64(z * (t ^ -0.5)))); elseif (t <= 4.8e+22) tmp = t_1; else tmp = Float64(-t); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = log(t) * a; tmp = 0.0; if (t <= 7.5e-208) tmp = t_1; elseif (t <= 6.2e-175) tmp = log((z * (y * (t ^ -0.5)))); elseif (t <= 7.2e-106) tmp = t_1; elseif (t <= 2.5e-43) tmp = log((y * (z * (t ^ -0.5)))); elseif (t <= 4.8e+22) tmp = t_1; else tmp = -t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]}, If[LessEqual[t, 7.5e-208], t$95$1, If[LessEqual[t, 6.2e-175], N[Log[N[(z * N[(y * N[Power[t, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 7.2e-106], t$95$1, If[LessEqual[t, 2.5e-43], N[Log[N[(y * N[(z * N[Power[t, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[t, 4.8e+22], t$95$1, (-t)]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log t \cdot a\\
\mathbf{if}\;t \leq 7.5 \cdot 10^{-208}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 6.2 \cdot 10^{-175}:\\
\;\;\;\;\log \left(z \cdot \left(y \cdot {t}^{-0.5}\right)\right)\\
\mathbf{elif}\;t \leq 7.2 \cdot 10^{-106}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.5 \cdot 10^{-43}:\\
\;\;\;\;\log \left(y \cdot \left(z \cdot {t}^{-0.5}\right)\right)\\
\mathbf{elif}\;t \leq 4.8 \cdot 10^{+22}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;-t\\
\end{array}
\end{array}
if t < 7.4999999999999999e-208 or 6.19999999999999997e-175 < t < 7.20000000000000025e-106 or 2.50000000000000009e-43 < t < 4.8e22Initial program 99.4%
associate-+l-99.4%
associate--l+99.4%
sub-neg99.4%
+-commutative99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
fma-undefine99.4%
sub-neg99.4%
+-commutative99.4%
distribute-neg-in99.4%
metadata-eval99.4%
metadata-eval99.4%
unsub-neg99.4%
Simplified99.4%
Taylor expanded in a around inf 55.6%
*-commutative55.6%
Simplified55.6%
if 7.4999999999999999e-208 < t < 6.19999999999999997e-175Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.7%
+-commutative99.7%
fma-define99.7%
sub-neg99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around 0 70.7%
associate--l+70.7%
remove-double-neg70.7%
log-rec70.7%
mul-1-neg70.7%
mul-1-neg70.7%
log-rec70.7%
remove-double-neg70.7%
sub-neg70.7%
metadata-eval70.7%
Simplified70.7%
Taylor expanded in t around 0 70.7%
Taylor expanded in a around 0 56.3%
+-commutative56.3%
*-commutative56.3%
Simplified56.3%
*-un-lft-identity56.3%
add-log-exp34.2%
sum-log34.6%
+-commutative34.6%
exp-sum34.6%
add-exp-log34.5%
pow-to-exp34.6%
Applied egg-rr34.6%
*-lft-identity34.6%
*-commutative34.6%
associate-*r*44.0%
Simplified44.0%
if 7.20000000000000025e-106 < t < 2.50000000000000009e-43Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.1%
+-commutative99.1%
fma-define99.1%
sub-neg99.1%
metadata-eval99.1%
Simplified99.1%
Taylor expanded in x around 0 59.6%
associate--l+59.6%
remove-double-neg59.6%
log-rec59.6%
mul-1-neg59.6%
mul-1-neg59.6%
log-rec59.6%
remove-double-neg59.6%
sub-neg59.6%
metadata-eval59.6%
Simplified59.6%
Taylor expanded in t around 0 59.6%
Taylor expanded in a around 0 40.1%
+-commutative40.1%
*-commutative40.1%
Simplified40.1%
*-un-lft-identity40.1%
add-log-exp40.1%
sum-log33.1%
+-commutative33.1%
exp-sum33.2%
add-exp-log33.2%
pow-to-exp33.2%
Applied egg-rr33.2%
*-lft-identity33.2%
Simplified33.2%
if 4.8e22 < 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-undefine99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in t around inf 78.9%
neg-mul-178.9%
Simplified78.9%
Final simplification64.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (* (log t) a)) (t_2 (log (* y (* z (pow t -0.5))))))
(if (<= t 2.5e-207)
t_1
(if (<= t 5.8e-175)
t_2
(if (<= t 9.2e-103)
t_1
(if (<= t 9.5e-43) t_2 (if (<= t 4.2e+21) t_1 (- t))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = log(t) * a;
double t_2 = log((y * (z * pow(t, -0.5))));
double tmp;
if (t <= 2.5e-207) {
tmp = t_1;
} else if (t <= 5.8e-175) {
tmp = t_2;
} else if (t <= 9.2e-103) {
tmp = t_1;
} else if (t <= 9.5e-43) {
tmp = t_2;
} else if (t <= 4.2e+21) {
tmp = t_1;
} 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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = log(t) * a
t_2 = log((y * (z * (t ** (-0.5d0)))))
if (t <= 2.5d-207) then
tmp = t_1
else if (t <= 5.8d-175) then
tmp = t_2
else if (t <= 9.2d-103) then
tmp = t_1
else if (t <= 9.5d-43) then
tmp = t_2
else if (t <= 4.2d+21) then
tmp = t_1
else
tmp = -t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = Math.log(t) * a;
double t_2 = Math.log((y * (z * Math.pow(t, -0.5))));
double tmp;
if (t <= 2.5e-207) {
tmp = t_1;
} else if (t <= 5.8e-175) {
tmp = t_2;
} else if (t <= 9.2e-103) {
tmp = t_1;
} else if (t <= 9.5e-43) {
tmp = t_2;
} else if (t <= 4.2e+21) {
tmp = t_1;
} else {
tmp = -t;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = math.log(t) * a t_2 = math.log((y * (z * math.pow(t, -0.5)))) tmp = 0 if t <= 2.5e-207: tmp = t_1 elif t <= 5.8e-175: tmp = t_2 elif t <= 9.2e-103: tmp = t_1 elif t <= 9.5e-43: tmp = t_2 elif t <= 4.2e+21: tmp = t_1 else: tmp = -t return tmp
function code(x, y, z, t, a) t_1 = Float64(log(t) * a) t_2 = log(Float64(y * Float64(z * (t ^ -0.5)))) tmp = 0.0 if (t <= 2.5e-207) tmp = t_1; elseif (t <= 5.8e-175) tmp = t_2; elseif (t <= 9.2e-103) tmp = t_1; elseif (t <= 9.5e-43) tmp = t_2; elseif (t <= 4.2e+21) tmp = t_1; else tmp = Float64(-t); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = log(t) * a; t_2 = log((y * (z * (t ^ -0.5)))); tmp = 0.0; if (t <= 2.5e-207) tmp = t_1; elseif (t <= 5.8e-175) tmp = t_2; elseif (t <= 9.2e-103) tmp = t_1; elseif (t <= 9.5e-43) tmp = t_2; elseif (t <= 4.2e+21) tmp = t_1; else tmp = -t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision]}, Block[{t$95$2 = N[Log[N[(y * N[(z * N[Power[t, -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t, 2.5e-207], t$95$1, If[LessEqual[t, 5.8e-175], t$95$2, If[LessEqual[t, 9.2e-103], t$95$1, If[LessEqual[t, 9.5e-43], t$95$2, If[LessEqual[t, 4.2e+21], t$95$1, (-t)]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log t \cdot a\\
t_2 := \log \left(y \cdot \left(z \cdot {t}^{-0.5}\right)\right)\\
\mathbf{if}\;t \leq 2.5 \cdot 10^{-207}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 5.8 \cdot 10^{-175}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq 9.2 \cdot 10^{-103}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 9.5 \cdot 10^{-43}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq 4.2 \cdot 10^{+21}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;-t\\
\end{array}
\end{array}
if t < 2.50000000000000007e-207 or 5.79999999999999998e-175 < t < 9.2000000000000003e-103 or 9.50000000000000044e-43 < t < 4.2e21Initial program 99.4%
associate-+l-99.4%
associate--l+99.4%
sub-neg99.4%
+-commutative99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
fma-undefine99.4%
sub-neg99.4%
+-commutative99.4%
distribute-neg-in99.4%
metadata-eval99.4%
metadata-eval99.4%
unsub-neg99.4%
Simplified99.4%
Taylor expanded in a around inf 55.1%
*-commutative55.1%
Simplified55.1%
if 2.50000000000000007e-207 < t < 5.79999999999999998e-175 or 9.2000000000000003e-103 < t < 9.50000000000000044e-43Initial program 99.3%
associate--l+99.3%
+-commutative99.3%
associate-+l+99.2%
+-commutative99.2%
fma-define99.2%
sub-neg99.2%
metadata-eval99.2%
Simplified99.2%
Taylor expanded in x around 0 61.7%
associate--l+61.7%
remove-double-neg61.7%
log-rec61.7%
mul-1-neg61.7%
mul-1-neg61.7%
log-rec61.7%
remove-double-neg61.7%
sub-neg61.7%
metadata-eval61.7%
Simplified61.7%
Taylor expanded in t around 0 61.7%
Taylor expanded in a around 0 43.2%
+-commutative43.2%
*-commutative43.2%
Simplified43.2%
*-un-lft-identity43.2%
add-log-exp39.5%
sum-log34.4%
+-commutative34.4%
exp-sum34.5%
add-exp-log34.5%
pow-to-exp34.5%
Applied egg-rr34.5%
*-lft-identity34.5%
Simplified34.5%
if 4.2e21 < 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-undefine99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in t around inf 78.9%
neg-mul-178.9%
Simplified78.9%
(FPCore (x y z t a)
:precision binary64
(if (<= t 5.2e+96)
(- (+ (log (* z (+ x y))) (* (log t) (- a 0.5))) t)
(if (<= t 1.85e+167)
(- t)
(if (<= t 1.12e+229)
(- (- (log (* y z)) (* (log t) (- 0.5 a))) t)
(+ (log (+ x y)) (- (log z) t))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 5.2e+96) {
tmp = (log((z * (x + y))) + (log(t) * (a - 0.5))) - t;
} else if (t <= 1.85e+167) {
tmp = -t;
} else if (t <= 1.12e+229) {
tmp = (log((y * z)) - (log(t) * (0.5 - a))) - t;
} else {
tmp = log((x + y)) + (log(z) - 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 <= 5.2d+96) then
tmp = (log((z * (x + y))) + (log(t) * (a - 0.5d0))) - t
else if (t <= 1.85d+167) then
tmp = -t
else if (t <= 1.12d+229) then
tmp = (log((y * z)) - (log(t) * (0.5d0 - a))) - t
else
tmp = log((x + y)) + (log(z) - 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 <= 5.2e+96) {
tmp = (Math.log((z * (x + y))) + (Math.log(t) * (a - 0.5))) - t;
} else if (t <= 1.85e+167) {
tmp = -t;
} else if (t <= 1.12e+229) {
tmp = (Math.log((y * z)) - (Math.log(t) * (0.5 - a))) - t;
} else {
tmp = Math.log((x + y)) + (Math.log(z) - t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 5.2e+96: tmp = (math.log((z * (x + y))) + (math.log(t) * (a - 0.5))) - t elif t <= 1.85e+167: tmp = -t elif t <= 1.12e+229: tmp = (math.log((y * z)) - (math.log(t) * (0.5 - a))) - t else: tmp = math.log((x + y)) + (math.log(z) - t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 5.2e+96) tmp = Float64(Float64(log(Float64(z * Float64(x + y))) + Float64(log(t) * Float64(a - 0.5))) - t); elseif (t <= 1.85e+167) tmp = Float64(-t); elseif (t <= 1.12e+229) tmp = Float64(Float64(log(Float64(y * z)) - Float64(log(t) * Float64(0.5 - a))) - t); else tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 5.2e+96) tmp = (log((z * (x + y))) + (log(t) * (a - 0.5))) - t; elseif (t <= 1.85e+167) tmp = -t; elseif (t <= 1.12e+229) tmp = (log((y * z)) - (log(t) * (0.5 - a))) - t; else tmp = log((x + y)) + (log(z) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 5.2e+96], N[(N[(N[Log[N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[(N[Log[t], $MachinePrecision] * N[(a - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], If[LessEqual[t, 1.85e+167], (-t), If[LessEqual[t, 1.12e+229], N[(N[(N[Log[N[(y * z), $MachinePrecision]], $MachinePrecision] - N[(N[Log[t], $MachinePrecision] * N[(0.5 - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 5.2 \cdot 10^{+96}:\\
\;\;\;\;\left(\log \left(z \cdot \left(x + y\right)\right) + \log t \cdot \left(a - 0.5\right)\right) - t\\
\mathbf{elif}\;t \leq 1.85 \cdot 10^{+167}:\\
\;\;\;\;-t\\
\mathbf{elif}\;t \leq 1.12 \cdot 10^{+229}:\\
\;\;\;\;\left(\log \left(y \cdot z\right) - \log t \cdot \left(0.5 - a\right)\right) - t\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\end{array}
\end{array}
if t < 5.2e96Initial program 99.4%
associate-+l-99.4%
associate--l+99.4%
sub-neg99.4%
+-commutative99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
fma-undefine99.4%
sub-neg99.4%
+-commutative99.4%
distribute-neg-in99.4%
metadata-eval99.4%
metadata-eval99.4%
unsub-neg99.4%
Simplified99.4%
associate-+r-99.4%
fma-undefine99.4%
associate--r+99.4%
sum-log75.5%
Applied egg-rr75.5%
if 5.2e96 < t < 1.85e167Initial program 100.0%
associate-+l-100.0%
associate--l+100.0%
sub-neg100.0%
+-commutative100.0%
*-commutative100.0%
distribute-rgt-neg-in100.0%
fma-undefine100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
Taylor expanded in t around inf 85.3%
neg-mul-185.3%
Simplified85.3%
if 1.85e167 < t < 1.12e229Initial program 99.8%
associate--l+99.8%
+-commutative99.8%
associate-+l+99.8%
+-commutative99.8%
fma-define99.8%
sub-neg99.8%
metadata-eval99.8%
Simplified99.8%
add-sqr-sqrt26.8%
pow226.8%
Applied egg-rr26.8%
Taylor expanded in x around 0 73.7%
if 1.12e229 < t Initial program 100.0%
associate-+l-100.0%
associate--l+100.0%
sub-neg100.0%
+-commutative100.0%
*-commutative100.0%
distribute-rgt-neg-in100.0%
fma-undefine100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
Taylor expanded in t around inf 96.8%
Final simplification80.2%
(FPCore (x y z t a) :precision binary64 (if (<= t 3.1e+95) (- (- (log (* y z)) (* (log t) (- 0.5 a))) t) (+ (log (+ x y)) (- (log z) t))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 3.1e+95) {
tmp = (log((y * z)) - (log(t) * (0.5 - a))) - t;
} else {
tmp = log((x + y)) + (log(z) - 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 <= 3.1d+95) then
tmp = (log((y * z)) - (log(t) * (0.5d0 - a))) - t
else
tmp = log((x + y)) + (log(z) - 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 <= 3.1e+95) {
tmp = (Math.log((y * z)) - (Math.log(t) * (0.5 - a))) - t;
} else {
tmp = Math.log((x + y)) + (Math.log(z) - t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 3.1e+95: tmp = (math.log((y * z)) - (math.log(t) * (0.5 - a))) - t else: tmp = math.log((x + y)) + (math.log(z) - t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 3.1e+95) tmp = Float64(Float64(log(Float64(y * z)) - Float64(log(t) * Float64(0.5 - a))) - t); else tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 3.1e+95) tmp = (log((y * z)) - (log(t) * (0.5 - a))) - t; else tmp = log((x + y)) + (log(z) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 3.1e+95], N[(N[(N[Log[N[(y * z), $MachinePrecision]], $MachinePrecision] - N[(N[Log[t], $MachinePrecision] * N[(0.5 - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 3.1 \cdot 10^{+95}:\\
\;\;\;\;\left(\log \left(y \cdot z\right) - \log t \cdot \left(0.5 - a\right)\right) - t\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\end{array}
\end{array}
if t < 3.1000000000000003e95Initial program 99.4%
associate--l+99.4%
+-commutative99.4%
associate-+l+99.4%
+-commutative99.4%
fma-define99.4%
sub-neg99.4%
metadata-eval99.4%
Simplified99.4%
add-sqr-sqrt55.0%
pow255.0%
Applied egg-rr40.6%
Taylor expanded in x around 0 49.9%
if 3.1000000000000003e95 < 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-undefine100.0%
sub-neg100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
Taylor expanded in t around inf 83.6%
Final simplification63.5%
(FPCore (x y z t a) :precision binary64 (if (<= t 3.5e+20) (+ (* (log t) (+ a -0.5)) (log (* z (+ x y)))) (+ (log (+ x y)) (- (log z) t))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 3.5e+20) {
tmp = (log(t) * (a + -0.5)) + log((z * (x + y)));
} else {
tmp = log((x + y)) + (log(z) - 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 <= 3.5d+20) then
tmp = (log(t) * (a + (-0.5d0))) + log((z * (x + y)))
else
tmp = log((x + y)) + (log(z) - 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 <= 3.5e+20) {
tmp = (Math.log(t) * (a + -0.5)) + Math.log((z * (x + y)));
} else {
tmp = Math.log((x + y)) + (Math.log(z) - t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 3.5e+20: tmp = (math.log(t) * (a + -0.5)) + math.log((z * (x + y))) else: tmp = math.log((x + y)) + (math.log(z) - t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 3.5e+20) tmp = Float64(Float64(log(t) * Float64(a + -0.5)) + log(Float64(z * Float64(x + y)))); else tmp = Float64(log(Float64(x + y)) + Float64(log(z) - t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 3.5e+20) tmp = (log(t) * (a + -0.5)) + log((z * (x + y))); else tmp = log((x + y)) + (log(z) - t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 3.5e+20], N[(N[(N[Log[t], $MachinePrecision] * N[(a + -0.5), $MachinePrecision]), $MachinePrecision] + N[Log[N[(z * N[(x + y), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[Log[N[(x + y), $MachinePrecision]], $MachinePrecision] + N[(N[Log[z], $MachinePrecision] - t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 3.5 \cdot 10^{+20}:\\
\;\;\;\;\log t \cdot \left(a + -0.5\right) + \log \left(z \cdot \left(x + y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + y\right) + \left(\log z - t\right)\\
\end{array}
\end{array}
if t < 3.5e20Initial program 99.4%
associate--l+99.4%
+-commutative99.4%
associate-+l+99.4%
+-commutative99.4%
fma-define99.3%
sub-neg99.3%
metadata-eval99.3%
Simplified99.3%
add-sqr-sqrt61.2%
pow261.2%
Applied egg-rr45.2%
Taylor expanded in t around 0 70.7%
+-commutative70.7%
sub-neg70.7%
metadata-eval70.7%
+-commutative70.7%
Simplified70.7%
if 3.5e20 < 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-undefine99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in t around inf 78.9%
Final simplification74.7%
(FPCore (x y z t a) :precision binary64 (if (<= t 5.2e+22) (* (log t) a) (- t)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 5.2e+22) {
tmp = log(t) * a;
} 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 <= 5.2d+22) then
tmp = log(t) * a
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 <= 5.2e+22) {
tmp = Math.log(t) * a;
} else {
tmp = -t;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 5.2e+22: tmp = math.log(t) * a else: tmp = -t return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 5.2e+22) tmp = Float64(log(t) * a); else tmp = Float64(-t); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 5.2e+22) tmp = log(t) * a; else tmp = -t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 5.2e+22], N[(N[Log[t], $MachinePrecision] * a), $MachinePrecision], (-t)]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 5.2 \cdot 10^{+22}:\\
\;\;\;\;\log t \cdot a\\
\mathbf{else}:\\
\;\;\;\;-t\\
\end{array}
\end{array}
if t < 5.2e22Initial program 99.4%
associate-+l-99.4%
associate--l+99.4%
sub-neg99.4%
+-commutative99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
fma-undefine99.4%
sub-neg99.4%
+-commutative99.4%
distribute-neg-in99.4%
metadata-eval99.4%
metadata-eval99.4%
unsub-neg99.4%
Simplified99.4%
Taylor expanded in a around inf 46.8%
*-commutative46.8%
Simplified46.8%
if 5.2e22 < 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-undefine99.9%
sub-neg99.9%
+-commutative99.9%
distribute-neg-in99.9%
metadata-eval99.9%
metadata-eval99.9%
unsub-neg99.9%
Simplified99.9%
Taylor expanded in t around inf 78.9%
neg-mul-178.9%
Simplified78.9%
(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.6%
associate-+l-99.6%
associate--l+99.6%
sub-neg99.6%
+-commutative99.6%
*-commutative99.6%
distribute-rgt-neg-in99.6%
fma-undefine99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
metadata-eval99.6%
metadata-eval99.6%
unsub-neg99.6%
Simplified99.6%
Taylor expanded in t around inf 40.7%
neg-mul-140.7%
Simplified40.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(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 2024096
(FPCore (x y z t a)
:name "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"
:precision binary64
:alt
(+ (log (+ x y)) (+ (- (log z) t) (* (- a 0.5) (log t))))
(+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))