multiway (version 1.0-7)

smpower: Symmetric Matrix Power

Description

Raise symmetric matrix to specified power. Default calculates symmetric square root.

Usage

smpower(X, power = 0.5, tol = NULL)

Value

Returns X raised to specified power.

Arguments

X

Symmetric real-valued matrix.

power

Power to apply to eigenvalues of X.

tol

Stability tolerance for eigenvalues.

Author

Nathaniel E. Helwig <[e-mail]>

Details

Basically returns tcrossprod(Y$vec%*%diag(Y$val^power),Y$vec) where Y = eigen(X,symmetric=TRUE).

Examples

##########   EXAMPLE   ##########

X <- crossprod(matrix(rnorm(2000),100,20))
Xsqrt <- smpower(X)         # square root
Xinv <- smpower(X,-1)       # inverse
Xisqrt <- smpower(X,-0.5)   # inverse square root