Average Error: 20.7 → 0.1
Time: 4.8s
Precision: binary64
\[\]
↓
\[\]
↓
double code(double x, double y, double z) {
	return ((double) (x + ((double) (((double) (y * ((double) (((double) (((double) (((double) (z * 0.0692910599291889)) + 0.4917317610505968)) * z)) + 0.279195317918525)))) / ((double) (((double) (((double) (z + 6.012459259764103)) * z)) + 3.350343815022304))))));
}

↓
double code(double x, double y, double z) {
	double VAR;
	if (((z <= -1.0595051241030401e+27) || !(z <= 135609.02453109337))) {
		VAR = ((double) (x + ((double) (((double) (0.0692910599291889 * y)) + ((double) (((double) (y / z)) * ((double) (0.07512208616047561 - ((double) (0.40462203869992125 / z))))))))));
	} else {
		VAR = ((double) (x + ((double) (((double) (y * ((double) (((double) (z * ((double) (((double) (z * 0.0692910599291889)) + 0.4917317610505968)))) + 0.279195317918525)))) / ((double) (((double) (z * ((double) (z + 6.012459259764103)))) + 3.350343815022304))))));
	}
	return VAR;
}

Error
Bits error versus x
Bits error versus y
Bits error versus z
Try it out
Your Program's Arguments
Results
Enter valid numbers for all inputs
Target
Original20.7
Target0.1
Herbie0.1
\[\]
Derivation
  1. Split input into 2 regimes
  2. if  z  < -1.0595051241030401e27 or 135609.02453109337 <  z 
    1. Initial program 42.6
      \[\]
    2. Simplified34.0
      \[\leadsto \]
    3. Taylor expanded around inf 0.0
      \[\leadsto \]
    4. Simplified0.0
      \[\leadsto \]
    if -1.0595051241030401e27 <  z  < 135609.02453109337
    1. Initial program 0.3
      \[\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.1
    \[\leadsto \]
Reproduce
herbie shell --seed 2020180 
(FPCore (x y z)
  :name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, B"
  :precision binary64

  :herbie-target
  (if (< z -8120153.652456675) (- (* (+ (/ 0.07512208616047561 z) 0.0692910599291889) y) (- (/ (* 0.40462203869992125 y) (* z z)) x)) (if (< z 6.576118972787377e+20) (+ x (* (* y (+ (* (+ (* z 0.0692910599291889) 0.4917317610505968) z) 0.279195317918525)) (/ 1.0 (+ (* (+ z 6.012459259764103) z) 3.350343815022304)))) (- (* (+ (/ 0.07512208616047561 z) 0.0692910599291889) y) (- (/ (* 0.40462203869992125 y) (* z z)) x))))

  (+ x (/ (* y (+ (* (+ (* z 0.0692910599291889) 0.4917317610505968) z) 0.279195317918525)) (+ (* (+ z 6.012459259764103) z) 3.350343815022304))))