When running irlba on a big.matrix as in the example from the vignette (on
Windows R 2.15.2 64 bit) it returns with an error (see below). irlba works fine
on a regular R matrix
> library(bigalgebra)
> library(irlba)
Loading required package: Matrix
Loading required package: lattice> matmul <- function(A, x, transpose=FALSE)
+ {
+ if(transpose) {return(t( t(x) %*% A))}# i.e., t(A) %*% x
+ return (A %*% x)
+ }> set.seed(1)
> A <- matrix(rnorm(100),10)
> a <- as.big.matrix(A)
> La <- irlba(a, nu=2, nv=2, matmul=matmul)
Error in W[, j] <- matmul(A, V[, j, drop = FALSE]) :
number of items to replace is not a multiple of replacement
length> LA <- irlba(A, nu=2, nv=2, matmul=matmul)
> LA$d
[1] 5.154081 4.816501
Any clues or is this a bug with the package integration?
Cheers,
Gavan
[[alternative HTML version deleted]]
