Sparse Cholesky Factorization?

Hi, I am wondering whether Dask currently supports parallelization for the Cholesky factorization of a sparse, banded matrix? I know that Dask supports sparse arrays but at first glance it doesn’t seem like there is a sparse Cholesky factorization akin to the old sksparse.cholmod.

I am trying to (efficiently) sample a high-dimensional Gaussian with a banded precision matrix Q, which requires me to factorize Q = L L^T and solve the linear system L^T x = b. If Dask won’t work for this, other suggestions appreciated!