Exception for ABCF->ab-angle b

Math

\[\frac{-\sqrt{\left(2 \cdot \left(\left({B}^{2} - \left(4 \cdot A\right) \cdot C\right) \cdot F\right)\right) \cdot \left(\left(A + C\right) - \sqrt{{\left(A - C\right)}^{2} + {B}^{2}}\right)}}{{B}^{2} - \left(4 \cdot A\right) \cdot C} \]

FPCore

(FPCore (A B C F)
 :precision binary64
 (/
  (-
   (sqrt
    (*
     (* 2.0 (* (- (pow B 2.0) (* (* 4.0 A) C)) F))
     (- (+ A C) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0)))))))
  (- (pow B 2.0) (* (* 4.0 A) C))))

C

double code(double A, double B, double C, double F) {
	return -sqrt(((2.0 * ((pow(B, 2.0) - ((4.0 * A) * C)) * F)) * ((A + C) - sqrt((pow((A - C), 2.0) + pow(B, 2.0)))))) / (pow(B, 2.0) - ((4.0 * A) * C));
}

Fortran

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

Java

public static double code(double A, double B, double C, double F) {
	return -Math.sqrt(((2.0 * ((Math.pow(B, 2.0) - ((4.0 * A) * C)) * F)) * ((A + C) - Math.sqrt((Math.pow((A - C), 2.0) + Math.pow(B, 2.0)))))) / (Math.pow(B, 2.0) - ((4.0 * A) * C));
}

Python

def code(A, B, C, F):
	return -math.sqrt(((2.0 * ((math.pow(B, 2.0) - ((4.0 * A) * C)) * F)) * ((A + C) - math.sqrt((math.pow((A - C), 2.0) + math.pow(B, 2.0)))))) / (math.pow(B, 2.0) - ((4.0 * A) * C))

Julia

function code(A, B, C, F)
	return Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(Float64((B ^ 2.0) - Float64(Float64(4.0 * A) * C)) * F)) * Float64(Float64(A + C) - sqrt(Float64((Float64(A - C) ^ 2.0) + (B ^ 2.0))))))) / Float64((B ^ 2.0) - Float64(Float64(4.0 * A) * C)))
end

MATLAB

function tmp = code(A, B, C, F)
	tmp = -sqrt(((2.0 * (((B ^ 2.0) - ((4.0 * A) * C)) * F)) * ((A + C) - sqrt((((A - C) ^ 2.0) + (B ^ 2.0)))))) / ((B ^ 2.0) - ((4.0 * A) * C));
end

Wolfram

code[A_, B_, C_, F_] := N[((-N[Sqrt[N[(N[(2.0 * N[(N[(N[Power[B, 2.0], $MachinePrecision] - N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]), $MachinePrecision] * F), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] - N[Sqrt[N[(N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[B, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / N[(N[Power[B, 2.0], $MachinePrecision] - N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2022150 
(FPCore (A B C F)
  :name "ABCF->ab-angle b"
  :precision binary64
  (/ (- (sqrt (* (* 2.0 (* (- (pow B 2.0) (* (* 4.0 A) C)) F)) (- (+ A C) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0))))))) (- (pow B 2.0) (* (* 4.0 A) C))))

Backtrace

application: required keyword argument not supplied procedure: batch-egg-rewrite required keyword: #:rules arguments...: '(#<rule erfc-erf_binary64> #<rule erf-erfc_binary64> #<rule erf-odd_binary64> #<rule prod-diff_binary64> #<rule hypot-udef_binary64> #<rule expm1-log1p-u_binary64> #<rule log1p-expm1-u_binary64> #<rule log1p-udef_binary64> #<rule expm1-udef_binary64> #... #:depths '(1 1 1 1) #:roots '((2) (2) (2) (2))		L	C
(unnamed)	/opt/racket-8.1/collects/racket/private/kw.rkt	1558	14
rewrite-expressions	/data/pavpan/nightlies/herbie/timeline/src/core/matcher.rkt	171	0
gen-rewrites!	/data/pavpan/nightlies/herbie/timeline/src/patch.rkt	144	0
patch-table-run	/data/pavpan/nightlies/herbie/timeline/src/patch.rkt	291	0
(unnamed)	/opt/racket-8.1/collects/racket/contract/private/arrow-val-first.rkt	486	18
run-iter!	/data/pavpan/nightlies/herbie/timeline/src/mainloop.rkt	246	0
mutate!	/data/pavpan/nightlies/herbie/timeline/src/mainloop.rkt	379	0
run-improve!	/data/pavpan/nightlies/herbie/timeline/src/mainloop.rkt	326	0
(unnamed)	/data/pavpan/nightlies/herbie/timeline/src/sandbox.rkt	66	6
(unnamed)	/opt/racket-8.1/share/pkgs/profile-lib/main.rkt	40	10
profile-thunk	/opt/racket-8.1/share/pkgs/profile-lib/main.rkt	9	0
(unnamed)	/opt/racket-8.1/collects/racket/engine.rkt	42	24