(FPCore (x y) :precision binary64 (/ x (* y y)))
(FPCore (x y) :precision binary64 (/ (/ x y) y))
double code(double x, double y) {
return x / (y * y);
}
double code(double x, double y) {
return (x / y) / y;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x / (y * y)
end function
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = (x / y) / y
end function
public static double code(double x, double y) {
return x / (y * y);
}
public static double code(double x, double y) {
return (x / y) / y;
}
def code(x, y): return x / (y * y)
def code(x, y): return (x / y) / y
function code(x, y) return Float64(x / Float64(y * y)) end
function code(x, y) return Float64(Float64(x / y) / y) end
function tmp = code(x, y) tmp = x / (y * y); end
function tmp = code(x, y) tmp = (x / y) / y; end
code[x_, y_] := N[(x / N[(y * y), $MachinePrecision]), $MachinePrecision]
code[x_, y_] := N[(N[(x / y), $MachinePrecision] / y), $MachinePrecision]
\frac{x}{y \cdot y}
\frac{\frac{x}{y}}{y}




Bits error versus x




Bits error versus y
Results
| Original | 10.5 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 10.5
Taylor expanded in x around 0 10.5
Simplified0.2
Final simplification0.2
herbie shell --seed 2022150
(FPCore (x y)
:name "Physics.ForceLayout:coulombForce from force-layout-0.4.0.2"
:precision binary64
:herbie-target
(/ (/ x y) y)
(/ x (* y y)))