Average Error: 1.0 → 0.0
Time: 5.0s
Precision: 64
\[\frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 - v \cdot v\right)\right) \cdot \sqrt{2 - 6 \cdot \left(v \cdot v\right)}}\]
↓
\[\frac{\frac{4}{3 \cdot \pi}}{\mathsf{fma}\left(v, v, 1\right) \cdot \left(\left(1 - v \cdot v\right) \cdot \sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}\right)} \cdot \left(\left(1 + v \cdot v\right) \cdot \sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}\right)\]
\frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 - v \cdot v\right)\right) \cdot \sqrt{2 - 6 \cdot \left(v \cdot v\right)}}
↓
\frac{\frac{4}{3 \cdot \pi}}{\mathsf{fma}\left(v, v, 1\right) \cdot \left(\left(1 - v \cdot v\right) \cdot \sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}\right)} \cdot \left(\left(1 + v \cdot v\right) \cdot \sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}\right)
double f(double v) {
        double r196435 = 4.0;
        double r196436 = 3.0;
        double r196437 = atan2(1.0, 0.0);
        double r196438 = r196436 * r196437;
        double r196439 = 1.0;
        double r196440 = v;
        double r196441 = r196440 * r196440;
        double r196442 = r196439 - r196441;
        double r196443 = r196438 * r196442;
        double r196444 = 2.0;
        double r196445 = 6.0;
        double r196446 = r196445 * r196441;
        double r196447 = r196444 - r196446;
        double r196448 = sqrt(r196447);
        double r196449 = r196443 * r196448;
        double r196450 = r196435 / r196449;
        return r196450;
}
↓
double f(double v) {
        double r196451 = 4.0;
        double r196452 = 3.0;
        double r196453 = atan2(1.0, 0.0);
        double r196454 = r196452 * r196453;
        double r196455 = r196451 / r196454;
        double r196456 = v;
        double r196457 = 1.0;
        double r196458 = fma(r196456, r196456, r196457);
        double r196459 = r196456 * r196456;
        double r196460 = r196457 - r196459;
        double r196461 = 2.0;
        double r196462 = 3.0;
        double r196463 = pow(r196461, r196462);
        double r196464 = 6.0;
        double r196465 = r196464 * r196459;
        double r196466 = pow(r196465, r196462);
        double r196467 = r196463 - r196466;
        double r196468 = sqrt(r196467);
        double r196469 = r196460 * r196468;
        double r196470 = r196458 * r196469;
        double r196471 = r196455 / r196470;
        double r196472 = r196457 + r196459;
        double r196473 = r196461 * r196461;
        double r196474 = r196465 * r196465;
        double r196475 = r196461 * r196465;
        double r196476 = r196474 + r196475;
        double r196477 = r196473 + r196476;
        double r196478 = sqrt(r196477);
        double r196479 = r196472 * r196478;
        double r196480 = r196471 * r196479;
        return r196480;
}

Error

Bits error versus v

Derivation

  1. Initial program 1.0

    \[\frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 - v \cdot v\right)\right) \cdot \sqrt{2 - 6 \cdot \left(v \cdot v\right)}}\]
  2. Using strategy rm

  3. Applied flip3--1.0

    \[\leadsto \frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 - v \cdot v\right)\right) \cdot \sqrt{\color{blue}{\frac{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}}}}\]

  4. Applied sqrt-div1.0

    \[\leadsto \frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 - v \cdot v\right)\right) \cdot \color{blue}{\frac{\sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}}{\sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}}}}\]

  5. Applied flip--1.0

    \[\leadsto \frac{4}{\left(\left(3 \cdot \pi\right) \cdot \color{blue}{\frac{1 \cdot 1 - \left(v \cdot v\right) \cdot \left(v \cdot v\right)}{1 + v \cdot v}}\right) \cdot \frac{\sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}}{\sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}}}\]

  6. Applied associate-*r/1.0

    \[\leadsto \frac{4}{\color{blue}{\frac{\left(3 \cdot \pi\right) \cdot \left(1 \cdot 1 - \left(v \cdot v\right) \cdot \left(v \cdot v\right)\right)}{1 + v \cdot v}} \cdot \frac{\sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}}{\sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}}}\]

  7. Applied frac-times1.0

    \[\leadsto \frac{4}{\color{blue}{\frac{\left(\left(3 \cdot \pi\right) \cdot \left(1 \cdot 1 - \left(v \cdot v\right) \cdot \left(v \cdot v\right)\right)\right) \cdot \sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}}{\left(1 + v \cdot v\right) \cdot \sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}}}}\]

  8. Applied associate-/r/1.0

    \[\leadsto \color{blue}{\frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 \cdot 1 - \left(v \cdot v\right) \cdot \left(v \cdot v\right)\right)\right) \cdot \sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}} \cdot \left(\left(1 + v \cdot v\right) \cdot \sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}\right)}\]

  9. Simplified0.0

    \[\leadsto \color{blue}{\frac{\frac{4}{3 \cdot \pi}}{\mathsf{fma}\left(v, v, 1\right) \cdot \left(\left(1 - v \cdot v\right) \cdot \sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}\right)}} \cdot \left(\left(1 + v \cdot v\right) \cdot \sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}\right)\]

  10. Final simplification0.0

    \[\leadsto \frac{\frac{4}{3 \cdot \pi}}{\mathsf{fma}\left(v, v, 1\right) \cdot \left(\left(1 - v \cdot v\right) \cdot \sqrt{{2}^{3} - {\left(6 \cdot \left(v \cdot v\right)\right)}^{3}}\right)} \cdot \left(\left(1 + v \cdot v\right) \cdot \sqrt{2 \cdot 2 + \left(\left(6 \cdot \left(v \cdot v\right)\right) \cdot \left(6 \cdot \left(v \cdot v\right)\right) + 2 \cdot \left(6 \cdot \left(v \cdot v\right)\right)\right)}\right)\]

Reproduce

herbie shell --seed 2019353 +o rules:numerics
(FPCore (v)
  :name "Falkner and Boettcher, Equation (22+)"
  :precision binary64
  (/ 4 (* (* (* 3 PI) (- 1 (* v v))) (sqrt (- 2 (* 6 (* v v)))))))