
(FPCore (x y) :precision binary64 (+ x (/ (- x y) 2.0)))
double code(double x, double y) {
return x + ((x - y) / 2.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + ((x - y) / 2.0d0)
end function
public static double code(double x, double y) {
return x + ((x - y) / 2.0);
}
def code(x, y): return x + ((x - y) / 2.0)
function code(x, y) return Float64(x + Float64(Float64(x - y) / 2.0)) end
function tmp = code(x, y) tmp = x + ((x - y) / 2.0); end
code[x_, y_] := N[(x + N[(N[(x - y), $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{x - y}{2}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (+ x (/ (- x y) 2.0)))
double code(double x, double y) {
return x + ((x - y) / 2.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + ((x - y) / 2.0d0)
end function
public static double code(double x, double y) {
return x + ((x - y) / 2.0);
}
def code(x, y): return x + ((x - y) / 2.0)
function code(x, y) return Float64(x + Float64(Float64(x - y) / 2.0)) end
function tmp = code(x, y) tmp = x + ((x - y) / 2.0); end
code[x_, y_] := N[(x + N[(N[(x - y), $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{x - y}{2}
\end{array}
(FPCore (x y) :precision binary64 (+ x (/ (- x y) 2.0)))
double code(double x, double y) {
return x + ((x - y) / 2.0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x + ((x - y) / 2.0d0)
end function
public static double code(double x, double y) {
return x + ((x - y) / 2.0);
}
def code(x, y): return x + ((x - y) / 2.0)
function code(x, y) return Float64(x + Float64(Float64(x - y) / 2.0)) end
function tmp = code(x, y) tmp = x + ((x - y) / 2.0); end
code[x_, y_] := N[(x + N[(N[(x - y), $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{x - y}{2}
\end{array}
Initial program 99.8%
(FPCore (x y) :precision binary64 (/ (- x y) 2.0))
double code(double x, double y) {
return (x - y) / 2.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x - y) / 2.0d0
end function
public static double code(double x, double y) {
return (x - y) / 2.0;
}
def code(x, y): return (x - y) / 2.0
function code(x, y) return Float64(Float64(x - y) / 2.0) end
function tmp = code(x, y) tmp = (x - y) / 2.0; end
code[x_, y_] := N[(N[(x - y), $MachinePrecision] / 2.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - y}{2}
\end{array}
Initial program 99.8%
Taylor expanded in x around 0
Applied rewrites59.3%
(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 + \left(x - y\right)
\end{array}
Initial program 99.8%
Taylor expanded in x around 0
Applied rewrites20.0%
(FPCore (x y) :precision binary64 (- x y))
double code(double x, double y) {
return x - y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x - y
end function
public static double code(double x, double y) {
return x - y;
}
def code(x, y): return x - y
function code(x, y) return Float64(x - y) end
function tmp = code(x, y) tmp = x - y; end
code[x_, y_] := N[(x - y), $MachinePrecision]
\begin{array}{l}
\\
x - y
\end{array}
Initial program 99.8%
Taylor expanded in x around 0
Applied rewrites59.3%
Taylor expanded in x around 0
Applied rewrites19.3%
(FPCore (x y) :precision binary64 (- (* 1.5 x) (* 0.5 y)))
double code(double x, double y) {
return (1.5 * x) - (0.5 * y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (1.5d0 * x) - (0.5d0 * y)
end function
public static double code(double x, double y) {
return (1.5 * x) - (0.5 * y);
}
def code(x, y): return (1.5 * x) - (0.5 * y)
function code(x, y) return Float64(Float64(1.5 * x) - Float64(0.5 * y)) end
function tmp = code(x, y) tmp = (1.5 * x) - (0.5 * y); end
code[x_, y_] := N[(N[(1.5 * x), $MachinePrecision] - N[(0.5 * y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1.5 \cdot x - 0.5 \cdot y
\end{array}
herbie shell --seed 2024321
(FPCore (x y)
:name "Graphics.Rendering.Chart.Axis.Types:hBufferRect from Chart-1.5.3"
:precision binary64
:pre (TRUE)
:alt
(! :herbie-platform default (- (* 3/2 x) (* 1/2 y)))
(+ x (/ (- x y) 2.0)))