Exception for Kahan's Unum-Targeted 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) + {\left({10}^{-300}\right)}^{\left(10000 \cdot \left(y + 1\right)\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) + {\left({10}^{-300}\right)}^{\left(10000 \cdot \left(y + 1\right)\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) + {\left({10}^{-300}\right)}^{\left(10000 \cdot \left(y + 1\right)\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))))))
       (pow (pow 10.0 -300.0) (* 10000.0 (+ 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))))))
       (pow (pow 10.0 -300.0) (* 10000.0 (+ 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))))))
     (pow (pow 10.0 -300.0) (* 10000.0 (+ 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)))))) + pow(pow(10.0, -300.0), (10000.0 * (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)))))) + pow(pow(10.0, -300.0), (10000.0 * (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)))))) + pow(pow(10.0, -300.0), (10000.0 * (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)))))) + ((10.0d0 ** (-300.0d0)) ** (10000.0d0 * (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)))))) + ((10.0d0 ** (-300.0d0)) ** (10000.0d0 * (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)))))) + ((10.0d0 ** (-300.0d0)) ** (10000.0d0 * (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)))))) + Math.pow(Math.pow(10.0, -300.0), (10000.0 * (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)))))) + Math.pow(Math.pow(10.0, -300.0), (10000.0 * (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)))))) + Math.pow(Math.pow(10.0, -300.0), (10000.0 * (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)))))) + math.pow(math.pow(10.0, -300.0), (10000.0 * (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)))))) + math.pow(math.pow(10.0, -300.0), (10000.0 * (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)))))) + math.pow(math.pow(10.0, -300.0), (10000.0 * (y + 1.0))))
	return tmp

Julia

function code(y)
	tmp = 0.0
	if (Float64(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)))))) + ((10.0 ^ -300.0) ^ Float64(10000.0 * Float64(y + 1.0)))) == 0.0)
		tmp = 1.0;
	else
		tmp = Float64(Float64(exp(Float64(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)))))) + ((10.0 ^ -300.0) ^ Float64(10000.0 * Float64(y + 1.0))))) - 1.0) / Float64(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)))))) + ((10.0 ^ -300.0) ^ Float64(10000.0 * Float64(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)))))) + ((10.0 ^ -300.0) ^ (10000.0 * (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)))))) + ((10.0 ^ -300.0) ^ (10000.0 * (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)))))) + ((10.0 ^ -300.0) ^ (10000.0 * (y + 1.0))));
	end
	tmp_2 = tmp;
end

Wolfram

code[y_] := If[Equal[N[(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] + N[Power[N[Power[10.0, -300.0], $MachinePrecision], N[(10000.0 * N[(y + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 0.0], 1.0, N[(N[(N[Exp[N[(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] + N[Power[N[Power[10.0, -300.0], $MachinePrecision], N[(10000.0 * N[(y + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] - 1.0), $MachinePrecision] / N[(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] + N[Power[N[Power[10.0, -300.0], $MachinePrecision], N[(10000.0 * N[(y + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]