Exception for b parameter of renormalized beta distribution

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

Specification

\[\left(0 < m \land 0 < v\right) \land v < 0.25\]

Math

\[\begin{array}{l} \\ \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right) \end{array} \]

FPCore

(FPCore (m v) :precision binary64 (* (- (/ (* m (- 1.0 m)) v) 1.0) (- 1.0 m)))

C

double code(double m, double v) {
	return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m);
}

Fortran

real(8) function code(m, v)
    real(8), intent (in) :: m
    real(8), intent (in) :: v
    code = (((m * (1.0d0 - m)) / v) - 1.0d0) * (1.0d0 - m)
end function

Java

public static double code(double m, double v) {
	return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m);
}

Python

def code(m, v):
	return (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m)

Julia

function code(m, v)
	return Float64(Float64(Float64(Float64(m * Float64(1.0 - m)) / v) - 1.0) * Float64(1.0 - m))
end

MATLAB

function tmp = code(m, v)
	tmp = (((m * (1.0 - m)) / v) - 1.0) * (1.0 - m);
end

Wolfram

code[m_, v_] := N[(N[(N[(N[(m * N[(1.0 - m), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision] - 1.0), $MachinePrecision] * N[(1.0 - m), $MachinePrecision]), $MachinePrecision]

Reproduce

herbie shell --seed 2024039 
(FPCore (m v)
  :name "b parameter of renormalized beta distribution"
  :precision binary64
  :pre (and (and (< 0.0 m) (< 0.0 v)) (< v 0.25))
  (* (- (/ (* m (- 1.0 m)) v) 1.0) (- 1.0 m)))

Please file a bug report with this information.

Backtrace

member: not a proper list in: 1		L	C
gather-multiplicative-terms	/home/nightlies/herbie/sound-rules/src/core/reduce.rkt	125	0
simplify-node*	/home/nightlies/herbie/sound-rules/src/core/reduce.rkt	84	0
hash-ref!	/usr/local/racket/collects/racket/private/more-scheme.rkt	377	2
hash-ref!	/usr/local/racket/collects/racket/private/more-scheme.rkt	377	2
loop	.../private/map.rkt	40	19
loop	.../private/map.rkt	40	19
f19	.../match/compiler.rkt	559	40
hash-ref!	/usr/local/racket/collects/racket/private/more-scheme.rkt	377	2
approximate	/home/nightlies/herbie/sound-rules/src/core/taylor.rkt	7	0
taylor-expr	/home/nightlies/herbie/sound-rules/src/patch.rkt	89	0
taylor-alt	/home/nightlies/herbie/sound-rules/src/patch.rkt	96	0
gen-series!	/home/nightlies/herbie/sound-rules/src/patch.rkt	109	0
patch-table-run	/home/nightlies/herbie/sound-rules/src/patch.rkt	257	0
(unnamed)	/usr/local/racket/collects/racket/contract/private/arrow-val-first.rkt	486	18
run-iter!	/home/nightlies/herbie/sound-rules/src/mainloop.rkt	288	0
mutate!	/home/nightlies/herbie/sound-rules/src/mainloop.rkt	340	0
run-improve!	/home/nightlies/herbie/sound-rules/src/mainloop.rkt	326	0
get-alternatives/report	/home/nightlies/herbie/sound-rules/src/sandbox.rkt	169	0
(unnamed)	/home/nightlies/herbie/sound-rules/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