我试图找到Sympy中某些矩阵的零空间的基础。
我有以下代码:
def get_null_space(vectors_list):
A = sympy.Matrix(vectors_list).T
return A.nullspace()
vectors_list
的每个元素都是SymPy.core.numbers
的元组,元素是1,-1,0。
矩阵形状为(384,120)。
运行代码时,使用所有内存(60GB),卡住并最终导致内核崩溃。
设置
os.environ["SYMPY_USE_CACHE"]="no"
和
os.environ["SYMPY_INT_TRACE"]="no"
没有做任何改变。
在数值上使用NumPy的svd
,零空间是零维的。
我使用了sympy的最新版本。如何以符号方式计算零空间?