Re: [R] add diagonal to matrix

Hi, I have another method: set.seed(48) n <- 4500 m1 <- matrix(sample(1:40,n*(n-1),replace=TRUE),ncol=n) m2 <- vector("numeric",length=n^2) system.time({sqn <- seq.int(n);m2[!seq_along(m2)%in%(sqn+(sqn-1)*n)] <- as.vector(m1);m3<-array(m2,dim=c(n,n))}) # user system elapsed # 1.70 0.19 1.92

Re: [R] add diagonal to matrix

Hi, I have another method: set.seed(48) n <- 4500 m1 <- matrix(sample(1:40,n*(n-1),replace=TRUE),ncol=n) m2 <- vector("numeric",length=n^2) system.time({m2[!m2%in%seq.int(n)^2] <- as.vector(m1);m3<-array(m2,dim=c(n,n))}) # user system elapsed # 1.43 0.28 1.87 Best Regards! Fechy - Ori

Re: [R] add diagonal to matrix (sorry for having a little bug in last email)

Hi, I have another method: set.seed(48) n <- 4500 m1 <- matrix(sample(1:40,n*(n-1),replace=TRUE),ncol=n) m2 <- vector("numeric",length=n^2) system.time({m2[!seq_along(m2)%in%seq.int(n)^2] <- as.vector(m1);m3<-array(m2,dim=c(n,n))}) # user system elapsed # 1.62 0.39 2.04 Best Regards! Fechy