/*Gauss*/

m=ones(200000,1);

ret=0;
starttime=date;

for i (1,20,1);

ret=ret+sumc(sumr(m));

endfor;

endtime=date;

print etstr( ethsec(starttime,endtime) );

/*OX*/

#include

main(){

decl tim,i,m,ret,k;

m = ones(200000,10);

ret =0;
tim = timer();

for (i=0;i<20;++i){

ret=ret+sumc(sumr(m));

}

print(”time lapsed: “, timespan(tim, timer()), “n”);

}
#R

# File name: loop_rowSums.r
my.loop <- 20

m.dim <- list(nrow = 200000, ncol = 10)

m <- matrix(1, nrow = m.dim$nrow, ncol = m.dim$ncol)

ret <- 0
start <- Sys.time()

for(k in 1 : my.loop){

ret <- ret + sum(rowSums(m))

}

Sys.time() - start

==========================

Speed:

OX(0.18)

Gauss(0.06)

R(0.66)

http://www.math.ncu.edu.tw/~chenwc/R_note/index.php?item=loop