\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)}1.5 \cdot \frac{{v}^{2}}{t \cdot \left(\sqrt{2} \cdot \left(\sqrt{1} \cdot \pi\right)\right)} + \left(1 \cdot \frac{\sqrt{1}}{t \cdot \left(\sqrt{2} \cdot \pi\right)} - \left(\left(1.5 \cdot \frac{{v}^{4}}{t \cdot \left(\sqrt{2} \cdot \left(\sqrt{1} \cdot \pi\right)\right)} + 1.125 \cdot \frac{{v}^{4}}{t \cdot \left(\sqrt{2} \cdot \left({\left(\sqrt{1}\right)}^{3} \cdot \pi\right)\right)}\right) + 4 \cdot \left(\frac{{v}^{2} \cdot \sqrt{1}}{t \cdot \left(\sqrt{2} \cdot \pi\right)} + \frac{{v}^{4} \cdot \sqrt{1}}{t \cdot \left(\sqrt{2} \cdot \pi\right)}\right)\right)\right)double f(double v, double t) {
double r280037 = 1.0;
double r280038 = 5.0;
double r280039 = v;
double r280040 = r280039 * r280039;
double r280041 = r280038 * r280040;
double r280042 = r280037 - r280041;
double r280043 = atan2(1.0, 0.0);
double r280044 = t;
double r280045 = r280043 * r280044;
double r280046 = 2.0;
double r280047 = 3.0;
double r280048 = r280047 * r280040;
double r280049 = r280037 - r280048;
double r280050 = r280046 * r280049;
double r280051 = sqrt(r280050);
double r280052 = r280045 * r280051;
double r280053 = r280037 - r280040;
double r280054 = r280052 * r280053;
double r280055 = r280042 / r280054;
return r280055;
}
double f(double v, double t) {
double r280056 = 1.5;
double r280057 = v;
double r280058 = 2.0;
double r280059 = pow(r280057, r280058);
double r280060 = t;
double r280061 = 2.0;
double r280062 = sqrt(r280061);
double r280063 = 1.0;
double r280064 = sqrt(r280063);
double r280065 = atan2(1.0, 0.0);
double r280066 = r280064 * r280065;
double r280067 = r280062 * r280066;
double r280068 = r280060 * r280067;
double r280069 = r280059 / r280068;
double r280070 = r280056 * r280069;
double r280071 = r280062 * r280065;
double r280072 = r280060 * r280071;
double r280073 = r280064 / r280072;
double r280074 = r280063 * r280073;
double r280075 = 4.0;
double r280076 = pow(r280057, r280075);
double r280077 = r280076 / r280068;
double r280078 = r280056 * r280077;
double r280079 = 1.125;
double r280080 = 3.0;
double r280081 = pow(r280064, r280080);
double r280082 = r280081 * r280065;
double r280083 = r280062 * r280082;
double r280084 = r280060 * r280083;
double r280085 = r280076 / r280084;
double r280086 = r280079 * r280085;
double r280087 = r280078 + r280086;
double r280088 = 4.0;
double r280089 = r280059 * r280064;
double r280090 = r280089 / r280072;
double r280091 = r280076 * r280064;
double r280092 = r280091 / r280072;
double r280093 = r280090 + r280092;
double r280094 = r280088 * r280093;
double r280095 = r280087 + r280094;
double r280096 = r280074 - r280095;
double r280097 = r280070 + r280096;
return r280097;
}



Bits error versus v



Bits error versus t
Results
Initial program 0.5
Taylor expanded around 0 0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2019322
(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)))))