Exception for a parameter of renormalized beta distribution

• 42.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 m \end{array} \]

FPCore

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

C

double code(double m, double v) {
	return (((m * (1.0 - m)) / v) - 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) * m
end function

Java

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

Python

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

Julia

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

MATLAB

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

Wolfram

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

Reproduce

herbie shell --seed 2024089 
(FPCore (m v)
  :name "a 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) m))

Please file a bug report with this information.

Backtrace

dict-ref: no value for key: '((*.f64 #s(literal -1 binary64) (/.f64 m v)) /.f64 m (neg.f64 v)) in: '((((*.f64 #s(literal -1 binary64) m) neg.f64 m) (*.f64 #s(literal -1 binary64) m) (Rewrite<= neg-mul-1_binary64 (neg.f64 m))) (((*.f64 m (-.f64 (/.f64 m v) #s(literal 1 binary64))) -.f64 (/.f64 (pow.f64 m #s(literal 2 binary64)) v) m) (*.f64 m (-.f64 (...		L	C
canonicalize-proof	/home/nightlies/herbie/simplify-patch/src/soundiness.rkt	30	0
hash-ref!	/usr/local/racket/collects/racket/private/more-scheme.rkt	377	2
add-soundiness-to	/home/nightlies/herbie/simplify-patch/src/soundiness.rkt	55	0
loop	.../private/map.rkt	40	19
loop	.../private/map.rkt	40	19
alt-map	/home/nightlies/herbie/simplify-patch/src/alternative.rkt	38	0
loop	.../private/map.rkt	40	19
alt-map	/home/nightlies/herbie/simplify-patch/src/alternative.rkt	38	0
add-soundiness	/home/nightlies/herbie/simplify-patch/src/soundiness.rkt	87	0
extract!	/home/nightlies/herbie/simplify-patch/src/mainloop.rkt	83	0
run-improve!	/home/nightlies/herbie/simplify-patch/src/mainloop.rkt	51	0
get-alternatives/report	/home/nightlies/herbie/simplify-patch/src/sandbox.rkt	171	0
(unnamed)	/home/nightlies/herbie/simplify-patch/src/sandbox.rkt	266	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