sla.chol2inv {RScaLAPACK}R Documentation

Invert a Matrix from the Cholesky Decomposition (Using ScaLAPACK)

Description

Invert a symmetric, positive definite, square matrix from its Choleski decomposition using ScaLAPACK.

Usage

sla.chol2inv(A, NPROWS, NPCOLS, MB, RFLAG, SPAWN)

Arguments

A a matrix. The upper triangle should hold the Choleski Decomposition to compute the inverse of.
NPROWS Number of Process Rows in the Process Grid.
NPCOLS Number of Process Cols in the Process Grid.
MB Block Size.
RFLAG Flag saying whether the Process Grid should be released after computation.
SPAWN Flag saying whether a new Process Grid should be spawned.

Details

If the number of processor rows and columns are both zero, one processor is used. If the number of processor rows is nonzero and the number of processor columns is zero, then the processor rows is taken to be a number of processors, and a grid is made accordingly.

Value

The inverse of the decomposed matrix is returned.

References

http://mathworld.wolfram.com/CholeskyDecomposition.html

http://www.netlib.org/scalapack/

See Also

chol2inv for the normal R version.

Examples

library(RScaLAPACK)
rnorm(100)->x
dim(x)=c(10,10)
# First, generate a symmetric, positive definite matrix
sla.chol2inv(x)->y
print(y)
# Invert
sla.solve(y)->y2
# Perform a Choleski decomposition
sla.chol(y2)->x


[Package RScaLAPACK version 0.5.1 Index]