Average Error: 26.9 → 12.9
Time: 26.1s
Precision: binary64
\[\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(\left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{h}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\left(\frac{1}{2} \cdot {\left(\frac{M}{2} \cdot \frac{D}{d}\right)}^{2}\right) \cdot \frac{\sqrt{h}}{\sqrt{\ell}}\right) \cdot \frac{\sqrt{h}}{\sqrt{\ell}}\right)\]
\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(\left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{h}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\left(\frac{1}{2} \cdot {\left(\frac{M}{2} \cdot \frac{D}{d}\right)}^{2}\right) \cdot \frac{\sqrt{h}}{\sqrt{\ell}}\right) \cdot \frac{\sqrt{h}}{\sqrt{\ell}}\right)
double code(double d, double h, double l, double M, double D) {
	return ((double) (((double) (((double) pow(((double) (d / h)), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (d / l)), ((double) (1.0 / 2.0)))))) * ((double) (1.0 - ((double) (((double) (((double) (1.0 / 2.0)) * ((double) pow(((double) (((double) (M * D)) / ((double) (2.0 * d)))), 2.0)))) * ((double) (h / l))))))));
}
↓
double code(double d, double h, double l, double M, double D) {
	return ((double) (((double) (((double) (((double) pow(d, ((double) (1.0 / 2.0)))) * ((double) pow(((double) (1.0 / h)), ((double) (1.0 / 2.0)))))) * ((double) (((double) pow(d, ((double) (1.0 / 2.0)))) * ((double) pow(((double) (1.0 / l)), ((double) (1.0 / 2.0)))))))) * ((double) (1.0 - ((double) (((double) (((double) (((double) (1.0 / 2.0)) * ((double) pow(((double) (((double) (M / 2.0)) * ((double) (D / d)))), 2.0)))) * ((double) (((double) sqrt(h)) / ((double) sqrt(l)))))) * ((double) (((double) sqrt(h)) / ((double) sqrt(l))))))))));
}

Error

Bits error versus d

Bits error versus h

Bits error versus l

Bits error versus M

Bits error versus D

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 26.9

    \[\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 div-inv26.9

    \[\leadsto \left({\color{blue}{\left(d \cdot \frac{1}{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 unpow-prod-down20.8

    \[\leadsto \left(\color{blue}{\left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{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)\]
  5. Using strategy rm

  6. Applied add-sqr-sqrt20.8

    \[\leadsto \left(\left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{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}{\color{blue}{\sqrt{\ell} \cdot \sqrt{\ell}}}\right)\]

  7. Applied add-sqr-sqrt20.8

    \[\leadsto \left(\left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{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{\color{blue}{\sqrt{h} \cdot \sqrt{h}}}{\sqrt{\ell} \cdot \sqrt{\ell}}\right)\]

  8. Applied times-frac20.8

    \[\leadsto \left(\left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{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 \color{blue}{\left(\frac{\sqrt{h}}{\sqrt{\ell}} \cdot \frac{\sqrt{h}}{\sqrt{\ell}}\right)}\right)\]

  9. Applied associate-*r*18.4

    \[\leadsto \left(\left({d}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{1}{h}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\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{\sqrt{h}}{\sqrt{\ell}}\right) \cdot \frac{\sqrt{h}}{\sqrt{\ell}}}\right)\]
  10. Using strategy rm

  11. Applied div-inv18.5

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

  12. Applied unpow-prod-down12.6

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

  14. Applied times-frac12.9

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

  15. Final simplification12.9

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

Reproduce

herbie shell --seed 2020156 
(FPCore (d h l M D)
  :name "Henrywood and Agarwal, Equation (12)"
  :precision binary64
  (* (* (pow (/ d h) (/ 1.0 2.0)) (pow (/ d l) (/ 1.0 2.0))) (- 1.0 (* (* (/ 1.0 2.0) (pow (/ (* M D) (* 2.0 d)) 2.0)) (/ h l)))))