
(FPCore (i) :precision binary64 (let* ((t_0 (* (* 2.0 i) (* 2.0 i)))) (/ (/ (* (* i i) (* i i)) t_0) (- t_0 1.0))))
double code(double i) {
double t_0 = (2.0 * i) * (2.0 * i);
return (((i * i) * (i * i)) / t_0) / (t_0 - 1.0);
}
real(8) function code(i)
real(8), intent (in) :: i
real(8) :: t_0
t_0 = (2.0d0 * i) * (2.0d0 * i)
code = (((i * i) * (i * i)) / t_0) / (t_0 - 1.0d0)
end function
public static double code(double i) {
double t_0 = (2.0 * i) * (2.0 * i);
return (((i * i) * (i * i)) / t_0) / (t_0 - 1.0);
}
def code(i): t_0 = (2.0 * i) * (2.0 * i) return (((i * i) * (i * i)) / t_0) / (t_0 - 1.0)
function code(i) t_0 = Float64(Float64(2.0 * i) * Float64(2.0 * i)) return Float64(Float64(Float64(Float64(i * i) * Float64(i * i)) / t_0) / Float64(t_0 - 1.0)) end
function tmp = code(i) t_0 = (2.0 * i) * (2.0 * i); tmp = (((i * i) * (i * i)) / t_0) / (t_0 - 1.0); end
code[i_] := Block[{t$95$0 = N[(N[(2.0 * i), $MachinePrecision] * N[(2.0 * i), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(i * i), $MachinePrecision] * N[(i * i), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 - 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(2 \cdot i\right) \cdot \left(2 \cdot i\right)\\
\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{t_0}}{t_0 - 1}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (i) :precision binary64 (let* ((t_0 (* (* 2.0 i) (* 2.0 i)))) (/ (/ (* (* i i) (* i i)) t_0) (- t_0 1.0))))
double code(double i) {
double t_0 = (2.0 * i) * (2.0 * i);
return (((i * i) * (i * i)) / t_0) / (t_0 - 1.0);
}
real(8) function code(i)
real(8), intent (in) :: i
real(8) :: t_0
t_0 = (2.0d0 * i) * (2.0d0 * i)
code = (((i * i) * (i * i)) / t_0) / (t_0 - 1.0d0)
end function
public static double code(double i) {
double t_0 = (2.0 * i) * (2.0 * i);
return (((i * i) * (i * i)) / t_0) / (t_0 - 1.0);
}
def code(i): t_0 = (2.0 * i) * (2.0 * i) return (((i * i) * (i * i)) / t_0) / (t_0 - 1.0)
function code(i) t_0 = Float64(Float64(2.0 * i) * Float64(2.0 * i)) return Float64(Float64(Float64(Float64(i * i) * Float64(i * i)) / t_0) / Float64(t_0 - 1.0)) end
function tmp = code(i) t_0 = (2.0 * i) * (2.0 * i); tmp = (((i * i) * (i * i)) / t_0) / (t_0 - 1.0); end
code[i_] := Block[{t$95$0 = N[(N[(2.0 * i), $MachinePrecision] * N[(2.0 * i), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(i * i), $MachinePrecision] * N[(i * i), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 - 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(2 \cdot i\right) \cdot \left(2 \cdot i\right)\\
\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{t_0}}{t_0 - 1}
\end{array}
\end{array}
(FPCore (i) :precision binary64 (if (<= i 0.5) (* (* i i) (- -0.25 (* i i))) (+ 0.0625 (/ 0.015625 (* i i)))))
double code(double i) {
double tmp;
if (i <= 0.5) {
tmp = (i * i) * (-0.25 - (i * i));
} else {
tmp = 0.0625 + (0.015625 / (i * i));
}
return tmp;
}
real(8) function code(i)
real(8), intent (in) :: i
real(8) :: tmp
if (i <= 0.5d0) then
tmp = (i * i) * ((-0.25d0) - (i * i))
else
tmp = 0.0625d0 + (0.015625d0 / (i * i))
end if
code = tmp
end function
public static double code(double i) {
double tmp;
if (i <= 0.5) {
tmp = (i * i) * (-0.25 - (i * i));
} else {
tmp = 0.0625 + (0.015625 / (i * i));
}
return tmp;
}
def code(i): tmp = 0 if i <= 0.5: tmp = (i * i) * (-0.25 - (i * i)) else: tmp = 0.0625 + (0.015625 / (i * i)) return tmp
function code(i) tmp = 0.0 if (i <= 0.5) tmp = Float64(Float64(i * i) * Float64(-0.25 - Float64(i * i))); else tmp = Float64(0.0625 + Float64(0.015625 / Float64(i * i))); end return tmp end
function tmp_2 = code(i) tmp = 0.0; if (i <= 0.5) tmp = (i * i) * (-0.25 - (i * i)); else tmp = 0.0625 + (0.015625 / (i * i)); end tmp_2 = tmp; end
code[i_] := If[LessEqual[i, 0.5], N[(N[(i * i), $MachinePrecision] * N[(-0.25 - N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0625 + N[(0.015625 / N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq 0.5:\\
\;\;\;\;\left(i \cdot i\right) \cdot \left(-0.25 - i \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;0.0625 + \frac{0.015625}{i \cdot i}\\
\end{array}
\end{array}
if i < 0.5Initial program 28.6%
associate-/l*56.8%
swap-sqr56.8%
metadata-eval56.8%
sub-neg56.8%
swap-sqr56.8%
metadata-eval56.8%
metadata-eval56.8%
Simplified56.8%
frac-2neg56.8%
div-inv56.8%
div-inv56.8%
clear-num56.8%
distribute-rgt-neg-in56.8%
*-commutative56.8%
associate-/r*56.8%
*-inverses100.0%
metadata-eval100.0%
metadata-eval100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
*-commutative100.0%
distribute-rgt-neg-in100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in i around 0 99.7%
unpow299.7%
*-commutative99.7%
associate-*r*99.7%
Simplified99.7%
Taylor expanded in i around 0 99.7%
+-commutative99.7%
unpow299.7%
fma-def99.7%
mul-1-neg99.7%
fma-neg99.7%
metadata-eval99.7%
pow-sqr99.7%
unpow299.7%
unpow299.7%
distribute-rgt-out--99.7%
Simplified99.7%
if 0.5 < i Initial program 20.1%
times-frac45.9%
associate-/l*45.9%
associate-/l*46.0%
associate-/l/46.0%
associate-/r/46.0%
associate-/l*46.0%
*-inverses46.0%
metadata-eval46.0%
associate-*l*46.0%
*-commutative46.0%
associate-/r*46.0%
Simplified100.0%
Taylor expanded in i around inf 99.5%
unpow299.5%
associate-*r/99.5%
metadata-eval99.5%
Simplified99.5%
Final simplification99.6%
(FPCore (i) :precision binary64 (if (<= i 0.5) (* i (* i -0.25)) (+ 0.0625 (/ 0.015625 (* i i)))))
double code(double i) {
double tmp;
if (i <= 0.5) {
tmp = i * (i * -0.25);
} else {
tmp = 0.0625 + (0.015625 / (i * i));
}
return tmp;
}
real(8) function code(i)
real(8), intent (in) :: i
real(8) :: tmp
if (i <= 0.5d0) then
tmp = i * (i * (-0.25d0))
else
tmp = 0.0625d0 + (0.015625d0 / (i * i))
end if
code = tmp
end function
public static double code(double i) {
double tmp;
if (i <= 0.5) {
tmp = i * (i * -0.25);
} else {
tmp = 0.0625 + (0.015625 / (i * i));
}
return tmp;
}
def code(i): tmp = 0 if i <= 0.5: tmp = i * (i * -0.25) else: tmp = 0.0625 + (0.015625 / (i * i)) return tmp
function code(i) tmp = 0.0 if (i <= 0.5) tmp = Float64(i * Float64(i * -0.25)); else tmp = Float64(0.0625 + Float64(0.015625 / Float64(i * i))); end return tmp end
function tmp_2 = code(i) tmp = 0.0; if (i <= 0.5) tmp = i * (i * -0.25); else tmp = 0.0625 + (0.015625 / (i * i)); end tmp_2 = tmp; end
code[i_] := If[LessEqual[i, 0.5], N[(i * N[(i * -0.25), $MachinePrecision]), $MachinePrecision], N[(0.0625 + N[(0.015625 / N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq 0.5:\\
\;\;\;\;i \cdot \left(i \cdot -0.25\right)\\
\mathbf{else}:\\
\;\;\;\;0.0625 + \frac{0.015625}{i \cdot i}\\
\end{array}
\end{array}
if i < 0.5Initial program 28.6%
times-frac99.9%
associate-/l*99.7%
associate-/l*99.8%
associate-/l/99.8%
associate-/r/98.6%
associate-/l*98.6%
*-inverses98.6%
metadata-eval98.6%
associate-*l*98.6%
*-commutative98.6%
associate-/r*98.8%
Simplified82.0%
Taylor expanded in i around 0 98.9%
unpow298.9%
*-commutative98.9%
associate-*l*98.9%
Simplified98.9%
if 0.5 < i Initial program 20.1%
times-frac45.9%
associate-/l*45.9%
associate-/l*46.0%
associate-/l/46.0%
associate-/r/46.0%
associate-/l*46.0%
*-inverses46.0%
metadata-eval46.0%
associate-*l*46.0%
*-commutative46.0%
associate-/r*46.0%
Simplified100.0%
Taylor expanded in i around inf 99.5%
unpow299.5%
associate-*r/99.5%
metadata-eval99.5%
Simplified99.5%
Final simplification99.2%
(FPCore (i) :precision binary64 (/ 0.25 (- 4.0 (/ 1.0 (* i i)))))
double code(double i) {
return 0.25 / (4.0 - (1.0 / (i * i)));
}
real(8) function code(i)
real(8), intent (in) :: i
code = 0.25d0 / (4.0d0 - (1.0d0 / (i * i)))
end function
public static double code(double i) {
return 0.25 / (4.0 - (1.0 / (i * i)));
}
def code(i): return 0.25 / (4.0 - (1.0 / (i * i)))
function code(i) return Float64(0.25 / Float64(4.0 - Float64(1.0 / Float64(i * i)))) end
function tmp = code(i) tmp = 0.25 / (4.0 - (1.0 / (i * i))); end
code[i_] := N[(0.25 / N[(4.0 - N[(1.0 / N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.25}{4 - \frac{1}{i \cdot i}}
\end{array}
Initial program 24.5%
times-frac73.7%
associate-/l*73.6%
associate-/l*73.7%
associate-/l/73.7%
associate-/r/73.1%
associate-/l*73.1%
*-inverses73.1%
metadata-eval73.1%
associate-*l*73.1%
*-commutative73.1%
associate-/r*73.2%
Simplified90.7%
Taylor expanded in i around 0 99.4%
unpow299.4%
Simplified99.4%
Final simplification99.4%
(FPCore (i) :precision binary64 (if (<= i 0.5) (* i (* i -0.25)) 0.0625))
double code(double i) {
double tmp;
if (i <= 0.5) {
tmp = i * (i * -0.25);
} else {
tmp = 0.0625;
}
return tmp;
}
real(8) function code(i)
real(8), intent (in) :: i
real(8) :: tmp
if (i <= 0.5d0) then
tmp = i * (i * (-0.25d0))
else
tmp = 0.0625d0
end if
code = tmp
end function
public static double code(double i) {
double tmp;
if (i <= 0.5) {
tmp = i * (i * -0.25);
} else {
tmp = 0.0625;
}
return tmp;
}
def code(i): tmp = 0 if i <= 0.5: tmp = i * (i * -0.25) else: tmp = 0.0625 return tmp
function code(i) tmp = 0.0 if (i <= 0.5) tmp = Float64(i * Float64(i * -0.25)); else tmp = 0.0625; end return tmp end
function tmp_2 = code(i) tmp = 0.0; if (i <= 0.5) tmp = i * (i * -0.25); else tmp = 0.0625; end tmp_2 = tmp; end
code[i_] := If[LessEqual[i, 0.5], N[(i * N[(i * -0.25), $MachinePrecision]), $MachinePrecision], 0.0625]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq 0.5:\\
\;\;\;\;i \cdot \left(i \cdot -0.25\right)\\
\mathbf{else}:\\
\;\;\;\;0.0625\\
\end{array}
\end{array}
if i < 0.5Initial program 28.6%
times-frac99.9%
associate-/l*99.7%
associate-/l*99.8%
associate-/l/99.8%
associate-/r/98.6%
associate-/l*98.6%
*-inverses98.6%
metadata-eval98.6%
associate-*l*98.6%
*-commutative98.6%
associate-/r*98.8%
Simplified82.0%
Taylor expanded in i around 0 98.9%
unpow298.9%
*-commutative98.9%
associate-*l*98.9%
Simplified98.9%
if 0.5 < i Initial program 20.1%
times-frac45.9%
associate-/l*45.9%
associate-/l*46.0%
associate-/l/46.0%
associate-/r/46.0%
associate-/l*46.0%
*-inverses46.0%
metadata-eval46.0%
associate-*l*46.0%
*-commutative46.0%
associate-/r*46.0%
Simplified100.0%
Taylor expanded in i around inf 99.4%
Final simplification99.2%
(FPCore (i) :precision binary64 0.0625)
double code(double i) {
return 0.0625;
}
real(8) function code(i)
real(8), intent (in) :: i
code = 0.0625d0
end function
public static double code(double i) {
return 0.0625;
}
def code(i): return 0.0625
function code(i) return 0.0625 end
function tmp = code(i) tmp = 0.0625; end
code[i_] := 0.0625
\begin{array}{l}
\\
0.0625
\end{array}
Initial program 24.5%
times-frac73.7%
associate-/l*73.6%
associate-/l*73.7%
associate-/l/73.7%
associate-/r/73.1%
associate-/l*73.1%
*-inverses73.1%
metadata-eval73.1%
associate-*l*73.1%
*-commutative73.1%
associate-/r*73.2%
Simplified90.7%
Taylor expanded in i around inf 49.5%
Final simplification49.5%
herbie shell --seed 2023201
(FPCore (i)
:name "Octave 3.8, jcobi/4, as called"
:precision binary64
:pre (> i 0.0)
(/ (/ (* (* i i) (* i i)) (* (* 2.0 i) (* 2.0 i))) (- (* (* 2.0 i) (* 2.0 i)) 1.0)))