
(FPCore (x) :precision binary64 (- (* 0.954929658551372 x) (* 0.12900613773279798 (* (* x x) x))))
double code(double x) {
return (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (0.954929658551372d0 * x) - (0.12900613773279798d0 * ((x * x) * x))
end function
public static double code(double x) {
return (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x));
}
def code(x): return (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x))
function code(x) return Float64(Float64(0.954929658551372 * x) - Float64(0.12900613773279798 * Float64(Float64(x * x) * x))) end
function tmp = code(x) tmp = (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x)); end
code[x_] := N[(N[(0.954929658551372 * x), $MachinePrecision] - N[(0.12900613773279798 * N[(N[(x * x), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.954929658551372 \cdot x - 0.12900613773279798 \cdot \left(\left(x \cdot x\right) \cdot x\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (- (* 0.954929658551372 x) (* 0.12900613773279798 (* (* x x) x))))
double code(double x) {
return (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (0.954929658551372d0 * x) - (0.12900613773279798d0 * ((x * x) * x))
end function
public static double code(double x) {
return (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x));
}
def code(x): return (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x))
function code(x) return Float64(Float64(0.954929658551372 * x) - Float64(0.12900613773279798 * Float64(Float64(x * x) * x))) end
function tmp = code(x) tmp = (0.954929658551372 * x) - (0.12900613773279798 * ((x * x) * x)); end
code[x_] := N[(N[(0.954929658551372 * x), $MachinePrecision] - N[(0.12900613773279798 * N[(N[(x * x), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.954929658551372 \cdot x - 0.12900613773279798 \cdot \left(\left(x \cdot x\right) \cdot x\right)
\end{array}
(FPCore (x) :precision binary64 (+ (* 0.954929658551372 x) (* (pow x 3.0) -0.12900613773279798)))
double code(double x) {
return (0.954929658551372 * x) + (pow(x, 3.0) * -0.12900613773279798);
}
real(8) function code(x)
real(8), intent (in) :: x
code = (0.954929658551372d0 * x) + ((x ** 3.0d0) * (-0.12900613773279798d0))
end function
public static double code(double x) {
return (0.954929658551372 * x) + (Math.pow(x, 3.0) * -0.12900613773279798);
}
def code(x): return (0.954929658551372 * x) + (math.pow(x, 3.0) * -0.12900613773279798)
function code(x) return Float64(Float64(0.954929658551372 * x) + Float64((x ^ 3.0) * -0.12900613773279798)) end
function tmp = code(x) tmp = (0.954929658551372 * x) + ((x ^ 3.0) * -0.12900613773279798); end
code[x_] := N[(N[(0.954929658551372 * x), $MachinePrecision] + N[(N[Power[x, 3.0], $MachinePrecision] * -0.12900613773279798), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.954929658551372 \cdot x + {x}^{3} \cdot -0.12900613773279798
\end{array}
Initial program 99.7%
sub-neg99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
unpow399.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (x) :precision binary64 (* x (- 0.954929658551372 (* 0.12900613773279798 (pow x 2.0)))))
double code(double x) {
return x * (0.954929658551372 - (0.12900613773279798 * pow(x, 2.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (0.954929658551372d0 - (0.12900613773279798d0 * (x ** 2.0d0)))
end function
public static double code(double x) {
return x * (0.954929658551372 - (0.12900613773279798 * Math.pow(x, 2.0)));
}
def code(x): return x * (0.954929658551372 - (0.12900613773279798 * math.pow(x, 2.0)))
function code(x) return Float64(x * Float64(0.954929658551372 - Float64(0.12900613773279798 * (x ^ 2.0)))) end
function tmp = code(x) tmp = x * (0.954929658551372 - (0.12900613773279798 * (x ^ 2.0))); end
code[x_] := N[(x * N[(0.954929658551372 - N[(0.12900613773279798 * N[Power[x, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(0.954929658551372 - 0.12900613773279798 \cdot {x}^{2}\right)
\end{array}
Initial program 99.7%
sub-neg99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
unpow399.8%
metadata-eval99.8%
Simplified99.8%
pow399.7%
*-commutative99.7%
metadata-eval99.7%
cancel-sign-sub-inv99.7%
associate-*r*99.8%
distribute-rgt-out--99.8%
pow299.8%
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (x) :precision binary64 (* x (fma (* x -0.12900613773279798) x 0.954929658551372)))
double code(double x) {
return x * fma((x * -0.12900613773279798), x, 0.954929658551372);
}
function code(x) return Float64(x * fma(Float64(x * -0.12900613773279798), x, 0.954929658551372)) end
code[x_] := N[(x * N[(N[(x * -0.12900613773279798), $MachinePrecision] * x + 0.954929658551372), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \mathsf{fma}\left(x \cdot -0.12900613773279798, x, 0.954929658551372\right)
\end{array}
Initial program 99.7%
sub-neg99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
unpow399.8%
metadata-eval99.8%
Simplified99.8%
+-commutative99.8%
pow399.7%
fma-define99.7%
pow399.8%
Applied egg-rr99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in x around 0 99.8%
*-commutative99.8%
*-commutative99.8%
cancel-sign-sub99.8%
distribute-lft-neg-in99.8%
*-commutative99.8%
distribute-lft-neg-in99.8%
metadata-eval99.8%
rem-square-sqrt0.0%
unpow20.0%
*-commutative0.0%
unpow30.0%
unpow20.0%
associate-*r*0.0%
distribute-rgt-out--0.0%
unpow20.0%
rem-square-sqrt99.8%
Simplified99.8%
unpow299.8%
associate-*r*99.8%
fma-neg99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (x) :precision binary64 (- (* 0.954929658551372 x) (* 0.12900613773279798 (* x (* x x)))))
double code(double x) {
return (0.954929658551372 * x) - (0.12900613773279798 * (x * (x * x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (0.954929658551372d0 * x) - (0.12900613773279798d0 * (x * (x * x)))
end function
public static double code(double x) {
return (0.954929658551372 * x) - (0.12900613773279798 * (x * (x * x)));
}
def code(x): return (0.954929658551372 * x) - (0.12900613773279798 * (x * (x * x)))
function code(x) return Float64(Float64(0.954929658551372 * x) - Float64(0.12900613773279798 * Float64(x * Float64(x * x)))) end
function tmp = code(x) tmp = (0.954929658551372 * x) - (0.12900613773279798 * (x * (x * x))); end
code[x_] := N[(N[(0.954929658551372 * x), $MachinePrecision] - N[(0.12900613773279798 * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.954929658551372 \cdot x - 0.12900613773279798 \cdot \left(x \cdot \left(x \cdot x\right)\right)
\end{array}
Initial program 99.7%
Final simplification99.7%
(FPCore (x) :precision binary64 (* 0.954929658551372 x))
double code(double x) {
return 0.954929658551372 * x;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.954929658551372d0 * x
end function
public static double code(double x) {
return 0.954929658551372 * x;
}
def code(x): return 0.954929658551372 * x
function code(x) return Float64(0.954929658551372 * x) end
function tmp = code(x) tmp = 0.954929658551372 * x; end
code[x_] := N[(0.954929658551372 * x), $MachinePrecision]
\begin{array}{l}
\\
0.954929658551372 \cdot x
\end{array}
Initial program 99.7%
sub-neg99.7%
*-commutative99.7%
distribute-rgt-neg-in99.7%
unpow399.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 47.7%
*-commutative47.7%
Simplified47.7%
Final simplification47.7%
herbie shell --seed 2024048
(FPCore (x)
:name "Rosa's Benchmark"
:precision binary64
(- (* 0.954929658551372 x) (* 0.12900613773279798 (* (* x x) x))))