| Alternative 1 | |
|---|---|
| Error | 52.0 |
| Cost | 256 |
\[\frac{-lo}{hi}
\]
Error
(FPCore (lo hi x) :precision binary64 (/ (- x lo) (- hi lo)))
(FPCore (lo hi x) :precision binary64 (/ (/ hi lo) (/ lo hi)))
double code(double lo, double hi, double x) {
return (x - lo) / (hi - lo);
}
double code(double lo, double hi, double x) {
return (hi / lo) / (lo / hi);
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = (x - lo) / (hi - lo)
end function
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = (hi / lo) / (lo / hi)
end function
public static double code(double lo, double hi, double x) {
return (x - lo) / (hi - lo);
}
public static double code(double lo, double hi, double x) {
return (hi / lo) / (lo / hi);
}
def code(lo, hi, x): return (x - lo) / (hi - lo)
def code(lo, hi, x): return (hi / lo) / (lo / hi)
function code(lo, hi, x) return Float64(Float64(x - lo) / Float64(hi - lo)) end
function code(lo, hi, x) return Float64(Float64(hi / lo) / Float64(lo / hi)) end
function tmp = code(lo, hi, x) tmp = (x - lo) / (hi - lo); end
function tmp = code(lo, hi, x) tmp = (hi / lo) / (lo / hi); end
code[lo_, hi_, x_] := N[(N[(x - lo), $MachinePrecision] / N[(hi - lo), $MachinePrecision]), $MachinePrecision]
code[lo_, hi_, x_] := N[(N[(hi / lo), $MachinePrecision] / N[(lo / hi), $MachinePrecision]), $MachinePrecision]
\frac{x - lo}{hi - lo}
\frac{\frac{hi}{lo}}{\frac{lo}{hi}}
Results
Initial program 62.0
Taylor expanded in lo around inf 64.0
Simplified51.9
Applied egg-rr51.9
Taylor expanded in hi around inf 64.0
Simplified51.6
Applied egg-rr51.6
Final simplification51.6
| Alternative 1 | |
|---|---|
| Error | 52.0 |
| Cost | 256 |
| Alternative 2 | |
|---|---|
| Error | 52.1 |
| Cost | 64 |
herbie shell --seed 2022298
(FPCore (lo hi x)
:name "(/ (- x lo) (- hi lo))"
:precision binary64
:pre (and (< lo -1e+308) (> hi 1e+308))
(/ (- x lo) (- hi lo)))