Average Error: 1.0 → 0.0
Time: 3.1s
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 r286248 = 4.0;
        double r286249 = 3.0;
        double r286250 = atan2(1.0, 0.0);
        double r286251 = r286249 * r286250;
        double r286252 = 1.0;
        double r286253 = v;
        double r286254 = r286253 * r286253;
        double r286255 = r286252 - r286254;
        double r286256 = r286251 * r286255;
        double r286257 = 2.0;
        double r286258 = 6.0;
        double r286259 = r286258 * r286254;
        double r286260 = r286257 - r286259;
        double r286261 = sqrt(r286260);
        double r286262 = r286256 * r286261;
        double r286263 = r286248 / r286262;
        return r286263;
}
↓
double f(double v) {
        double r286264 = 4.0;
        double r286265 = 3.0;
        double r286266 = atan2(1.0, 0.0);
        double r286267 = r286265 * r286266;
        double r286268 = r286264 / r286267;
        double r286269 = v;
        double r286270 = 1.0;
        double r286271 = fma(r286269, r286269, r286270);
        double r286272 = r286269 * r286269;
        double r286273 = r286270 - r286272;
        double r286274 = 2.0;
        double r286275 = 3.0;
        double r286276 = pow(r286274, r286275);
        double r286277 = 6.0;
        double r286278 = r286277 * r286272;
        double r286279 = pow(r286278, r286275);
        double r286280 = r286276 - r286279;
        double r286281 = sqrt(r286280);
        double r286282 = r286273 * r286281;
        double r286283 = r286271 * r286282;
        double r286284 = r286268 / r286283;
        double r286285 = r286270 + r286272;
        double r286286 = r286274 * r286274;
        double r286287 = r286278 * r286278;
        double r286288 = r286274 * r286278;
        double r286289 = r286287 + r286288;
        double r286290 = r286286 + r286289;
        double r286291 = sqrt(r286290);
        double r286292 = r286285 * r286291;
        double r286293 = r286284 * r286292;
        return r286293;
}

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 2020039 +o rules:numerics
(FPCore (v)
  :name "Falkner and Boettcher, Equation (22+)"
  :precision binary64
  (/ 4 (* (* (* 3 PI) (- 1 (* v v))) (sqrt (- 2 (* 6 (* v v)))))))