\frac{1 - 5 \cdot \left(v \cdot v\right)}{\left(\left(\pi \cdot t\right) \cdot \sqrt{2 \cdot \left(1 - 3 \cdot \left(v \cdot v\right)\right)}\right) \cdot \left(1 - v \cdot v\right)}\frac{1 - 5 \cdot \left(v \cdot v\right)}{\left(\left(\pi \cdot \left(t \cdot \left(\sqrt[3]{\sqrt{2 \cdot \left(1 - 3 \cdot \left(v \cdot v\right)\right)}} \cdot \sqrt[3]{\sqrt{2 \cdot \left(1 - 3 \cdot \left(v \cdot v\right)\right)}}\right)\right)\right) \cdot \sqrt[3]{\sqrt{2 \cdot \left(1 - 3 \cdot \left(v \cdot v\right)\right)}}\right) \cdot \left(1 - v \cdot v\right)}double f(double v, double t) {
double r264226 = 1.0;
double r264227 = 5.0;
double r264228 = v;
double r264229 = r264228 * r264228;
double r264230 = r264227 * r264229;
double r264231 = r264226 - r264230;
double r264232 = atan2(1.0, 0.0);
double r264233 = t;
double r264234 = r264232 * r264233;
double r264235 = 2.0;
double r264236 = 3.0;
double r264237 = r264236 * r264229;
double r264238 = r264226 - r264237;
double r264239 = r264235 * r264238;
double r264240 = sqrt(r264239);
double r264241 = r264234 * r264240;
double r264242 = r264226 - r264229;
double r264243 = r264241 * r264242;
double r264244 = r264231 / r264243;
return r264244;
}
double f(double v, double t) {
double r264245 = 1.0;
double r264246 = 5.0;
double r264247 = v;
double r264248 = r264247 * r264247;
double r264249 = r264246 * r264248;
double r264250 = r264245 - r264249;
double r264251 = atan2(1.0, 0.0);
double r264252 = t;
double r264253 = 2.0;
double r264254 = 3.0;
double r264255 = r264254 * r264248;
double r264256 = r264245 - r264255;
double r264257 = r264253 * r264256;
double r264258 = sqrt(r264257);
double r264259 = cbrt(r264258);
double r264260 = r264259 * r264259;
double r264261 = r264252 * r264260;
double r264262 = r264251 * r264261;
double r264263 = r264262 * r264259;
double r264264 = r264245 - r264248;
double r264265 = r264263 * r264264;
double r264266 = r264250 / r264265;
return r264266;
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
rmApplied add-cube-cbrt0.5
Applied associate-*r*0.5
rmApplied associate-*l*0.5
Final simplification0.5
herbie shell --seed 2019347 +o rules:numerics
(FPCore (v t)
:name "Falkner and Boettcher, Equation (20:1,3)"
:precision binary64
(/ (- 1 (* 5 (* v v))) (* (* (* PI t) (sqrt (* 2 (- 1 (* 3 (* v v)))))) (- 1 (* v v)))))