Average Error: 16.6 → 13.6
Time: 12.1s
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 ((t <= -9.173386750578478e-299)) {
		VAR = ((double) (((double) (x + ((double) (((double) (y / t)) * z)))) / ((double) (a + ((double) (1.0 + ((double) (((double) (y / t)) * b))))))));
	} else {
		double VAR_1;
		if ((t <= 1.8495783259794747e-30)) {
			VAR_1 = ((double) (((double) (x + ((double) (((double) (y * z)) / t)))) / ((double) (((double) (a + 1.0)) + ((double) (((double) (y * b)) * ((double) (1.0 / t))))))));
		} else {
			VAR_1 = ((double) (((double) (x + ((double) (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.6
Target13.4
Herbie13.6
\[\]
Derivation
  1. Split input into 3 regimes
  2. if  t  < -9.1733867505784784e-299
    1. Initial program 16.5
      \[\]
    2. Using strategy rm
    3. Applied add-cube-cbrt16.6
      \[\leadsto \]
    4. Applied associate-/r*16.6
      \[\leadsto \]
    5. Simplified15.9
      \[\leadsto \]
    6. Using strategy rm
    7. Applied div-inv15.9
      \[\leadsto \]
    8. Simplified15.9
      \[\leadsto \]
    9. Using strategy rm
    10. Applied associate-*r/15.8
      \[\leadsto \]
    11. Simplified14.2
      \[\leadsto \]
    if -9.1733867505784784e-299 <  t  < 1.8495783259794747e-30
    1. Initial program 23.3
      \[\]
    2. Using strategy rm
    3. Applied div-inv23.4
      \[\leadsto \]
    if 1.8495783259794747e-30 <  t 
    1. Initial program 11.6
      \[\]
    2. Simplified4.8
      \[\leadsto \]
  3. Recombined 3 regimes into one program.
  4. Final simplification13.6
    \[\leadsto \]
Reproduce
herbie shell --seed 2020179 
(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))))