mkl_?coogemv

Computes matrix-vector product of a sparse general matrix stored in the coordinate format with one-based indexing.

Syntax

Fortran:

call mkl_scoogemv(transa, m, val, rowind, colind, nnz, x, y)

call mkl_dcoogemv(transa, m, val, rowind, colind, nnz, x, y)

call mkl_ccoogemv(transa, m, val, rowind, colind, nnz, x, y)

call mkl_zcoogemv(transa, m, val, rowind, colind, nnz, x, y)

C:

mkl_scoogemv(&transa, &m, val, rowind, colind, &nnz, x, y);

mkl_dcoogemv(&transa, &m, val, rowind, colind, &nnz, x, y);

mkl_ccoogemv(&transa, &m, val, rowind, colind, &nnz, x, y);

mkl_zcoogemv(&transa, &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_?coogemv routine performs a matrix-vector operation defined as y := A*x

or y := A'*x,

where:

x and y are vectors,

A is an m-by-m sparse square matrix in the coordinate format, A' is the transpose of A.

Note iconNote

This routine supports only one-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.

transa

CHARACTER*1. Specifies the operation.

If transa = 'N' or 'n', then the matrix-vector product is computed as y := A*x

If transa = 'T' or 't' or 'C' or 'c', then the matrix-vector product is computed as y := A'*x,

m

INTEGER. Number of rows of the matrix A.

val

REAL for mkl_scoogemv.

DOUBLE PRECISION for mkl_dcoogemv.

COMPLEX for mkl_ccoogemv.

DOUBLE COMPLEX for mkl_zcoogemv.

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_scoogemv.

DOUBLE PRECISION for mkl_dcoogemv.

COMPLEX for mkl_ccoogemv.

DOUBLE COMPLEX for mkl_zcoogemv.

Array, DIMENSION is m.

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

Output Parameters

y

REAL for mkl_scoogemv.

DOUBLE PRECISION for mkl_dcoogemv.

COMPLEX for mkl_ccoogemv.

DOUBLE COMPLEX for mkl_zcoogemv.

Array, DIMENSION at least m.

On exit, the array y must contain the vector y.

Interfaces

FORTRAN 77:

SUBROUTINE mkl_scoogemv(transa, m, val, rowind, colind, nnz, x, y)
  CHARACTER*1   transa
  INTEGER       m, nnz
  INTEGER       rowind(*), colind(*)
  REAL            val(*), x(*), y(*)

SUBROUTINE mkl_dcoogemv(transa, m, val, rowind, colind, nnz, x, y)
  CHARACTER*1   transa
  INTEGER       m, nnz
  INTEGER       rowind(*), colind(*)
  DOUBLE PRECISION      val(*), x(*), y(*)

SUBROUTINE mkl_ccoogemv(transa, m, val, rowind, colind, nnz, x, y)
  CHARACTER*1   transa
  INTEGER       m, nnz
  INTEGER       rowind(*), colind(*)
  COMPLEX     val(*), x(*), y(*)

SUBROUTINE mkl_zcoogemv(transa, m, val, rowind, colind, nnz, x, y)
  CHARACTER*1   transa
  INTEGER       m, nnz
  INTEGER       rowind(*), colind(*)
  DOUBLE COMPLEX        val(*), x(*), y(*)

C:

void mkl_scoogemv(char *transa, int *m, float  *val, int *rowind,
int *colind, int *nnz, float *x, float *y);

void mkl_dcoogemv(char *transa, int *m, double  *val, int *rowind,
int *colind, int *nnz, double *x, double *y);

void mkl_ccoogemv(char *transa, int *m, MKL_Complex8  *val, int *rowind,
int *colind, int *nnz, MKL_Complex8 *x, MKL_Complex8 *y);

void mkl_zcoogemv(char *transa, int *m, MKL_Complex16  *val, int *rowind, 
int *colind, int *nnz, MKL_Complex16 *x, MKL_Complex16 *y);

Submit feedback on this help topic

Copyright © 1994 - 2011, Intel Corporation. All rights reserved.