Question

以下代码：

#include <iostream>
#include <armadillo>

using namespace std;
using namespace arma;

int main()
{
        sp_mat A = speye<sp_mat>(5,5);
        rowvec s1 = max(A,0);

        return 0;
}

给出以下编译时错误：

benchmark.cpp: In function ‘int main()’:
benchmark.cpp:11:21: error: conversion from ‘arma::enable_if2<true, const arma::SpOp<arma::SpMat<double>, arma::spop_max> >::result {aka const arma::SpOp<arma::SpMat<double>, arma::spop_max>}’ to non-scalar type ‘arma::rowvec {aka arma::Row<double>}’ requested
  rowvec s1 = max(A,0);
                     ^
make: *** [all] Error 1

对于稀疏矩阵的min，sum和其他操作也是如此，而它们对于密集矩阵非常有效。我在这里做错了吗？

Answer 1

稀疏矩阵的最大运算将产生稀疏矩阵。

将您的代码更改为：

sp_mat A = speye<sp_mat>(5,5);
sp_mat s1 = max(A,0);

Answer 2

为了计算每一行的总和，我使用了矩阵乘法：

sp_mat A = sprandu<sp_mat>(5, 5);
mat sumRows = A * ones(A.n_cols, 1);

在c ++上的armadillo中，稀疏矩阵上的sum（<sp_mat>，<dim>）不起作用</dim> </sp_mat>

2 个答案: