Exception for Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2

• 36.0% of points produce a very large (infinite) output. You may want to add a precondition.

Specification

Math

\[\begin{array}{l} \\ \left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))

C

double code(double x, double y, double z, double t, double a, double b) {
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}

Python

def code(x, y, z, t, a, b):
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)

Julia

function code(x, y, z, t, a, b)
	return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b))
end

MATLAB

function tmp = code(x, y, z, t, a, b)
	tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]

Reproduce

herbie shell --seed 2024277 
(FPCore (x y z t a b)
  :name "Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2"
  :precision binary64
  (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))

Please file a bug report with this information.

Backtrace

rival-apply: contract violation expected: (or/c ^TYPE? boolean?) given: 0 in: an element of the 2nd argument of (-> rival-machine? (vectorof (or/c ^TYPE? boolean?)) (vectorof any/c)) contract from: <pkgs>/rival/main.rkt blaming: <pkgs>/herbie/core/rival.rkt (assuming the contract is correct) at: <pkgs>/rival/main.rkt:95:24		L	C
raise-blame-error	/usr/share/racket/collects/racket/contract/private/blame.rkt	346	0
loop	.../private/map.rkt	40	19
vector-map/update	/usr/share/racket/collects/racket/vector.rkt	62	0
vector-map	/usr/share/racket/collects/racket/vector.rkt	106	0
rival-apply	/home/nightlies/.local/share/racket/8.10/pkgs/rival/eval/main.rkt	65	0
(unnamed)	/usr/share/racket/collects/racket/contract/private/arrow-val-first.rkt	486	18
(unnamed)	/home/nightlies/herbie/zane-concrete-values/src/core/rival.rkt	76	4
real-apply	/home/nightlies/herbie/zane-concrete-values/src/core/rival.rkt	67	0
(unnamed)	.../private/arrow-val-first.rkt	486	18
<eval-prog-real>	/home/nightlies/herbie/zane-concrete-values/src/core/sampling.rkt	136	2
compute-local-errors	/home/nightlies/herbie/zane-concrete-values/src/core/localize.rkt	216	0
actual-errors	/home/nightlies/herbie/zane-concrete-values/src/core/explain.rkt	35	0
generate-timelines	/home/nightlies/herbie/zane-concrete-values/src/core/explain.rkt	498	0
explain!	/home/nightlies/herbie/zane-concrete-values/src/core/mainloop.rkt	323	0
run-improve!	/home/nightlies/herbie/zane-concrete-values/src/core/mainloop.rkt	45	0
get-alternatives/report	/home/nightlies/herbie/zane-concrete-values/src/api/sandbox.rkt	186	0
(unnamed)	/home/nightlies/herbie/zane-concrete-values/src/api/sandbox.rkt	282	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