C ++中的矩阵/向量与R math独立库？

Question

所有，

我一直在使用C ++中的R math独立库。 我非常希望能够生成随机数并使用R中我熟悉的分布函数。我的问题是：是否可以使用R中可用的矩阵运算（乘法，转置，逆，Chol等）一个独立的图书馆？ 我在Rmath.h中没有看到它们。

如果可以在独立的C ++代码中使用矩阵运算，那么R API似乎成为完美的开源计算引擎。 人们是否以这种方式使用它？

非常感谢您的指导和建议！

Answer 1

大多数（如果不是全部）你提到的东西都是由R链接的LAPACK或BLAS Fortran代码提供的，而不是R为自己提供新代码的东西。

如果您对使用C ++和R感兴趣，请查看Dirk Eddelbuettel和Romain Francois的Rcpp软件包 。 Dirk还编写了RcppArmadillo软件包作为Armadillo C ++线性代数库的接口，可以执行您提到的矩阵运算。

这是否有用将取决于您是否想要编写在R外部访问的单独C ++代码，或者在R中访问您访问的C ++代码.Rcpp（大大）促进后者。 如果是以前的情况，请直接看看Armadillo。

Answer 2

我认为原始问题从错误的一端开始。 R中没有C ++，R本身也没有C ++ API。

因此，如果要访问R矩阵函数，则必须通过R API - 这是一个C层而非 C ++。

Gavin在他的回答（和评论）和Josh中都是正确的：你可能想要像Armadillo这样的高级矩阵代数。 Armadillo（就像相关的库一样）最终进入BLAS和LAPACK ---就像R一样。[Doug Bates通常提醒我们，与旋转分解有一个重要的区别; 这在fastLm()实现和示例中有所涉及。 ]

最后，Rcpp可以帮助解决这两种情况：

• 是否要通过调用您编写的C ++代码来扩展R，其中Rcpp可以轻松地来回传递对象，或者

• 是否要使用RInside将R嵌入到C ++应用程序中，因为Rcpp再次提供了粘合剂

rcpp-devel列表是进行更深入讨论和示例的好地方。