
(FPCore (x y) :precision binary64 (* x (+ y 1.0)))
double code(double x, double y) {
return x * (y + 1.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x * (y + 1.0d0)
end function
public static double code(double x, double y) {
return x * (y + 1.0);
}
def code(x, y): return x * (y + 1.0)
function code(x, y) return Float64(x * Float64(y + 1.0)) end
function tmp = code(x, y) tmp = x * (y + 1.0); end
code[x_, y_] := N[(x * N[(y + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(y + 1\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 2 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (* x (+ y 1.0)))
double code(double x, double y) {
return x * (y + 1.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x * (y + 1.0d0)
end function
public static double code(double x, double y) {
return x * (y + 1.0);
}
def code(x, y): return x * (y + 1.0)
function code(x, y) return Float64(x * Float64(y + 1.0)) end
function tmp = code(x, y) tmp = x * (y + 1.0); end
code[x_, y_] := N[(x * N[(y + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(y + 1\right)
\end{array}
(FPCore (x y) :precision binary64 (* x (+ y 1.0)))
double code(double x, double y) {
return x * (y + 1.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x * (y + 1.0d0)
end function
public static double code(double x, double y) {
return x * (y + 1.0);
}
def code(x, y): return x * (y + 1.0)
function code(x, y) return Float64(x * Float64(y + 1.0)) end
function tmp = code(x, y) tmp = x * (y + 1.0); end
code[x_, y_] := N[(x * N[(y + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(y + 1\right)
\end{array}
Initial program 100.0%
(FPCore (x y) :precision binary64 (+ y 1.0))
double code(double x, double y) {
return y + 1.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = y + 1.0d0
end function
public static double code(double x, double y) {
return y + 1.0;
}
def code(x, y): return y + 1.0
function code(x, y) return Float64(y + 1.0) end
function tmp = code(x, y) tmp = y + 1.0; end
code[x_, y_] := N[(y + 1.0), $MachinePrecision]
\begin{array}{l}
\\
y + 1
\end{array}
Initial program 100.0%
Taylor expanded in y around 0
Applied rewrites3.2%
(FPCore (x y) :precision binary64 (+ x (* x y)))
double code(double x, double y) {
return x + (x * y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + (x * y)
end function
public static double code(double x, double y) {
return x + (x * y);
}
def code(x, y): return x + (x * y)
function code(x, y) return Float64(x + Float64(x * y)) end
function tmp = code(x, y) tmp = x + (x * y); end
code[x_, y_] := N[(x + N[(x * y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + x \cdot y
\end{array}
herbie shell --seed 2024321
(FPCore (x y)
:name "Data.Colour.RGBSpace.HSL:hsl from colour-2.3.3, B"
:precision binary64
:pre (TRUE)
:alt
(! :herbie-platform default (+ x (* x y)))
(* x (+ y 1.0)))