fix bounds error when upper range is too small

Author: simonbyrne

src/bivariate.jl

@@ -35,7 +35,7 @@ function tabulate(data::(RealVector, RealVector), midpoints::(Range, Range))
     for (x, y) in zip(xdata,ydata)
         kx, ky = searchsortedfirst(xmid,x), searchsortedfirst(ymid,y)
         jx, jy = kx-1, ky-1
-        if 1 <= jx <= nx && 1 <= jy <= ny
+        if 1 <= jx <= nx-1 && 1 <= jy <= ny-1
             grid[jx,jy] += (xmid[kx]-x)*(ymid[ky]-y)*ainc
             grid[kx,jy] += (x-xmid[jx])*(ymid[ky]-y)*ainc
             grid[jx,ky] += (xmid[kx]-x)*(y-ymid[jy])*ainc

src/univariate.jl

@@ -80,7 +80,7 @@ function tabulate(data::RealVector, midpoints::Range)
     for x in data
         k = searchsortedfirst(midpoints,x)
         j = k-1
-        if 1 <= j <= npoints
+        if 1 <= j <= npoints-1
             grid[j] += (midpoints[k]-x)*ainc
             grid[k] += (x-midpoints[j])*ainc
         end