\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -6.3372341303149501 \cdot 10^{130}:\\
\;\;\;\;\frac{1}{\sqrt{\log 10}} \cdot \log \left({\left(-1 \cdot re\right)}^{\left(\frac{1}{\sqrt{\log 10}}\right)}\right)\\
\mathbf{elif}\;re \le 1.3249567362868981 \cdot 10^{86}:\\
\;\;\;\;\frac{1}{\sqrt{\log 10}} \cdot \log \left({\left(\sqrt{\sqrt{re \cdot re + im \cdot im}}\right)}^{\left(\frac{1}{\sqrt{\log 10}}\right)} \cdot {\left(\sqrt{\sqrt{re \cdot re + im \cdot im}}\right)}^{\left(\frac{1}{\sqrt{\log 10}}\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\sqrt{\log 10}} \cdot \log \left({\left(\frac{1}{re}\right)}^{\left(-\sqrt{\frac{1}{\log 10}}\right)}\right)\\
\end{array}double f(double re, double im) {
double r47865 = re;
double r47866 = r47865 * r47865;
double r47867 = im;
double r47868 = r47867 * r47867;
double r47869 = r47866 + r47868;
double r47870 = sqrt(r47869);
double r47871 = log(r47870);
double r47872 = 10.0;
double r47873 = log(r47872);
double r47874 = r47871 / r47873;
return r47874;
}
double f(double re, double im) {
double r47875 = re;
double r47876 = -6.33723413031495e+130;
bool r47877 = r47875 <= r47876;
double r47878 = 1.0;
double r47879 = 10.0;
double r47880 = log(r47879);
double r47881 = sqrt(r47880);
double r47882 = r47878 / r47881;
double r47883 = -1.0;
double r47884 = r47883 * r47875;
double r47885 = pow(r47884, r47882);
double r47886 = log(r47885);
double r47887 = r47882 * r47886;
double r47888 = 1.3249567362868981e+86;
bool r47889 = r47875 <= r47888;
double r47890 = r47875 * r47875;
double r47891 = im;
double r47892 = r47891 * r47891;
double r47893 = r47890 + r47892;
double r47894 = sqrt(r47893);
double r47895 = sqrt(r47894);
double r47896 = pow(r47895, r47882);
double r47897 = r47896 * r47896;
double r47898 = log(r47897);
double r47899 = r47882 * r47898;
double r47900 = r47878 / r47875;
double r47901 = r47878 / r47880;
double r47902 = sqrt(r47901);
double r47903 = -r47902;
double r47904 = pow(r47900, r47903);
double r47905 = log(r47904);
double r47906 = r47882 * r47905;
double r47907 = r47889 ? r47899 : r47906;
double r47908 = r47877 ? r47887 : r47907;
return r47908;
}



Bits error versus re



Bits error versus im
Results
if re < -6.33723413031495e+130Initial program 58.7
rmApplied add-sqr-sqrt58.7
Applied pow158.7
Applied log-pow58.7
Applied times-frac58.7
rmApplied add-log-exp58.7
Simplified58.7
Taylor expanded around -inf 7.6
if -6.33723413031495e+130 < re < 1.3249567362868981e+86Initial program 21.7
rmApplied add-sqr-sqrt21.7
Applied pow121.7
Applied log-pow21.7
Applied times-frac21.6
rmApplied add-log-exp21.6
Simplified21.5
rmApplied add-sqr-sqrt21.5
Applied sqrt-prod21.5
Applied unpow-prod-down21.5
if 1.3249567362868981e+86 < re Initial program 50.1
rmApplied add-sqr-sqrt50.1
Applied pow150.1
Applied log-pow50.1
Applied times-frac50.1
rmApplied add-log-exp50.1
Simplified50.1
Taylor expanded around inf 9.5
Simplified9.4
Final simplification17.3
herbie shell --seed 2020027
(FPCore (re im)
:name "math.log10 on complex, real part"
:precision binary64
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))