Average Error: 62.0 → 51.9
Time: 3.7s
Precision: binary64
\[lo < -1 \cdot 10^{+308} \land hi > 10^{+308}\]
\[\frac{x - lo}{hi - lo} \]
↓
\[\frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(\left(-\sqrt[3]{lo} \cdot \sqrt[3]{lo}\right) \cdot \frac{\sqrt[3]{lo}}{hi}\right) - \frac{lo}{hi}\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]
\frac{x - lo}{hi - lo}
↓
\frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(\left(-\sqrt[3]{lo} \cdot \sqrt[3]{lo}\right) \cdot \frac{\sqrt[3]{lo}}{hi}\right) - \frac{lo}{hi}\right) - {\left(\frac{lo}{hi}\right)}^{3}\right)

(FPCore (lo hi x) :precision binary64 (/ (- x lo) (- hi lo)))
↓
(FPCore (lo hi x)
 :precision binary64
 (+
  (/ x hi)
  (-
   (fma
    (/ x hi)
    (* (/ lo hi) (/ lo hi))
    (-
     (* (/ lo hi) (* (- (* (cbrt lo) (cbrt lo))) (/ (cbrt lo) hi)))
     (/ lo hi)))
   (pow (/ lo hi) 3.0))))
double code(double lo, double hi, double x) {
	return (x - lo) / (hi - lo);
}
↓
double code(double lo, double hi, double x) {
	return (x / hi) + (fma((x / hi), ((lo / hi) * (lo / hi)), (((lo / hi) * (-(cbrt(lo) * cbrt(lo)) * (cbrt(lo) / hi))) - (lo / hi))) - pow((lo / hi), 3.0));
}

Error

Bits error versus lo

Bits error versus hi

Bits error versus x

Derivation

  1. Initial program 62.0

    \[\frac{x - lo}{hi - lo} \]

  2. Taylor expanded in hi around inf 64.0

    \[\leadsto \color{blue}{\left(\frac{x}{hi} + \left(\frac{{lo}^{2} \cdot x}{{hi}^{3}} + \frac{lo \cdot x}{{hi}^{2}}\right)\right) - \left(\frac{{lo}^{3}}{{hi}^{3}} + \left(\frac{lo}{hi} + \frac{{lo}^{2}}{{hi}^{2}}\right)\right)} \]

  3. Simplified51.9

    \[\leadsto \color{blue}{\frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \left(\frac{lo}{hi} + 1\right) \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(-1 - \frac{lo}{hi}\right)\right) - {\left(\frac{lo}{hi}\right)}^{3}\right)} \]

  4. Taylor expanded in lo around inf 51.9

    \[\leadsto \frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \color{blue}{\frac{lo}{hi}} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(-1 - \frac{lo}{hi}\right)\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]

  5. Applied *-un-lft-identity_binary6451.9

    \[\leadsto \frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(-1 - \frac{lo}{\color{blue}{1 \cdot hi}}\right)\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]

  6. Applied add-cube-cbrt_binary6451.9

    \[\leadsto \frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(-1 - \frac{\color{blue}{\left(\sqrt[3]{lo} \cdot \sqrt[3]{lo}\right) \cdot \sqrt[3]{lo}}}{1 \cdot hi}\right)\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]

  7. Applied times-frac_binary6451.9

    \[\leadsto \frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(-1 - \color{blue}{\frac{\sqrt[3]{lo} \cdot \sqrt[3]{lo}}{1} \cdot \frac{\sqrt[3]{lo}}{hi}}\right)\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]

  8. Applied cancel-sign-sub-inv_binary6451.9

    \[\leadsto \frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \color{blue}{\left(-1 + \left(-\frac{\sqrt[3]{lo} \cdot \sqrt[3]{lo}}{1}\right) \cdot \frac{\sqrt[3]{lo}}{hi}\right)}\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]

  9. Applied distribute-rgt-in_binary6451.9

    \[\leadsto \frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \color{blue}{-1 \cdot \frac{lo}{hi} + \left(\left(-\frac{\sqrt[3]{lo} \cdot \sqrt[3]{lo}}{1}\right) \cdot \frac{\sqrt[3]{lo}}{hi}\right) \cdot \frac{lo}{hi}}\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]

  10. Final simplification51.9

    \[\leadsto \frac{x}{hi} + \left(\mathsf{fma}\left(\frac{x}{hi}, \frac{lo}{hi} \cdot \frac{lo}{hi}, \frac{lo}{hi} \cdot \left(\left(-\sqrt[3]{lo} \cdot \sqrt[3]{lo}\right) \cdot \frac{\sqrt[3]{lo}}{hi}\right) - \frac{lo}{hi}\right) - {\left(\frac{lo}{hi}\right)}^{3}\right) \]

Reproduce

herbie shell --seed 2021344 
(FPCore (lo hi x)
  :name "(/ (- x lo) (- hi lo))"
  :precision binary64
  :pre (and (< lo -1e+308) (> hi 1e+308))
  (/ (- x lo) (- hi lo)))