Как ускорить вычисления корреляционной матрицы .corr()?

Строю корреляционную матрицу большого датасета

# featureslst датафрейм признаков
# features_work список имен из 10000 признаков 
# создание корреляционной матрицы
corr_matrix = features_work[featureslst].corr()

Итог 35 минут вычислений с загрузкой процессора 4%. Как ускорить вычисления?