I was using scidbpy’s mean() function on some 2-d arrays of different sizes, and here are some interesting results. Numpy and SciDB on [100,100] are pretty much the same; [1000,1000] starts showing some discrepancy; [4094,2046] is completely messed up (see below). Any ideas?

a.shape = (100, 100)

Numpy: a.mean = 286.45

SciDB: mean = 286.451162834

a.shape = (1000, 1000)

Numpy: a.mean = 287.57

Numpy: a_new.mean = 287.57

SciDB: mean = 287.490968998

a.shape = (4094, 2046)

Numpy: a.mean = 273.766

SciDB: mean = 291.174422006