Python Scipy.sparse RuntimeError: коэффициент ровно единичный ⇐ Python

[Anonymous]

Я пытался найти факторизацию нижнего-верхнего (LU) для разреженной матрицы (размер 159990 на 159990). Ниже я прикрепляю три строки кода Python, которые я использовал для LU-декомпозиции

Код: Выделить всё

from scipy.sparse import csc_matrix, linalg as sla
interior_stiff = CSC_matrix(159990 by 159990)
LU = sla.splu(interior_stiff, options=dict(SymmetricMode=True))

Свойства этой матрицы включают:

• Симметричность
• Диагональные элементы являются самыми большими среди соответствующих столбцов матрицы.
• Несингулярная матрица

При запуске кода я получил следующую ошибку в Spyder Графический интерфейс,

Код: Выделить всё

  File "", line 1, in 
    runfile('C:/temp/new_LU/Run_1_using_LU.py', wdir='C:/temp/new_LU')
  File "C:\ProgramData\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 866, in runfile
    execfile(filename, namespace)
  File "C:\ProgramData\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 102, in execfile
    exec(compile(f.read(), filename, 'exec'), namespace)
  File "C:/temp/new_LU/Run_1_using_LU.py", line 574, in 
    LU= sla.splu(interior_stiff,options=dict(SymmetricMode=True))
  File "C:\ProgramData\Anaconda3\lib\site-packages\scipy\sparse\linalg\dsolve\linsolve.py", line 261, in splu
    ilu=False, options=_options)
RuntimeError: Factor is exactly singular

Может ли кто-нибудь предложить возможное решение этой проблемы?
Любая помощь будет оценена по достоинству.

Подробнее здесь: https://stackoverflow.com/questions/456 ... y-singular