sla.solve {RScaLAPACK} | R Documentation |
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.
sla.solve(A, B=NULL, NPROWS=0, NPCOLS=0, MB=16,RFLAG=1, SPAWN=1 )
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. |
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.
Nagiza Samatova (samatovan@ornl.gov), Guruprasad Kora (koragh@ornl.gov), Srikanth Yoginath (yoginathsb@ornl.gov), David Bauer (bauerda@ornl.gov)
http://www.netlib.org/scalapack/
solve
for the normal solving function.
sla.qr
the SCALAPACK version of the R function qr
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)