Average Error: 32.5 → 18.5
Time: 6.3s
Precision: binary64
\[\]
↓
\[\]
↓
double code(double re, double im, double base) {
	return ((double) (((double) (((double) (((double) log(((double) sqrt(((double) (((double) (re * re)) + ((double) (im * im)))))))) * ((double) log(base)))) + ((double) (((double) atan2(im, re)) * 0.0)))) / ((double) (((double) (((double) log(base)) * ((double) log(base)))) + ((double) (0.0 * 0.0))))));
}

↓
double code(double re, double im, double base) {
	double VAR;
	if ((re <= -5.852465244811857e+36)) {
		VAR = ((double) (((double) (((double) log(1.0)) - ((double) log(((double) (-1.0 / re)))))) / ((double) log(base))));
	} else {
		double VAR_1;
		if ((re <= -9.484903639874046e-162)) {
			VAR_1 = ((double) (((double) (((double) (((double) (((double) log(base)) * ((double) log(((double) sqrt(((double) (((double) (re * re)) + ((double) (im * im)))))))))) + ((double) (((double) atan2(im, re)) * 0.0)))) / ((double) sqrt(((double) (((double) pow(((double) log(base)), 2.0)) + ((double) (0.0 * 0.0)))))))) / ((double) sqrt(((double) (((double) (0.0 * 0.0)) + ((double) (((double) log(base)) * ((double) log(base))))))))));
		} else {
			double VAR_2;
			if ((re <= -8.579711569628729e-271)) {
				VAR_2 = ((double) (((double) (((double) log(1.0)) + ((double) log(im)))) / ((double) log(base))));
			} else {
				double VAR_3;
				if ((re <= 2.370492566939671e+98)) {
					VAR_3 = ((double) (((double) (((double) (((double) (((double) log(base)) * ((double) log(((double) sqrt(((double) (((double) (re * re)) + ((double) (im * im)))))))))) + ((double) (((double) atan2(im, re)) * 0.0)))) / ((double) sqrt(((double) (((double) pow(((double) log(base)), 2.0)) + ((double) (0.0 * 0.0)))))))) / ((double) sqrt(((double) (((double) (0.0 * 0.0)) + ((double) (((double) log(base)) * ((double) log(base))))))))));
				} else {
					VAR_3 = ((double) (((double) log(re)) / ((double) log(base))));
				}
				VAR_2 = VAR_3;
			}
			VAR_1 = VAR_2;
		}
		VAR = VAR_1;
	}
	return VAR;
}

Error
Bits error versus re
Bits error versus im
Bits error versus base
Try it out
Your Program's Arguments
Results
Enter valid numbers for all inputs
Derivation
  1. Split input into 4 regimes
  2. if  re  < -5.8524652448118575e36
    1. Initial program 43.0
      \[\]
    2. Taylor expanded around -inf 64.0
      \[\leadsto \]
    3. Simplified12.4
      \[\leadsto \]
    if -5.8524652448118575e36 <  re  < -9.48490363987404599e-162 or -8.57971156962872867e-271 <  re  < 2.370492566939671e98
    1. Initial program 21.5
      \[\]
    2. Using strategy rm
    3. Applied add-sqr-sqrt21.5
      \[\leadsto \]
    4. Applied associate-/r*21.4
      \[\leadsto \]
    5. Simplified21.4
      \[\leadsto \]
    if -9.48490363987404599e-162 <  re  < -8.57971156962872867e-271
    1. Initial program 32.2
      \[\]
    2. Taylor expanded around 0 36.1
      \[\leadsto \]
    3. Simplified36.1
      \[\leadsto \]
    if 2.370492566939671e98 <  re 
    1. Initial program 52.0
      \[\]
    2. Taylor expanded around inf 8.9
      \[\leadsto \]
    3. Simplified8.9
      \[\leadsto \]
  3. Recombined 4 regimes into one program.
  4. Final simplification18.5
    \[\leadsto \]
Reproduce
herbie shell --seed 2020192 
(FPCore (re im base)
  :name "math.log/2 on complex, real part"
  :precision binary64
  (/ (+ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (atan2 im re) 0.0)) (+ (* (log base) (log base)) (* 0.0 0.0))))