\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -3.1879146296582742 \cdot 10^{94}:\\
\;\;\;\;\frac{\frac{1}{2}}{\sqrt{\log 10}} \cdot \left(-2 \cdot \left(\log \left(\frac{-1}{re}\right) \cdot \sqrt{\frac{1}{\log 10}}\right)\right)\\
\mathbf{elif}\;re \le 1.0437409384331536 \cdot 10^{133}:\\
\;\;\;\;\frac{\frac{1}{2}}{\sqrt{\log 10}} \cdot \log \left({\left(re \cdot re + im \cdot im\right)}^{\left(\frac{1}{\sqrt{\log 10}}\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\log re}{\log 10}\\
\end{array}double code(double re, double im) {
return (log(sqrt(((re * re) + (im * im)))) / log(10.0));
}
double code(double re, double im) {
double VAR;
if ((re <= -3.187914629658274e+94)) {
VAR = ((0.5 / sqrt(log(10.0))) * (-2.0 * (log((-1.0 / re)) * sqrt((1.0 / log(10.0))))));
} else {
double VAR_1;
if ((re <= 1.0437409384331536e+133)) {
VAR_1 = ((0.5 / sqrt(log(10.0))) * log(pow(((re * re) + (im * im)), (1.0 / sqrt(log(10.0))))));
} else {
VAR_1 = (log(re) / log(10.0));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus re



Bits error versus im
Results
if re < -3.187914629658274e+94Initial program 50.8
rmApplied add-sqr-sqrt50.8
Applied pow1/250.8
Applied log-pow50.8
Applied times-frac50.8
Taylor expanded around -inf 9.4
if -3.187914629658274e+94 < re < 1.0437409384331536e+133Initial program 21.7
rmApplied add-sqr-sqrt21.7
Applied pow1/221.7
Applied log-pow21.7
Applied times-frac21.7
rmApplied add-log-exp21.7
Simplified21.5
if 1.0437409384331536e+133 < re Initial program 59.6
Taylor expanded around inf 7.7
Final simplification17.6
herbie shell --seed 2020105
(FPCore (re im)
:name "math.log10 on complex, real part"
:precision binary64
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))