Exception for VandenBroeck and Keller, Equation (23)

Specification

Math

\[\begin{array}{l} \\ \left(-x \cdot \frac{1}{\tan B}\right) + \frac{F}{\sin B} \cdot {\left(\left(F \cdot F + 2\right) + 2 \cdot x\right)}^{\left(-\frac{1}{2}\right)} \end{array} \]

FPCore

(FPCore (F B x)
 :precision binary64
 (+
  (- (* x (/ 1.0 (tan B))))
  (* (/ F (sin B)) (pow (+ (+ (* F F) 2.0) (* 2.0 x)) (- (/ 1.0 2.0))))))

C

double code(double F, double B, double x) {
	return -(x * (1.0 / tan(B))) + ((F / sin(B)) * pow((((F * F) + 2.0) + (2.0 * x)), -(1.0 / 2.0)));
}

Fortran

real(8) function code(f, b, x)
    real(8), intent (in) :: f
    real(8), intent (in) :: b
    real(8), intent (in) :: x
    code = -(x * (1.0d0 / tan(b))) + ((f / sin(b)) * ((((f * f) + 2.0d0) + (2.0d0 * x)) ** -(1.0d0 / 2.0d0)))
end function

Java

public static double code(double F, double B, double x) {
	return -(x * (1.0 / Math.tan(B))) + ((F / Math.sin(B)) * Math.pow((((F * F) + 2.0) + (2.0 * x)), -(1.0 / 2.0)));
}

Python

def code(F, B, x):
	return -(x * (1.0 / math.tan(B))) + ((F / math.sin(B)) * math.pow((((F * F) + 2.0) + (2.0 * x)), -(1.0 / 2.0)))

Julia

function code(F, B, x)
	return Float64(Float64(-Float64(x * Float64(1.0 / tan(B)))) + Float64(Float64(F / sin(B)) * (Float64(Float64(Float64(F * F) + 2.0) + Float64(2.0 * x)) ^ Float64(-Float64(1.0 / 2.0)))))
end

MATLAB

function tmp = code(F, B, x)
	tmp = -(x * (1.0 / tan(B))) + ((F / sin(B)) * ((((F * F) + 2.0) + (2.0 * x)) ^ -(1.0 / 2.0)));
end

Wolfram

code[F_, B_, x_] := N[((-N[(x * N[(1.0 / N[Tan[B], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]) + N[(N[(F / N[Sin[B], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[(N[(F * F), $MachinePrecision] + 2.0), $MachinePrecision] + N[(2.0 * x), $MachinePrecision]), $MachinePrecision], (-N[(1.0 / 2.0), $MachinePrecision])], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Reproduce

herbie shell --seed 2024035 
(FPCore (F B x)
  :name "VandenBroeck and Keller, Equation (23)"
  :precision binary64
  (+ (- (* x (/ 1.0 (tan B)))) (* (/ F (sin B)) (pow (+ (+ (* F F) 2.0) (* 2.0 x)) (- (/ 1.0 2.0))))))

Please file a bug report with this information.

Backtrace

list-ref: index too large for list index: 256 in: '(#f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t #t...		L	C
valid-splitindices?	/home/nightlies/herbie/zane-regimes-2/src/core/regimes.rkt	170	0
(unnamed)	.../private/arr-i.rkt	961	8
option-on-expr	/home/nightlies/herbie/zane-regimes-2/src/core/regimes.rkt	106	0
infer-splitpoints	/home/nightlies/herbie/zane-regimes-2/src/core/regimes.rkt	39	0
loop	/home/nightlies/herbie/zane-regimes-2/src/core/regimes.rkt	27	2
extract!	/home/nightlies/herbie/zane-regimes-2/src/mainloop.rkt	365	0
run-improve!	/home/nightlies/herbie/zane-regimes-2/src/mainloop.rkt	339	0
get-alternatives/report	/home/nightlies/herbie/zane-regimes-2/src/sandbox.rkt	169	0
(unnamed)	/home/nightlies/herbie/zane-regimes-2/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