Snakemake and dask incompatibility

Snakefile

import numpy as np
import dask.array as da


rule all:
    output: "out.hdf5"
    run:
        n = 1000
        s = 1001

        X = np.random.rand(s, n)
        X -= X.mean(axis=0)[None,:]
        da.from_array(X, chunks = (100,100))\
            .to_hdf5(output[0], "/out")

test.py

import numpy as np
import dask.array as da


n = 1000
s = 1001

X = np.random.rand(s, n)
X -= X.mean(axis=0)[None,:]
da.from_array(X, chunks = (100,100))\
    .to_hdf5("out.hdf5", "/out")