When the dask library is installed as a dependency (through the result_caching module, see here), calls to scipy.stats.pearsonr become significantly slower. This slowdown notably impacts performance during ceiling calculations, which rely on result_caching. Removing or bypassing dask in this workflow returns pearsonr to its expected speed. We need to investigate why dask is causing this overhead and how to mitigate it without sacrificing functionality in result_caching.