矩阵指数与Armadillo

Question

我目前正在用我自己的C++/Mex代码和Matlab开发，但我的项目越来越大，我正在考虑切换到一个适当的线性代数库。我读过关于Armadillo的一些很好的东西，但是我找不到我的项目所需要的一些基本功能。

我理解Armadillo链接到LAPACK和BLAS库，但我在Armdaillo API和LAPACK函数中找不到矩阵指数函数。

有人能告诉我是否有一个附加的计算矩阵指数的Armadillo吗？如果是这样的话，一个简短的示例代码将是非常感谢的。

Answer 1

这已作为expmat添加到最新版本中，请参阅http://arma.sourceforge.net/docs.html#expmat。

Answer 2

矩阵指数是Matlab所拥有的。因此，八度实现了它。因此，其他自由软件项目研究了八度拥有什么，并通过借用这个实现重新实现了它。

我通过R包和[医]鲤鱼包(我是该包的合著者)在RcppArmadillo和[医]鲤鱼上工作了很多。在最近的一篇文章中，我需要expm()，并从R包外接那里借用它供[医]鲤鱼使用。

代码是这样的：

arma::mat expm(arma::mat x) {
    arma::mat z(x.n_rows, x.n_cols);
    (*expmat)(x.begin(), x.n_rows, z.begin(), Ward_2);
    return z;
}

但是，它当然取决于这样一个事实，即我从expmat获得指向R包外接的函数指针。包含枚举就在吉特布的完整文件Ward_2。

Answer 3

这是John Burkardt的c/c++实现的一个港口，它是Armadillo的19种可疑之一.

https://gist.github.com/tesch1/0c03e43885cd66eceabe