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




Bits error versus x




Bits error versus y
Results
| Original | 15.0 |
|---|---|
| Target | 0.0 |
| Herbie | 0.0 |
Initial program 15.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2022150
(FPCore (x y)
:name "Linear.Projection:inversePerspective from linear-1.19.1.3, B"
:precision binary64
:herbie-target
(- (/ 0.5 y) (/ 0.5 x))
(/ (- x y) (* (* x 2.0) y)))