Scipy稀疏矩阵求解输出dgstrf info 10200且结果全NaN的问题

1. dgstrf info 10200的含义 #

该信息来自Scipy稀疏求解器依赖的SuperLU库：

• dgstrf是SuperLU中负责双精度稀疏矩阵LU分解的核心函数
• info 10200表示待求解的稀疏矩阵AMatrixT是奇异矩阵（不可逆），无法完成LU分解以求解线性方程组，因此返回的foo全为NaN。
• 结合你的dist维度为(10200,1)，说明方程组AMatrixT * x = dist要么无解，要么有无穷多解，本质原因是矩阵的秩不足。

2. 如何阻止该输出 #

该信息是SuperLU库直接输出到标准错误流的，属于底层库打印内容，可通过以下两种方式抑制：

方法1：重定向错误流 #

使用Python标准库的contextlib.redirect_stderr将底层错误输出导向空设备：

from contextlib import redirect_stderr
import os
from scipy.sparse.linalg import solve as spsolve  # 注：原代码导入路径有误，正确路径为scipy.sparse.linalg

with redirect_stderr(open(os.devnull, 'w')):
    foo = spsolve(AMatrixT, dist)

方法2：切换求解器（推荐） #

Scipy的spsolve支持通过use_umfpack参数切换到UMFPACK求解器，该求解器默认不会输出这类底层信息：

from scipy.sparse.linalg import solve as spsolve

foo = spsolve(AMatrixT, dist, use_umfpack=True)

注：若环境未编译UMFPACK组件，可 fallback 到方法1。

额外建议 #

矩阵奇异说明你的建模或矩阵构造可能存在问题，建议：

• 检查AMatrixT是否存在线性相关的行/列
• 验证问题逻辑是否存在约束不足、重复约束等情况
• 用scipy.sparse.linalg.norm计算矩阵条件数，或scipy.sparse.linalg.eigs求解最小特征值，进一步验证矩阵奇异性

内容的提问来源于stack exchange，提问作者FooBar