Average Error: 16.2 → 13.1
Time: 8.6s
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 <= -3.126118079743165e-109)) {
		VAR = ((double) (((double) (x + ((double) (((double) (((double) cbrt(z)) / ((double) cbrt(t)))) * ((double) (y * ((double) (((double) (((double) cbrt(z)) / ((double) cbrt(t)))) * ((double) (((double) cbrt(z)) / ((double) cbrt(t)))))))))))) / ((double) (a + ((double) (1.0 + ((double) (y * ((double) (b / t))))))))));
	} else {
		double VAR_1;
		if ((y <= 9.475764007628411e+79)) {
			VAR_1 = ((double) (((double) (x + ((double) (((double) (z * ((double) (y / ((double) pow(((double) cbrt(t)), 2.0)))))) * ((double) (1.0 / ((double) cbrt(t)))))))) / ((double) (a + ((double) (1.0 + ((double) (((double) (((double) cbrt(b)) / ((double) cbrt(t)))) * ((double) (y * ((double) (((double) (((double) cbrt(b)) / ((double) cbrt(t)))) * ((double) (((double) cbrt(b)) / ((double) cbrt(t))))))))))))))));
		} else {
			double VAR_2;
			if ((y <= 1.0014304840356311e+244)) {
				VAR_2 = ((double) (1.0 / ((double) (((double) (a + ((double) (1.0 + ((double) (y * ((double) (b / t)))))))) / ((double) (x + ((double) (y * ((double) (z / t))))))))));
			} else {
				VAR_2 = ((double) (((double) (x + ((double) (((double) (y * z)) * ((double) (1.0 / t)))))) / ((double) (a + ((double) (1.0 + ((double) (y * ((double) (b / t))))))))));
			}
			VAR_1 = VAR_2;
		}
		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.2
Target13.2
Herbie13.1
\[\]
Derivation
  1. Split input into 4 regimes
  2. if  y  < -3.12611807974316519e-109
    1. Initial program 23.3
      \[\]
    2. Simplified18.8
      \[\leadsto \]
    3. Using strategy rm
    4. Applied add-cube-cbrt19.0
      \[\leadsto \]
    5. Applied add-cube-cbrt19.1
      \[\leadsto \]
    6. Applied times-frac19.1
      \[\leadsto \]
    7. Applied associate-*r*18.3
      \[\leadsto \]
    8. Simplified18.3
      \[\leadsto \]
    if -3.12611807974316519e-109 <  y  < 9.4757640076284113e79
    1. Initial program 5.0
      \[\]
    2. Simplified11.3
      \[\leadsto \]
    3. Using strategy rm
    4. Applied add-cube-cbrt11.4
      \[\leadsto \]
    5. Applied add-cube-cbrt11.4
      \[\leadsto \]
    6. Applied times-frac11.4
      \[\leadsto \]
    7. Applied associate-*r*8.9
      \[\leadsto \]
    8. Simplified8.9
      \[\leadsto \]
    9. Using strategy rm
    10. Applied add-cube-cbrt9.0
      \[\leadsto \]
    11. Applied *-un-lft-identity9.0
      \[\leadsto \]
    12. Applied times-frac9.0
      \[\leadsto \]
    13. Applied associate-*r*6.4
      \[\leadsto \]
    14. Simplified6.4
      \[\leadsto \]
    15. Using strategy rm
    16. Applied div-inv6.5
      \[\leadsto \]
    17. Applied associate-*r*5.1
      \[\leadsto \]
    18. Simplified5.1
      \[\leadsto \]
    if 9.4757640076284113e79 <  y  < 1.00143048403563111e244
    1. Initial program 30.3
      \[\]
    2. Simplified20.2
      \[\leadsto \]
    3. Using strategy rm
    4. Applied clear-num20.5
      \[\leadsto \]
    if 1.00143048403563111e244 <  y 
    1. Initial program 43.0
      \[\]
    2. Simplified33.7
      \[\leadsto \]
    3. Using strategy rm
    4. Applied div-inv33.7
      \[\leadsto \]
    5. Applied associate-*r*39.0
      \[\leadsto \]
  3. Recombined 4 regimes into one program.
  4. Final simplification13.1
    \[\leadsto \]
Reproduce
herbie shell --seed 2020191 
(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))))