double code(double t, double l, double Om, double Omc) {
return ((double) asin(((double) sqrt(((double) (((double) (1.0 - ((double) pow(((double) (Om / Omc)), 2.0)))) / ((double) (1.0 + ((double) (2.0 * ((double) pow(((double) (t / l)), 2.0))))))))))));
}
double code(double t, double l, double Om, double Omc) {
return ((double) asin(((double) sqrt(((double) (((double) (1.0 - ((double) pow(((double) (Om / Omc)), 2.0)))) / ((double) (1.0 + ((double) (2.0 * ((double) pow(((double) (t / l)), 2.0))))))))))));
}



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus Omc
Results
Initial program 10.6
Final simplification10.6
herbie shell --seed 2020192
(FPCore (t l Om Omc)
:name "Toniolo and Linder, Equation (2)"
:precision binary64
(asin (sqrt (/ (- 1.0 (pow (/ Om Omc) 2.0)) (+ 1.0 (* 2.0 (pow (/ t l) 2.0)))))))