On 8 April 2011 at 14:18, [email protected] wrote:
| Hmmm. I think you just need "i<n" instead of "i<=n" in your for construct.
==:-)
Indeed. I shouldn't mail before coffee. Here is the fixed version:
library(Rcpp)
library(inline)
src <- '
Rcpp::NumericMatrix xA(A);
int n=xA.nrow();
NumericVector dist(n);
for(int i=0;i<n;i++){dist[i]=1.0;}
return dist;'
cxxfun <- cxxfunction(sig=signature(A="matrix"),body=src,plugin="Rcpp")
## Tests.
doWork <- function(n) {
A <- matrix(0, n, n)
print(sum(dist <- cxxfun(A)))
NULL
}
doWork(20)
doWork(200)
doWork(2000)
Dirk
--
Dirk Eddelbuettel | [email protected] | http://dirk.eddelbuettel.com
_______________________________________________
Rcpp-devel mailing list
[email protected]
https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel
