Armadillo之计算矩阵的行列式（determinant）

计算矩阵的行列式很简单，用det方法或是log_det方法

1 det(A)

  若是A不是方阵的（square），将抛出std::logic_error异常

例：

mat m = "3,2,4;1,-2,3;2,3,2;";
double d = det(m);
cout << d << endl;

运行结果是-3

2 log_det(value, sign，A)

    文档里推荐当矩阵A比较大时，使用本函数来代替det函数（估计会加快计算速度）

    det(A)=exp(value)*sign

    若是A不是方阵的（square），将抛出std::logic_error异常

例：

mat m = "3,2,4;1,-2,3;2,3,2;";
double value, sign;
log_det(value, sign, m);
cout << exp(value)*sign << endl;

运行结果是-3

    咱们知道判断一个方阵（square matrix）是不是奇异（singular）的能够看此方阵的行列式（determinant）是否为0，这对于精确的数值计算而言是没有问题的，但咱们的计算机使用的是有限数字系统，老是存在偏差的，因此一个矩阵即便是奇异，它的行列式用计算机计算出来的结果也有可能不为0，但老是比较接近0。但通常而言,det(A)的值并非判断A是否接近奇异的好的判断标准。