
(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 6 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 (+ (+ alpha beta) (* i 2.0)))
(t_1 (* t_0 t_0))
(t_2 (+ i (+ alpha beta)))
(t_3 (* i t_2))
(t_4 (pow (fma i 2.0 (+ alpha beta)) 2.0)))
(if (<= (/ (/ (* t_3 (+ t_3 (* alpha beta))) t_1) (+ t_1 -1.0)) INFINITY)
(* (/ t_3 (+ t_4 -1.0)) (/ (fma i t_2 (* alpha beta)) t_4))
(-
(+ 0.0625 (* 0.0625 (/ (+ (* alpha 2.0) (* beta 2.0)) i)))
(* 0.125 (/ (+ alpha beta) i))))))
double code(double alpha, double beta, double i) {
double t_0 = (alpha + beta) + (i * 2.0);
double t_1 = t_0 * t_0;
double t_2 = i + (alpha + beta);
double t_3 = i * t_2;
double t_4 = pow(fma(i, 2.0, (alpha + beta)), 2.0);
double tmp;
if ((((t_3 * (t_3 + (alpha * beta))) / t_1) / (t_1 + -1.0)) <= ((double) INFINITY)) {
tmp = (t_3 / (t_4 + -1.0)) * (fma(i, t_2, (alpha * beta)) / t_4);
} else {
tmp = (0.0625 + (0.0625 * (((alpha * 2.0) + (beta * 2.0)) / i))) - (0.125 * ((alpha + beta) / i));
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(alpha + beta) + Float64(i * 2.0)) t_1 = Float64(t_0 * t_0) t_2 = Float64(i + Float64(alpha + beta)) t_3 = Float64(i * t_2) t_4 = fma(i, 2.0, Float64(alpha + beta)) ^ 2.0 tmp = 0.0 if (Float64(Float64(Float64(t_3 * Float64(t_3 + Float64(alpha * beta))) / t_1) / Float64(t_1 + -1.0)) <= Inf) tmp = Float64(Float64(t_3 / Float64(t_4 + -1.0)) * Float64(fma(i, t_2, Float64(alpha * beta)) / t_4)); else tmp = Float64(Float64(0.0625 + Float64(0.0625 * Float64(Float64(Float64(alpha * 2.0) + Float64(beta * 2.0)) / i))) - Float64(0.125 * Float64(Float64(alpha + beta) / i))); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(alpha + beta), $MachinePrecision] + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(i + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(i * t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[Power[N[(i * 2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[N[(N[(N[(t$95$3 * N[(t$95$3 + N[(alpha * beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] / N[(t$95$1 + -1.0), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(t$95$3 / N[(t$95$4 + -1.0), $MachinePrecision]), $MachinePrecision] * N[(N[(i * t$95$2 + N[(alpha * beta), $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision]), $MachinePrecision], N[(N[(0.0625 + N[(0.0625 * N[(N[(N[(alpha * 2.0), $MachinePrecision] + N[(beta * 2.0), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.125 * N[(N[(alpha + beta), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\alpha + \beta\right) + i \cdot 2\\
t_1 := t\_0 \cdot t\_0\\
t_2 := i + \left(\alpha + \beta\right)\\
t_3 := i \cdot t\_2\\
t_4 := {\left(\mathsf{fma}\left(i, 2, \alpha + \beta\right)\right)}^{2}\\
\mathbf{if}\;\frac{\frac{t\_3 \cdot \left(t\_3 + \alpha \cdot \beta\right)}{t\_1}}{t\_1 + -1} \leq \infty:\\
\;\;\;\;\frac{t\_3}{t\_4 + -1} \cdot \frac{\mathsf{fma}\left(i, t\_2, \alpha \cdot \beta\right)}{t\_4}\\
\mathbf{else}:\\
\;\;\;\;\left(0.0625 + 0.0625 \cdot \frac{\alpha \cdot 2 + \beta \cdot 2}{i}\right) - 0.125 \cdot \frac{\alpha + \beta}{i}\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (*.f64 i (+.f64 (+.f64 alpha beta) i)) (+.f64 (*.f64 beta alpha) (*.f64 i (+.f64 (+.f64 alpha beta) i)))) (*.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)))) (-.f64 (*.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) #s(literal 1 binary64))) < +inf.0Initial program 53.1%
associate-/l/48.3%
Simplified48.3%
Applied egg-rr99.8%
if +inf.0 < (/.f64 (/.f64 (*.f64 (*.f64 i (+.f64 (+.f64 alpha beta) i)) (+.f64 (*.f64 beta alpha) (*.f64 i (+.f64 (+.f64 alpha beta) i)))) (*.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)))) (-.f64 (*.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) #s(literal 1 binary64))) Initial program 0.0%
associate-/l/0.0%
associate-/l*4.5%
+-commutative4.5%
+-commutative4.5%
+-commutative4.5%
associate-+l+4.5%
+-commutative4.5%
associate-*l*4.5%
Simplified4.5%
Taylor expanded in i around inf 73.9%
Final simplification82.5%
(FPCore (alpha beta i)
:precision binary64
(if (<= beta 9e+208)
(-
(+ 0.0625 (* 0.0625 (/ (+ (* alpha 2.0) (* beta 2.0)) i)))
(* 0.125 (/ (+ alpha beta) i)))
(* i (/ (/ (+ i alpha) beta) beta))))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 9e+208) {
tmp = (0.0625 + (0.0625 * (((alpha * 2.0) + (beta * 2.0)) / i))) - (0.125 * ((alpha + beta) / i));
} else {
tmp = i * (((i + alpha) / beta) / 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 <= 9d+208) then
tmp = (0.0625d0 + (0.0625d0 * (((alpha * 2.0d0) + (beta * 2.0d0)) / i))) - (0.125d0 * ((alpha + beta) / i))
else
tmp = i * (((i + alpha) / beta) / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 9e+208) {
tmp = (0.0625 + (0.0625 * (((alpha * 2.0) + (beta * 2.0)) / i))) - (0.125 * ((alpha + beta) / i));
} else {
tmp = i * (((i + alpha) / beta) / beta);
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 9e+208: tmp = (0.0625 + (0.0625 * (((alpha * 2.0) + (beta * 2.0)) / i))) - (0.125 * ((alpha + beta) / i)) else: tmp = i * (((i + alpha) / beta) / beta) return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 9e+208) tmp = Float64(Float64(0.0625 + Float64(0.0625 * Float64(Float64(Float64(alpha * 2.0) + Float64(beta * 2.0)) / i))) - Float64(0.125 * Float64(Float64(alpha + beta) / i))); else tmp = Float64(i * Float64(Float64(Float64(i + alpha) / beta) / beta)); end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 9e+208) tmp = (0.0625 + (0.0625 * (((alpha * 2.0) + (beta * 2.0)) / i))) - (0.125 * ((alpha + beta) / i)); else tmp = i * (((i + alpha) / beta) / beta); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 9e+208], N[(N[(0.0625 + N[(0.0625 * N[(N[(N[(alpha * 2.0), $MachinePrecision] + N[(beta * 2.0), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.125 * N[(N[(alpha + beta), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(i * N[(N[(N[(i + alpha), $MachinePrecision] / beta), $MachinePrecision] / beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 9 \cdot 10^{+208}:\\
\;\;\;\;\left(0.0625 + 0.0625 \cdot \frac{\alpha \cdot 2 + \beta \cdot 2}{i}\right) - 0.125 \cdot \frac{\alpha + \beta}{i}\\
\mathbf{else}:\\
\;\;\;\;i \cdot \frac{\frac{i + \alpha}{\beta}}{\beta}\\
\end{array}
\end{array}
if beta < 9.00000000000000029e208Initial program 19.4%
associate-/l/17.7%
associate-/l*20.6%
+-commutative20.6%
+-commutative20.6%
+-commutative20.6%
associate-+l+20.6%
+-commutative20.6%
associate-*l*20.5%
Simplified20.5%
Taylor expanded in i around inf 83.5%
if 9.00000000000000029e208 < beta Initial program 0.0%
associate-/l/0.0%
associate-/l*8.7%
+-commutative8.7%
+-commutative8.7%
+-commutative8.7%
associate-+l+8.7%
+-commutative8.7%
associate-*l*8.7%
Simplified8.7%
Taylor expanded in beta around inf 21.9%
associate-/l*24.5%
Simplified24.5%
add-exp-log24.5%
div-inv24.5%
pow-flip24.5%
metadata-eval24.5%
Applied egg-rr24.5%
rem-exp-log24.5%
*-commutative24.5%
sqr-pow24.5%
associate-*l*44.8%
metadata-eval44.8%
unpow-144.8%
metadata-eval44.8%
unpow-144.8%
Applied egg-rr44.8%
associate-*l/44.7%
*-un-lft-identity44.7%
associate-*l/44.7%
*-un-lft-identity44.7%
+-commutative44.7%
Applied egg-rr44.7%
Final simplification79.9%
(FPCore (alpha beta i) :precision binary64 (if (<= beta 8.2e+208) 0.0625 (* i (/ (/ (+ i alpha) beta) beta))))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 8.2e+208) {
tmp = 0.0625;
} else {
tmp = i * (((i + alpha) / beta) / 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 <= 8.2d+208) then
tmp = 0.0625d0
else
tmp = i * (((i + alpha) / beta) / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 8.2e+208) {
tmp = 0.0625;
} else {
tmp = i * (((i + alpha) / beta) / beta);
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 8.2e+208: tmp = 0.0625 else: tmp = i * (((i + alpha) / beta) / beta) return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 8.2e+208) tmp = 0.0625; else tmp = Float64(i * Float64(Float64(Float64(i + alpha) / beta) / beta)); end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 8.2e+208) tmp = 0.0625; else tmp = i * (((i + alpha) / beta) / beta); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 8.2e+208], 0.0625, N[(i * N[(N[(N[(i + alpha), $MachinePrecision] / beta), $MachinePrecision] / beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 8.2 \cdot 10^{+208}:\\
\;\;\;\;0.0625\\
\mathbf{else}:\\
\;\;\;\;i \cdot \frac{\frac{i + \alpha}{\beta}}{\beta}\\
\end{array}
\end{array}
if beta < 8.1999999999999996e208Initial program 19.4%
associate-/l/17.7%
associate-/l*20.6%
+-commutative20.6%
+-commutative20.6%
+-commutative20.6%
associate-+l+20.6%
+-commutative20.6%
associate-*l*20.5%
Simplified20.5%
Taylor expanded in i around inf 80.9%
if 8.1999999999999996e208 < beta Initial program 0.0%
associate-/l/0.0%
associate-/l*8.7%
+-commutative8.7%
+-commutative8.7%
+-commutative8.7%
associate-+l+8.7%
+-commutative8.7%
associate-*l*8.7%
Simplified8.7%
Taylor expanded in beta around inf 21.9%
associate-/l*24.5%
Simplified24.5%
add-exp-log24.5%
div-inv24.5%
pow-flip24.5%
metadata-eval24.5%
Applied egg-rr24.5%
rem-exp-log24.5%
*-commutative24.5%
sqr-pow24.5%
associate-*l*44.8%
metadata-eval44.8%
unpow-144.8%
metadata-eval44.8%
unpow-144.8%
Applied egg-rr44.8%
associate-*l/44.7%
*-un-lft-identity44.7%
associate-*l/44.7%
*-un-lft-identity44.7%
+-commutative44.7%
Applied egg-rr44.7%
(FPCore (alpha beta i) :precision binary64 (if (<= beta 3.3e+208) 0.0625 (* i (* (/ 1.0 beta) (/ i beta)))))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 3.3e+208) {
tmp = 0.0625;
} else {
tmp = i * ((1.0 / 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 <= 3.3d+208) then
tmp = 0.0625d0
else
tmp = i * ((1.0d0 / beta) * (i / beta))
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 3.3e+208) {
tmp = 0.0625;
} else {
tmp = i * ((1.0 / beta) * (i / beta));
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 3.3e+208: tmp = 0.0625 else: tmp = i * ((1.0 / beta) * (i / beta)) return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 3.3e+208) tmp = 0.0625; else tmp = Float64(i * Float64(Float64(1.0 / beta) * Float64(i / beta))); end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 3.3e+208) tmp = 0.0625; else tmp = i * ((1.0 / beta) * (i / beta)); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 3.3e+208], 0.0625, N[(i * N[(N[(1.0 / beta), $MachinePrecision] * N[(i / beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 3.3 \cdot 10^{+208}:\\
\;\;\;\;0.0625\\
\mathbf{else}:\\
\;\;\;\;i \cdot \left(\frac{1}{\beta} \cdot \frac{i}{\beta}\right)\\
\end{array}
\end{array}
if beta < 3.3e208Initial program 19.4%
associate-/l/17.7%
associate-/l*20.6%
+-commutative20.6%
+-commutative20.6%
+-commutative20.6%
associate-+l+20.6%
+-commutative20.6%
associate-*l*20.5%
Simplified20.5%
Taylor expanded in i around inf 80.9%
if 3.3e208 < beta Initial program 0.0%
associate-/l/0.0%
associate-/l*8.7%
+-commutative8.7%
+-commutative8.7%
+-commutative8.7%
associate-+l+8.7%
+-commutative8.7%
associate-*l*8.7%
Simplified8.7%
Taylor expanded in beta around inf 21.9%
associate-/l*24.5%
Simplified24.5%
add-exp-log24.5%
div-inv24.5%
pow-flip24.5%
metadata-eval24.5%
Applied egg-rr24.5%
rem-exp-log24.5%
*-commutative24.5%
sqr-pow24.5%
associate-*l*44.8%
metadata-eval44.8%
unpow-144.8%
metadata-eval44.8%
unpow-144.8%
Applied egg-rr44.8%
Taylor expanded in alpha around 0 44.7%
(FPCore (alpha beta i) :precision binary64 (if (<= beta 9.5e+209) 0.0625 0.0))
double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 9.5e+209) {
tmp = 0.0625;
} else {
tmp = 0.0;
}
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 <= 9.5d+209) then
tmp = 0.0625d0
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (beta <= 9.5e+209) {
tmp = 0.0625;
} else {
tmp = 0.0;
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if beta <= 9.5e+209: tmp = 0.0625 else: tmp = 0.0 return tmp
function code(alpha, beta, i) tmp = 0.0 if (beta <= 9.5e+209) tmp = 0.0625; else tmp = 0.0; end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (beta <= 9.5e+209) tmp = 0.0625; else tmp = 0.0; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[beta, 9.5e+209], 0.0625, 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 9.5 \cdot 10^{+209}:\\
\;\;\;\;0.0625\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if beta < 9.50000000000000069e209Initial program 19.4%
associate-/l/17.7%
associate-/l*20.6%
+-commutative20.6%
+-commutative20.6%
+-commutative20.6%
associate-+l+20.6%
+-commutative20.6%
associate-*l*20.5%
Simplified20.5%
Taylor expanded in i around inf 80.9%
if 9.50000000000000069e209 < beta Initial program 0.0%
associate-/l/0.0%
associate-/l*8.7%
+-commutative8.7%
+-commutative8.7%
+-commutative8.7%
associate-+l+8.7%
+-commutative8.7%
associate-*l*8.7%
Simplified8.7%
Taylor expanded in i around inf 32.9%
Taylor expanded in i around 0 24.5%
distribute-lft-in24.5%
associate-*r*24.5%
metadata-eval24.5%
div-sub24.5%
+-inverses24.5%
Simplified24.5%
(FPCore (alpha beta i) :precision binary64 0.0)
double code(double alpha, double beta, double i) {
return 0.0;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
code = 0.0d0
end function
public static double code(double alpha, double beta, double i) {
return 0.0;
}
def code(alpha, beta, i): return 0.0
function code(alpha, beta, i) return 0.0 end
function tmp = code(alpha, beta, i) tmp = 0.0; end
code[alpha_, beta_, i_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 17.6%
associate-/l/16.0%
associate-/l*19.5%
+-commutative19.5%
+-commutative19.5%
+-commutative19.5%
associate-+l+19.5%
+-commutative19.5%
associate-*l*19.4%
Simplified19.4%
Taylor expanded in i around inf 78.8%
Taylor expanded in i around 0 7.7%
distribute-lft-in7.7%
associate-*r*7.7%
metadata-eval7.7%
div-sub7.7%
+-inverses7.7%
Simplified7.7%
herbie shell --seed 2024156
(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)))