Exception for Kahan's Monster

\[\]

Math

\[\]

C

double code(double y) {
	double VAR;
	if ((((double) (((double) (((double) fabs(((double) (y - ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))) - ((double) (1.0 / ((double) (y + ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))))) * ((double) (((double) fabs(((double) (y - ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))) - ((double) (1.0 / ((double) (y + ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))))))) == 0.0)) {
		VAR = 1.0;
	} else {
		VAR = ((double) (((double) (((double) exp(((double) (((double) (((double) fabs(((double) (y - ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))) - ((double) (1.0 / ((double) (y + ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))))) * ((double) (((double) fabs(((double) (y - ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))) - ((double) (1.0 / ((double) (y + ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))))))))) - 1.0)) / ((double) (((double) (((double) fabs(((double) (y - ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))) - ((double) (1.0 / ((double) (y + ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))))) * ((double) (((double) fabs(((double) (y - ((double) sqrt(((double) (((double) (y * y)) + 1.0)))))))) - ((double) (1.0 / ((double) (y + ((double) sqrt(((double) (((double) (y * y)) + 1.0))))))))))))));
	}
	return VAR;
}

Reproduce

herbie shell --seed 2020191 
(FPCore (y)
  :name "Kahan's Monster"
  :precision binary64
  :pre (<= 1.0 y 9999.0)
  (if (== (* (- (fabs (- y (sqrt (+ (* y y) 1.0)))) (/ 1.0 (+ y (sqrt (+ (* y y) 1.0))))) (- (fabs (- y (sqrt (+ (* y y) 1.0)))) (/ 1.0 (+ y (sqrt (+ (* y y) 1.0)))))) 0.0) 1.0 (/ (- (exp (* (- (fabs (- y (sqrt (+ (* y y) 1.0)))) (/ 1.0 (+ y (sqrt (+ (* y y) 1.0))))) (- (fabs (- y (sqrt (+ (* y y) 1.0)))) (/ 1.0 (+ y (sqrt (+ (* y y) 1.0))))))) 1.0) (* (- (fabs (- y (sqrt (+ (* y y) 1.0)))) (/ 1.0 (+ y (sqrt (+ (* y y) 1.0))))) (- (fabs (- y (sqrt (+ (* y y) 1.0)))) (/ 1.0 (+ y (sqrt (+ (* y y) 1.0)))))))))