SuiteSparse(4.5.1) 的 SPQR - 调用 cholmod_allocate_triplet 总是 returns NULL

Question

我正在尝试使用 SuiteSparse SPQR 求解线性方程组 x = A\b；我的 A 矩阵是稀疏的，它是一个矩形矩阵，所以我选择 SPQR 来解决这个问题。我使用 MS Visual Studio 2012 在 Windows 7 x64 上使用 https://github.com/jlblancoc/suitesparse-metis-for-windows 提供的那些构建了 SuiteSparse。

为了测试功能，我修改了 spqr_example 项目以在转换为稀疏矩阵之前分配三元组，而不是最初从标准输入读取输入来创建稀疏矩阵。我输入一个小的 A 和 b 矩阵进行测试。程序编译成功。我调试了程序，发现我对 cholmod_allocate_triplet() 的调用失败了，因为在这个函数的声明中它有下面的代码：

RETURN_IF_NULL_COMMON (NULL) ;

这总是 return false（即使我的 common 成功启动）。

我不想明确地更改这一行，因为我可能在某处犯了错误，或者我忘记做一些我必须做的事情，因为我是新来使用这个库的。

任何人都可以给我一些关于如何正确制作我的程序运行的建议吗？我下面的代码是根据提供的 spqr_example 修改的。非常感谢。

#include <iostream>
#include "SuiteSparseQR.hpp"

int main (int argc, char **argv)
{
cholmod_common Common, *cc ;
cholmod_sparse *A ;
cholmod_dense *X, *B, *Residual ;
double rnorm, one [2] = {1,0}, minusone [2] = {-1,0} ;
int mtype ;
// start CHOLMOD
cc = &Common ;
cholmod_l_start (cc) ;

// load A
//A = (cholmod_sparse *) cholmod_l_read_matrix (stdin, 1, &mtype, cc) ;

        // A = [ 1  0  0  0; 
        //      -1  1  0  0; ...
        //       0 -1  1  0; ...
        //       0  0 -1  1; ...
        //       0  0  0 -1];
        int row[] = {0, 1, 1, 2, 2, 3, 3, 4};
        int col[] = {0, 0, 1, 1, 2, 2, 3, 3};
        double val[] = {1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0};
        int numEq = 5;
        int numElement = 8;
        int numSol = 4;
        double b[] = {5.0, -5.0, 2.0, 1.0, 0.0};

cholmod_triplet* triplet = cholmod_allocate_triplet(5,4,5*4,0,CHOLMOD_REAL,cc);
int * triplet_i = (int *)(triplet->i);
int * triplet_j = (int *)(triplet->j);
double * triplet_x = (double *)(triplet->x);
for (int ne=0; ne<numElement; ne++)
{
    triplet_i[triplet->nnz] = row[ne];
    triplet_j[triplet->nnz] = col[ne];
    triplet_x[triplet->nnz] = val[ne];

    triplet->nnz++;
}

// Convert triplet to sparse matrix
A = cholmod_triplet_to_sparse(triplet, numElement, cc);
cholmod_free_triplet(&triplet, cc);

// B = ones (size (A,1),1)
//B = cholmod_l_ones (A->nrow, 1, A->xtype, cc) ;
B = cholmod_l_zeros(numEq, 1, CHOLMOD_REAL, cc);
for (int ne=0; ne<numEq; ne++)
{
    ((double *)(B->x))[ne] = val[ne];
}

// X = A\B
X = SuiteSparseQR<double>(A,B,cc);
//X = SuiteSparseQR <double> (A, B, cc) ;

// Print out the result
double *sol = static_cast<double *>(malloc(sizeof(X->x)));
sol = (double *)(X->x);
for (int r=0; r<numSol; r++)
{
    std::cout << "x[" << r << "] = " << sol << std::endl;
    sol++;
}

///// END HERE

// rnorm = norm (B-A*X)
Residual = cholmod_l_copy_dense (B, cc) ;
cholmod_l_sdmult (A, 0, minusone, one, X, Residual, cc) ;
rnorm = cholmod_l_norm_dense (Residual, 2, cc) ;
printf ("2-norm of residual: %8.1e\n", rnorm) ;
printf ("rank %ld\n", cc->SPQR_istat [4]) ;
// free everything and finish CHOLMOD
cholmod_l_free_dense (&Residual, cc) ;
cholmod_l_free_sparse (&A, cc) ;
cholmod_l_free_dense (&X, cc) ;
cholmod_l_free_dense (&B, cc) ;
cholmod_l_finish (cc) ;
return (0) ;
}

Answer 1

我终于找到了为什么我的程序在下面一行之后崩溃了

cholmod_triplet* 三元组 = cholmod_allocate_triplet(5,4,5*4,0,CHOLMOD_REAL,cc);

作为 cholmod_allocate_triplet() 内部调用 RETURN_IF_NULL_COMMON (NULL) 的结果，return 为假。

原因是我启动进程调用

cholmod_l_start (抄送);

这是 cholmod_start() 的 long int 版本。

为了解决这个问题，我必须调用 cholmod_l_allocate_triplet() 而不是 cholmod_allocate_triplet() 并将所有其他函数更改为使用 cholmod_l 而不是仅调用 cholmod_

SuiteSparse(4.5.1) 的 SPQR - 调用 cholmod_allocate_triplet 总是 returns NULL

SuiteSparse(4.5.1)'s SPQR - calling to cholmod_allocate_triplet always returns NULL

c++

matlab

sparse-matrix

suitesparse