
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (* i (+ (+ alpha beta) i)))
(t_1 (+ (+ alpha beta) (* 2.0 i)))
(t_2 (* t_1 t_1)))
(/ (/ (* t_0 (+ (* beta alpha) t_0)) t_2) (- t_2 1.0))))
double code(double alpha, double beta, double i) {
double t_0 = i * ((alpha + beta) + i);
double t_1 = (alpha + beta) + (2.0 * i);
double t_2 = t_1 * t_1;
return ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0);
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
t_0 = i * ((alpha + beta) + i)
t_1 = (alpha + beta) + (2.0d0 * i)
t_2 = t_1 * t_1
code = ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0d0)
end function
public static double code(double alpha, double beta, double i) {
double t_0 = i * ((alpha + beta) + i);
double t_1 = (alpha + beta) + (2.0 * i);
double t_2 = t_1 * t_1;
return ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0);
}
def code(alpha, beta, i): t_0 = i * ((alpha + beta) + i) t_1 = (alpha + beta) + (2.0 * i) t_2 = t_1 * t_1 return ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0)
function code(alpha, beta, i) t_0 = Float64(i * Float64(Float64(alpha + beta) + i)) t_1 = Float64(Float64(alpha + beta) + Float64(2.0 * i)) t_2 = Float64(t_1 * t_1) return Float64(Float64(Float64(t_0 * Float64(Float64(beta * alpha) + t_0)) / t_2) / Float64(t_2 - 1.0)) end
function tmp = code(alpha, beta, i) t_0 = i * ((alpha + beta) + i); t_1 = (alpha + beta) + (2.0 * i); t_2 = t_1 * t_1; tmp = ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0); end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(i * N[(N[(alpha + beta), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * t$95$1), $MachinePrecision]}, N[(N[(N[(t$95$0 * N[(N[(beta * alpha), $MachinePrecision] + t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision] / N[(t$95$2 - 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot \left(\left(\alpha + \beta\right) + i\right)\\
t_1 := \left(\alpha + \beta\right) + 2 \cdot i\\
t_2 := t\_1 \cdot t\_1\\
\frac{\frac{t\_0 \cdot \left(\beta \cdot \alpha + t\_0\right)}{t\_2}}{t\_2 - 1}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (* i (+ (+ alpha beta) i)))
(t_1 (+ (+ alpha beta) (* 2.0 i)))
(t_2 (* t_1 t_1)))
(/ (/ (* t_0 (+ (* beta alpha) t_0)) t_2) (- t_2 1.0))))
double code(double alpha, double beta, double i) {
double t_0 = i * ((alpha + beta) + i);
double t_1 = (alpha + beta) + (2.0 * i);
double t_2 = t_1 * t_1;
return ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0);
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
t_0 = i * ((alpha + beta) + i)
t_1 = (alpha + beta) + (2.0d0 * i)
t_2 = t_1 * t_1
code = ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0d0)
end function
public static double code(double alpha, double beta, double i) {
double t_0 = i * ((alpha + beta) + i);
double t_1 = (alpha + beta) + (2.0 * i);
double t_2 = t_1 * t_1;
return ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0);
}
def code(alpha, beta, i): t_0 = i * ((alpha + beta) + i) t_1 = (alpha + beta) + (2.0 * i) t_2 = t_1 * t_1 return ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0)
function code(alpha, beta, i) t_0 = Float64(i * Float64(Float64(alpha + beta) + i)) t_1 = Float64(Float64(alpha + beta) + Float64(2.0 * i)) t_2 = Float64(t_1 * t_1) return Float64(Float64(Float64(t_0 * Float64(Float64(beta * alpha) + t_0)) / t_2) / Float64(t_2 - 1.0)) end
function tmp = code(alpha, beta, i) t_0 = i * ((alpha + beta) + i); t_1 = (alpha + beta) + (2.0 * i); t_2 = t_1 * t_1; tmp = ((t_0 * ((beta * alpha) + t_0)) / t_2) / (t_2 - 1.0); end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(i * N[(N[(alpha + beta), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * t$95$1), $MachinePrecision]}, N[(N[(N[(t$95$0 * N[(N[(beta * alpha), $MachinePrecision] + t$95$0), $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision] / N[(t$95$2 - 1.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot \left(\left(\alpha + \beta\right) + i\right)\\
t_1 := \left(\alpha + \beta\right) + 2 \cdot i\\
t_2 := t\_1 \cdot t\_1\\
\frac{\frac{t\_0 \cdot \left(\beta \cdot \alpha + t\_0\right)}{t\_2}}{t\_2 - 1}
\end{array}
\end{array}
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ beta (* i 2.0)))
(t_1 (+ alpha t_0))
(t_2 (* i (+ i (+ alpha beta)))))
(if (<= i 3.1e+144)
(/
(/ (/ (+ (* alpha beta) t_2) (/ t_1 (/ t_2 t_1))) (+ alpha (+ t_0 1.0)))
(+ alpha (+ t_0 -1.0)))
0.0625)))
double code(double alpha, double beta, double i) {
double t_0 = beta + (i * 2.0);
double t_1 = alpha + t_0;
double t_2 = i * (i + (alpha + beta));
double tmp;
if (i <= 3.1e+144) {
tmp = ((((alpha * beta) + t_2) / (t_1 / (t_2 / t_1))) / (alpha + (t_0 + 1.0))) / (alpha + (t_0 + -1.0));
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = beta + (i * 2.0d0)
t_1 = alpha + t_0
t_2 = i * (i + (alpha + beta))
if (i <= 3.1d+144) then
tmp = ((((alpha * beta) + t_2) / (t_1 / (t_2 / t_1))) / (alpha + (t_0 + 1.0d0))) / (alpha + (t_0 + (-1.0d0)))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = beta + (i * 2.0);
double t_1 = alpha + t_0;
double t_2 = i * (i + (alpha + beta));
double tmp;
if (i <= 3.1e+144) {
tmp = ((((alpha * beta) + t_2) / (t_1 / (t_2 / t_1))) / (alpha + (t_0 + 1.0))) / (alpha + (t_0 + -1.0));
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): t_0 = beta + (i * 2.0) t_1 = alpha + t_0 t_2 = i * (i + (alpha + beta)) tmp = 0 if i <= 3.1e+144: tmp = ((((alpha * beta) + t_2) / (t_1 / (t_2 / t_1))) / (alpha + (t_0 + 1.0))) / (alpha + (t_0 + -1.0)) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) t_0 = Float64(beta + Float64(i * 2.0)) t_1 = Float64(alpha + t_0) t_2 = Float64(i * Float64(i + Float64(alpha + beta))) tmp = 0.0 if (i <= 3.1e+144) tmp = Float64(Float64(Float64(Float64(Float64(alpha * beta) + t_2) / Float64(t_1 / Float64(t_2 / t_1))) / Float64(alpha + Float64(t_0 + 1.0))) / Float64(alpha + Float64(t_0 + -1.0))); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = beta + (i * 2.0); t_1 = alpha + t_0; t_2 = i * (i + (alpha + beta)); tmp = 0.0; if (i <= 3.1e+144) tmp = ((((alpha * beta) + t_2) / (t_1 / (t_2 / t_1))) / (alpha + (t_0 + 1.0))) / (alpha + (t_0 + -1.0)); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(alpha + t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(i * N[(i + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[i, 3.1e+144], N[(N[(N[(N[(N[(alpha * beta), $MachinePrecision] + t$95$2), $MachinePrecision] / N[(t$95$1 / N[(t$95$2 / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(alpha + N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(alpha + N[(t$95$0 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0625]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \beta + i \cdot 2\\
t_1 := \alpha + t\_0\\
t_2 := i \cdot \left(i + \left(\alpha + \beta\right)\right)\\
\mathbf{if}\;i \leq 3.1 \cdot 10^{+144}:\\
\;\;\;\;\frac{\frac{\frac{\alpha \cdot \beta + t\_2}{\frac{t\_1}{\frac{t\_2}{t\_1}}}}{\alpha + \left(t\_0 + 1\right)}}{\alpha + \left(t\_0 + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
\end{array}
if i < 3.1000000000000002e144Initial program 41.9%
/-lowering-/.f64N/A
Simplified83.2%
Applied egg-rr88.4%
if 3.1000000000000002e144 < i Initial program 0.1%
/-lowering-/.f64N/A
Simplified1.7%
Taylor expanded in i around inf
Simplified89.0%
Final simplification88.7%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ beta (* i 2.0)))
(t_1 (* i (+ i (+ alpha beta))))
(t_2 (+ alpha t_0)))
(if (<= i 6.4e+144)
(*
(/ (/ t_1 t_2) (+ alpha (+ t_0 1.0)))
(/ (/ (+ (* alpha beta) t_1) t_2) (+ alpha (+ t_0 -1.0))))
0.0625)))
double code(double alpha, double beta, double i) {
double t_0 = beta + (i * 2.0);
double t_1 = i * (i + (alpha + beta));
double t_2 = alpha + t_0;
double tmp;
if (i <= 6.4e+144) {
tmp = ((t_1 / t_2) / (alpha + (t_0 + 1.0))) * ((((alpha * beta) + t_1) / t_2) / (alpha + (t_0 + -1.0)));
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = beta + (i * 2.0d0)
t_1 = i * (i + (alpha + beta))
t_2 = alpha + t_0
if (i <= 6.4d+144) then
tmp = ((t_1 / t_2) / (alpha + (t_0 + 1.0d0))) * ((((alpha * beta) + t_1) / t_2) / (alpha + (t_0 + (-1.0d0))))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = beta + (i * 2.0);
double t_1 = i * (i + (alpha + beta));
double t_2 = alpha + t_0;
double tmp;
if (i <= 6.4e+144) {
tmp = ((t_1 / t_2) / (alpha + (t_0 + 1.0))) * ((((alpha * beta) + t_1) / t_2) / (alpha + (t_0 + -1.0)));
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): t_0 = beta + (i * 2.0) t_1 = i * (i + (alpha + beta)) t_2 = alpha + t_0 tmp = 0 if i <= 6.4e+144: tmp = ((t_1 / t_2) / (alpha + (t_0 + 1.0))) * ((((alpha * beta) + t_1) / t_2) / (alpha + (t_0 + -1.0))) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) t_0 = Float64(beta + Float64(i * 2.0)) t_1 = Float64(i * Float64(i + Float64(alpha + beta))) t_2 = Float64(alpha + t_0) tmp = 0.0 if (i <= 6.4e+144) tmp = Float64(Float64(Float64(t_1 / t_2) / Float64(alpha + Float64(t_0 + 1.0))) * Float64(Float64(Float64(Float64(alpha * beta) + t_1) / t_2) / Float64(alpha + Float64(t_0 + -1.0)))); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = beta + (i * 2.0); t_1 = i * (i + (alpha + beta)); t_2 = alpha + t_0; tmp = 0.0; if (i <= 6.4e+144) tmp = ((t_1 / t_2) / (alpha + (t_0 + 1.0))) * ((((alpha * beta) + t_1) / t_2) / (alpha + (t_0 + -1.0))); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(i * N[(i + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(alpha + t$95$0), $MachinePrecision]}, If[LessEqual[i, 6.4e+144], N[(N[(N[(t$95$1 / t$95$2), $MachinePrecision] / N[(alpha + N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(alpha * beta), $MachinePrecision] + t$95$1), $MachinePrecision] / t$95$2), $MachinePrecision] / N[(alpha + N[(t$95$0 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0625]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \beta + i \cdot 2\\
t_1 := i \cdot \left(i + \left(\alpha + \beta\right)\right)\\
t_2 := \alpha + t\_0\\
\mathbf{if}\;i \leq 6.4 \cdot 10^{+144}:\\
\;\;\;\;\frac{\frac{t\_1}{t\_2}}{\alpha + \left(t\_0 + 1\right)} \cdot \frac{\frac{\alpha \cdot \beta + t\_1}{t\_2}}{\alpha + \left(t\_0 + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
\end{array}
if i < 6.4000000000000002e144Initial program 41.9%
/-lowering-/.f64N/A
Simplified83.2%
Applied egg-rr88.3%
if 6.4000000000000002e144 < i Initial program 0.1%
/-lowering-/.f64N/A
Simplified1.7%
Taylor expanded in i around inf
Simplified89.0%
Final simplification88.6%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ beta (* i 2.0))) (t_1 (+ alpha t_0)))
(if (<= i 4.9e+144)
(/
(*
(+ (* alpha beta) (* i (+ i (+ alpha beta))))
(/ (* i (+ i beta)) (* t_0 t_0)))
(+ (* t_1 t_1) -1.0))
0.0625)))
double code(double alpha, double beta, double i) {
double t_0 = beta + (i * 2.0);
double t_1 = alpha + t_0;
double tmp;
if (i <= 4.9e+144) {
tmp = (((alpha * beta) + (i * (i + (alpha + beta)))) * ((i * (i + beta)) / (t_0 * t_0))) / ((t_1 * t_1) + -1.0);
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = beta + (i * 2.0d0)
t_1 = alpha + t_0
if (i <= 4.9d+144) then
tmp = (((alpha * beta) + (i * (i + (alpha + beta)))) * ((i * (i + beta)) / (t_0 * t_0))) / ((t_1 * t_1) + (-1.0d0))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = beta + (i * 2.0);
double t_1 = alpha + t_0;
double tmp;
if (i <= 4.9e+144) {
tmp = (((alpha * beta) + (i * (i + (alpha + beta)))) * ((i * (i + beta)) / (t_0 * t_0))) / ((t_1 * t_1) + -1.0);
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): t_0 = beta + (i * 2.0) t_1 = alpha + t_0 tmp = 0 if i <= 4.9e+144: tmp = (((alpha * beta) + (i * (i + (alpha + beta)))) * ((i * (i + beta)) / (t_0 * t_0))) / ((t_1 * t_1) + -1.0) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) t_0 = Float64(beta + Float64(i * 2.0)) t_1 = Float64(alpha + t_0) tmp = 0.0 if (i <= 4.9e+144) tmp = Float64(Float64(Float64(Float64(alpha * beta) + Float64(i * Float64(i + Float64(alpha + beta)))) * Float64(Float64(i * Float64(i + beta)) / Float64(t_0 * t_0))) / Float64(Float64(t_1 * t_1) + -1.0)); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = beta + (i * 2.0); t_1 = alpha + t_0; tmp = 0.0; if (i <= 4.9e+144) tmp = (((alpha * beta) + (i * (i + (alpha + beta)))) * ((i * (i + beta)) / (t_0 * t_0))) / ((t_1 * t_1) + -1.0); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(alpha + t$95$0), $MachinePrecision]}, If[LessEqual[i, 4.9e+144], N[(N[(N[(N[(alpha * beta), $MachinePrecision] + N[(i * N[(i + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(i * N[(i + beta), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$1 * t$95$1), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision], 0.0625]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \beta + i \cdot 2\\
t_1 := \alpha + t\_0\\
\mathbf{if}\;i \leq 4.9 \cdot 10^{+144}:\\
\;\;\;\;\frac{\left(\alpha \cdot \beta + i \cdot \left(i + \left(\alpha + \beta\right)\right)\right) \cdot \frac{i \cdot \left(i + \beta\right)}{t\_0 \cdot t\_0}}{t\_1 \cdot t\_1 + -1}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
\end{array}
if i < 4.9e144Initial program 41.9%
/-lowering-/.f64N/A
Simplified83.2%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.4%
Simplified75.4%
if 4.9e144 < i Initial program 0.1%
/-lowering-/.f64N/A
Simplified1.7%
Taylor expanded in i around inf
Simplified89.0%
Final simplification81.7%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (* i (+ i beta))) (t_1 (+ beta (* i 2.0))) (t_2 (+ alpha t_1)))
(if (<= i 7.8e+144)
(/ (* t_0 (/ t_0 (* t_1 t_1))) (+ (* t_2 t_2) -1.0))
0.0625)))
double code(double alpha, double beta, double i) {
double t_0 = i * (i + beta);
double t_1 = beta + (i * 2.0);
double t_2 = alpha + t_1;
double tmp;
if (i <= 7.8e+144) {
tmp = (t_0 * (t_0 / (t_1 * t_1))) / ((t_2 * t_2) + -1.0);
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = i * (i + beta)
t_1 = beta + (i * 2.0d0)
t_2 = alpha + t_1
if (i <= 7.8d+144) then
tmp = (t_0 * (t_0 / (t_1 * t_1))) / ((t_2 * t_2) + (-1.0d0))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = i * (i + beta);
double t_1 = beta + (i * 2.0);
double t_2 = alpha + t_1;
double tmp;
if (i <= 7.8e+144) {
tmp = (t_0 * (t_0 / (t_1 * t_1))) / ((t_2 * t_2) + -1.0);
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): t_0 = i * (i + beta) t_1 = beta + (i * 2.0) t_2 = alpha + t_1 tmp = 0 if i <= 7.8e+144: tmp = (t_0 * (t_0 / (t_1 * t_1))) / ((t_2 * t_2) + -1.0) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) t_0 = Float64(i * Float64(i + beta)) t_1 = Float64(beta + Float64(i * 2.0)) t_2 = Float64(alpha + t_1) tmp = 0.0 if (i <= 7.8e+144) tmp = Float64(Float64(t_0 * Float64(t_0 / Float64(t_1 * t_1))) / Float64(Float64(t_2 * t_2) + -1.0)); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = i * (i + beta); t_1 = beta + (i * 2.0); t_2 = alpha + t_1; tmp = 0.0; if (i <= 7.8e+144) tmp = (t_0 * (t_0 / (t_1 * t_1))) / ((t_2 * t_2) + -1.0); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(i * N[(i + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(alpha + t$95$1), $MachinePrecision]}, If[LessEqual[i, 7.8e+144], N[(N[(t$95$0 * N[(t$95$0 / N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$2 * t$95$2), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision], 0.0625]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot \left(i + \beta\right)\\
t_1 := \beta + i \cdot 2\\
t_2 := \alpha + t\_1\\
\mathbf{if}\;i \leq 7.8 \cdot 10^{+144}:\\
\;\;\;\;\frac{t\_0 \cdot \frac{t\_0}{t\_1 \cdot t\_1}}{t\_2 \cdot t\_2 + -1}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
\end{array}
if i < 7.80000000000000036e144Initial program 41.9%
/-lowering-/.f64N/A
Simplified83.2%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.4%
Simplified75.4%
Taylor expanded in alpha around 0
*-lowering-*.f64N/A
+-lowering-+.f6478.5%
Simplified78.5%
if 7.80000000000000036e144 < i Initial program 0.1%
/-lowering-/.f64N/A
Simplified1.7%
Taylor expanded in i around inf
Simplified89.0%
Final simplification83.3%
(FPCore (alpha beta i) :precision binary64 (let* ((t_0 (* i (+ i beta))) (t_1 (+ beta (* i 2.0))) (t_2 (* t_1 t_1))) (if (<= i 1.9e+144) (/ (/ t_0 (/ t_2 t_0)) (+ t_2 -1.0)) 0.0625)))
double code(double alpha, double beta, double i) {
double t_0 = i * (i + beta);
double t_1 = beta + (i * 2.0);
double t_2 = t_1 * t_1;
double tmp;
if (i <= 1.9e+144) {
tmp = (t_0 / (t_2 / t_0)) / (t_2 + -1.0);
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = i * (i + beta)
t_1 = beta + (i * 2.0d0)
t_2 = t_1 * t_1
if (i <= 1.9d+144) then
tmp = (t_0 / (t_2 / t_0)) / (t_2 + (-1.0d0))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = i * (i + beta);
double t_1 = beta + (i * 2.0);
double t_2 = t_1 * t_1;
double tmp;
if (i <= 1.9e+144) {
tmp = (t_0 / (t_2 / t_0)) / (t_2 + -1.0);
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): t_0 = i * (i + beta) t_1 = beta + (i * 2.0) t_2 = t_1 * t_1 tmp = 0 if i <= 1.9e+144: tmp = (t_0 / (t_2 / t_0)) / (t_2 + -1.0) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) t_0 = Float64(i * Float64(i + beta)) t_1 = Float64(beta + Float64(i * 2.0)) t_2 = Float64(t_1 * t_1) tmp = 0.0 if (i <= 1.9e+144) tmp = Float64(Float64(t_0 / Float64(t_2 / t_0)) / Float64(t_2 + -1.0)); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = i * (i + beta); t_1 = beta + (i * 2.0); t_2 = t_1 * t_1; tmp = 0.0; if (i <= 1.9e+144) tmp = (t_0 / (t_2 / t_0)) / (t_2 + -1.0); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(i * N[(i + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * t$95$1), $MachinePrecision]}, If[LessEqual[i, 1.9e+144], N[(N[(t$95$0 / N[(t$95$2 / t$95$0), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 + -1.0), $MachinePrecision]), $MachinePrecision], 0.0625]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot \left(i + \beta\right)\\
t_1 := \beta + i \cdot 2\\
t_2 := t\_1 \cdot t\_1\\
\mathbf{if}\;i \leq 1.9 \cdot 10^{+144}:\\
\;\;\;\;\frac{\frac{t\_0}{\frac{t\_2}{t\_0}}}{t\_2 + -1}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
\end{array}
if i < 1.90000000000000013e144Initial program 41.9%
/-lowering-/.f64N/A
Simplified83.2%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.4%
Simplified75.4%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified35.0%
unswap-sqrN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
times-fracN/A
pow2N/A
*-commutativeN/A
pow2N/A
*-lowering-*.f64N/A
Applied egg-rr73.3%
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr73.3%
if 1.90000000000000013e144 < i Initial program 0.1%
/-lowering-/.f64N/A
Simplified1.7%
Taylor expanded in i around inf
Simplified89.0%
Final simplification80.5%
(FPCore (alpha beta i) :precision binary64 (let* ((t_0 (* i (+ i beta))) (t_1 (+ beta (* i 2.0))) (t_2 (* t_1 t_1))) (if (<= i 6e+144) (* (/ t_0 t_2) (/ t_0 (+ t_2 -1.0))) 0.0625)))
double code(double alpha, double beta, double i) {
double t_0 = i * (i + beta);
double t_1 = beta + (i * 2.0);
double t_2 = t_1 * t_1;
double tmp;
if (i <= 6e+144) {
tmp = (t_0 / t_2) * (t_0 / (t_2 + -1.0));
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = i * (i + beta)
t_1 = beta + (i * 2.0d0)
t_2 = t_1 * t_1
if (i <= 6d+144) then
tmp = (t_0 / t_2) * (t_0 / (t_2 + (-1.0d0)))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = i * (i + beta);
double t_1 = beta + (i * 2.0);
double t_2 = t_1 * t_1;
double tmp;
if (i <= 6e+144) {
tmp = (t_0 / t_2) * (t_0 / (t_2 + -1.0));
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): t_0 = i * (i + beta) t_1 = beta + (i * 2.0) t_2 = t_1 * t_1 tmp = 0 if i <= 6e+144: tmp = (t_0 / t_2) * (t_0 / (t_2 + -1.0)) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) t_0 = Float64(i * Float64(i + beta)) t_1 = Float64(beta + Float64(i * 2.0)) t_2 = Float64(t_1 * t_1) tmp = 0.0 if (i <= 6e+144) tmp = Float64(Float64(t_0 / t_2) * Float64(t_0 / Float64(t_2 + -1.0))); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = i * (i + beta); t_1 = beta + (i * 2.0); t_2 = t_1 * t_1; tmp = 0.0; if (i <= 6e+144) tmp = (t_0 / t_2) * (t_0 / (t_2 + -1.0)); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(i * N[(i + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * t$95$1), $MachinePrecision]}, If[LessEqual[i, 6e+144], N[(N[(t$95$0 / t$95$2), $MachinePrecision] * N[(t$95$0 / N[(t$95$2 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0625]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot \left(i + \beta\right)\\
t_1 := \beta + i \cdot 2\\
t_2 := t\_1 \cdot t\_1\\
\mathbf{if}\;i \leq 6 \cdot 10^{+144}:\\
\;\;\;\;\frac{t\_0}{t\_2} \cdot \frac{t\_0}{t\_2 + -1}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
\end{array}
if i < 5.9999999999999998e144Initial program 41.9%
/-lowering-/.f64N/A
Simplified83.2%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.4%
Simplified75.4%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified35.0%
unswap-sqrN/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
times-fracN/A
pow2N/A
*-commutativeN/A
pow2N/A
*-lowering-*.f64N/A
Applied egg-rr73.3%
if 5.9999999999999998e144 < i Initial program 0.1%
/-lowering-/.f64N/A
Simplified1.7%
Taylor expanded in i around inf
Simplified89.0%
Final simplification80.5%
(FPCore (alpha beta i) :precision binary64 (if (<= beta 2.15e+201) (+ 0.0625 (/ 0.015625 (* i i))) (* (/ (+ i alpha) beta) (/ i beta))))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 2.15e+201) {
tmp = 0.0625 + (0.015625 / (i * i));
} else {
tmp = ((i + alpha) / beta) * (i / beta);
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: tmp
if (beta <= 2.15d+201) then
tmp = 0.0625d0 + (0.015625d0 / (i * i))
else
tmp = ((i + alpha) / beta) * (i / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 2.15e+201) {
tmp = 0.0625 + (0.015625 / (i * i));
} else {
tmp = ((i + alpha) / beta) * (i / beta);
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 2.15e+201: tmp = 0.0625 + (0.015625 / (i * i)) else: tmp = ((i + alpha) / beta) * (i / beta) return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 2.15e+201) tmp = Float64(0.0625 + Float64(0.015625 / Float64(i * i))); else tmp = Float64(Float64(Float64(i + alpha) / beta) * Float64(i / beta)); end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 2.15e+201) tmp = 0.0625 + (0.015625 / (i * i)); else tmp = ((i + alpha) / beta) * (i / beta); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 2.15e+201], N[(0.0625 + N[(0.015625 / N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(i + alpha), $MachinePrecision] / beta), $MachinePrecision] * N[(i / beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2.15 \cdot 10^{+201}:\\
\;\;\;\;0.0625 + \frac{0.015625}{i \cdot i}\\
\mathbf{else}:\\
\;\;\;\;\frac{i + \alpha}{\beta} \cdot \frac{i}{\beta}\\
\end{array}
\end{array}
if beta < 2.14999999999999995e201Initial program 24.8%
/-lowering-/.f64N/A
Simplified48.3%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6443.7%
Simplified43.7%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified20.7%
Taylor expanded in i around inf
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified76.6%
Taylor expanded in beta around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6479.9%
Simplified79.9%
if 2.14999999999999995e201 < beta Initial program 0.0%
/-lowering-/.f64N/A
Simplified19.5%
Taylor expanded in beta around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6426.0%
Simplified26.0%
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f6482.3%
Applied egg-rr82.3%
Final simplification80.1%
(FPCore (alpha beta i) :precision binary64 (if (<= beta 1.85e+201) (+ 0.0625 (/ 0.015625 (* i i))) (* (/ i beta) (/ i beta))))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 1.85e+201) {
tmp = 0.0625 + (0.015625 / (i * i));
} else {
tmp = (i / beta) * (i / beta);
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: tmp
if (beta <= 1.85d+201) then
tmp = 0.0625d0 + (0.015625d0 / (i * i))
else
tmp = (i / beta) * (i / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 1.85e+201) {
tmp = 0.0625 + (0.015625 / (i * i));
} else {
tmp = (i / beta) * (i / beta);
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 1.85e+201: tmp = 0.0625 + (0.015625 / (i * i)) else: tmp = (i / beta) * (i / beta) return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 1.85e+201) tmp = Float64(0.0625 + Float64(0.015625 / Float64(i * i))); else tmp = Float64(Float64(i / beta) * Float64(i / beta)); end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 1.85e+201) tmp = 0.0625 + (0.015625 / (i * i)); else tmp = (i / beta) * (i / beta); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 1.85e+201], N[(0.0625 + N[(0.015625 / N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(i / beta), $MachinePrecision] * N[(i / beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 1.85 \cdot 10^{+201}:\\
\;\;\;\;0.0625 + \frac{0.015625}{i \cdot i}\\
\mathbf{else}:\\
\;\;\;\;\frac{i}{\beta} \cdot \frac{i}{\beta}\\
\end{array}
\end{array}
if beta < 1.8499999999999999e201Initial program 24.8%
/-lowering-/.f64N/A
Simplified48.3%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6443.7%
Simplified43.7%
Taylor expanded in alpha around 0
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
Simplified20.7%
Taylor expanded in i around inf
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
Simplified76.6%
Taylor expanded in beta around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6479.9%
Simplified79.9%
if 1.8499999999999999e201 < beta Initial program 0.0%
/-lowering-/.f64N/A
Simplified19.5%
Taylor expanded in beta around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6426.0%
Simplified26.0%
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f6482.3%
Applied egg-rr82.3%
Taylor expanded in i around inf
/-lowering-/.f6473.9%
Simplified73.9%
Final simplification79.4%
(FPCore (alpha beta i) :precision binary64 (if (<= beta 1.82e+201) 0.0625 (* (/ i beta) (/ i beta))))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 1.82e+201) {
tmp = 0.0625;
} else {
tmp = (i / beta) * (i / beta);
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: tmp
if (beta <= 1.82d+201) then
tmp = 0.0625d0
else
tmp = (i / beta) * (i / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 1.82e+201) {
tmp = 0.0625;
} else {
tmp = (i / beta) * (i / beta);
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 1.82e+201: tmp = 0.0625 else: tmp = (i / beta) * (i / beta) return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 1.82e+201) tmp = 0.0625; else tmp = Float64(Float64(i / beta) * Float64(i / beta)); end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 1.82e+201) tmp = 0.0625; else tmp = (i / beta) * (i / beta); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 1.82e+201], 0.0625, N[(N[(i / beta), $MachinePrecision] * N[(i / beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 1.82 \cdot 10^{+201}:\\
\;\;\;\;0.0625\\
\mathbf{else}:\\
\;\;\;\;\frac{i}{\beta} \cdot \frac{i}{\beta}\\
\end{array}
\end{array}
if beta < 1.82e201Initial program 24.8%
/-lowering-/.f64N/A
Simplified48.3%
Taylor expanded in i around inf
Simplified79.5%
if 1.82e201 < beta Initial program 0.0%
/-lowering-/.f64N/A
Simplified19.5%
Taylor expanded in beta around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6426.0%
Simplified26.0%
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f6482.3%
Applied egg-rr82.3%
Taylor expanded in i around inf
/-lowering-/.f6473.9%
Simplified73.9%
(FPCore (alpha beta i) :precision binary64 (if (<= beta 4.8e+253) 0.0625 (* (/ i beta) (/ alpha beta))))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 4.8e+253) {
tmp = 0.0625;
} else {
tmp = (i / beta) * (alpha / beta);
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: tmp
if (beta <= 4.8d+253) then
tmp = 0.0625d0
else
tmp = (i / beta) * (alpha / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 4.8e+253) {
tmp = 0.0625;
} else {
tmp = (i / beta) * (alpha / beta);
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 4.8e+253: tmp = 0.0625 else: tmp = (i / beta) * (alpha / beta) return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 4.8e+253) tmp = 0.0625; else tmp = Float64(Float64(i / beta) * Float64(alpha / beta)); end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 4.8e+253) tmp = 0.0625; else tmp = (i / beta) * (alpha / beta); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 4.8e+253], 0.0625, N[(N[(i / beta), $MachinePrecision] * N[(alpha / beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 4.8 \cdot 10^{+253}:\\
\;\;\;\;0.0625\\
\mathbf{else}:\\
\;\;\;\;\frac{i}{\beta} \cdot \frac{\alpha}{\beta}\\
\end{array}
\end{array}
if beta < 4.79999999999999982e253Initial program 23.8%
/-lowering-/.f64N/A
Simplified46.8%
Taylor expanded in i around inf
Simplified77.3%
if 4.79999999999999982e253 < beta Initial program 0.0%
/-lowering-/.f64N/A
Simplified27.3%
Taylor expanded in beta around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6438.5%
Simplified38.5%
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f6486.8%
Applied egg-rr86.8%
Taylor expanded in i around 0
/-lowering-/.f6451.1%
Simplified51.1%
Final simplification76.2%
(FPCore (alpha beta i) :precision binary64 0.0625)
double code(double alpha, double beta, double i) {
return 0.0625;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
code = 0.0625d0
end function
public static double code(double alpha, double beta, double i) {
return 0.0625;
}
def code(alpha, beta, i): return 0.0625
function code(alpha, beta, i) return 0.0625 end
function tmp = code(alpha, beta, i) tmp = 0.0625; end
code[alpha_, beta_, i_] := 0.0625
\begin{array}{l}
\\
0.0625
\end{array}
Initial program 22.8%
/-lowering-/.f64N/A
Simplified45.9%
Taylor expanded in i around inf
Simplified74.5%
herbie shell --seed 2024161
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:precision binary64
:pre (and (and (> alpha -1.0) (> beta -1.0)) (> i 1.0))
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i)))) (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0)))