为什么这两个函数给出不同的结果？

Question

我使用 RcppArmadillo 包定义了两个函数并将它们保存到文件 cxxFuns.cpp 中。 f01 和 f02 之间的唯一区别是 V(0, 0):

的位置

#include <RcppArmadillo.h>
using namespace Rcpp;

// [[Rcpp::export]]
arma::mat f01 (arma::mat x) {
    unsigned int p0 = x.n_cols ;
    unsigned int iR = 0, iC = 0 ;
    arma::mat V(3, 3) ; V.fill(NA_REAL) ;
    for (iR = 0; iR < p0; iR++) { 
        V(0, 0) = arma::sum(x.col(iR) % x.col(iR)) ;
        for (iC = iR+1; iC < p0; iC++) {
            ;
        }
    }
    
    return V ;
}

// [[Rcpp::export]]
arma::mat f02 (arma::mat x) {
    unsigned int p0 = x.n_cols ;
    unsigned int iR = 0, iC = 0 ;
    arma::mat V(3, 3) ; V.fill(NA_REAL) ;
    for (iR = 0; iR < p0; iR++) { 
        for (iC = iR+1; iC < p0; iC++) {
            V(0, 0) = arma::sum(x.col(iR) % x.col(iR)) ;
        }
    }
    
    return V ;
}

据我了解，函数 f01 和 f02 应该给出相同的结果。但是，测试未显示预期结果。

rm(list=ls())
set.seed(2020)
Rcpp::sourceCpp('cxxFuns.cpp')
x <- matrix(rnorm(100*10), 10)
(egg01 <- f01(x))
         [,1] [,2] [,3]
[1,] 12.78607   NA   NA
[2,]       NA   NA   NA
[3,]       NA   NA   NA

(egg02 <- f02(x))
         [,1] [,2] [,3]
[1,] 14.80855   NA   NA
[2,]       NA   NA   NA
[3,]       NA   NA   NA

发生了什么事？

Answer 1

第一个块的最后一次执行是在iR = p0-1.

时

第二块的最后一次执行是在iC=p0-1.
时由于 iC 以 iR+1 开始，您最后一次执行是针对 iR=p0-2.

在调试器下计算 V(0,0) 或运行之前打印 iR 的值应该会立即明确这一点。

为什么这两个函数给出不同的结果？

Why do the two functions give different results?

c++

r

rcpp

rcpparmadillo