
(FPCore (x) :precision binary64 (* (* x x) (- 3.0 (* x 2.0))))
double code(double x) {
return (x * x) * (3.0 - (x * 2.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (x * x) * (3.0d0 - (x * 2.0d0))
end function
public static double code(double x) {
return (x * x) * (3.0 - (x * 2.0));
}
def code(x): return (x * x) * (3.0 - (x * 2.0))
function code(x) return Float64(Float64(x * x) * Float64(3.0 - Float64(x * 2.0))) end
function tmp = code(x) tmp = (x * x) * (3.0 - (x * 2.0)); end
code[x_] := N[(N[(x * x), $MachinePrecision] * N[(3.0 - N[(x * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot x\right) \cdot \left(3 - x \cdot 2\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (* (* x x) (- 3.0 (* x 2.0))))
double code(double x) {
return (x * x) * (3.0 - (x * 2.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (x * x) * (3.0d0 - (x * 2.0d0))
end function
public static double code(double x) {
return (x * x) * (3.0 - (x * 2.0));
}
def code(x): return (x * x) * (3.0 - (x * 2.0))
function code(x) return Float64(Float64(x * x) * Float64(3.0 - Float64(x * 2.0))) end
function tmp = code(x) tmp = (x * x) * (3.0 - (x * 2.0)); end
code[x_] := N[(N[(x * x), $MachinePrecision] * N[(3.0 - N[(x * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot x\right) \cdot \left(3 - x \cdot 2\right)
\end{array}
(FPCore (x) :precision binary64 (* x (* x (+ 3.0 (* x -2.0)))))
double code(double x) {
return x * (x * (3.0 + (x * -2.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (x * (3.0d0 + (x * (-2.0d0))))
end function
public static double code(double x) {
return x * (x * (3.0 + (x * -2.0)));
}
def code(x): return x * (x * (3.0 + (x * -2.0)))
function code(x) return Float64(x * Float64(x * Float64(3.0 + Float64(x * -2.0)))) end
function tmp = code(x) tmp = x * (x * (3.0 + (x * -2.0))); end
code[x_] := N[(x * N[(x * N[(3.0 + N[(x * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(x \cdot \left(3 + x \cdot -2\right)\right)
\end{array}
Initial program 99.8%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x) :precision binary64 (if (or (<= x -1.5) (not (<= x 1.5))) (* x (* x (* x -2.0))) (* x (* x 3.0))))
double code(double x) {
double tmp;
if ((x <= -1.5) || !(x <= 1.5)) {
tmp = x * (x * (x * -2.0));
} else {
tmp = x * (x * 3.0);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if ((x <= (-1.5d0)) .or. (.not. (x <= 1.5d0))) then
tmp = x * (x * (x * (-2.0d0)))
else
tmp = x * (x * 3.0d0)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if ((x <= -1.5) || !(x <= 1.5)) {
tmp = x * (x * (x * -2.0));
} else {
tmp = x * (x * 3.0);
}
return tmp;
}
def code(x): tmp = 0 if (x <= -1.5) or not (x <= 1.5): tmp = x * (x * (x * -2.0)) else: tmp = x * (x * 3.0) return tmp
function code(x) tmp = 0.0 if ((x <= -1.5) || !(x <= 1.5)) tmp = Float64(x * Float64(x * Float64(x * -2.0))); else tmp = Float64(x * Float64(x * 3.0)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if ((x <= -1.5) || ~((x <= 1.5))) tmp = x * (x * (x * -2.0)); else tmp = x * (x * 3.0); end tmp_2 = tmp; end
code[x_] := If[Or[LessEqual[x, -1.5], N[Not[LessEqual[x, 1.5]], $MachinePrecision]], N[(x * N[(x * N[(x * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * 3.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.5 \lor \neg \left(x \leq 1.5\right):\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot -2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot 3\right)\\
\end{array}
\end{array}
if x < -1.5 or 1.5 < x Initial program 99.8%
associate-*l*99.8%
sub-neg99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 96.6%
*-commutative96.6%
unpow296.6%
associate-*r*96.6%
Simplified96.6%
if -1.5 < x < 1.5Initial program 99.8%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 98.8%
Final simplification97.6%
(FPCore (x) :precision binary64 (if (or (<= x -1.12) (not (<= x 1.5))) (* x (* x (* x -2.0))) (/ x (+ (/ 0.3333333333333333 x) 0.2222222222222222))))
double code(double x) {
double tmp;
if ((x <= -1.12) || !(x <= 1.5)) {
tmp = x * (x * (x * -2.0));
} else {
tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if ((x <= (-1.12d0)) .or. (.not. (x <= 1.5d0))) then
tmp = x * (x * (x * (-2.0d0)))
else
tmp = x / ((0.3333333333333333d0 / x) + 0.2222222222222222d0)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if ((x <= -1.12) || !(x <= 1.5)) {
tmp = x * (x * (x * -2.0));
} else {
tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222);
}
return tmp;
}
def code(x): tmp = 0 if (x <= -1.12) or not (x <= 1.5): tmp = x * (x * (x * -2.0)) else: tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222) return tmp
function code(x) tmp = 0.0 if ((x <= -1.12) || !(x <= 1.5)) tmp = Float64(x * Float64(x * Float64(x * -2.0))); else tmp = Float64(x / Float64(Float64(0.3333333333333333 / x) + 0.2222222222222222)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if ((x <= -1.12) || ~((x <= 1.5))) tmp = x * (x * (x * -2.0)); else tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222); end tmp_2 = tmp; end
code[x_] := If[Or[LessEqual[x, -1.12], N[Not[LessEqual[x, 1.5]], $MachinePrecision]], N[(x * N[(x * N[(x * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(N[(0.3333333333333333 / x), $MachinePrecision] + 0.2222222222222222), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.12 \lor \neg \left(x \leq 1.5\right):\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot -2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{0.3333333333333333}{x} + 0.2222222222222222}\\
\end{array}
\end{array}
if x < -1.1200000000000001 or 1.5 < x Initial program 99.8%
associate-*l*99.8%
sub-neg99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 96.6%
*-commutative96.6%
unpow296.6%
associate-*r*96.6%
Simplified96.6%
if -1.1200000000000001 < x < 1.5Initial program 99.8%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
distribute-lft-in99.9%
+-commutative99.9%
Applied egg-rr99.9%
distribute-lft-out99.9%
fma-udef99.9%
associate-*r*99.8%
remove-double-div99.8%
/-rgt-identity99.8%
associate-/l*99.7%
frac-2neg99.7%
clear-num99.7%
frac-times99.8%
fma-udef99.8%
distribute-neg-in99.8%
+-commutative99.8%
metadata-eval99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
associate-/l*99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 99.3%
+-commutative99.3%
associate-*r/99.4%
metadata-eval99.4%
Simplified99.4%
Final simplification97.9%
(FPCore (x)
:precision binary64
(if (<= x -1.12)
(/ x (/ -0.5 (* x x)))
(if (<= x 1.5)
(/ x (+ (/ 0.3333333333333333 x) 0.2222222222222222))
(* x (* x (* x -2.0))))))
double code(double x) {
double tmp;
if (x <= -1.12) {
tmp = x / (-0.5 / (x * x));
} else if (x <= 1.5) {
tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222);
} else {
tmp = x * (x * (x * -2.0));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-1.12d0)) then
tmp = x / ((-0.5d0) / (x * x))
else if (x <= 1.5d0) then
tmp = x / ((0.3333333333333333d0 / x) + 0.2222222222222222d0)
else
tmp = x * (x * (x * (-2.0d0)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -1.12) {
tmp = x / (-0.5 / (x * x));
} else if (x <= 1.5) {
tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222);
} else {
tmp = x * (x * (x * -2.0));
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.12: tmp = x / (-0.5 / (x * x)) elif x <= 1.5: tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222) else: tmp = x * (x * (x * -2.0)) return tmp
function code(x) tmp = 0.0 if (x <= -1.12) tmp = Float64(x / Float64(-0.5 / Float64(x * x))); elseif (x <= 1.5) tmp = Float64(x / Float64(Float64(0.3333333333333333 / x) + 0.2222222222222222)); else tmp = Float64(x * Float64(x * Float64(x * -2.0))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.12) tmp = x / (-0.5 / (x * x)); elseif (x <= 1.5) tmp = x / ((0.3333333333333333 / x) + 0.2222222222222222); else tmp = x * (x * (x * -2.0)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -1.12], N[(x / N[(-0.5 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.5], N[(x / N[(N[(0.3333333333333333 / x), $MachinePrecision] + 0.2222222222222222), $MachinePrecision]), $MachinePrecision], N[(x * N[(x * N[(x * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.12:\\
\;\;\;\;\frac{x}{\frac{-0.5}{x \cdot x}}\\
\mathbf{elif}\;x \leq 1.5:\\
\;\;\;\;\frac{x}{\frac{0.3333333333333333}{x} + 0.2222222222222222}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(x \cdot \left(x \cdot -2\right)\right)\\
\end{array}
\end{array}
if x < -1.1200000000000001Initial program 99.9%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
distribute-lft-in99.9%
+-commutative99.9%
Applied egg-rr99.9%
distribute-lft-out99.9%
fma-udef99.9%
associate-*r*99.9%
remove-double-div99.9%
/-rgt-identity99.9%
associate-/l*99.8%
frac-2neg99.8%
clear-num99.9%
frac-times99.8%
fma-udef99.8%
distribute-neg-in99.8%
+-commutative99.8%
metadata-eval99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
associate-/l*100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 97.9%
unpow297.9%
Simplified97.9%
if -1.1200000000000001 < x < 1.5Initial program 99.8%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
distribute-lft-in99.9%
+-commutative99.9%
Applied egg-rr99.9%
distribute-lft-out99.9%
fma-udef99.9%
associate-*r*99.8%
remove-double-div99.8%
/-rgt-identity99.8%
associate-/l*99.7%
frac-2neg99.7%
clear-num99.7%
frac-times99.8%
fma-udef99.8%
distribute-neg-in99.8%
+-commutative99.8%
metadata-eval99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
associate-/l*99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 99.3%
+-commutative99.3%
associate-*r/99.4%
metadata-eval99.4%
Simplified99.4%
if 1.5 < x Initial program 99.8%
associate-*l*99.8%
sub-neg99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 95.5%
*-commutative95.5%
unpow295.5%
associate-*r*95.5%
Simplified95.5%
Final simplification97.9%
(FPCore (x) :precision binary64 (* 3.0 (* x x)))
double code(double x) {
return 3.0 * (x * x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = 3.0d0 * (x * x)
end function
public static double code(double x) {
return 3.0 * (x * x);
}
def code(x): return 3.0 * (x * x)
function code(x) return Float64(3.0 * Float64(x * x)) end
function tmp = code(x) tmp = 3.0 * (x * x); end
code[x_] := N[(3.0 * N[(x * x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
3 \cdot \left(x \cdot x\right)
\end{array}
Initial program 99.8%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 58.5%
unpow258.5%
Simplified58.5%
Final simplification58.5%
(FPCore (x) :precision binary64 (* x (* x 3.0)))
double code(double x) {
return x * (x * 3.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (x * 3.0d0)
end function
public static double code(double x) {
return x * (x * 3.0);
}
def code(x): return x * (x * 3.0)
function code(x) return Float64(x * Float64(x * 3.0)) end
function tmp = code(x) tmp = x * (x * 3.0); end
code[x_] := N[(x * N[(x * 3.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(x \cdot 3\right)
\end{array}
Initial program 99.8%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 58.5%
Final simplification58.5%
(FPCore (x) :precision binary64 (* x 4.5))
double code(double x) {
return x * 4.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * 4.5d0
end function
public static double code(double x) {
return x * 4.5;
}
def code(x): return x * 4.5
function code(x) return Float64(x * 4.5) end
function tmp = code(x) tmp = x * 4.5; end
code[x_] := N[(x * 4.5), $MachinePrecision]
\begin{array}{l}
\\
x \cdot 4.5
\end{array}
Initial program 99.8%
associate-*l*99.9%
sub-neg99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
distribute-lft-in99.8%
+-commutative99.8%
Applied egg-rr99.8%
distribute-lft-out99.9%
fma-udef99.9%
associate-*r*99.8%
remove-double-div99.8%
/-rgt-identity99.8%
associate-/l*99.7%
frac-2neg99.7%
clear-num99.7%
frac-times99.8%
fma-udef99.8%
distribute-neg-in99.8%
+-commutative99.8%
metadata-eval99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
metadata-eval99.8%
Applied egg-rr99.8%
*-commutative99.8%
associate-/l*99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 46.0%
+-commutative46.0%
associate-*r/46.1%
metadata-eval46.1%
Simplified46.1%
Taylor expanded in x around inf 3.1%
*-commutative3.1%
Simplified3.1%
Final simplification3.1%
(FPCore (x) :precision binary64 (* x (* x (- 3.0 (* x 2.0)))))
double code(double x) {
return x * (x * (3.0 - (x * 2.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (x * (3.0d0 - (x * 2.0d0)))
end function
public static double code(double x) {
return x * (x * (3.0 - (x * 2.0)));
}
def code(x): return x * (x * (3.0 - (x * 2.0)))
function code(x) return Float64(x * Float64(x * Float64(3.0 - Float64(x * 2.0)))) end
function tmp = code(x) tmp = x * (x * (3.0 - (x * 2.0))); end
code[x_] := N[(x * N[(x * N[(3.0 - N[(x * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(x \cdot \left(3 - x \cdot 2\right)\right)
\end{array}
herbie shell --seed 2023297
(FPCore (x)
:name "Data.Spline.Key:interpolateKeys from smoothie-0.4.0.2"
:precision binary64
:herbie-target
(* x (* x (- 3.0 (* x 2.0))))
(* (* x x) (- 3.0 (* x 2.0))))