Bryan O'Sullivan avatar Bryan O'Sullivan committed 016af22

incompleteBeta: check x for NaN (bug 8)

Comments (0)

Files changed (1)

Numeric/SpecFunctions.hs

                 -> Double -- ^ /x/, must lie in [0,1] range
                 -> Double
 incompleteBeta_ beta p q x
-  | p <= 0 || q <= 0 = error "p <= 0 || q <= 0"
-  | x <  0 || x >  1 = error "x <  0 || x >  1"
-  | x == 0 || x == 1 = x
+  | p <= 0 || q <= 0            = error "p <= 0 || q <= 0"
+  | x <  0 || x >  1 || isNaN x = error "x out of [0,1] range"
+  | x == 0 || x == 1            = x
   | p >= (p+q) * x   = incompleteBetaWorker beta p q x
   | otherwise        = 1 - incompleteBetaWorker beta q p (1 - x)
 
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.