numsamples = int(np.prod(shp))