Timing BigDataAlgorithms::rsvd against base::svd

BDAvSVD.R

library(TENxBrainData)
library(ggplot2)
library(tidyr)

tenx <- TENxBrainData()

ncols <- c(1,2,5,10,20,50,100,200,500,1000)

svd <- bda <- numeric(length = length(ncols))

for(i in seq_along(ncols)) {
  ncol = ncols[i]
  message(ncol)
  svd[i] <- system.time( svd1 <- svd(counts(tenx[,1:ncol])) )['elapsed']
  message(ncol)
  bda[i] <- system.time( svd2 <- BigDataAlgorithms::rsvd(counts(tenx[,1:ncol])) )['elapsed']
}

res <- data.frame(ncols, svd, bda) %>%
  tidyr::gather(key = 'algorithm', value = time, -ncols)

ggplot(res, aes(x = ncols, 
                y = time, 
                colour = algorithm)) + 
  geom_line()