\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}}{\left(1 \cdot 1 - {v}^{3} \cdot v\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)double f(double v) {
double r130291 = 4.0;
double r130292 = 3.0;
double r130293 = atan2(1.0, 0.0);
double r130294 = r130292 * r130293;
double r130295 = 1.0;
double r130296 = v;
double r130297 = r130296 * r130296;
double r130298 = r130295 - r130297;
double r130299 = r130294 * r130298;
double r130300 = 2.0;
double r130301 = 6.0;
double r130302 = r130301 * r130297;
double r130303 = r130300 - r130302;
double r130304 = sqrt(r130303);
double r130305 = r130299 * r130304;
double r130306 = r130291 / r130305;
return r130306;
}
double f(double v) {
double r130307 = 4.0;
double r130308 = 3.0;
double r130309 = atan2(1.0, 0.0);
double r130310 = r130308 * r130309;
double r130311 = r130307 / r130310;
double r130312 = 1.0;
double r130313 = r130312 * r130312;
double r130314 = v;
double r130315 = 3.0;
double r130316 = pow(r130314, r130315);
double r130317 = r130316 * r130314;
double r130318 = r130313 - r130317;
double r130319 = 2.0;
double r130320 = pow(r130319, r130315);
double r130321 = 6.0;
double r130322 = r130314 * r130314;
double r130323 = r130321 * r130322;
double r130324 = pow(r130323, r130315);
double r130325 = r130320 - r130324;
double r130326 = sqrt(r130325);
double r130327 = r130318 * r130326;
double r130328 = r130311 / r130327;
double r130329 = r130312 + r130322;
double r130330 = r130319 * r130319;
double r130331 = r130323 * r130323;
double r130332 = r130319 * r130323;
double r130333 = r130331 + r130332;
double r130334 = r130330 + r130333;
double r130335 = sqrt(r130334);
double r130336 = r130329 * r130335;
double r130337 = r130328 * r130336;
return r130337;
}



Bits error versus v
Results
Initial program 1.0
rmApplied flip3--1.0
Applied sqrt-div1.0
Applied flip--1.0
Applied associate-*r/1.0
Applied frac-times1.0
Applied associate-/r/1.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019323 +o rules:numerics
(FPCore (v)
:name "Falkner and Boettcher, Equation (22+)"
:precision binary64
(/ 4 (* (* (* 3 PI) (- 1 (* v v))) (sqrt (- 2 (* 6 (* v v)))))))