Average Error: 16.5 → 14.2
Time: 11.0s
Precision: binary64
\[\]
↓
\[\]
↓
double code(double x, double y, double z, double t, double a, double b) {
	return ((double) (((double) (x + ((double) (((double) (y * z)) / t)))) / ((double) (((double) (a + 1.0)) + ((double) (((double) (y * b)) / t))))));
}

↓
double code(double x, double y, double z, double t, double a, double b) {
	double VAR;
	if ((y <= -8.796357687792906e-230)) {
		VAR = ((double) (((double) (x + ((double) (((double) (((double) cbrt(y)) * ((double) (((double) cbrt(y)) / t)))) * ((double) (((double) cbrt(y)) * z)))))) * ((double) (1.0 / ((double) (a + ((double) (1.0 + ((double) (y * ((double) (b / t))))))))))));
	} else {
		double VAR_1;
		if ((y <= 1.649048289790591e-41)) {
			VAR_1 = ((double) (((double) (x + ((double) (((double) (y * z)) / t)))) / ((double) (((double) (a + 1.0)) + ((double) (((double) (y * b)) / t))))));
		} else {
			VAR_1 = ((double) (((double) (x + ((double) (((double) sqrt(y)) * ((double) (((double) sqrt(y)) * ((double) (z / t)))))))) / ((double) (a + ((double) (1.0 + ((double) (y * ((double) (b / t))))))))));
		}
		VAR = VAR_1;
	}
	return VAR;
}

Error
Bits error versus x
Bits error versus y
Bits error versus z
Bits error versus t
Bits error versus a
Bits error versus b
Try it out
Your Program's Arguments
Results
Enter valid numbers for all inputs
Target
Original16.5
Target13.1
Herbie14.2
\[\]
Derivation
  1. Split input into 3 regimes
  2. if  y  < -8.7963576877929056e-230
    1. Initial program 18.3
      \[\]
    2. Simplified16.6
      \[\leadsto \]
    3. Using strategy rm
    4. Applied div-inv16.6
      \[\leadsto \]
    5. Using strategy rm
    6. Applied associate-*r/18.1
      \[\leadsto \]
    7. Using strategy rm
    8. Applied associate-/l*16.5
      \[\leadsto \]
    9. Using strategy rm
    10. Applied div-inv16.5
      \[\leadsto \]
    11. Applied add-cube-cbrt16.7
      \[\leadsto \]
    12. Applied times-frac17.5
      \[\leadsto \]
    13. Simplified17.5
      \[\leadsto \]
    14. Simplified17.4
      \[\leadsto \]
    if -8.7963576877929056e-230 <  y  < 1.649048289790591e-41
    1. Initial program 3.2
      \[\]
    if 1.649048289790591e-41 <  y 
    1. Initial program 27.5
      \[\]
    2. Simplified20.3
      \[\leadsto \]
    3. Using strategy rm
    4. Applied add-sqr-sqrt20.4
      \[\leadsto \]
    5. Applied associate-*l*20.4
      \[\leadsto \]
    6. Simplified20.4
      \[\leadsto \]
  3. Recombined 3 regimes into one program.
  4. Final simplification14.2
    \[\leadsto \]
Reproduce
herbie shell --seed 2020181 
(FPCore (x y z t a b)
  :name "Diagrams.Solve.Tridiagonal:solveCyclicTriDiagonal from diagrams-solve-0.1, B"
  :precision binary64

  :herbie-target
  (if (< t -1.3659085366310088e-271) (* 1.0 (* (+ x (* (/ y t) z)) (/ 1.0 (+ (+ a 1.0) (* (/ y t) b))))) (if (< t 3.036967103737246e-130) (/ z b) (* 1.0 (* (+ x (* (/ y t) z)) (/ 1.0 (+ (+ a 1.0) (* (/ y t) b)))))))

  (/ (+ x (/ (* y z) t)) (+ (+ a 1.0) (/ (* y b) t))))