Hi guys,<div><br></div><div>I'm trying to use Rcpp to implement an 'apply'-type function for matrices, whereby I apply a function to each row or column of a matrix. I'm testing it here with the 'mean' just to see how well we can do vs. the highly optimized rowMeans, colMeans functions, but of course any other function taking a vector and returning a scalar might fit here.</div>

<div><br></div><div>I'm wondering if I can do better. Please see the gist in the link here; you can sourceCpp it in R.<br><div><br></div><div><a href="https://gist.github.com/4561281" target="_blank">https://gist.github.com/4561281</a></div>

</div><div><br></div><div>I try to limit the amount of copying as much as possible with NumericMatrix::Column and NumericMatrix::Row.</div><div><br></div><div>Note that the Rcpp solution over columns is just as fast as colMeans, but over rows it's a fair bit slower relative to rowMeans. Is there any way I could improve this?</div>
<div><br></div><div>I plan to submit an expanded version of this to the Rcpp gallery so any advice is appreciated!</div>
<div><br></div><div>Thanks,</div><div>-Kevin</div>