Привет, ребята!!!
Я хочу вычислить обобщенное собственное разложение формы:
Lf = лямбда Af
с помощью функции scipy.sparse.linalg.eigs, но получите эту ошибку:
/usr/local/lib/python2.7/dist-packages/scipy/linalg/decomp_lu.py:61: RuntimeWarning: диагональное число 65 равно нулю. Сингулярная матрица. RuntimeWarning) ** При входе в DLASCL параметр номер 4 имел недопустимое значение
Я передаю три аргумента, диагональную матрицу, положительную полуопределенную (PSD) матрицу и числовое значение K (первые K собственных значений). Функция eigs Matlab хорошо работает с теми же входными параметрами, но в SciPy, как я понял, для вычислений с PSD мне также нужно указать параметр sigma.
Итак, мой вопрос: есть ли способ избежать установки параметра сигмы, как в MatLab, или, если нет, как подобрать значение сигмы?
С нетерпением жду советов или подсказок... Заранее спасибо!