Python Scipy.sparse运行时错误:因子完全是单数

Python Scipy.sparse运行时错误:因子完全是单数,python,Python,我试图找到稀疏矩阵(大小159990×159990)的LU分解。下面是我在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)) 该矩阵的性质包括: 对称的 对角线元素是对应矩阵列中最大的元素 非奇

我试图找到稀疏矩阵(大小159990×159990)的LU分解。下面是我在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 GUI中遇到以下错误:

      File "<ipython-input-1-7243a2294501>", line 1, in <module>
    
        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 <module>
    
        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
    
    文件“”,第1行,在
    运行文件('C:/temp/new\u LU/Run\u 1\u使用\u LU.py',wdir='C:/temp/new\u LU')
    文件“C:\ProgramData\Anaconda3\lib\site packages\spyder\utils\site\site customize.py”,第866行,在runfile中
    execfile(文件名、命名空间)
    文件“C:\ProgramData\Anaconda3\lib\site packages\spyder\utils\site\sitecustomize.py”,第102行,在execfile中
    exec(编译(f.read(),文件名,'exec'),命名空间)
    文件“C:/temp/new_LU/Run_1_using_LU.py”,第574行,在
    LU=sla.splu(interior_stiff,options=dict(SymmetricMode=True))
    splu中的文件“C:\ProgramData\Anaconda3\lib\site packages\scipy\sparse\linalg\dsolve\linsolve.py”,第261行
    ilu=假,选项=_选项)
    RuntimeError:因子完全是单数
    
    有人能提出一个可能的解决方案吗

    任何帮助都将不胜感激

    谢谢,


    Paul Thomas可能与此相关?谢谢您的建议,但哪个=LM选项仅适用于eigh(),而不适用于splu()。因此,这在这方面无法帮助我。使用“”解决问题