\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{\frac{\sqrt[3]{\frac{\mathsf{fma}\left(5 \cdot \left(-v\right), v, 1\right)}{{\left(\sqrt{\mathsf{fma}\left(-v, v \cdot 3, 1\right) \cdot 2}\right)}^{3}}}}{\mathsf{fma}\left(v, -v, 1\right)} \cdot \frac{\sqrt[3]{\mathsf{fma}\left(5 \cdot \left(-v\right), v, 1\right)} \cdot \sqrt[3]{\mathsf{fma}\left(5 \cdot \left(-v\right), v, 1\right)}}{\pi}}{t}double f(double v, double t) {
double r292092 = 1.0;
double r292093 = 5.0;
double r292094 = v;
double r292095 = r292094 * r292094;
double r292096 = r292093 * r292095;
double r292097 = r292092 - r292096;
double r292098 = atan2(1.0, 0.0);
double r292099 = t;
double r292100 = r292098 * r292099;
double r292101 = 2.0;
double r292102 = 3.0;
double r292103 = r292102 * r292095;
double r292104 = r292092 - r292103;
double r292105 = r292101 * r292104;
double r292106 = sqrt(r292105);
double r292107 = r292100 * r292106;
double r292108 = r292092 - r292095;
double r292109 = r292107 * r292108;
double r292110 = r292097 / r292109;
return r292110;
}
double f(double v, double t) {
double r292111 = 5.0;
double r292112 = v;
double r292113 = -r292112;
double r292114 = r292111 * r292113;
double r292115 = 1.0;
double r292116 = fma(r292114, r292112, r292115);
double r292117 = 3.0;
double r292118 = r292112 * r292117;
double r292119 = fma(r292113, r292118, r292115);
double r292120 = 2.0;
double r292121 = r292119 * r292120;
double r292122 = sqrt(r292121);
double r292123 = 3.0;
double r292124 = pow(r292122, r292123);
double r292125 = r292116 / r292124;
double r292126 = cbrt(r292125);
double r292127 = fma(r292112, r292113, r292115);
double r292128 = r292126 / r292127;
double r292129 = cbrt(r292116);
double r292130 = r292129 * r292129;
double r292131 = atan2(1.0, 0.0);
double r292132 = r292130 / r292131;
double r292133 = r292128 * r292132;
double r292134 = t;
double r292135 = r292133 / r292134;
return r292135;
}



Bits error versus v



Bits error versus t
Initial program 0.4
Simplified0.4
rmApplied *-un-lft-identity0.4
Applied add-cube-cbrt0.4
Applied times-frac0.4
Applied times-frac0.5
Simplified0.4
Simplified0.4
rmApplied associate-*r/0.3
Applied associate-*l/0.1
Simplified0.1
rmApplied add-cbrt-cube0.9
Applied cbrt-undiv0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2019179 +o rules:numerics
(FPCore (v t)
:name "Falkner and Boettcher, Equation (20:1,3)"
(/ (- 1.0 (* 5.0 (* v v))) (* (* (* PI t) (sqrt (* 2.0 (- 1.0 (* 3.0 (* v v)))))) (- 1.0 (* v v)))))