(FPCore (lo hi x) :precision binary64 (/ (- x lo) (- hi lo)))
(FPCore (lo hi x) :precision binary64 (sqrt (pow (/ lo hi) 6.0)))
double code(double lo, double hi, double x) {
return (x - lo) / (hi - lo);
}
double code(double lo, double hi, double x) {
return sqrt(pow((lo / hi), 6.0));
}
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 = sqrt(((lo / hi) ** 6.0d0))
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 Math.sqrt(Math.pow((lo / hi), 6.0));
}
def code(lo, hi, x): return (x - lo) / (hi - lo)
def code(lo, hi, x): return math.sqrt(math.pow((lo / hi), 6.0))
function code(lo, hi, x) return Float64(Float64(x - lo) / Float64(hi - lo)) end
function code(lo, hi, x) return sqrt((Float64(lo / hi) ^ 6.0)) end
function tmp = code(lo, hi, x) tmp = (x - lo) / (hi - lo); end
function tmp = code(lo, hi, x) tmp = sqrt(((lo / hi) ^ 6.0)); end
code[lo_, hi_, x_] := N[(N[(x - lo), $MachinePrecision] / N[(hi - lo), $MachinePrecision]), $MachinePrecision]
code[lo_, hi_, x_] := N[Sqrt[N[Power[N[(lo / hi), $MachinePrecision], 6.0], $MachinePrecision]], $MachinePrecision]
\frac{x - lo}{hi - lo}
\sqrt{{\left(\frac{lo}{hi}\right)}^{6}}
Results
Initial program 62.0
Taylor expanded in hi around inf 64.0
Simplified51.9
Applied egg-rr51.9
Taylor expanded in lo around -inf 51.6
Simplified51.6
Applied egg-rr51.6
Final simplification51.6
herbie shell --seed 2022192
(FPCore (lo hi x)
:name "(/ (- x lo) (- hi lo))"
:precision binary64
:pre (and (< lo -1e+308) (> hi 1e+308))
(/ (- x lo) (- hi lo)))