\frac{\left(-t1\right) \cdot v}{\left(t1 + u\right) \cdot \left(t1 + u\right)}\frac{\sqrt[3]{-t1} \cdot \sqrt[3]{-t1}}{\sqrt[3]{t1 + u} \cdot \sqrt[3]{t1 + u}} \cdot \left(\frac{\sqrt[3]{-t1}}{\sqrt[3]{t1 + u}} \cdot \frac{v}{t1 + u}\right)double code(double u, double v, double t1) {
return ((double) (((double) (((double) -(t1)) * v)) / ((double) (((double) (t1 + u)) * ((double) (t1 + u))))));
}
double code(double u, double v, double t1) {
return ((double) (((double) (((double) (((double) cbrt(((double) -(t1)))) * ((double) cbrt(((double) -(t1)))))) / ((double) (((double) cbrt(((double) (t1 + u)))) * ((double) cbrt(((double) (t1 + u)))))))) * ((double) (((double) (((double) cbrt(((double) -(t1)))) / ((double) cbrt(((double) (t1 + u)))))) * ((double) (v / ((double) (t1 + u))))))));
}



Bits error versus u



Bits error versus v



Bits error versus t1
Results
Initial program 19.1
rmApplied times-frac1.5
rmApplied add-cube-cbrt2.2
Applied add-cube-cbrt1.8
Applied times-frac1.8
Applied associate-*l*1.2
Final simplification1.2
herbie shell --seed 2020121 +o rules:numerics
(FPCore (u v t1)
:name "Rosa's DopplerBench"
:precision binary64
(/ (* (- t1) v) (* (+ t1 u) (+ t1 u))))