sla.solve {RScaLAPACK}R Documentation

Solve a System of Equations using Scalapack functions

Description

This generic function solves the equation a %*% x = b for x using the ScaLAPACK library, where b can be either a vector or a matrix.

Usage


sla.solve(A, B=NULL, NPROWS=0, NPCOLS=0, MB=16,RFLAG=1, SPAWN=1 )

Arguments

A

A square numeric matrix containing the coefficients of the linear system.

B

A numeric vector or matrix giving the RHS of the linear system.

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.

Author(s)

Nagiza Samatova (samatovan@ornl.gov), Guruprasad Kora (koragh@ornl.gov), Srikanth Yoginath (yoginathsb@ornl.gov), David Bauer (bauerda@ornl.gov)

References

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

See Also

solve for the normal solving function.

sla.qr the SCALAPACK version of the R function qr

Examples

library(RScaLAPACK)
rnorm(10000)->x
dim(x)=c(100,100)
c(1,0*(1:99))->y
dim(y)=c(100,1)
sla.solve(x,y, NPROWS=2, NPCOLS=2)

[Package RScaLAPACK version 0.6.1 Index]