\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 r24223 = im;
double r24224 = re;
double r24225 = atan2(r24223, r24224);
double r24226 = 10.0;
double r24227 = log(r24226);
double r24228 = r24225 / r24227;
return r24228;
}
double f(double re, double im) {
double r24229 = 1.0;
double r24230 = 10.0;
double r24231 = log(r24230);
double r24232 = sqrt(r24231);
double r24233 = r24229 / r24232;
double r24234 = im;
double r24235 = re;
double r24236 = atan2(r24234, r24235);
double r24237 = r24236 * r24233;
double r24238 = r24233 * r24237;
double r24239 = expm1(r24238);
double r24240 = log1p(r24239);
return r24240;
}



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