mkl_sparse_?

Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

Download PDF

ID 766686

Date 12/16/2022

Version

Public

A newer version of this document is available. Customers should click here to go to the newest version.

Visible to Intel only — GUID: GUID-9E3880E2-E7FE-4636-B312-6A67551584CB

View Details

mkl_sparse_?_syprd

Computes the symmetric triple product of a sparse matrix and a dense matrix and stores the result as a dense matrix.

Syntax

stat = mkl_sparse_s_syprd (operation, A, B, denselayoutB, ldb, alpha, beta, C, denselayoutC, ldc)

stat = mkl_sparse_d_syprd (operation, A, B, denselayoutB, ldb, alpha, beta, C, denselayoutC, ldc)

stat = mkl_sparse_c_syprd (operation, A, B, denselayoutB, ldb, alpha, beta, C, denselayoutC, ldc)

stat = mkl_sparse_z_syprd (operation, A, B, denselayoutB, ldb, alpha, beta, C, denselayoutC, ldc)

Include Files

mkl_spblas.f90

Description

The mkl_sparse_?_syprd routine performs a multiplication of three sparse matrices that results in a symmetric or Hermitian matrix, C.

C:=alpha*A*B*op(A) + beta*C

C:=alpha*op(A)*B*A + beta*C

depending on the matrix modifier operation. Here A is a sparse matrix, B and C are dense and symmetric (or Hermitian) matrices.

op is the transpose (real precision) or conjugate transpose (complex precision) operator.

NOTE:

This routine is not supported for sparse matrices in COO or CSC formats. It supports only CSR and BSR formats. In addition, this routine supports only the sorted CSR and sorted BSR formats for the input matrix. If the data is unsorted, call the mkl_sparse_order routine before either mkl_sparse_sypr or mkl_sparse_?_syprd.

Input Parameters

operation

C_INT.

Specifies operation on the input sparse matrix.

SPARSE_OPERATION_NON_TRANSPOSE

Non-transpose case.

C:=alpha*A*B*(A^T)+beta*C for real precision.

C:=alpha*A*B*(A^H)+beta*C for complex precision.

SPARSE_OPERATION_TRANSPOSE

Transpose case. This is not supported for complex matrices.

C:=alpha*(A^T)*B*A +beta*C

SPARSE_OPERATION_CONJUGATE_TRANSPOSE

Conjugate transpose case. This is not supported for real matrices.

C:=alpha*(A^H)*B*A+beta*C

SPARSE_MATRIX_T.

Handle which contains the sparse matrix A.

SPARSE_MATRIX_T.

Input dense matrix. Only the upper triangular part of the matrix is used for computation.

denselayoutB

C_INT.

Structure that describes the storage scheme for the dense matrix.

SPARSE_LAYOUT_COLUMN_MAJOR	Store elements in a column-major layout.
SPARSE_LAYOUT_ROW_MAJOR	Store elements in a row-major layout.

ldb

SPARSE_MATRIX_T.

Leading dimension of matrix B.

alpha

Scalar parameter.

mkl_sparse_s_syprd	C_FLOAT
mkl_sparse_d_syprd	C_DOUBLE
mkl_sparse_c_syprd	C_FLOAT_COMPLEX
mkl_sparse_z_syprd	C_DOUBLE_COMPLEX

beta

Scalar parameter.

mkl_sparse_s_syprd	C_FLOAT
mkl_sparse_d_syprd	C_DOUBLE
mkl_sparse_c_syprd	C_FLOAT_COMPLEX
mkl_sparse_z_syprd	C_DOUBLE_COMPLEX

NOTE:

Since the upper triangular part of matrix C is the only portion that is processed, set real values of alpha and beta in the complex case to obtain the Hermitian matrix.

denselayoutC

C_INT.

Structure that describes the storage scheme for the dense matrix.

SPARSE_LAYOUT_COLUMN_MAJOR	Store elements in a column-major layout.
SPARSE_LAYOUT_ROW_MAJOR	Store elements in a row-major layout.

ldc

C_INT.

Leading dimension of matrix C.

Output Parameters

SPARSE_MATRIX_T.

Handle which contains the resulting dense matrix. Only the upper-triangular part of the matrix is computed.

stat

INTEGER

Value indicating whether the operation was successful, or the reason why it failed.

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.

Parent topic: Inspector-Executor Sparse BLAS Execution 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_sparse_?_syprd