By default uses a random initialization vector that you can make deterministic using set.seed or override by specifying opts = list(initvec = ...).

eig_decomp(M, n_eigs, sym, ..., opts = list())

Arguments

M

A matrix (e.g. from the Matrix package) or a function (see eigs).

n_eigs

Number of eigenvectors to return.

sym

defunct and ignored.

...

Passed to eigs.

opts

Passed to eigs.

Value

see eigs.

Examples

eig_decomp(cbind(c(1,0,-1), c(0,1,0), c(-1,0,1)), 2)
#> $values #> [1] 2 1 #> #> $vectors #> [,1] [,2] #> [1,] 7.071068e-01 1.665335e-16 #> [2,] -1.387779e-16 -1.000000e+00 #> [3,] -7.071068e-01 -3.608225e-16 #> #> $nconv #> [1] 2 #> #> $niter #> [1] 1 #> #> $nops #> [1] 3 #>