Closed pneuvial closed 7 years ago
Reproducible example:
W <- matrix(1:6, 3, 2)
L <- 5
p <- ncol(W)
Dm1 <- -Matrix::bandSparse(L, L, k=0:(L-1))
X <- kronecker(Dm1, W)
Can be fixed (e.g.) by converting W to a sparse matrix:
W <- matrix(1:6, 3, 2)
W <- as(W, "sparseMatrix")
L <- 5
p <- ncol(W)
Dm1 <- -Matrix::bandSparse(L, L, k=0:(L-1), giveCsparse = FALSE)
X <- kronecker(Dm1, W)
Fixed likewise in [develop cb8a136]