mkl_?csrmultd

Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

Download PDF

ID 766686

Date 10/31/2024

Version

Public

Visible to Intel only — GUID: GUID-D6F15254-4106-4112-9B11-F19DA5C8DDB1

View Details

mkl_?csrmultd

Computes product of two sparse matrices stored in the CSR format (3-array variation) with one-based indexing. The result is stored in the dense matrix (deprecated).

Syntax

call mkl_scsrmultd(trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

call mkl_dcsrmultd(trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

call mkl_ccsrmultd(trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

call mkl_zcsrmultd(trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

Include Files

mkl.fi

Description

This routine is deprecated. Use mkl_sparse_?_spmmdfrom the Intel® oneAPI Math Kernel Library (oneMKL) Inspector-executor Sparse BLAS interface instead.

The mkl_?csrmultd routine performs a matrix-matrix operation defined as

C := op(A)*B

where:

A, B are the sparse matrices in the CSR format (3-array variation), C is dense matrix;

op(A) is one of op(A) = A, or op(A) =A^T, or op(A) = A^H .

The routine works correctly if and only if the column indices in sparse matrix representations of matrices A and B are arranged in the increasing order for each row.

NOTE:

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.

trans

CHARACTER*1. Specifies the operation.

If trans = 'N' or 'n', then C := A*B

If trans = 'T' or 't' or 'C' or 'c', then C := A^T*B.

m

INTEGER. Number of rows of the matrix A.

n

INTEGER. Number of columns of the matrix A.

k

INTEGER. Number of columns of the matrix B.

a

REAL for mkl_scsrmultd.

DOUBLE PRECISION for mkl_dcsrmultd.

COMPLEX for mkl_ccsrmultd.

DOUBLE COMPLEX for mkl_zcsrmultd.

Array containing non-zero elements of the matrix A. Its length is equal to the number of non-zero elements in the matrix A. Refer to values array description in Sparse Matrix Storage Formats for more details.

ja

INTEGER. Array containing the column indices for each non-zero element of the matrix A. For each row the column indices must be arranged in the increasing order.

The length of this array is equal to the length of the array a. Refer to columns array description in Sparse Matrix Storage Formats for more details.

ia

INTEGER. Array of length m + 1 when trans = 'N' or 'n', or n + 1 otherwise.

This array contains indices of elements in the array a, such that ia(i) is the index in the array a of the first non-zero element from the row i. The value of the last element ia(m + 1) or ia(n + 1) is equal to the number of non-zero elements of the matrix A plus one. Refer to rowIndex array description in Sparse Matrix Storage Formats for more details.

b

REAL for mkl_scsrmultd.

DOUBLE PRECISION for mkl_dcsrmultd.

COMPLEX for mkl_ccsrmultd.

DOUBLE COMPLEX for mkl_zcsrmultd.

Array containing non-zero elements of the matrix B. Its length is equal to the number of non-zero elements in the matrix B. Refer to values array description in Sparse Matrix Storage Formats for more details.

jb

INTEGER. Array containing the column indices for each non-zero element of the matrix B. For each row the column indices must be arranged in the increasing order.

The length of this array is equal to the length of the array b. Refer to columns array description in Sparse Matrix Storage Formats for more details.

ib

INTEGER. Array of length m + 1.

This array contains indices of elements in the array b, such that ib(i) is the index in the array b of the first non-zero element from the row i. The value of the last element ib(m + 1) is equal to the number of non-zero elements of the matrix B plus one. Refer to rowIndex array description in Sparse Matrix Storage Formats for more details.

Output Parameters

c

REAL for mkl_scsrmultd.

DOUBLE PRECISION for mkl_dcsrmultd.

COMPLEX for mkl_ccsrmultd.

DOUBLE COMPLEX for mkl_zcsrmultd.

Array containing non-zero elements of the result matrix C.

ldc

INTEGER. Specifies the leading dimension of the dense matrix C as declared in the calling (sub)program. Must be at least max(m, 1) when trans = 'N' or 'n', or max(1, n) otherwise.

Interfaces

FORTRAN 77:

SUBROUTINE mkl_scsrmultd( trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

   CHARACTER*1  trans

   INTEGER       m, n, k, ldc

   INTEGER       ja(*), jb(*), ia(*), ib(*)

   REAL            a(*), b(*), c(ldc, *)

SUBROUTINE mkl_dcsrmultd( trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

   CHARACTER*1  trans

   INTEGER       m, n, k, ldc

   INTEGER       ja(*), jb(*), ia(*), ib(*)

   DOUBLE PRECISION         a(*), b(*), c(ldc, *)

SUBROUTINE mkl_ccsrmultd( trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

   CHARACTER*1  trans

   INTEGER       m, n, k, ldc

   INTEGER       ja(*), jb(*), ia(*), ib(*)

   COMPLEX     a(*), b(*), c(ldc, *)

SUBROUTINE mkl_zcsrmultd( trans, m, n, k, a, ja, ia, b, jb, ib, c, ldc)

   CHARACTER*1  trans

   INTEGER       m, n, k, ldc

   INTEGER       ja(*), jb(*), ia(*), ib(*)

   DOUBLE COMPLEX      a(*), b(*), c(ldc, *)

Parent topic: Sparse BLAS Level 2 and Level 3 Routines.

Select Your Language

Using Intel.com Search

Quick Links

Recent Searches

Advanced Search

Only search in

Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

mkl_?csrmultd

Syntax

Include Files

Description

Input Parameters

Output Parameters

Interfaces

FORTRAN 77: