\sqrt{1 - x \cdot x}\mathsf{log1p}\left(\mathsf{expm1}\left(\sqrt{1 - x \cdot x}\right)\right)double code(double x) {
return ((double) sqrt(((double) (1.0 - ((double) (x * x))))));
}
double code(double x) {
return ((double) log1p(((double) expm1(((double) sqrt(((double) (1.0 - ((double) (x * x))))))))));
}



Bits error versus x
Results
Initial program 0.0
rmApplied log1p-expm1-u0.0
Final simplification0.0
herbie shell --seed 2020120 +o rules:numerics
(FPCore (x)
:name "Diagrams.TwoD.Ellipse:ellipse from diagrams-lib-1.3.0.3"
:precision binary64
(sqrt (- 1 (* x x))))