我正在处理需要进行元素级矩阵乘法的代码。由于代码涉及到一些昂贵的循环操作,因此我尝试在Rcpp中实现这个功能。我对Rcpp相对陌生,可能会有所遗漏,但是我目前无法使元素级矩阵乘法正常工作。
// [[Rcpp::export]]
NumericMatrix multMat(NumericMatrix m1, NumericMatrix m2) {
NumericMatrix multMatrix = m1 * m2 // How can this be implemented ?
}
我可能错过了一些非常微不足道的东西,并想问是否有任何方法可以做到这一点(除了使用循环迭代每个元素并相乘之外)。
提前感谢。