\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -1.1845577271365163 \cdot 10^{+79}:\\
\;\;\;\;\frac{\frac{1}{2}}{\sqrt{\log 10}} \cdot \left(\left(\sqrt{\frac{1}{\log 10}} \cdot \log \left(\frac{-1}{re}\right)\right) \cdot -2\right)\\
\mathbf{elif}\;re \le 1.1144412990449351 \cdot 10^{-247}:\\
\;\;\;\;\left(\frac{\sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10}}} \cdot \frac{\log \left(re \cdot re + im \cdot im\right)}{\sqrt{\log 10}}\right) \cdot \frac{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10} \cdot \sqrt[3]{\log 10}}}\\
\mathbf{elif}\;re \le 6.513147769961864 \cdot 10^{-175}:\\
\;\;\;\;\left(2 \cdot \left(\log im \cdot \sqrt{\frac{1}{\log 10}}\right)\right) \cdot \frac{\frac{1}{2}}{\sqrt{\log 10}}\\
\mathbf{elif}\;re \le 1.0899421385131229 \cdot 10^{+149}:\\
\;\;\;\;\left(\frac{\sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10}}} \cdot \frac{\log \left(re \cdot re + im \cdot im\right)}{\sqrt{\log 10}}\right) \cdot \frac{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10} \cdot \sqrt[3]{\log 10}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{2}}{\sqrt{\log 10}} \cdot \left(\left(2 \cdot \log re\right) \cdot \sqrt{\frac{1}{\log 10}}\right)\\
\end{array}double f(double re, double im) {
double r739624 = re;
double r739625 = r739624 * r739624;
double r739626 = im;
double r739627 = r739626 * r739626;
double r739628 = r739625 + r739627;
double r739629 = sqrt(r739628);
double r739630 = log(r739629);
double r739631 = 10.0;
double r739632 = log(r739631);
double r739633 = r739630 / r739632;
return r739633;
}
double f(double re, double im) {
double r739634 = re;
double r739635 = -1.1845577271365163e+79;
bool r739636 = r739634 <= r739635;
double r739637 = 0.5;
double r739638 = 10.0;
double r739639 = log(r739638);
double r739640 = sqrt(r739639);
double r739641 = r739637 / r739640;
double r739642 = 1.0;
double r739643 = r739642 / r739639;
double r739644 = sqrt(r739643);
double r739645 = -1.0;
double r739646 = r739645 / r739634;
double r739647 = log(r739646);
double r739648 = r739644 * r739647;
double r739649 = -2.0;
double r739650 = r739648 * r739649;
double r739651 = r739641 * r739650;
double r739652 = 1.1144412990449351e-247;
bool r739653 = r739634 <= r739652;
double r739654 = cbrt(r739637);
double r739655 = cbrt(r739639);
double r739656 = sqrt(r739655);
double r739657 = r739654 / r739656;
double r739658 = r739634 * r739634;
double r739659 = im;
double r739660 = r739659 * r739659;
double r739661 = r739658 + r739660;
double r739662 = log(r739661);
double r739663 = r739662 / r739640;
double r739664 = r739657 * r739663;
double r739665 = r739654 * r739654;
double r739666 = r739655 * r739655;
double r739667 = sqrt(r739666);
double r739668 = r739665 / r739667;
double r739669 = r739664 * r739668;
double r739670 = 6.513147769961864e-175;
bool r739671 = r739634 <= r739670;
double r739672 = 2.0;
double r739673 = log(r739659);
double r739674 = r739673 * r739644;
double r739675 = r739672 * r739674;
double r739676 = r739675 * r739641;
double r739677 = 1.0899421385131229e+149;
bool r739678 = r739634 <= r739677;
double r739679 = log(r739634);
double r739680 = r739672 * r739679;
double r739681 = r739680 * r739644;
double r739682 = r739641 * r739681;
double r739683 = r739678 ? r739669 : r739682;
double r739684 = r739671 ? r739676 : r739683;
double r739685 = r739653 ? r739669 : r739684;
double r739686 = r739636 ? r739651 : r739685;
return r739686;
}



Bits error versus re



Bits error versus im
Results
if re < -1.1845577271365163e+79Initial program 47.7
rmApplied add-sqr-sqrt47.7
Applied pow1/247.7
Applied log-pow47.7
Applied times-frac47.7
Taylor expanded around -inf 9.8
if -1.1845577271365163e+79 < re < 1.1144412990449351e-247 or 6.513147769961864e-175 < re < 1.0899421385131229e+149Initial program 20.5
rmApplied add-sqr-sqrt20.5
Applied pow1/220.5
Applied log-pow20.5
Applied times-frac20.5
rmApplied add-cube-cbrt21.0
Applied sqrt-prod21.0
Applied add-cube-cbrt20.5
Applied times-frac20.5
Applied associate-*l*20.4
if 1.1144412990449351e-247 < re < 6.513147769961864e-175Initial program 33.6
rmApplied add-sqr-sqrt33.6
Applied pow1/233.6
Applied log-pow33.6
Applied times-frac33.6
Taylor expanded around 0 38.3
if 1.0899421385131229e+149 < re Initial program 60.9
rmApplied add-sqr-sqrt60.9
Applied pow1/260.9
Applied log-pow60.9
Applied times-frac60.9
Taylor expanded around inf 6.6
Simplified6.6
Final simplification17.7
herbie shell --seed 2019133
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))