Average Error: 25.3 → 10.8
Time: 4.8m
Precision: 64
\[\left({\left(\frac{d}{h}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
↓
\[\left(1 - \frac{\frac{\frac{M \cdot D}{2 \cdot d}}{\sqrt{2}}}{\frac{\ell}{\frac{\frac{M \cdot D}{2 \cdot d}}{\sqrt{2}}}} \cdot h\right) \cdot \left(\left({\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\frac{1}{2}} \cdot \left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right|\right) \cdot \left({\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\frac{1}{2}} \cdot \left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right|\right)\right)\]
double f(double d, double h, double l, double M, double D) {
        double r47098254 = d;
        double r47098255 = h;
        double r47098256 = r47098254 / r47098255;
        double r47098257 = 1.0;
        double r47098258 = 2.0;
        double r47098259 = r47098257 / r47098258;
        double r47098260 = pow(r47098256, r47098259);
        double r47098261 = l;
        double r47098262 = r47098254 / r47098261;
        double r47098263 = pow(r47098262, r47098259);
        double r47098264 = r47098260 * r47098263;
        double r47098265 = M;
        double r47098266 = D;
        double r47098267 = r47098265 * r47098266;
        double r47098268 = r47098258 * r47098254;
        double r47098269 = r47098267 / r47098268;
        double r47098270 = pow(r47098269, r47098258);
        double r47098271 = r47098259 * r47098270;
        double r47098272 = r47098255 / r47098261;
        double r47098273 = r47098271 * r47098272;
        double r47098274 = r47098257 - r47098273;
        double r47098275 = r47098264 * r47098274;
        return r47098275;
}
↓
double f(double d, double h, double l, double M, double D) {
        double r47098276 = 1.0;
        double r47098277 = M;
        double r47098278 = D;
        double r47098279 = r47098277 * r47098278;
        double r47098280 = 2.0;
        double r47098281 = d;
        double r47098282 = r47098280 * r47098281;
        double r47098283 = r47098279 / r47098282;
        double r47098284 = sqrt(r47098280);
        double r47098285 = r47098283 / r47098284;
        double r47098286 = l;
        double r47098287 = r47098286 / r47098285;
        double r47098288 = r47098285 / r47098287;
        double r47098289 = h;
        double r47098290 = r47098288 * r47098289;
        double r47098291 = r47098276 - r47098290;
        double r47098292 = cbrt(r47098281);
        double r47098293 = cbrt(r47098286);
        double r47098294 = r47098292 / r47098293;
        double r47098295 = 0.5;
        double r47098296 = pow(r47098294, r47098295);
        double r47098297 = fabs(r47098294);
        double r47098298 = r47098296 * r47098297;
        double r47098299 = cbrt(r47098289);
        double r47098300 = r47098292 / r47098299;
        double r47098301 = pow(r47098300, r47098295);
        double r47098302 = fabs(r47098300);
        double r47098303 = r47098301 * r47098302;
        double r47098304 = r47098298 * r47098303;
        double r47098305 = r47098291 * r47098304;
        return r47098305;
}


\left({\left(\frac{d}{h}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)
↓
\left(1 - \frac{\frac{\frac{M \cdot D}{2 \cdot d}}{\sqrt{2}}}{\frac{\ell}{\frac{\frac{M \cdot D}{2 \cdot d}}{\sqrt{2}}}} \cdot h\right) \cdot \left(\left({\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\frac{1}{2}} \cdot \left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right|\right) \cdot \left({\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\frac{1}{2}} \cdot \left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right|\right)\right)

Error

Bits error versus d

Bits error versus h

Bits error versus l

Bits error versus M

Bits error versus D

Derivation

  1. Initial program 25.3

    \[\left({\left(\frac{d}{h}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
  2. Using strategy rm

  3. Applied add-cube-cbrt25.5

    \[\leadsto \left({\left(\frac{d}{\color{blue}{\left(\sqrt[3]{h} \cdot \sqrt[3]{h}\right) \cdot \sqrt[3]{h}}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  4. Applied add-cube-cbrt25.6

    \[\leadsto \left({\left(\frac{\color{blue}{\left(\sqrt[3]{d} \cdot \sqrt[3]{d}\right) \cdot \sqrt[3]{d}}}{\left(\sqrt[3]{h} \cdot \sqrt[3]{h}\right) \cdot \sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  5. Applied times-frac25.6

    \[\leadsto \left({\color{blue}{\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{h} \cdot \sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  6. Applied unpow-prod-down20.6

    \[\leadsto \left(\color{blue}{\left({\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{h} \cdot \sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  7. Simplified20.1

    \[\leadsto \left(\left(\color{blue}{\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right|} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
  8. Using strategy rm

  9. Applied add-cube-cbrt20.2

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot {\left(\frac{d}{\color{blue}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  10. Applied add-cube-cbrt20.3

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot {\left(\frac{\color{blue}{\left(\sqrt[3]{d} \cdot \sqrt[3]{d}\right) \cdot \sqrt[3]{d}}}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  11. Applied times-frac20.3

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot {\color{blue}{\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  12. Applied unpow-prod-down16.3

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \color{blue}{\left({\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

  13. Simplified16.1

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\color{blue}{\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right|} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
  14. Using strategy rm

  15. Applied pow116.1

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \color{blue}{{\left(\frac{h}{\ell}\right)}^{1}}\right)\]

  16. Applied pow116.1

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \color{blue}{{\left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right)}^{1}} \cdot {\left(\frac{h}{\ell}\right)}^{1}\right)\]

  17. Applied pow-prod-down16.1

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \color{blue}{{\left(\left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)}^{1}}\right)\]

  18. Simplified13.1

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - {\color{blue}{\left(\frac{\frac{\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}}{2}}{\ell} \cdot h\right)}}^{1}\right)\]
  19. Using strategy rm

  20. Applied add-sqr-sqrt13.1

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - {\left(\frac{\frac{\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}}{\color{blue}{\sqrt{2} \cdot \sqrt{2}}}}{\ell} \cdot h\right)}^{1}\right)\]

  21. Applied times-frac13.1

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - {\left(\frac{\color{blue}{\frac{\frac{M \cdot D}{d \cdot 2}}{\sqrt{2}} \cdot \frac{\frac{M \cdot D}{d \cdot 2}}{\sqrt{2}}}}{\ell} \cdot h\right)}^{1}\right)\]

  22. Applied associate-/l*10.8

    \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right| \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - {\left(\color{blue}{\frac{\frac{\frac{M \cdot D}{d \cdot 2}}{\sqrt{2}}}{\frac{\ell}{\frac{\frac{M \cdot D}{d \cdot 2}}{\sqrt{2}}}}} \cdot h\right)}^{1}\right)\]

  23. Final simplification10.8

    \[\leadsto \left(1 - \frac{\frac{\frac{M \cdot D}{2 \cdot d}}{\sqrt{2}}}{\frac{\ell}{\frac{\frac{M \cdot D}{2 \cdot d}}{\sqrt{2}}}} \cdot h\right) \cdot \left(\left({\left(\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right)}^{\frac{1}{2}} \cdot \left|\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}\right|\right) \cdot \left({\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\frac{1}{2}} \cdot \left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right|\right)\right)\]

Reproduce

herbie shell --seed 2019101 +o rules:numerics
(FPCore (d h l M D)
  :name "Henrywood and Agarwal, Equation (12)"
  (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))))