C++armadillo princom()分段错误

Question

我正在运行如下代码：

#define ARMA_64BIT_WORD
#include <iostream>
#include </usr/include/armadillo>
using namespace arma;
typedef Mat<float> fmat;
typedef Col<float> fvec;
int main (int argc, char** argv)
{
    fmat A;
    A.load("matrix.txt");
    fmat coeff;
    fmat score;
    fvec latent;
    fvec tsquared;
    princomp(coeff, score, latent, tsquared, A);
 return 0;
 }

matrix.txt包含2乘46341（2行，46341列）矩阵。它运行良好，但当我将矩阵增加到2乘46342（2行，46342列）时，发生分段故障（堆芯转储）。

列的数量有限制吗？或者，我写错什么了吗？

我在Ubuntu 12.04.4 LTS（GNU/Linux 3.2.0-40-virtual x86_64）中运行，Menmory:16GB

谢谢！

Answer 1

princom（）将计算协方差矩阵（X'*X），在您的情况下，协方差矩阵的大小为46342x43642。它的大小为16.001 GB，而您只有16 GB。因此失败了。

如果使用调试符号（或不使用-DNDEBUG或类似符号）进行编译，则应显示实际的内存分配失败。

我不能推荐任何其他解决方案，因为我不知道你想解决什么问题。