Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

ID 766686
Date 7/13/2023
Public

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

Document Table of Contents

mkl_sparse_syrk

Computes the product of sparse matrix with its transpose (or conjugate transpose) and stores the result in a newly allocated sparse matrix.

Syntax

stat = mkl_sparse_syrk (operation, A, C)

Include Files

  • mkl_spblas.f90

Description

The mkl_sparse_syrk routine performs a sparse matrix-matrix operation which results in a sparse matrix C that is either Symmetric (real) or Hermitian (complex):

C := A*op(A)
where op(*) is the transpose for real matrices and conjugate transpose for complex matrices OR
C := op(A)*A
depending on the matrix modifier op which can be the transpose for real matrices or conjugate transpose for complex matrices.

Here, A and C are sparse matrices.

NOTE:
This routine is not supported for sparse matrices in COO or CSC formats. It supports only CSR and BSR formats. Additionally, this routine supports only the sorted CSR and sorted BSR formats for the input matrix. If data is unsorted, call the mkl_sparse_order routine before either mkl_sparse_syrk or mkl_sparse_?_syrkd.

Input Parameters

operation

Specifies the operation op() on input matrix .

SPARSE_OPERATION_NON_TRANSPOSE, Non-transpose,C := A*op(A) where op(*) is the transpose for real matrices and conjugate transpose for complex matrices

SPARSE_OPERATION_TRANSPOSE, Transpose,C := (AT)*Afor real matrix A

SPARSE_OPERATION_CONJUGATE_TRANSPOSE, Conjugate transpose, C := (AH)*A for complex matrix A.

A

SPARSE_MATRIX_T.

Handle which contains the sparse matrix A.

Output Parameters

C

SPARSE_MATRIX_T.

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

stat

INTEGER

Value indicating whether the operation was successful or not, and 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.