Rcpp: my distance matrix program is slower than the function in package
问题 I would like to calculate the pairwise euclidean distance matrix. I wrote Rcpp programs by the suggestion of Dirk Eddelbuettel as follows NumericMatrix calcPWD1 (NumericMatrix x){ int outrows = x.nrow(); double d; NumericMatrix out(outrows,outrows); for (int i = 0 ; i < outrows - 1; i++){ for (int j = i + 1 ; j < outrows ; j ++){ NumericVector v1= x.row(i); NumericVector v2= x.row(j); NumericVector v3=v1-v2; d = sqrt(sum(pow(v3,2))); out(j,i)=d; out(i,j)=d; } } return (out) ; } But I find my