|
kusano |
7d535a |
|
|
kusano |
7d535a |
<meta content="text/html;charset=UTF-8" http-equiv="Content-Type">
|
|
kusano |
7d535a |
<title>SuperLU: SRC/csp_blas2.c File Reference</title>
|
|
kusano |
7d535a |
<link href="doxygen.css" rel="stylesheet" type="text/css">
|
|
kusano |
7d535a |
<link href="tabs.css" rel="stylesheet" type="text/css">
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Main Page
|
|
kusano |
7d535a |
Data Structures
|
|
kusano |
7d535a |
Files
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
SRC/csp_blas2.c File ReferenceSparse BLAS 2, using some dense BLAS 2 operations. More...
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
#include "slu_cdefs.h"
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Functions
|
|
kusano |
7d535a |
void cusolve (int, int, complex *, complex *)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Solves a dense upper triangular system.
|
|
kusano |
7d535a |
void clsolve (int, int, complex *, complex *)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Solves a dense UNIT lower triangular system.
|
|
kusano |
7d535a |
void cmatvec (int, int, int, complex *, complex *, complex *)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Performs a dense matrix-vector multiply: Mxvec = Mxvec + M * vec.
|
|
kusano |
7d535a |
int sp_ctrsv (char *uplo, char *trans, char *diag, SuperMatrix *L, SuperMatrix *U, complex *x, SuperLUStat_t *stat, int *info)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Solves one of the systems of equations A*x = b, or A'*x = b.
|
|
kusano |
7d535a |
int sp_cgemv (char *trans, complex alpha, SuperMatrix *A, complex *x, int incx, complex beta, complex *y, int incy)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Performs one of the matrix-vector operations y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y.
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Detailed Description
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
-- SuperLU routine (version 3.0) --
|
|
kusano |
7d535a |
Univ. of California Berkeley, Xerox Palo Alto Research Center,
|
|
kusano |
7d535a |
and Lawrence Berkeley National Lab.
|
|
kusano |
7d535a |
October 15, 2003
|
|
kusano |
7d535a |
Function Documentation
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
void clsolve
|
|
kusano |
7d535a |
(
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
ldm,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
ncol,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
M,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
rhs
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
The unit lower triangular matrix is stored in a 2D array M(1:nrow,1:ncol). The solution will be returned in the rhs vector.
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
void cmatvec
|
|
kusano |
7d535a |
(
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
ldm,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
nrow,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
ncol,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
M,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
vec,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
Mxvec
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
The input matrix is M(1:nrow,1:ncol); The product is returned in Mxvec[].
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
void cusolve
|
|
kusano |
7d535a |
(
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
ldm,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
ncol,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
M,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
rhs
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
The upper triangular matrix is stored in a 2-dim array M(1:ldm,1:ncol). The solution will be returned in the rhs vector.
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int sp_cgemv
|
|
kusano |
7d535a |
(
|
|
kusano |
7d535a |
char *
|
|
kusano |
7d535a |
trans,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex
|
|
kusano |
7d535a |
alpha,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
SuperMatrix *
|
|
kusano |
7d535a |
A,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
x,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
incx,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex
|
|
kusano |
7d535a |
beta,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
y,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int
|
|
kusano |
7d535a |
incy
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Purpose
|
|
kusano |
7d535a |
=======
|
|
kusano |
7d535a |
sp_cgemv() performs one of the matrix-vector operations
|
|
kusano |
7d535a |
y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y,
|
|
kusano |
7d535a |
where alpha and beta are scalars, x and y are vectors and A is a
|
|
kusano |
7d535a |
sparse A->nrow by A->ncol matrix.
|
|
kusano |
7d535a |
Parameters
|
|
kusano |
7d535a |
==========
|
|
kusano |
7d535a |
TRANS - (input) char*
|
|
kusano |
7d535a |
On entry, TRANS specifies the operation to be performed as
|
|
kusano |
7d535a |
follows:
|
|
kusano |
7d535a |
TRANS = 'N' or 'n' y := alpha*A*x + beta*y.
|
|
kusano |
7d535a |
TRANS = 'T' or 't' y := alpha*A'*x + beta*y.
|
|
kusano |
7d535a |
TRANS = 'C' or 'c' y := alpha*A'*x + beta*y.
|
|
kusano |
7d535a |
ALPHA - (input) complex
|
|
kusano |
7d535a |
On entry, ALPHA specifies the scalar alpha.
|
|
kusano |
7d535a |
A - (input) SuperMatrix*
|
|
kusano |
7d535a |
Before entry, the leading m by n part of the array A must
|
|
kusano |
7d535a |
contain the matrix of coefficients.
|
|
kusano |
7d535a |
X - (input) complex*, array of DIMENSION at least
|
|
kusano |
7d535a |
( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n'
|
|
kusano |
7d535a |
and at least
|
|
kusano |
7d535a |
( 1 + ( m - 1 )*abs( INCX ) ) otherwise.
|
|
kusano |
7d535a |
Before entry, the incremented array X must contain the
|
|
kusano |
7d535a |
vector x.
|
|
kusano |
7d535a |
INCX - (input) int
|
|
kusano |
7d535a |
On entry, INCX specifies the increment for the elements of
|
|
kusano |
7d535a |
X. INCX must not be zero.
|
|
kusano |
7d535a |
BETA - (input) complex
|
|
kusano |
7d535a |
On entry, BETA specifies the scalar beta. When BETA is
|
|
kusano |
7d535a |
supplied as zero then Y need not be set on input.
|
|
kusano |
7d535a |
Y - (output) complex*, array of DIMENSION at least
|
|
kusano |
7d535a |
( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n'
|
|
kusano |
7d535a |
and at least
|
|
kusano |
7d535a |
( 1 + ( n - 1 )*abs( INCY ) ) otherwise.
|
|
kusano |
7d535a |
Before entry with BETA non-zero, the incremented array Y
|
|
kusano |
7d535a |
must contain the vector y. On exit, Y is overwritten by the
|
|
kusano |
7d535a |
updated vector y.
|
|
kusano |
7d535a |
INCY - (input) int
|
|
kusano |
7d535a |
On entry, INCY specifies the increment for the elements of
|
|
kusano |
7d535a |
Y. INCY must not be zero.
|
|
kusano |
7d535a |
==== Sparse Level 2 Blas routine.
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int sp_ctrsv
|
|
kusano |
7d535a |
(
|
|
kusano |
7d535a |
char *
|
|
kusano |
7d535a |
uplo,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
char *
|
|
kusano |
7d535a |
trans,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
char *
|
|
kusano |
7d535a |
diag,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
SuperMatrix *
|
|
kusano |
7d535a |
L,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
SuperMatrix *
|
|
kusano |
7d535a |
U,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
complex *
|
|
kusano |
7d535a |
x,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
SuperLUStat_t *
|
|
kusano |
7d535a |
stat,
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
int *
|
|
kusano |
7d535a |
info
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
)
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
Purpose
|
|
kusano |
7d535a |
=======
|
|
kusano |
7d535a |
sp_ctrsv() solves one of the systems of equations
|
|
kusano |
7d535a |
A*x = b, or A'*x = b,
|
|
kusano |
7d535a |
where b and x are n element vectors and A is a sparse unit , or
|
|
kusano |
7d535a |
non-unit, upper or lower triangular matrix.
|
|
kusano |
7d535a |
No test for singularity or near-singularity is included in this
|
|
kusano |
7d535a |
routine. Such tests must be performed before calling this routine.
|
|
kusano |
7d535a |
Parameters
|
|
kusano |
7d535a |
==========
|
|
kusano |
7d535a |
uplo - (input) char*
|
|
kusano |
7d535a |
On entry, uplo specifies whether the matrix is an upper or
|
|
kusano |
7d535a |
lower triangular matrix as follows:
|
|
kusano |
7d535a |
uplo = 'U' or 'u' A is an upper triangular matrix.
|
|
kusano |
7d535a |
uplo = 'L' or 'l' A is a lower triangular matrix.
|
|
kusano |
7d535a |
trans - (input) char*
|
|
kusano |
7d535a |
On entry, trans specifies the equations to be solved as
|
|
kusano |
7d535a |
follows:
|
|
kusano |
7d535a |
trans = 'N' or 'n' A*x = b.
|
|
kusano |
7d535a |
trans = 'T' or 't' A'*x = b.
|
|
kusano |
7d535a |
trans = 'C' or 'c' A^H*x = b.
|
|
kusano |
7d535a |
diag - (input) char*
|
|
kusano |
7d535a |
On entry, diag specifies whether or not A is unit
|
|
kusano |
7d535a |
triangular as follows:
|
|
kusano |
7d535a |
diag = 'U' or 'u' A is assumed to be unit triangular.
|
|
kusano |
7d535a |
diag = 'N' or 'n' A is not assumed to be unit
|
|
kusano |
7d535a |
triangular.
|
|
kusano |
7d535a |
L - (input) SuperMatrix*
|
|
kusano |
7d535a |
The factor L from the factorization Pr*A*Pc=L*U. Use
|
|
kusano |
7d535a |
compressed row subscripts storage for supernodes,
|
|
kusano |
7d535a |
i.e., L has types: Stype = SC, Dtype = SLU_C, Mtype = TRLU.
|
|
kusano |
7d535a |
U - (input) SuperMatrix*
|
|
kusano |
7d535a |
The factor U from the factorization Pr*A*Pc=L*U.
|
|
kusano |
7d535a |
U has types: Stype = NC, Dtype = SLU_C, Mtype = TRU.
|
|
kusano |
7d535a |
x - (input/output) complex*
|
|
kusano |
7d535a |
Before entry, the incremented array X must contain the n
|
|
kusano |
7d535a |
element right-hand side vector b. On exit, X is overwritten
|
|
kusano |
7d535a |
with the solution vector x.
|
|
kusano |
7d535a |
info - (output) int*
|
|
kusano |
7d535a |
If *info = -i, the i-th argument had an illegal value.
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
<address style="text-align: right;"><small>Generated on Mon Nov 22 10:23:47 2010 for SuperLU by </small></address>
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
 1.5.5
|
|
kusano |
7d535a |
|
|
kusano |
7d535a |
|