Плотное обновление Холецкого в Python

Может ли кто-нибудь указать мне на библиотеку / код, позволяющий мне выполнять низкоуровневые обновления декомпозиции Холецкого в python (numpy)? Matlab предлагает эту функциональность в виде функции под названием cholupdate. LINPACK также имеет эту функцию, но (насколько мне известно) еще не был перенесен в LAPACK и, следовательно, недоступен, например, в scipy.

Я узнал, что scikits.sparse предлагает аналогичную функцию на основе CHOLMOD, но мои матрицы плотные.

Есть ли какой-нибудь код для python с функцией cholupdate, совместимый с numpy?

Спасибо!