\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 r36009 = re;
double r36010 = r36009 * r36009;
double r36011 = im;
double r36012 = r36011 * r36011;
double r36013 = r36010 + r36012;
double r36014 = sqrt(r36013);
double r36015 = log(r36014);
double r36016 = 10.0;
double r36017 = log(r36016);
double r36018 = r36015 / r36017;
return r36018;
}
double f(double re, double im) {
double r36019 = 1.0;
double r36020 = 10.0;
double r36021 = log(r36020);
double r36022 = sqrt(r36021);
double r36023 = r36019 / r36022;
double r36024 = re;
double r36025 = im;
double r36026 = hypot(r36024, r36025);
double r36027 = pow(r36026, r36023);
double r36028 = log(r36027);
double r36029 = r36023 * r36028;
return r36029;
}



Bits error versus re



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