\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\frac{1}{\sqrt{\log 10}} \cdot \log \left({\left(\mathsf{hypot}\left(re, im\right)\right)}^{\left(\frac{1}{\sqrt{\log 10}}\right)}\right)double f(double re, double im) {
double r90362 = re;
double r90363 = r90362 * r90362;
double r90364 = im;
double r90365 = r90364 * r90364;
double r90366 = r90363 + r90365;
double r90367 = sqrt(r90366);
double r90368 = log(r90367);
double r90369 = 10.0;
double r90370 = log(r90369);
double r90371 = r90368 / r90370;
return r90371;
}
double f(double re, double im) {
double r90372 = 1.0;
double r90373 = 10.0;
double r90374 = log(r90373);
double r90375 = sqrt(r90374);
double r90376 = r90372 / r90375;
double r90377 = re;
double r90378 = im;
double r90379 = hypot(r90377, r90378);
double r90380 = pow(r90379, r90376);
double r90381 = log(r90380);
double r90382 = r90376 * r90381;
return r90382;
}



Bits error versus re



Bits error versus im
Results
Initial program 32.0
rmApplied hypot-def0.6
rmApplied add-sqr-sqrt0.6
Applied pow10.6
Applied log-pow0.6
Applied times-frac0.5
rmApplied add-log-exp0.5
Simplified0.3
Final simplification0.3
herbie shell --seed 2019362 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, real part"
:precision binary64
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))