Exception for Kahan's Monster

• could not determine a ground truth

  1. y = 12.309990923302163

\[1 \leq y \land y \leq 9999\]

Math

\[\begin{array}{l} \mathbf{if}\;\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) = 0:\\ \;\;\;\;1\\ \mathbf{else}:\\ \;\;\;\;\frac{e^{\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right)} - 1}{\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right)}\\ \end{array} \]

FPCore

(FPCore (y)
 :precision binary64
 (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)))))))))

C

double code(double y) {
	double tmp;
	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) {
		tmp = 1.0;
	} else {
		tmp = (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))))));
	}
	return tmp;
}

Fortran

real(8) function code(y)
    real(8), intent (in) :: y
    real(8) :: tmp
    if (((abs((y - sqrt(((y * y) + 1.0d0)))) - (1.0d0 / (y + sqrt(((y * y) + 1.0d0))))) * (abs((y - sqrt(((y * y) + 1.0d0)))) - (1.0d0 / (y + sqrt(((y * y) + 1.0d0)))))) == 0.0d0) then
        tmp = 1.0d0
    else
        tmp = (exp(((abs((y - sqrt(((y * y) + 1.0d0)))) - (1.0d0 / (y + sqrt(((y * y) + 1.0d0))))) * (abs((y - sqrt(((y * y) + 1.0d0)))) - (1.0d0 / (y + sqrt(((y * y) + 1.0d0))))))) - 1.0d0) / ((abs((y - sqrt(((y * y) + 1.0d0)))) - (1.0d0 / (y + sqrt(((y * y) + 1.0d0))))) * (abs((y - sqrt(((y * y) + 1.0d0)))) - (1.0d0 / (y + sqrt(((y * y) + 1.0d0))))))
    end if
    code = tmp
end function

Java

public static double code(double y) {
	double tmp;
	if (((Math.abs((y - Math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + Math.sqrt(((y * y) + 1.0))))) * (Math.abs((y - Math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + Math.sqrt(((y * y) + 1.0)))))) == 0.0) {
		tmp = 1.0;
	} else {
		tmp = (Math.exp(((Math.abs((y - Math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + Math.sqrt(((y * y) + 1.0))))) * (Math.abs((y - Math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + Math.sqrt(((y * y) + 1.0))))))) - 1.0) / ((Math.abs((y - Math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + Math.sqrt(((y * y) + 1.0))))) * (Math.abs((y - Math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + Math.sqrt(((y * y) + 1.0))))));
	}
	return tmp;
}

Python

def code(y):
	tmp = 0
	if ((math.fabs((y - math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + math.sqrt(((y * y) + 1.0))))) * (math.fabs((y - math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + math.sqrt(((y * y) + 1.0)))))) == 0.0:
		tmp = 1.0
	else:
		tmp = (math.exp(((math.fabs((y - math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + math.sqrt(((y * y) + 1.0))))) * (math.fabs((y - math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + math.sqrt(((y * y) + 1.0))))))) - 1.0) / ((math.fabs((y - math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + math.sqrt(((y * y) + 1.0))))) * (math.fabs((y - math.sqrt(((y * y) + 1.0)))) - (1.0 / (y + math.sqrt(((y * y) + 1.0))))))
	return tmp

Julia

function code(y)
	tmp = 0.0
	if (Float64(Float64(abs(Float64(y - sqrt(Float64(Float64(y * y) + 1.0)))) - Float64(1.0 / Float64(y + sqrt(Float64(Float64(y * y) + 1.0))))) * Float64(abs(Float64(y - sqrt(Float64(Float64(y * y) + 1.0)))) - Float64(1.0 / Float64(y + sqrt(Float64(Float64(y * y) + 1.0)))))) == 0.0)
		tmp = 1.0;
	else
		tmp = Float64(Float64(exp(Float64(Float64(abs(Float64(y - sqrt(Float64(Float64(y * y) + 1.0)))) - Float64(1.0 / Float64(y + sqrt(Float64(Float64(y * y) + 1.0))))) * Float64(abs(Float64(y - sqrt(Float64(Float64(y * y) + 1.0)))) - Float64(1.0 / Float64(y + sqrt(Float64(Float64(y * y) + 1.0))))))) - 1.0) / Float64(Float64(abs(Float64(y - sqrt(Float64(Float64(y * y) + 1.0)))) - Float64(1.0 / Float64(y + sqrt(Float64(Float64(y * y) + 1.0))))) * Float64(abs(Float64(y - sqrt(Float64(Float64(y * y) + 1.0)))) - Float64(1.0 / Float64(y + sqrt(Float64(Float64(y * y) + 1.0)))))));
	end
	return tmp
end

MATLAB

function tmp_2 = code(y)
	tmp = 0.0;
	if (((abs((y - sqrt(((y * y) + 1.0)))) - (1.0 / (y + sqrt(((y * y) + 1.0))))) * (abs((y - sqrt(((y * y) + 1.0)))) - (1.0 / (y + sqrt(((y * y) + 1.0)))))) == 0.0)
		tmp = 1.0;
	else
		tmp = (exp(((abs((y - sqrt(((y * y) + 1.0)))) - (1.0 / (y + sqrt(((y * y) + 1.0))))) * (abs((y - sqrt(((y * y) + 1.0)))) - (1.0 / (y + sqrt(((y * y) + 1.0))))))) - 1.0) / ((abs((y - sqrt(((y * y) + 1.0)))) - (1.0 / (y + sqrt(((y * y) + 1.0))))) * (abs((y - sqrt(((y * y) + 1.0)))) - (1.0 / (y + sqrt(((y * y) + 1.0))))));
	end
	tmp_2 = tmp;
end

Wolfram

code[y_] := If[Equal[N[(N[(N[Abs[N[(y - N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(1.0 / N[(y + N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Abs[N[(y - N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(1.0 / N[(y + N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.0], 1.0, N[(N[(N[Exp[N[(N[(N[Abs[N[(y - N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(1.0 / N[(y + N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Abs[N[(y - N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(1.0 / N[(y + N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - 1.0), $MachinePrecision] / N[(N[(N[Abs[N[(y - N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(1.0 / N[(y + N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Abs[N[(y - N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - N[(1.0 / N[(y + N[Sqrt[N[(N[(y * y), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]