\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 r563233 = im;
double r563234 = re;
double r563235 = atan2(r563233, r563234);
double r563236 = 10.0;
double r563237 = log(r563236);
double r563238 = r563235 / r563237;
return r563238;
}
double f(double re, double im) {
double r563239 = 1.0;
double r563240 = 10.0;
double r563241 = log(r563240);
double r563242 = sqrt(r563241);
double r563243 = r563239 / r563242;
double r563244 = im;
double r563245 = re;
double r563246 = atan2(r563244, r563245);
double r563247 = r563246 * r563243;
double r563248 = r563243 * r563247;
double r563249 = expm1(r563248);
double r563250 = log1p(r563249);
return r563250;
}



Bits error versus re



Bits error versus im
Results
Initial program 0.8
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 2019156 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, imaginary part"
(/ (atan2 im re) (log 10)))