C++ 尖点稀疏库给出奇怪的结果

我尝试在CUSP中执行以下操作：

A=[
    1,1,0,0;
    2,2,2,0;
    0,3,3,3;
    0,0,4,4];
B=[1,1,1,1]';
disp(mldivide(A,B));

那是

X=[0.9167,0.0833,-0.5000,0.7500]

另一方面，我从CUSP得到了一个奇怪的答案

#include <cusp/dia_matrix.h>
#include <cusp/krylov/cg.h>
#include <cusp/print.h>
int main()
{
    cusp::dia_matrix<int,float,cusp::host_memory> A(4,4,10,3);
    A.diagonal_offsets[0] = -1;
    A.diagonal_offsets[1] =  0;
    A.diagonal_offsets[2] =  1;
    for (int i = 0;i <3;i++)
    {
        for (int q = 0 ;q < A.num_cols;q++)
        {
            A.values(q,i)=q+1;
        }
    }
    //copy
    cusp::dia_matrix<int,float,cusp::device_memory> AA = A;
    cusp::array1d<float,cusp::device_memory> BB(A.num_rows,1);
    cusp::array1d<float,cusp::device_memory> XX(A.num_rows,0);
    cusp::print(AA);
    cusp::print(XX);
    cusp::print(BB);
    cusp::krylov::cg(AA,XX,BB);\
    cusp::print(XX);
    return 0;
}

#包括
#包括
#包括
int main（）
{
尖点：直径矩阵A（4,4,10,3）；
A.对角线_偏移[0]=-1；
A.对角线_偏移量[1]=0；
A.对角线_偏移量[2]=1；
for（int i=0；i仅在对称正定矩阵中有效。您的矩阵不是对称的。这就是为什么它不是（也不能）对称的原因产生一个有效的解决方案。要么使用一个适当的、条件良好的SPD矩阵，要么使用不同的数值方法。
这相当令人失望。你知道解决一个非对称但仍然是对角的系统的方法吗？同样，如果系统需要对称，为什么数据结构允许不对称？biggstab或GMRES是最好的方法用于非SPD稀疏系统的常用Krylov子空间方法。CUSP中的数据结构和数值方法在很大程度上是独立的。CUSP基本上是一组GPU优化数据结构和一个高性能稀疏点积，并在该点积的基础上构建了一系列算法。它是（与所有数值库一样）一个先决条件是，在你尝试使用它们之前，你必须理解基本的数学知识。买主要当心，所有这些……我认为在比克斯塔布中，它仍然必须是积极的？
sparse matrix <4, 4> with 10 entries
              0              0              1
              0              1              1
              1              0              2
              1              1              2
              1              2              2
              2              1              3
              2              2              3
              2              3              3
              3              2              4
              3              3              4
array1d <4>
             0
             0
             0
             0
array1d <4>
             1
             1
             1
             1
array1d <4>
      -39.9938
       -53.436
       87.9025
      -30.1429