\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)}{\frac{\pi \cdot \left(t \cdot \sqrt{2 \cdot \left({1}^{3} - {\left(3 \cdot \left(v \cdot v\right)\right)}^{3}\right)}\right)}{\sqrt{1 \cdot 1 + \left(\left(3 \cdot \left(v \cdot v\right)\right) \cdot \left(3 \cdot \left(v \cdot v\right)\right) + 1 \cdot \left(3 \cdot \left(v \cdot v\right)\right)\right)}} \cdot \left(1 - v \cdot v\right)}double f(double v, double t) {
double r255196 = 1.0;
double r255197 = 5.0;
double r255198 = v;
double r255199 = r255198 * r255198;
double r255200 = r255197 * r255199;
double r255201 = r255196 - r255200;
double r255202 = atan2(1.0, 0.0);
double r255203 = t;
double r255204 = r255202 * r255203;
double r255205 = 2.0;
double r255206 = 3.0;
double r255207 = r255206 * r255199;
double r255208 = r255196 - r255207;
double r255209 = r255205 * r255208;
double r255210 = sqrt(r255209);
double r255211 = r255204 * r255210;
double r255212 = r255196 - r255199;
double r255213 = r255211 * r255212;
double r255214 = r255201 / r255213;
return r255214;
}
double f(double v, double t) {
double r255215 = 1.0;
double r255216 = 5.0;
double r255217 = v;
double r255218 = r255217 * r255217;
double r255219 = r255216 * r255218;
double r255220 = r255215 - r255219;
double r255221 = atan2(1.0, 0.0);
double r255222 = t;
double r255223 = 2.0;
double r255224 = 3.0;
double r255225 = pow(r255215, r255224);
double r255226 = 3.0;
double r255227 = r255226 * r255218;
double r255228 = pow(r255227, r255224);
double r255229 = r255225 - r255228;
double r255230 = r255223 * r255229;
double r255231 = sqrt(r255230);
double r255232 = r255222 * r255231;
double r255233 = r255221 * r255232;
double r255234 = r255215 * r255215;
double r255235 = r255227 * r255227;
double r255236 = r255215 * r255227;
double r255237 = r255235 + r255236;
double r255238 = r255234 + r255237;
double r255239 = sqrt(r255238);
double r255240 = r255233 / r255239;
double r255241 = r255215 - r255218;
double r255242 = r255240 * r255241;
double r255243 = r255220 / r255242;
return r255243;
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
rmApplied associate-*l*0.4
rmApplied flip3--0.4
Applied associate-*r/0.4
Applied sqrt-div0.4
Applied associate-*r/0.4
Applied associate-*r/0.4
Final simplification0.4
herbie shell --seed 2020034 +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)))))