Developer Reference for Intel® oneAPI Math Kernel Library for Fortran
mkl_sparse_?_sp2md
Computes the product of two sparse matrices (support operations on both matrices) and stores the result as a dense matrix.
Syntax
stat = mkl_sparse_s_sp2md (opA, descrA, A, opB, descrB, B, alpha, beta, C, layout, ldc )
stat = mkl_sparse_d_sp2md (opA, descrA, A, opB, descrB, B, alpha, beta, C, layout, ldc )
stat = mkl_sparse_c_sp2md (opA, descrA, A, opB, descrB, B, alpha, beta, C, layout, ldc )
stat = mkl_sparse_z_sp2md (opA, descrA, A, opB, descrB, B, alpha, beta, C, layout, ldc )Include Files
mkl_spblas.f90
Description
The mkl_sparse_?_sp2md routine performs a matrix-matrix operation:
\[C \leftarrow \alpha \cdot opA(A) \cdot opB(B) + \beta\cdot C\]
where \(A\) and \(B\) are sparse matrices, opA() is a matrix modifier for matrix \(A\) , opB() is a matrix modifier for matrix \(B\), and \(C\) is a dense matrix, \(\alpha\) and \(\beta\) are scalars.
(SPARSE_INDEX_BASE_ZERO, SPARSE_LAYOUT_ROW_MAJOR)
(SPARSE_INDEX_BASE_ONE, SPARSE_LAYOUT_COLUMN_MAJOR)
Input Parameters
opA
C_INT.
Specifies operation on input matrix \(A\).
opA |
Description |
|---|---|
SPARSE_OPERATION_NON_TRANSPOSE |
\(opA(A) = A\) |
SPARSE_OPERATION_TRANSPOSE |
\(opA(A) = A^{T}\) |
SPARSE_OPERATION_CONJUGATE_TRANSPOSE |
\(opA(A) = A^{H}\) |
descrA
MATRIX_DESCR .
Descriptor specifying sparse matrix properties.
sparse_matrix_type_t type |
Specifies the type of a sparse matrix:
|
sparse_fill_mode_t mode |
Specifies the triangular matrix part for symmetric, Hermitian, triangular, and block-triangular matrices:
|
sparse_diag_type_t diag |
Specifies diagonal type for non-general matrices:
|
A
SPARSE_MATRIX_T.
Handle which contains the sparse matrix \(A\) .
opB
C_INT.
Specifies operation on input matrix \(B\).
opB |
Description |
|---|---|
SPARSE_OPERATION_NON_TRANSPOSE |
\(opB(B) = B\) |
SPARSE_OPERATION_TRANSPOSE |
\(opB(B) = B^{T}\) |
SPARSE_OPERATION_CONJUGATE_TRANSPOSE |
\(opB(B) = B^{H}\) |
descrB
MATRIX_DESCR .
Descriptor specifying sparse matrix properties.
sparse_matrix_type_t type |
Specifies the type of a sparse matrix:
|
sparse_fill_mode_t mode |
Specifies the triangular matrix part for symmetric, Hermitian, triangular, and block-triangular matrices:
|
sparse_diag_type_t diag |
Specifies diagonal type for non-general matrices:
|
B
SPARSE_MATRIX_T.
Handle which contains the sparse matrix \(B\) .
alpha
C_FLOAT |
for mkl_sparse_s_sp2md |
C_DOUBLE |
for mkl_sparse_d_sp2md |
C_FLOAT_COMPLEX |
for mkl_sparse_c_sp2md |
C_DOUBLE_COMPLEX |
for mkl_sparse_z_sp2md |
Specifies the scalar, \(\alpha\) .
beta
C_FLOAT |
for mkl_sparse_s_sp2md |
C_DOUBLE |
for mkl_sparse_d_sp2md |
C_FLOAT_COMPLEX |
for mkl_sparse_c_sp2md |
C_DOUBLE_COMPLEX |
for mkl_sparse_z_sp2md |
Specifies the scalar, \(\beta\) .
layout
C_INT .
Describes the storage scheme for the dense matrix \(C\):
SPARSE_LAYOUT_COLUMN_MAJOR |
Storage of elements uses column major layout. |
SPARSE_LAYOUT_ROW_MAJOR |
Storage of elements uses row major layout. |
ldc
C_INT .
Leading dimension of matrix \(C\) .
Output Parameters
C
C_FLOAT |
for mkl_sparse_s_sp2md |
C_DOUBLE |
for mkl_sparse_d_sp2md |
C_FLOAT_COMPLEX |
for mkl_sparse_c_sp2md |
C_DOUBLE_COMPLEX |
for mkl_sparse_z_sp2md |
The resulting dense matrix.
stat
INTEGER
Value indicating whether the operation was successful, and if not, why:
SPARSE_STATUS_SUCCESS |
The operation was successful. |
SPARSE_STATUS_NOT_INITIALIZED |
The routine encountered an empty handle or matrix array. |
SPARSE_STATUS_ALLOC_FAILED |
Internal memory allocation failed. |
SPARSE_STATUS_INVALID_VALUE |
The input parameters contain an invalid value. |
SPARSE_STATUS_EXECUTION_FAILED |
Execution failed. |
SPARSE_STATUS_INTERNAL_ERROR |
An error in algorithm implementation occurred. |
SPARSE_STATUS_NOT_SUPPORTED |
The requested operation is not supported. |
Return Values
See Output Parameter stat.