\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -1.2199670469480383 \cdot 10^{95}:\\
\;\;\;\;\frac{\sqrt{\frac{1}{2}}}{\sqrt{1}} \cdot \left(\frac{\sqrt{\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)\right)\\
\mathbf{elif}\;re \le 3.8849079706542796 \cdot 10^{45}:\\
\;\;\;\;\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{\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)\\
\end{array}double f(double re, double im) {
double r55846 = re;
double r55847 = r55846 * r55846;
double r55848 = im;
double r55849 = r55848 * r55848;
double r55850 = r55847 + r55849;
double r55851 = sqrt(r55850);
double r55852 = log(r55851);
double r55853 = 10.0;
double r55854 = log(r55853);
double r55855 = r55852 / r55854;
return r55855;
}
double f(double re, double im) {
double r55856 = re;
double r55857 = -1.2199670469480383e+95;
bool r55858 = r55856 <= r55857;
double r55859 = 0.5;
double r55860 = sqrt(r55859);
double r55861 = 1.0;
double r55862 = sqrt(r55861);
double r55863 = r55860 / r55862;
double r55864 = 10.0;
double r55865 = log(r55864);
double r55866 = sqrt(r55865);
double r55867 = r55860 / r55866;
double r55868 = -2.0;
double r55869 = -1.0;
double r55870 = r55869 / r55856;
double r55871 = log(r55870);
double r55872 = r55861 / r55865;
double r55873 = sqrt(r55872);
double r55874 = r55871 * r55873;
double r55875 = r55868 * r55874;
double r55876 = r55867 * r55875;
double r55877 = r55863 * r55876;
double r55878 = 3.8849079706542796e+45;
bool r55879 = r55856 <= r55878;
double r55880 = r55859 / r55866;
double r55881 = r55856 * r55856;
double r55882 = im;
double r55883 = r55882 * r55882;
double r55884 = r55881 + r55883;
double r55885 = r55861 / r55866;
double r55886 = pow(r55884, r55885);
double r55887 = log(r55886);
double r55888 = r55880 * r55887;
double r55889 = r55861 / r55856;
double r55890 = log(r55889);
double r55891 = r55890 * r55873;
double r55892 = r55868 * r55891;
double r55893 = r55880 * r55892;
double r55894 = r55879 ? r55888 : r55893;
double r55895 = r55858 ? r55877 : r55894;
return r55895;
}



Bits error versus re



Bits error versus im
Results
if re < -1.2199670469480383e+95Initial program 50.4
rmApplied add-sqr-sqrt50.4
Applied pow1/250.4
Applied log-pow50.4
Applied times-frac50.4
rmApplied pow150.4
Applied log-pow50.4
Applied sqrt-prod50.4
Applied add-sqr-sqrt50.4
Applied times-frac50.4
Applied associate-*l*50.4
Taylor expanded around -inf 9.2
if -1.2199670469480383e+95 < re < 3.8849079706542796e+45Initial program 21.9
rmApplied add-sqr-sqrt21.9
Applied pow1/221.9
Applied log-pow21.9
Applied times-frac21.8
rmApplied add-log-exp21.8
Simplified21.7
if 3.8849079706542796e+45 < re Initial program 44.8
rmApplied add-sqr-sqrt44.8
Applied pow1/244.8
Applied log-pow44.8
Applied times-frac44.8
Taylor expanded around inf 11.8
Final simplification17.4
herbie shell --seed 2020065
(FPCore (re im)
:name "math.log10 on complex, real part"
:precision binary64
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))