R kronecker-methods “Matrix”包中函數“kronecker()”的方法

說明

計算繼承自 "Matrix" 的對象的克羅內克積。

為了保持稀疏性，我們對待0 * NA作為0，不作為NA通常在R(並且用於base函數kronecker)。

方法

克羅內克

signature(X = "Matrix", Y = "ANY") ......

克羅內克

signature(X = "ANY", Y = "Matrix") ......

克羅內克

signature(X = "diagonalMatrix", Y = "ANY") ......

克羅內克

signature(X = "sparseMatrix", Y = "ANY") ......

克羅內克

signature(X = "TsparseMatrix", Y = "TsparseMatrix")…………

克羅內克

signature(X = "dgTMatrix", Y = "dgTMatrix") ......

克羅內克

signature(X = "dtTMatrix", Y = "dtTMatrix") ......

克羅內克

signature(X = "indMatrix", Y = "indMatrix") ......

例子

(t1 <- spMatrix(5,4, x= c(3,2,-7,11), i= 1:4, j=4:1)) #  5 x  4
(t2 <- kronecker(Diagonal(3, 2:4), t1))               # 15 x 12

## should also work with special-cased logical matrices
l3 <- upper.tri(matrix(,3,3))
M <- Matrix(l3)
(N <- as(M, "nsparseMatrix")) # "ntCMatrix" (upper triangular)
N2 <- as(N, "generalMatrix")  # (lost "t"riangularity)
MM <- kronecker(M,M)
NN <- kronecker(N,N) # "dtTMatrix" i.e. did keep
NN2 <- kronecker(N2,N2)
stopifnot(identical(NN,MM),
          is(NN2, "sparseMatrix"), all(NN2 == NN),
          is(NN, "triangularMatrix"))