
(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)
use fmin_fmax_functions
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}
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}
Herbie found 5 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)
use fmin_fmax_functions
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}
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}
(FPCore (alpha beta i) :precision binary64 (if (<= i 3.4e-19) (/ (/ (* alpha 0.0) beta) beta) (- (+ 0.0625 (* 0.0625 (/ (fma 2.0 alpha (* 2.0 beta)) i))) (* 0.125 (/ (+ alpha beta) i)))))
double code(double alpha, double beta, double i) {
double tmp;
if (i <= 3.4e-19) {
tmp = ((alpha * 0.0) / beta) / beta;
} else {
tmp = (0.0625 + (0.0625 * (fma(2.0, alpha, (2.0 * beta)) / i))) - (0.125 * ((alpha + beta) / i));
}
return tmp;
}
function code(alpha, beta, i) tmp = 0.0 if (i <= 3.4e-19) tmp = Float64(Float64(Float64(alpha * 0.0) / beta) / beta); else tmp = Float64(Float64(0.0625 + Float64(0.0625 * Float64(fma(2.0, alpha, Float64(2.0 * beta)) / i))) - Float64(0.125 * Float64(Float64(alpha + beta) / i))); end return tmp end
code[alpha_, beta_, i_] := If[LessEqual[i, 3.4e-19], N[(N[(N[(alpha * 0.0), $MachinePrecision] / beta), $MachinePrecision] / beta), $MachinePrecision], N[(N[(0.0625 + N[(0.0625 * N[(N[(2.0 * alpha + N[(2.0 * beta), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.125 * N[(N[(alpha + beta), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\mathbf{if}\;i \leq 3.4 \cdot 10^{-19}:\\
\;\;\;\;\frac{\frac{\alpha \cdot 0}{\beta}}{\beta}\\
\mathbf{else}:\\
\;\;\;\;\left(0.0625 + 0.0625 \cdot \frac{\mathsf{fma}\left(2, \alpha, 2 \cdot \beta\right)}{i}\right) - 0.125 \cdot \frac{\alpha + \beta}{i}\\
\end{array}
if i < 3.4000000000000002e-19Initial program 30.1%
Taylor expanded in i around 0
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower--.f64N/A
lower-pow.f64N/A
lower-+.f6453.0%
Applied rewrites53.0%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-*.f64N/A
lower-pow.f6431.3%
Applied rewrites31.3%
lift-/.f64N/A
lift-pow.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6431.3%
Applied rewrites31.3%
Taylor expanded in undef-var around zero
Applied rewrites87.1%
if 3.4000000000000002e-19 < i Initial program 30.1%
Taylor expanded in i around 0
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower--.f64N/A
lower-pow.f64N/A
lower-+.f6453.0%
Applied rewrites53.0%
Taylor expanded in i around inf
lower--.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-fma.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-+.f6450.5%
Applied rewrites50.5%
(FPCore (alpha beta i) :precision binary64 (if (<= i 5.2e+235) (/ (/ (* alpha 0.0) beta) beta) 0.0625))
double code(double alpha, double beta, double i) {
double tmp;
if (i <= 5.2e+235) {
tmp = ((alpha * 0.0) / beta) / beta;
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
use fmin_fmax_functions
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: tmp
if (i <= 5.2d+235) then
tmp = ((alpha * 0.0d0) / beta) / beta
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (i <= 5.2e+235) {
tmp = ((alpha * 0.0) / beta) / beta;
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if i <= 5.2e+235: tmp = ((alpha * 0.0) / beta) / beta else: tmp = 0.0625 return tmp
function code(alpha, beta, i) tmp = 0.0 if (i <= 5.2e+235) tmp = Float64(Float64(Float64(alpha * 0.0) / beta) / beta); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (i <= 5.2e+235) tmp = ((alpha * 0.0) / beta) / beta; else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[i, 5.2e+235], N[(N[(N[(alpha * 0.0), $MachinePrecision] / beta), $MachinePrecision] / beta), $MachinePrecision], 0.0625]
\begin{array}{l}
\mathbf{if}\;i \leq 5.2 \cdot 10^{+235}:\\
\;\;\;\;\frac{\frac{\alpha \cdot 0}{\beta}}{\beta}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
if i < 5.1999999999999996e235Initial program 30.1%
Taylor expanded in i around 0
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower--.f64N/A
lower-pow.f64N/A
lower-+.f6453.0%
Applied rewrites53.0%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-*.f64N/A
lower-pow.f6431.3%
Applied rewrites31.3%
lift-/.f64N/A
lift-pow.f64N/A
unpow2N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6431.3%
Applied rewrites31.3%
Taylor expanded in undef-var around zero
Applied rewrites87.1%
if 5.1999999999999996e235 < i Initial program 30.1%
Taylor expanded in i around inf
Applied rewrites13.0%
(FPCore (alpha beta i) :precision binary64 (if (<= i 65000000.0) (* 0.0 (/ alpha (* beta beta))) 0.0625))
double code(double alpha, double beta, double i) {
double tmp;
if (i <= 65000000.0) {
tmp = 0.0 * (alpha / (beta * beta));
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
use fmin_fmax_functions
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: tmp
if (i <= 65000000.0d0) then
tmp = 0.0d0 * (alpha / (beta * beta))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (i <= 65000000.0) {
tmp = 0.0 * (alpha / (beta * beta));
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if i <= 65000000.0: tmp = 0.0 * (alpha / (beta * beta)) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) tmp = 0.0 if (i <= 65000000.0) tmp = Float64(0.0 * Float64(alpha / Float64(beta * beta))); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (i <= 65000000.0) tmp = 0.0 * (alpha / (beta * beta)); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[i, 65000000.0], N[(0.0 * N[(alpha / N[(beta * beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0625]
\begin{array}{l}
\mathbf{if}\;i \leq 65000000:\\
\;\;\;\;0 \cdot \frac{\alpha}{\beta \cdot \beta}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
if i < 6.5e7Initial program 30.1%
Taylor expanded in i around 0
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower--.f64N/A
lower-pow.f64N/A
lower-+.f6453.0%
Applied rewrites53.0%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-*.f64N/A
lower-pow.f6431.3%
Applied rewrites31.3%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6434.0%
lift-pow.f64N/A
unpow2N/A
lower-*.f6434.0%
Applied rewrites34.0%
Taylor expanded in undef-var around zero
Applied rewrites60.8%
if 6.5e7 < i Initial program 30.1%
Taylor expanded in i around inf
Applied rewrites13.0%
(FPCore (alpha beta i) :precision binary64 (if (<= i 5.5e+25) (* i (/ alpha (* beta beta))) 0.0625))
double code(double alpha, double beta, double i) {
double tmp;
if (i <= 5.5e+25) {
tmp = i * (alpha / (beta * beta));
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(alpha, beta, i)
use fmin_fmax_functions
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: tmp
if (i <= 5.5d+25) then
tmp = i * (alpha / (beta * beta))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double tmp;
if (i <= 5.5e+25) {
tmp = i * (alpha / (beta * beta));
} else {
tmp = 0.0625;
}
return tmp;
}
def code(alpha, beta, i): tmp = 0 if i <= 5.5e+25: tmp = i * (alpha / (beta * beta)) else: tmp = 0.0625 return tmp
function code(alpha, beta, i) tmp = 0.0 if (i <= 5.5e+25) tmp = Float64(i * Float64(alpha / Float64(beta * beta))); else tmp = 0.0625; end return tmp end
function tmp_2 = code(alpha, beta, i) tmp = 0.0; if (i <= 5.5e+25) tmp = i * (alpha / (beta * beta)); else tmp = 0.0625; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := If[LessEqual[i, 5.5e+25], N[(i * N[(alpha / N[(beta * beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0625]
\begin{array}{l}
\mathbf{if}\;i \leq 5.5 \cdot 10^{+25}:\\
\;\;\;\;i \cdot \frac{\alpha}{\beta \cdot \beta}\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
if i < 5.5000000000000002e25Initial program 30.1%
Taylor expanded in i around 0
lower-/.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower--.f64N/A
lower-pow.f64N/A
lower-+.f6453.0%
Applied rewrites53.0%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-*.f64N/A
lower-pow.f6431.3%
Applied rewrites31.3%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6434.0%
lift-pow.f64N/A
unpow2N/A
lower-*.f6434.0%
Applied rewrites34.0%
if 5.5000000000000002e25 < i Initial program 30.1%
Taylor expanded in i around inf
Applied rewrites13.0%
(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)
use fmin_fmax_functions
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
0.0625
Initial program 30.1%
Taylor expanded in i around inf
Applied rewrites13.0%
herbie shell --seed 2025313 -o setup:search
(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)))