\frac{\tan^{-1}_* \frac{im}{re}}{\log 10}\mathsf{log1p}\left(\mathsf{expm1}\left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\right)\right)\right)double f(double re, double im) {
double r43788 = im;
double r43789 = re;
double r43790 = atan2(r43788, r43789);
double r43791 = 10.0;
double r43792 = log(r43791);
double r43793 = r43790 / r43792;
return r43793;
}
double f(double re, double im) {
double r43794 = 1.0;
double r43795 = 10.0;
double r43796 = log(r43795);
double r43797 = sqrt(r43796);
double r43798 = r43794 / r43797;
double r43799 = sqrt(r43798);
double r43800 = im;
double r43801 = re;
double r43802 = atan2(r43800, r43801);
double r43803 = r43794 / r43796;
double r43804 = sqrt(r43803);
double r43805 = r43802 * r43804;
double r43806 = r43799 * r43805;
double r43807 = r43799 * r43806;
double r43808 = expm1(r43807);
double r43809 = log1p(r43808);
return r43809;
}



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
Taylor expanded around 0 0.7
rmApplied add-sqr-sqrt0.7
Applied associate-*l*0.7
Final simplification0.7
herbie shell --seed 2019354 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, imaginary part"
:precision binary64
(/ (atan2 im re) (log 10)))