Exception for quadm (p42, negative)

Math

\[\frac{\left(-b\right) - \sqrt{b \cdot b - 4 \cdot \left(a \cdot c\right)}}{2 \cdot a} \]

FPCore

(FPCore (a b c)
 :precision binary64
 (/ (- (- b) (sqrt (- (* b b) (* 4.0 (* a c))))) (* 2.0 a)))

C

double code(double a, double b, double c) {
	return (-b - sqrt(((b * b) - (4.0 * (a * c))))) / (2.0 * a);
}

Fortran

real(8) function code(a, b, c)
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8), intent (in) :: c
    code = (-b - sqrt(((b * b) - (4.0d0 * (a * c))))) / (2.0d0 * a)
end function

Java

public static double code(double a, double b, double c) {
	return (-b - Math.sqrt(((b * b) - (4.0 * (a * c))))) / (2.0 * a);
}

Python

def code(a, b, c):
	return (-b - math.sqrt(((b * b) - (4.0 * (a * c))))) / (2.0 * a)

Julia

function code(a, b, c)
	return Float64(Float64(Float64(-b) - sqrt(Float64(Float64(b * b) - Float64(4.0 * Float64(a * c))))) / Float64(2.0 * a))
end

MATLAB

function tmp = code(a, b, c)
	tmp = (-b - sqrt(((b * b) - (4.0 * (a * c))))) / (2.0 * a);
end

Wolfram

code[a_, b_, c_] := N[(N[((-b) - N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(4.0 * N[(a * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2023096 
(FPCore (a b c)
  :name "quadm (p42, negative)"
  :precision binary64

  :herbie-target
  (if (< b 0.0) (/ c (* a (/ (+ (- b) (sqrt (- (* b b) (* 4.0 (* a c))))) (* 2.0 a)))) (/ (- (- b) (sqrt (- (* b b) (* 4.0 (* a c))))) (* 2.0 a)))

  (/ (- (- b) (sqrt (- (* b b) (* 4.0 (* a c))))) (* 2.0 a)))

Backtrace

match-define: no matching clause for '((x -.f64 (exp.f64 (log1p.f64 (sqrt.f64 (fma.f64 c (*.f64 a -4) (*.f64 b b))))) 1))		L	C
program-body	/home/nightlies/herbie/remove-change/src/programs.rkt	21	0
(unnamed)	/usr/share/racket/collects/racket/contract/private/arrow-higher-order.rkt	379	33
gen-rewrites!	/home/nightlies/herbie/remove-change/src/patch.rkt	145	0
patch-table-run	/home/nightlies/herbie/remove-change/src/patch.rkt	296	0
(unnamed)	/usr/share/racket/collects/racket/contract/private/arrow-val-first.rkt	486	18
run-iter!	/home/nightlies/herbie/remove-change/src/mainloop.rkt	267	0
mutate!	/home/nightlies/herbie/remove-change/src/mainloop.rkt	377	0
run-improve!	/home/nightlies/herbie/remove-change/src/mainloop.rkt	325	0
run-herbie	/home/nightlies/herbie/remove-change/src/sandbox.rkt	57	0
(unnamed)	/home/nightlies/herbie/remove-change/src/sandbox.rkt	133	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