c – 如何在反转奇异矩阵时使Armadillo函数不打印错误？

一个伙伴和我正在研究一个R包,并使用RcppArmadillo包来处理一些较重的矩阵代数.到目前为止它已经非常甜蜜,但我们在矩阵求逆方面遇到了一些问题.简而言之,程序正在搜索特定类型的矩阵,并且必须检查在循环的每次迭代中是否存在更新的矩阵的逆(逆本身也是必需的).现在我们正在使用函数inv(A,B),它返回一个布尔值,指示矩阵B是否可逆(如果不是,A设置为0x0矩阵,否则A = inv(B)).如果此函数没有打印错误,那将是很好的,因为返回的布尔值给循环提供正确进行所需的信息.并且似乎只是打印错误而不是“抛出”,如下面的程序所示：

#include <iostream>
#include <armadillo>

using namespace std;
using namespace arma;

int main(int argc,char** argv)
{
    mat A = randu<mat>(5,5);
    mat B = zeros<mat>(5,5);

    inv(A,B);

    cout << A << "error printed but not fatal" << endl;

    A = inv(B);

    cout << A << "never make it this far" << endl;

    return 0;
}

导致：

Johns-MacBook-Pro:test johnsherrill$g++ armaExample.cpp -o example -O2 -larmadillo
Johns-MacBook-Pro:test johnsherrill$./example

error: inv(): matrix appears to be singular

[matrix size: 0x0]
error printed but not fatal

error: inv(): matrix appears to be singular

terminate called after throwing an instance of 'std::runtime_error'
  what():  inv(): matrix appears to be singular
Abort trap: 6

有没有办法解决这个问题,而无需先单独检查B是否可逆？这种类型的错误也印在R中.

#define ARMA_DONT_PRINT_ERRORS
#include <armadillo>   // or #include <RcppArmadillo.h> if you're using Rcpp