\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -3.962076038462479 \cdot 10^{+99}:\\
\;\;\;\;\frac{\log \left(-re\right)}{\log 10}\\
\mathbf{elif}\;re \le 9.628393082979933 \cdot 10^{+82}:\\
\;\;\;\;\frac{1}{\sqrt{\log 10}} \cdot \left(\left(\frac{\sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt{\log 10}}} \cdot \log \left(re \cdot re + im \cdot im\right)\right) \cdot \frac{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt{\log 10}}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\sqrt{\log 10}} \cdot \left(\log re \cdot \sqrt{\frac{1}{\log 10}}\right)\\
\end{array}double f(double re, double im) {
double r1243497 = re;
double r1243498 = r1243497 * r1243497;
double r1243499 = im;
double r1243500 = r1243499 * r1243499;
double r1243501 = r1243498 + r1243500;
double r1243502 = sqrt(r1243501);
double r1243503 = log(r1243502);
double r1243504 = 10.0;
double r1243505 = log(r1243504);
double r1243506 = r1243503 / r1243505;
return r1243506;
}
double f(double re, double im) {
double r1243507 = re;
double r1243508 = -3.962076038462479e+99;
bool r1243509 = r1243507 <= r1243508;
double r1243510 = -r1243507;
double r1243511 = log(r1243510);
double r1243512 = 10.0;
double r1243513 = log(r1243512);
double r1243514 = r1243511 / r1243513;
double r1243515 = 9.628393082979933e+82;
bool r1243516 = r1243507 <= r1243515;
double r1243517 = 1.0;
double r1243518 = sqrt(r1243513);
double r1243519 = r1243517 / r1243518;
double r1243520 = 0.5;
double r1243521 = cbrt(r1243520);
double r1243522 = sqrt(r1243518);
double r1243523 = r1243521 / r1243522;
double r1243524 = r1243507 * r1243507;
double r1243525 = im;
double r1243526 = r1243525 * r1243525;
double r1243527 = r1243524 + r1243526;
double r1243528 = log(r1243527);
double r1243529 = r1243523 * r1243528;
double r1243530 = r1243521 * r1243521;
double r1243531 = r1243530 / r1243522;
double r1243532 = r1243529 * r1243531;
double r1243533 = r1243519 * r1243532;
double r1243534 = log(r1243507);
double r1243535 = r1243517 / r1243513;
double r1243536 = sqrt(r1243535);
double r1243537 = r1243534 * r1243536;
double r1243538 = r1243519 * r1243537;
double r1243539 = r1243516 ? r1243533 : r1243538;
double r1243540 = r1243509 ? r1243514 : r1243539;
return r1243540;
}



Bits error versus re



Bits error versus im
Results
if re < -3.962076038462479e+99Initial program 49.7
Taylor expanded around -inf 8.2
Simplified8.2
if -3.962076038462479e+99 < re < 9.628393082979933e+82Initial program 21.7
rmApplied add-sqr-sqrt21.7
Applied pow1/221.7
Applied log-pow21.7
Applied times-frac21.7
rmApplied div-inv21.6
Applied associate-*r*21.6
rmApplied add-sqr-sqrt21.6
Applied sqrt-prod22.0
Applied add-cube-cbrt21.6
Applied times-frac21.6
Applied associate-*l*21.6
if 9.628393082979933e+82 < re Initial program 46.5
rmApplied add-sqr-sqrt46.5
Applied pow1/246.5
Applied log-pow46.5
Applied times-frac46.5
rmApplied div-inv46.5
Applied associate-*r*46.5
Taylor expanded around inf 10.1
Simplified10.1
Final simplification17.2
herbie shell --seed 2019168
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))