\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -1.300812438992646141617859246198844532718 \cdot 10^{100}:\\
\;\;\;\;\left(\left(\log \left(\frac{-1}{re}\right) \cdot -2\right) \cdot \frac{\frac{1}{2}}{\sqrt{\log 10}}\right) \cdot \frac{1}{\sqrt{\log 10}}\\
\mathbf{elif}\;re \le 5.714402214507161350041984173167312711037 \cdot 10^{91}:\\
\;\;\;\;\frac{1}{\sqrt{\log 10}} \cdot \left(\log \left(im \cdot im + re \cdot re\right) \cdot \frac{\frac{1}{2}}{\sqrt{\log 10}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\sqrt{\log 10}} \cdot \left(\sqrt{\frac{1}{\log 10}} \cdot \log re\right)\\
\end{array}double f(double re, double im) {
double r1408314 = re;
double r1408315 = r1408314 * r1408314;
double r1408316 = im;
double r1408317 = r1408316 * r1408316;
double r1408318 = r1408315 + r1408317;
double r1408319 = sqrt(r1408318);
double r1408320 = log(r1408319);
double r1408321 = 10.0;
double r1408322 = log(r1408321);
double r1408323 = r1408320 / r1408322;
return r1408323;
}
double f(double re, double im) {
double r1408324 = re;
double r1408325 = -1.3008124389926461e+100;
bool r1408326 = r1408324 <= r1408325;
double r1408327 = -1.0;
double r1408328 = r1408327 / r1408324;
double r1408329 = log(r1408328);
double r1408330 = -2.0;
double r1408331 = r1408329 * r1408330;
double r1408332 = 0.5;
double r1408333 = 10.0;
double r1408334 = log(r1408333);
double r1408335 = sqrt(r1408334);
double r1408336 = r1408332 / r1408335;
double r1408337 = r1408331 * r1408336;
double r1408338 = 1.0;
double r1408339 = r1408338 / r1408335;
double r1408340 = r1408337 * r1408339;
double r1408341 = 5.714402214507161e+91;
bool r1408342 = r1408324 <= r1408341;
double r1408343 = im;
double r1408344 = r1408343 * r1408343;
double r1408345 = r1408324 * r1408324;
double r1408346 = r1408344 + r1408345;
double r1408347 = log(r1408346);
double r1408348 = r1408347 * r1408336;
double r1408349 = r1408339 * r1408348;
double r1408350 = r1408338 / r1408334;
double r1408351 = sqrt(r1408350);
double r1408352 = log(r1408324);
double r1408353 = r1408351 * r1408352;
double r1408354 = r1408339 * r1408353;
double r1408355 = r1408342 ? r1408349 : r1408354;
double r1408356 = r1408326 ? r1408340 : r1408355;
return r1408356;
}



Bits error versus re



Bits error versus im
Results
if re < -1.3008124389926461e+100Initial program 51.9
rmApplied add-sqr-sqrt51.9
Applied pow1/251.9
Applied log-pow51.9
Applied times-frac51.9
rmApplied div-inv51.9
Applied associate-*r*51.9
Taylor expanded around -inf 9.2
Simplified9.2
if -1.3008124389926461e+100 < re < 5.714402214507161e+91Initial program 21.9
rmApplied add-sqr-sqrt21.9
Applied pow1/221.9
Applied log-pow21.9
Applied times-frac21.9
rmApplied div-inv21.8
Applied associate-*r*21.8
if 5.714402214507161e+91 < re Initial program 49.9
rmApplied add-sqr-sqrt49.9
Applied pow1/249.9
Applied log-pow49.9
Applied times-frac49.9
rmApplied div-inv49.9
Applied associate-*r*49.9
Taylor expanded around inf 8.9
Simplified8.9
Final simplification17.3
herbie shell --seed 2019171
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10.0)))