mkl_cspblas_?cootrsv

Syntax

Fortran:

call mkl_cspblas_scootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

call mkl_cspblas_dcootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

call mkl_cspblas_ccootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

call mkl_cspblas_zcootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

C:

mkl_cspblas_scootrsv(&uplo, &transa, &diag, &m, val, rowind, colind, &nnz, x, y);

mkl_cspblas_dcootrsv(&uplo, &transa, &diag, &m, val, rowind, colind, &nnz, x, y);

mkl_cspblas_ccootrsv(&uplo, &transa, &diag, &m, val, rowind, colind, &nnz, x, y);

mkl_cspblas_zcootrsv(&uplo, &transa, &diag, &m, val, rowind, colind, &nnz, x, y);

Include Files

The FORTRAN 77 interfaces are specified in the mkl_spblas.fi include file and the C interfaces are specified in the mkl_spblas.h include file.

Description

The mkl_cspblas_?cootrsv routine solves a system of linear equations with matrix-vector operations for a sparse matrix stored in the coordinate format with zero-based indexing: A*y = x

or A'*y = x,

where:

x and y are vectors,

A is a sparse upper or lower triangular matrix with unit or non-unit main diagonal, A' is the transpose of A.

Note

This routine supports only zero-based indexing of the input arrays.

Input Parameters

Parameter descriptions are common for all implemented interfaces with the exception of data types that refer here to the FORTRAN 77 standard types. Data types specific to the different interfaces are described in the section "Interfaces" below.

uplo

CHARACTER*1. Specifies whether the upper or low triangle of the matrix A is considered.

If uplo = 'U' or 'u', then the upper triangle of the matrix A is used.

If uplo = 'L' or 'l', then the low triangle of the matrix A is used.

transa

CHARACTER*1. Specifies the system of linear equations.

If transa = 'N' or 'n', then A*y = x

If transa = 'T' or 't' or 'C' or 'c', then A'*y = x,

diag

CHARACTER*1. Specifies whether A is unit triangular.

If diag = 'U' or 'u', then A is unit triangular.

If diag = 'N' or 'n', then A is not unit triangular.

m

INTEGER. Number of rows of the matrix A.

val

REAL for mkl_cspblas_scootrsv.

DOUBLE PRECISION for mkl_cspblas_dcootrsv.

COMPLEX for mkl_cspblas_ccootrsv.

DOUBLE COMPLEX for mkl_cspblas_zcootrsv.

Array of length nnz, contains non-zero elements of the matrix A in the arbitrary order.

Refer to values array description in Coordinate Format for more details.

rowind

INTEGER. Array of length nnz, contains the row indices for each non-zero element of the matrix A.

Refer to rows array description in Coordinate Format for more details.

colind

INTEGER. Array of length nnz, contains the column indices for each non-zero element of the matrix A. Refer to columns array description in Coordinate Format for more details.

nnz

INTEGER. Specifies the number of non-zero element of the matrix A.

Refer to nnz description in Coordinate Format for more details.

x

REAL for mkl_cspblas_scootrsv.

DOUBLE PRECISION for mkl_cspblas_dcootrsv.

COMPLEX for mkl_cspblas_ccootrsv.

DOUBLE COMPLEX for mkl_cspblas_zcootrsv.

Array, DIMENSION is m.

On entry, the array x must contain the vector x.

Output Parameters

y

REAL for mkl_cspblas_scootrsv.

DOUBLE PRECISION for mkl_cspblas_dcootrsv.

COMPLEX for mkl_cspblas_ccootrsv.

DOUBLE COMPLEX for mkl_cspblas_zcootrsv.

Array, DIMENSION at least m.

Contains the vector y.

Interfaces

SUBROUTINE mkl_cspblas_scootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

  CHARACTER*1   uplo, transa, diag

  INTEGER       m, nnz

  INTEGER       rowind(*), colind(*)

  REAL           val(*), x(*), y(*)

SUBROUTINE mkl_cspblas_dcootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

  CHARACTER*1   uplo, transa, diag

  INTEGER       m, nnz

  INTEGER       rowind(*), colind(*)

  DOUBLE PRECISION      val(*), x(*), y(*)

SUBROUTINE mkl_cspblas_ccootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

  CHARACTER*1   uplo, transa, diag

  INTEGER       m, nnz

  INTEGER       rowind(*), colind(*)

  COMPLEX        val(*), x(*), y(*)

SUBROUTINE mkl_cspblas_zcootrsv(uplo, transa, diag, m, val, rowind, colind, nnz, x, y)

  CHARACTER*1   uplo, transa, diag

  INTEGER       m, nnz

  INTEGER       rowind(*), colind(*)

  DOUBLE COMPLEX      val(*), x(*), y(*)

void mkl_cspblas_scootrsv(char *uplo, char *transa, char *diag, int *m,

float *val, int *rowind, int *colind, int *nnz, float *x, float *y);

void mkl_cspblas_dcootrsv(char *uplo, char *transa, char *diag, int *m,

double *val, int *rowind, int *colind, int *nnz, double *x, double *y);

void mkl_cspblas_ccootrsv(char *uplo, char *transa, char *diag, int *m,

MKL_Complex8 *val, int *rowind, int *colind, int *nnz, MKL_Complex8 *x, MKL_Complex8 *y);

void mkl_cspblas_zcootrsv(char *uplo, char *transa, char *diag, int *m,

MKL_Complex16 *val, int *rowind, int *colind, int *nnz, MKL_Complex16 *x, MKL_Complex16 *y);