\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 r42301 = im;
double r42302 = re;
double r42303 = atan2(r42301, r42302);
double r42304 = 10.0;
double r42305 = log(r42304);
double r42306 = r42303 / r42305;
return r42306;
}
double f(double re, double im) {
double r42307 = 1.0;
double r42308 = 10.0;
double r42309 = log(r42308);
double r42310 = sqrt(r42309);
double r42311 = r42307 / r42310;
double r42312 = im;
double r42313 = re;
double r42314 = atan2(r42312, r42313);
double r42315 = r42314 * r42311;
double r42316 = r42311 * r42315;
double r42317 = expm1(r42316);
double r42318 = log1p(r42317);
return r42318;
}



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 2019303 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, imaginary part"
:precision binary64
(/ (atan2 im re) (log 10)))