Exception for quad2p (problem 3.2.1, positive)

Specification

Math

\[\begin{array}{l} \\ \frac{\left(-b\_2\right) + \sqrt{b\_2 \cdot b\_2 - a \cdot c}}{a} \end{array} \]

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

herbie shell --seed 2024039 
(FPCore (a b_2 c)
  :name "quad2p (problem 3.2.1, positive)"
  :precision binary64
  :herbie-expected 10

  :herbie-target
  (if (< b_2 0.0) (/ (- (if (== (copysign a c) a) (* (sqrt (- (fabs b_2) (* (sqrt (fabs a)) (sqrt (fabs c))))) (sqrt (+ (fabs b_2) (* (sqrt (fabs a)) (sqrt (fabs c)))))) (hypot b_2 (* (sqrt (fabs a)) (sqrt (fabs c))))) b_2) a) (/ (- c) (+ b_2 (if (== (copysign a c) a) (* (sqrt (- (fabs b_2) (* (sqrt (fabs a)) (sqrt (fabs c))))) (sqrt (+ (fabs b_2) (* (sqrt (fabs a)) (sqrt (fabs c)))))) (hypot b_2 (* (sqrt (fabs a)) (sqrt (fabs c))))))))

  (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))

Please file a bug report with this information.

Backtrace

member: not a proper list in: 2		L	C
gather-multiplicative-terms	/home/nightlies/herbie/sound-rules/src/core/reduce.rkt	125	0
simplify-node*	/home/nightlies/herbie/sound-rules/src/core/reduce.rkt	84	0
hash-ref!	/usr/local/racket/collects/racket/private/more-scheme.rkt	377	2
hash-ref!	/usr/local/racket/collects/racket/private/more-scheme.rkt	377	2
taylor-sqrt	/home/nightlies/herbie/sound-rules/src/core/taylor.rkt	295	0
hash-ref!	/usr/local/racket/collects/racket/private/more-scheme.rkt	377	2
approximate	/home/nightlies/herbie/sound-rules/src/core/taylor.rkt	7	0
taylor-expr	/home/nightlies/herbie/sound-rules/src/patch.rkt	89	0
taylor-alt	/home/nightlies/herbie/sound-rules/src/patch.rkt	96	0
gen-series!	/home/nightlies/herbie/sound-rules/src/patch.rkt	109	0
patch-table-run	/home/nightlies/herbie/sound-rules/src/patch.rkt	257	0
(unnamed)	/usr/local/racket/collects/racket/contract/private/arrow-val-first.rkt	486	18
run-iter!	/home/nightlies/herbie/sound-rules/src/mainloop.rkt	288	0
mutate!	/home/nightlies/herbie/sound-rules/src/mainloop.rkt	340	0
run-improve!	/home/nightlies/herbie/sound-rules/src/mainloop.rkt	326	0
get-alternatives/report	/home/nightlies/herbie/sound-rules/src/sandbox.rkt	169	0
(unnamed)	/home/nightlies/herbie/sound-rules/src/sandbox.rkt	262	6
(unnamed)	/usr/local/racket/share/pkgs/profile-lib/main.rkt	40	10
profile-thunk	/usr/local/racket/share/pkgs/profile-lib/main.rkt	9	0
(unnamed)	/usr/local/racket/collects/racket/engine.rkt	42	24