\frac{\tan^{-1}_* \frac{im}{re}}{\log 10}\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{1}{\sqrt{\log 10}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \frac{1}{\sqrt{\log 10}}\right)\right)\right)double f(double re, double im) {
double r45667 = im;
double r45668 = re;
double r45669 = atan2(r45667, r45668);
double r45670 = 10.0;
double r45671 = log(r45670);
double r45672 = r45669 / r45671;
return r45672;
}
double f(double re, double im) {
double r45673 = 1.0;
double r45674 = 10.0;
double r45675 = log(r45674);
double r45676 = sqrt(r45675);
double r45677 = r45673 / r45676;
double r45678 = im;
double r45679 = re;
double r45680 = atan2(r45678, r45679);
double r45681 = r45680 * r45677;
double r45682 = r45677 * r45681;
double r45683 = expm1(r45682);
double r45684 = log1p(r45683);
return r45684;
}



Bits error versus re



Bits error versus im
Results
Initial program 0.9
rmApplied log1p-expm1-u0.7
rmApplied add-sqr-sqrt0.7
Applied *-un-lft-identity0.7
Applied times-frac0.7
rmApplied div-inv0.7
Final simplification0.7
herbie shell --seed 2019325 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, imaginary part"
:precision binary64
(/ (atan2 im re) (log 10)))