Average Error: 32.2 → 18.5
Time: 6.7s
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 <= -3.2792927917367154e+72)) {
		VAR = ((double) (((double) (((double) log(1.0)) - ((double) log(((double) (-1.0 / re)))))) / ((double) log(base))));
	} else {
		double VAR_1;
		if ((re <= 1.3644360600461893e-269)) {
			VAR_1 = ((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) (((double) cbrt(((double) pow(((double) log(base)), 6.0)))) + ((double) (0.0 * 0.0))))));
		} else {
			double VAR_2;
			if ((re <= 3.4160082491612644e-169)) {
				VAR_2 = ((double) (((double) (((double) log(1.0)) + ((double) log(im)))) / ((double) log(base))));
			} else {
				double VAR_3;
				if ((re <= 7.175104160147713e+17)) {
					VAR_3 = ((double) cbrt(((double) pow(((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) (((double) (0.0 * 0.0)) + ((double) pow(((double) log(base)), 2.0)))))), 3.0))));
				} 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 5 regimes
  2. if re < -3.2792927917367154e72

    1. Initial program 47.5

      \[\]
    2. Taylor expanded around -inf 64.0

      \[\leadsto \]
    3. Simplified10.2

      \[\leadsto \]

    if -3.2792927917367154e72 < re < 1.36443606004618934e-269

    1. Initial program 23.6

      \[\]
    2. Using strategy rm
    3. Applied add-cbrt-cube23.8

      \[\leadsto \]
    4. Applied add-cbrt-cube23.9

      \[\leadsto \]
    5. Applied cbrt-unprod23.8

      \[\leadsto \]
    6. Simplified23.7

      \[\leadsto \]

    if 1.36443606004618934e-269 < re < 3.41600824916126442e-169

    1. Initial program 34.7

      \[\]
    2. Taylor expanded around 0 35.7

      \[\leadsto \]
    3. Simplified35.7

      \[\leadsto \]

    if 3.41600824916126442e-169 < re < 717510416014771328

    1. Initial program 16.7

      \[\]
    2. Using strategy rm
    3. Applied add-cbrt-cube16.9

      \[\leadsto \]
    4. Applied add-cbrt-cube17.0

      \[\leadsto \]
    5. Applied cbrt-undiv16.9

      \[\leadsto \]
    6. Simplified16.9

      \[\leadsto \]

    if 717510416014771328 < re

    1. Initial program 42.3

      \[\]
    2. Taylor expanded around inf 12.7

      \[\leadsto \]
    3. Simplified12.7

      \[\leadsto \]
  3. Recombined 5 regimes into one program.
  4. Final simplification18.5

    \[\leadsto \]

Reproduce

herbie shell --seed 2020191 
(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))))