相关疑难解决方法(0)

我有这个问题,这需要解决X在AX=B.A是15000 x 15000的顺序,是稀疏和对称的.B是15000 X 7500并且不稀疏.解决X的最快方法是什么？

我可以想到两种方式.

1. 最简单的方式, X = A\B

2. 使用for循环,

   invA = A\speye(size(A))
   for i = 1:size(B,2)
       X(:,i) = invA*B(:,i);
   end
   

   Run Code Online (Sandbox Code Playgroud)

有没有比上面两个更好的方法？如果没有,我提到的两者之间哪一个最好？

是否有任何库或免费提供的代码可以完全在GPU上计算小（6x6），双精度矩阵的行列式？