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




Bits error versus x




Bits error versus y
Results
| Original | 8.5 |
|---|---|
| Target | 0.0 |
| Herbie | 0.0 |
Initial program 8.5
Simplified0.0
Applied egg-rr0.1
Applied egg-rr0.0
Final simplification0.0
herbie shell --seed 2022151
(FPCore (x y)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, B"
:precision binary64
:herbie-target
(if (< y -3693.8482788297247) (- (/ x (* y y)) (- (/ x y) x)) (if (< y 6799310503.41891) (/ (* x y) (+ y 1.0)) (- (/ x (* y y)) (- (/ x y) x))))
(/ (* x y) (+ y 1.0)))