mkl_?bsrsv
mkl_?bsrsv
Solves a system of linear equations for a sparse matrix in the BSR format (deprecated).
Syntax
void
mkl_sbsrsv
(
const
char
*transa
,
const
MKL_INT
*m
,
const
MKL_INT
*lb
,
const
float
*alpha
,
const
char
*matdescra
,
const
float
*val
,
const
MKL_INT
*indx
,
const
MKL_INT
*pntrb
,
const
MKL_INT
*pntre
,
const
float
*x
,
float
*y
);
void
mkl_dbsrsv
(
const
char
*transa
,
const
MKL_INT
*m
,
const
MKL_INT
*lb
,
const
double
*alpha
,
const
char
*matdescra
,
const
double
*val
,
const
MKL_INT
*indx
,
const
MKL_INT
*pntrb
,
const
MKL_INT
*pntre
,
const
double
*x
,
double
*y
);
void
mkl_cbsrsv
(
const
char
*transa
,
const
MKL_INT
*m
,
const
MKL_INT
*lb
,
const
MKL_Complex8
*alpha
,
const
char
*matdescra
,
const
MKL_Complex8
*val
,
const
MKL_INT
*indx
,
const
MKL_INT
*pntrb
,
const
MKL_INT
*pntre
,
const
MKL_Complex8
*x
,
MKL_Complex8
*y
);
void
mkl_zbsrsv
(
const
char
*transa
,
const
MKL_INT
*m
,
const
MKL_INT
*lb
,
const
MKL_Complex16
*alpha
,
const
char
*matdescra
,
const
MKL_Complex16
*val
,
const
MKL_INT
*indx
,
const
MKL_INT
*pntrb
,
const
MKL_INT
*pntre
,
const
MKL_Complex16
*x
,
MKL_Complex16
*y
);
Include Files
- mkl.h
Description
This routine is deprecated. Use mkl_sparse_?_trsvfrom the Inspector-executor Sparse BLAS interface instead.
Intel® oneAPI Math Kernel Library
The
mkl_?bsrsv
routine solves a system of linear equations with matrix-vector operations for a sparse matrix in the BSR format: y := alpha*inv(A)*x
or
y := alpha*inv(AT)* x,
where:
alpha
is scalar, x
and y
are vectors, A
is a sparse upper or lower triangular matrix with unit or non-unit main diagonal, A
T
is the transpose of A
.This routine supports a BSR format both with one-based indexing and zero-based indexing.
Input Parameters
- transa
- Specifies the operation.Ifortransa='N''n', theny:=alpha*inv(A)*xIfortransa='T''t'or'C'or'c', then,y:=alpha*inv(AT)*x
- m
- Number of block columns of the matrixA.
- lb
- Size of the block in the matrixA.
- alpha
- Specifies the scalaralpha.
- matdescra
- Array of six elements, specifies properties of the matrix used for operation. Only first four array elements are used, their possible values are given inTable “Possible Values of the Parameter. Possible combinations of element values of this parameter are given inmatdescra(descra)”Table “Possible Combinations of Element Values of the Parameter.matdescra”
- val
- Array containing elements of non-zero blocks of the matrixA. Its length is equal to the number of non-zero blocks in the matrixAmultiplied by.lb*lbRefer to thevaluesarray description in BSR Format for more details.The non-zero elements of the given row of the matrix must be stored in the same order as they appear in the row (from left to right).No diagonal element can be omitted from a sparse storage if the solver is called with the non-unit indicator.
- indx
- For one-based indexing, array containing the column indices plus one for each non-zero element of the matrixA. For zero-based indexing, array containing the column indices for each non-zero element of the matrixA.Its length is equal to the number of non-zero blocks in the matrixA.Refer to thecolumnsarray description in BSR Format for more details.
- pntrb
- Array of lengthm.This array contains row indices, such thatis the first index of block rowpntrb[i] -pntrb[0]iin the arrayindxRefer topointerBarray description in BSR Format for more details.
- pntre
- Array of lengthm.For one-based indexing this array contains row indices, such thatis the last index of block rowpntre[i] -pntrb[1]iin the arrayindx.For zero-based indexing this array contains row indices, such thatis the last index of block rowpntre[i] -pntrb[0] - 1iin the arrayindx.Refer topointerEarray description in BSR Format for more details.
- x
- Array, size at least(.m*lb)On entry, the arrayxmust contain the vectorx. The elements are accessed with unit increment.
- y
- Array, size at least(.m*lb)On entry, the arrayymust contain the vectory. The elements are accessed with unit increment.
Output Parameters
- y
- Contains solution vectorx.