Exception for quad2p (problem 3.2.1, positive)

Math

\[\frac{\left(-b_2\right) + \sqrt{b_2 \cdot b_2 - a \cdot c}}{a} \]

FPCore

(FPCore (a b_2 c)
 :precision binary64
 (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))

C

double code(double a, double b_2, double c) {
	return (-b_2 + sqrt(((b_2 * b_2) - (a * c)))) / a;
}

Fortran

real(8) function code(a, b_2, c)
    real(8), intent (in) :: a
    real(8), intent (in) :: b_2
    real(8), intent (in) :: c
    code = (-b_2 + sqrt(((b_2 * b_2) - (a * c)))) / a
end function

Java

public static double code(double a, double b_2, double c) {
	return (-b_2 + Math.sqrt(((b_2 * b_2) - (a * c)))) / a;
}

Python

def code(a, b_2, c):
	return (-b_2 + math.sqrt(((b_2 * b_2) - (a * c)))) / a

Julia

function code(a, b_2, c)
	return Float64(Float64(Float64(-b_2) + sqrt(Float64(Float64(b_2 * b_2) - Float64(a * c)))) / a)
end

MATLAB

function tmp = code(a, b_2, c)
	tmp = (-b_2 + sqrt(((b_2 * b_2) - (a * c)))) / a;
end

Wolfram

code[a_, b$95$2_, c_] := N[(N[((-b$95$2) + N[Sqrt[N[(N[(b$95$2 * b$95$2), $MachinePrecision] - N[(a * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2023241 
(FPCore (a b_2 c)
  :name "quad2p (problem 3.2.1, positive)"
  :precision binary64
  (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))

Backtrace

inexact?: contract violation expected: number? given: (ival (bf "-2.098578716467387692404358116883839070632e323228496") (bf "2.098578716467387692404358116883839070632e323228496"))		L	C
real->bigfloat	/usr/share/racket/pkgs/math-lib/math/private/bigfloat/mpfr.rkt	470	0
...row-higher-order.rkt:375:33	/usr/share/racket/collects/racket/contract/private/arrow-higher-order.rkt	375	33
arg->precision	/home/nightlies/herbie/arb-support/src/programs.rkt	138	11
f	/home/nightlies/herbie/arb-support/src/programs.rkt	190	2
(unnamed)	/home/nightlies/herbie/arb-support/src/ground-truth.rkt	23	2
search-step	/home/nightlies/herbie/arb-support/src/searchreals.rkt	37	0
find-intervals	/home/nightlies/herbie/arb-support/src/searchreals.rkt	87	0
make-sampler	/home/nightlies/herbie/arb-support/src/sampling.rkt	94	0
sample-points	/home/nightlies/herbie/arb-support/src/ground-truth.rkt	64	0
setup-context!	/home/nightlies/herbie/arb-support/src/mainloop.rkt	293	0
get-alternatives/report	/home/nightlies/herbie/arb-support/src/sandbox.rkt	148	0
(unnamed)	/home/nightlies/herbie/arb-support/src/sandbox.rkt	247	6
(unnamed)	/usr/share/racket/pkgs/profile-lib/main.rkt	40	10
profile-thunk	/usr/share/racket/pkgs/profile-lib/main.rkt	9	0
(unnamed)	/usr/share/racket/collects/racket/engine.rkt	42	24