可能是因为SQP算法必须提取经常约束矩阵的行子集。这可以更快的完全比稀疏矩阵(见下文)。另外,即使你的约束是稀疏的,没有理由认为算法的二次规划子问题。所以这是可疑的好处在稀疏约束独自的形式。
一个= sprand (4000、4000、4000/4e5) + speye (4000);
房颤=全(A);
抽搐
为i = 1:尺寸(1)
=(我:);
结束
toc
抽搐
为i = 1:尺寸(1)
= Af(我:);
结束
toc
最后,即使你的问题是稀疏的一切,任何矩阵反演中必须经常会更快完成完整形式比稀疏小规模问题,这sqp预计问题。的例子,
b =兰德(4000 1);
抽搐
一个\ b;
toc
抽搐
房颤\ b;
toc
