Actually, there does seem to be something fishy going on here. I find that I'm able to execute temp.rolling(x=100).construct('window').mean('window').compute() successfully but not temp.rolling(x=100).mean().compute(), even though that should mostly be equivalent to the former.

Thank you so much for pointing it out. I tried the rollling.construct and it worked! I also tried it on other netcdf files and it sure solved the problem. Thank you so much for your help!

If this is caused by Dask's scheduler and there is no quick fix yet, do you think mention the rolling.construct in the Xarray document as the recommended usage would worth doing? It can help newbies like me a lot.

Cheers, Joey

da.zeros((5000, 50000), chu

Tried but same error.

python import numpy as np import xarray as xr import dask.array as da temp= xr.DataArray(da.zeros((5000, 50000), chunks=(-1,100)),dims=("x","y")) temp.rolling(x=100).mean() Like I said, I have also saved to nc file and read it from disk (as below), but still same error.

python import numpy as np import xarray as xr import dask.array as da temp= xr.DataArray(da.zeros((5000, 50000), chunks=(-1,100)),dims=("x","y")) temp.to_netcdf("temp.nc") temp.close() test = xr.open_dataarray("temp.nc",chunks={"y":100,}) test.rolling(x=100).mean()

Did you try converting np.zeros((5000, 50000) to use dask.array.zeros instead? The former will allocate 2 GB of data within each chunk

Thank you for your suggestion. Tried as you suggested, still with same error.

```python import numpy as np import xarray as xr import dask.array as da

from dask.distributed import Client

temp= xr.DataArray(da.zeros((5000, 50000)),dims=("x","y")).chunk({"y":100, }) temp.rolling(x=100).mean() ```

I have also tried saving the array to nc file and read it after that. Still rolling gives same error (with or without bottleneck and different chunks). Even though it says memory error, it doesn't consume too much memory.

Have you tried adding more chunking, e.g., along the x dimension? That’s that usual recommendation if you’re running out of memory.

Hi Shoyer,

Thanks for your reply and help. However, I have tried various chunks along each and both dimension (like 200 on x dimension, 100 on y dimension; or larger chunks like 2000 on y dimension), it doesn't work.

In both a ubuntu machine with 100 Gb memory and a local windows10 machine, it simply crashed in couple of seconds. Even though it says memory error, the code does not use much memory at all. Also even with the one dimension setup, the temp.data shows that each chunk only takes 4 mb memory (which makes me think it might be too small and then used larger chunks). I also used a new conda environment with clean install of just the necessary libraries, and the problem is still there.

Here is the neat new environment under which I tried again but gives the same errors,

Output of xr.show_versions()

By the way, the above code seems to work ok with previous 0.12.1 version Xarray and bottleneck.

Cheers, Joey