Exception for Octave 3.8, jcobi/2

\[\left(\alpha > -1 \land \beta > -1\right) \land i > 0\]

Math

\[\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2} \]

FPCore

(FPCore (alpha beta i)
 :precision binary64
 (/
  (+
   (/
    (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2.0 i)))
    (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))
   1.0)
  2.0))

C

double code(double alpha, double beta, double i) {
	return (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0;
}

Fortran

real(8) function code(alpha, beta, i)
    real(8), intent (in) :: alpha
    real(8), intent (in) :: beta
    real(8), intent (in) :: i
    code = (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0d0 * i))) / (((alpha + beta) + (2.0d0 * i)) + 2.0d0)) + 1.0d0) / 2.0d0
end function

Java

public static double code(double alpha, double beta, double i) {
	return (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0;
}

Python

def code(alpha, beta, i):
	return (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0

Julia

function code(alpha, beta, i)
	return Float64(Float64(Float64(Float64(Float64(Float64(alpha + beta) * Float64(beta - alpha)) / Float64(Float64(alpha + beta) + Float64(2.0 * i))) / Float64(Float64(Float64(alpha + beta) + Float64(2.0 * i)) + 2.0)) + 1.0) / 2.0)
end

MATLAB

function tmp = code(alpha, beta, i)
	tmp = (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0;
end

Wolfram

code[alpha_, beta_, i_] := N[(N[(N[(N[(N[(N[(alpha + beta), $MachinePrecision] * N[(beta - alpha), $MachinePrecision]), $MachinePrecision] / N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2023241 
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/2"
  :precision binary64
  :pre (and (and (> alpha -1.0) (> beta -1.0)) (> i 0.0))
  (/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0) 2.0))

Backtrace

match: no matching clause for 'ival		L	C
compile-progs	/home/nightlies/herbie/arb-support/src/programs.rkt	122	0
make-search-func	/home/nightlies/herbie/arb-support/src/ground-truth.rkt	20	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