
(FPCore (x y) :precision binary64 (- (* x x) (* y y)))
double code(double x, double y) {
return (x * x) - (y * y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x * x) - (y * y)
end function
public static double code(double x, double y) {
return (x * x) - (y * y);
}
def code(x, y): return (x * x) - (y * y)
function code(x, y) return Float64(Float64(x * x) - Float64(y * y)) end
function tmp = code(x, y) tmp = (x * x) - (y * y); end
code[x_, y_] := N[(N[(x * x), $MachinePrecision] - N[(y * y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot x - y \cdot y
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (- (* x x) (* y y)))
double code(double x, double y) {
return (x * x) - (y * y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x * x) - (y * y)
end function
public static double code(double x, double y) {
return (x * x) - (y * y);
}
def code(x, y): return (x * x) - (y * y)
function code(x, y) return Float64(Float64(x * x) - Float64(y * y)) end
function tmp = code(x, y) tmp = (x * x) - (y * y); end
code[x_, y_] := N[(N[(x * x), $MachinePrecision] - N[(y * y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot x - y \cdot y
\end{array}
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (if (<= (* x_m x_m) 3e+289) (fma x_m x_m (* y_m (- y_m))) (* x_m (+ x_m (* y_m -2.0)))))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
double tmp;
if ((x_m * x_m) <= 3e+289) {
tmp = fma(x_m, x_m, (y_m * -y_m));
} else {
tmp = x_m * (x_m + (y_m * -2.0));
}
return tmp;
}
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) tmp = 0.0 if (Float64(x_m * x_m) <= 3e+289) tmp = fma(x_m, x_m, Float64(y_m * Float64(-y_m))); else tmp = Float64(x_m * Float64(x_m + Float64(y_m * -2.0))); end return tmp end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := If[LessEqual[N[(x$95$m * x$95$m), $MachinePrecision], 3e+289], N[(x$95$m * x$95$m + N[(y$95$m * (-y$95$m)), $MachinePrecision]), $MachinePrecision], N[(x$95$m * N[(x$95$m + N[(y$95$m * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
\begin{array}{l}
\mathbf{if}\;x_m \cdot x_m \leq 3 \cdot 10^{+289}:\\
\;\;\;\;\mathsf{fma}\left(x_m, x_m, y_m \cdot \left(-y_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x_m \cdot \left(x_m + y_m \cdot -2\right)\\
\end{array}
\end{array}
if (*.f64 x x) < 3.0000000000000002e289Initial program 100.0%
sqr-neg100.0%
cancel-sign-sub100.0%
fma-def100.0%
Simplified100.0%
if 3.0000000000000002e289 < (*.f64 x x) Initial program 78.6%
difference-of-squares100.0%
add-sqr-sqrt40.0%
sqrt-prod91.4%
sqr-neg91.4%
sqrt-unprod54.3%
add-sqr-sqrt90.0%
sub-neg90.0%
pow190.0%
pow190.0%
pow-prod-up90.0%
add-sqr-sqrt48.6%
add-sqr-sqrt17.1%
difference-of-squares17.1%
metadata-eval17.1%
unpow-prod-down17.1%
Applied egg-rr17.1%
unpow217.1%
unpow217.1%
unswap-sqr17.1%
difference-of-squares17.1%
unpow1/217.1%
unpow1/217.1%
pow-sqr17.1%
metadata-eval17.1%
unpow117.1%
unpow1/217.1%
unpow1/217.1%
pow-sqr17.1%
metadata-eval17.1%
unpow117.1%
difference-of-squares17.1%
unpow1/217.1%
unpow1/217.1%
pow-sqr35.7%
metadata-eval35.7%
unpow135.7%
Simplified90.0%
Taylor expanded in x around inf 78.6%
*-commutative78.6%
associate-*l*78.6%
unpow278.6%
distribute-lft-out97.1%
Simplified97.1%
Final simplification99.2%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (if (<= (* x_m x_m) 3e+289) (- (* x_m x_m) (* y_m y_m)) (* x_m (+ x_m (* y_m -2.0)))))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
double tmp;
if ((x_m * x_m) <= 3e+289) {
tmp = (x_m * x_m) - (y_m * y_m);
} else {
tmp = x_m * (x_m + (y_m * -2.0));
}
return tmp;
}
x_m = abs(x)
y_m = abs(y)
real(8) function code(x_m, y_m)
real(8), intent (in) :: x_m
real(8), intent (in) :: y_m
real(8) :: tmp
if ((x_m * x_m) <= 3d+289) then
tmp = (x_m * x_m) - (y_m * y_m)
else
tmp = x_m * (x_m + (y_m * (-2.0d0)))
end if
code = tmp
end function
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
double tmp;
if ((x_m * x_m) <= 3e+289) {
tmp = (x_m * x_m) - (y_m * y_m);
} else {
tmp = x_m * (x_m + (y_m * -2.0));
}
return tmp;
}
x_m = math.fabs(x) y_m = math.fabs(y) def code(x_m, y_m): tmp = 0 if (x_m * x_m) <= 3e+289: tmp = (x_m * x_m) - (y_m * y_m) else: tmp = x_m * (x_m + (y_m * -2.0)) return tmp
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) tmp = 0.0 if (Float64(x_m * x_m) <= 3e+289) tmp = Float64(Float64(x_m * x_m) - Float64(y_m * y_m)); else tmp = Float64(x_m * Float64(x_m + Float64(y_m * -2.0))); end return tmp end
x_m = abs(x); y_m = abs(y); function tmp_2 = code(x_m, y_m) tmp = 0.0; if ((x_m * x_m) <= 3e+289) tmp = (x_m * x_m) - (y_m * y_m); else tmp = x_m * (x_m + (y_m * -2.0)); end tmp_2 = tmp; end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := If[LessEqual[N[(x$95$m * x$95$m), $MachinePrecision], 3e+289], N[(N[(x$95$m * x$95$m), $MachinePrecision] - N[(y$95$m * y$95$m), $MachinePrecision]), $MachinePrecision], N[(x$95$m * N[(x$95$m + N[(y$95$m * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
\begin{array}{l}
\mathbf{if}\;x_m \cdot x_m \leq 3 \cdot 10^{+289}:\\
\;\;\;\;x_m \cdot x_m - y_m \cdot y_m\\
\mathbf{else}:\\
\;\;\;\;x_m \cdot \left(x_m + y_m \cdot -2\right)\\
\end{array}
\end{array}
if (*.f64 x x) < 3.0000000000000002e289Initial program 100.0%
if 3.0000000000000002e289 < (*.f64 x x) Initial program 78.6%
difference-of-squares100.0%
add-sqr-sqrt40.0%
sqrt-prod91.4%
sqr-neg91.4%
sqrt-unprod54.3%
add-sqr-sqrt90.0%
sub-neg90.0%
pow190.0%
pow190.0%
pow-prod-up90.0%
add-sqr-sqrt48.6%
add-sqr-sqrt17.1%
difference-of-squares17.1%
metadata-eval17.1%
unpow-prod-down17.1%
Applied egg-rr17.1%
unpow217.1%
unpow217.1%
unswap-sqr17.1%
difference-of-squares17.1%
unpow1/217.1%
unpow1/217.1%
pow-sqr17.1%
metadata-eval17.1%
unpow117.1%
unpow1/217.1%
unpow1/217.1%
pow-sqr17.1%
metadata-eval17.1%
unpow117.1%
difference-of-squares17.1%
unpow1/217.1%
unpow1/217.1%
pow-sqr35.7%
metadata-eval35.7%
unpow135.7%
Simplified90.0%
Taylor expanded in x around inf 78.6%
*-commutative78.6%
associate-*l*78.6%
unpow278.6%
distribute-lft-out97.1%
Simplified97.1%
Final simplification99.2%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (* x_m (+ x_m (* y_m -2.0))))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
return x_m * (x_m + (y_m * -2.0));
}
x_m = abs(x)
y_m = abs(y)
real(8) function code(x_m, y_m)
real(8), intent (in) :: x_m
real(8), intent (in) :: y_m
code = x_m * (x_m + (y_m * (-2.0d0)))
end function
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
return x_m * (x_m + (y_m * -2.0));
}
x_m = math.fabs(x) y_m = math.fabs(y) def code(x_m, y_m): return x_m * (x_m + (y_m * -2.0))
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) return Float64(x_m * Float64(x_m + Float64(y_m * -2.0))) end
x_m = abs(x); y_m = abs(y); function tmp = code(x_m, y_m) tmp = x_m * (x_m + (y_m * -2.0)); end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := N[(x$95$m * N[(x$95$m + N[(y$95$m * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
x_m \cdot \left(x_m + y_m \cdot -2\right)
\end{array}
Initial program 94.1%
difference-of-squares100.0%
add-sqr-sqrt48.4%
sqrt-prod76.6%
sqr-neg76.6%
sqrt-unprod28.9%
add-sqr-sqrt52.1%
sub-neg52.1%
pow152.1%
pow152.1%
pow-prod-up52.1%
add-sqr-sqrt27.6%
add-sqr-sqrt12.1%
difference-of-squares12.1%
metadata-eval12.1%
unpow-prod-down12.1%
Applied egg-rr12.1%
unpow212.1%
unpow212.1%
unswap-sqr12.1%
difference-of-squares12.1%
unpow1/212.1%
unpow1/212.1%
pow-sqr12.1%
metadata-eval12.1%
unpow112.1%
unpow1/212.1%
unpow1/212.1%
pow-sqr12.1%
metadata-eval12.1%
unpow112.1%
difference-of-squares12.1%
unpow1/212.1%
unpow1/212.1%
pow-sqr23.2%
metadata-eval23.2%
unpow123.2%
Simplified52.1%
Taylor expanded in x around inf 50.9%
*-commutative50.9%
associate-*l*50.9%
unpow250.9%
distribute-lft-out55.9%
Simplified55.9%
Final simplification55.9%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (* y_m (* x_m -2.0)))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
return y_m * (x_m * -2.0);
}
x_m = abs(x)
y_m = abs(y)
real(8) function code(x_m, y_m)
real(8), intent (in) :: x_m
real(8), intent (in) :: y_m
code = y_m * (x_m * (-2.0d0))
end function
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
return y_m * (x_m * -2.0);
}
x_m = math.fabs(x) y_m = math.fabs(y) def code(x_m, y_m): return y_m * (x_m * -2.0)
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) return Float64(y_m * Float64(x_m * -2.0)) end
x_m = abs(x); y_m = abs(y); function tmp = code(x_m, y_m) tmp = y_m * (x_m * -2.0); end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := N[(y$95$m * N[(x$95$m * -2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
y_m \cdot \left(x_m \cdot -2\right)
\end{array}
Initial program 94.1%
difference-of-squares100.0%
add-sqr-sqrt48.4%
sqrt-prod76.6%
sqr-neg76.6%
sqrt-unprod28.9%
add-sqr-sqrt52.1%
sub-neg52.1%
pow152.1%
pow152.1%
pow-prod-up52.1%
add-sqr-sqrt27.6%
add-sqr-sqrt12.1%
difference-of-squares12.1%
metadata-eval12.1%
unpow-prod-down12.1%
Applied egg-rr12.1%
unpow212.1%
unpow212.1%
unswap-sqr12.1%
difference-of-squares12.1%
unpow1/212.1%
unpow1/212.1%
pow-sqr12.1%
metadata-eval12.1%
unpow112.1%
unpow1/212.1%
unpow1/212.1%
pow-sqr12.1%
metadata-eval12.1%
unpow112.1%
difference-of-squares12.1%
unpow1/212.1%
unpow1/212.1%
pow-sqr23.2%
metadata-eval23.2%
unpow123.2%
Simplified52.1%
Taylor expanded in x around inf 50.9%
*-commutative50.9%
associate-*l*50.9%
unpow250.9%
distribute-lft-out55.9%
Simplified55.9%
Taylor expanded in x around 0 15.6%
associate-*r*15.6%
*-commutative15.6%
Simplified15.6%
Final simplification15.6%
herbie shell --seed 2024024
(FPCore (x y)
:name "Examples.Basics.BasicTests:f2 from sbv-4.4"
:precision binary64
(- (* x x) (* y y)))