Visible to Intel only — GUID: GUID-C09039C6-8ABD-403A-9980-3BEBF0FF3573
Visible to Intel only — GUID: GUID-C09039C6-8ABD-403A-9980-3BEBF0FF3573
mkl_sparse_sypr
Computes the symmetric product of three sparse matrices and stores the result in a newly allocated sparse matrix.
Syntax
stat = mkl_sparse_sypr (operation, A, B, descrB, C, request)
Include Files
- mkl_spblas.f90
Description
The mkl_sparse_sypr routine performs a multiplication of three sparse matrices that results in a symmetric or Hermitian matrix, C.
C:=A*B*opA(A)or
C:=opA(A)*B*Adepending on the matrix modifier operation.
Here, A, B, and C are sparse matrices, where A has a general structure while B and C are symmetric (for real data types) or Hermitian (for complex data types) matrices. opA is the transpose (real data types) or conjugate transpose (complex data types) operator.
This routine is not supported for sparse matrices in COO or CSC formats. This routine supports only CSR and BSR formats. In addition, it 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 matrices.
|
||||||||||
A |
SPARSE_MATRIX_T. Handle which contains the sparse matrix A. |
||||||||||
B |
SPARSE_MATRIX_T. Handle which contains the sparse matrix B. |
||||||||||
descrB |
MATRIX_DESCR. Structure specifying properties of the sparse matrix. sparse_matrix_type_t type specifies the type of a sparse matrix
sparse_fill_mode_t mode specifies the triangular matrix part.
sparse_diag_type_t diag specifies the type of diagonal.
NOTE:
This routine also supports C=AAT,H with these parameters: descrB.type=SPARSE_MATRIX_TYPE_DIAGONAL descrB.diag=SPARSE_DIAG_UNIT In this case, you do not need to allocate structure B. Use the routine as a 2-stage version of mkl_sparse_syrk. |
||||||||||
request |
INTEGER. Use this routine to specify if the computations should be performed in a single step or using the two-stage algorithm. See Two-stage Algorithm for Inspector-executor Sparse BLAS Routines for more information.
|
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 the reason why it failed.
|