Interface for matrix vector, andmatrix matrix products

Propose an interface for MatrixVector and MAtrixMatrix products.

src/z_spm.h

@@ -41,6 +41,69 @@ spm_complex64_t *z_spm2dense( const spmatrix_t *spm );
 /**
  * Matrix-Vector and matrix-matrix product routines
  */
+int spm_zgemv_csc( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+int spm_zgemv_csr( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+int spm_zgemv_ijv( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+
+int spm_zhemv_csc( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+int spm_zhemv_csr( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+int spm_zhemv_ijv( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+
+int spm_zsymv_csc( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+int spm_zsymv_csr( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+int spm_zsymv_ijv( const spm_trans_t trans, const spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, int incx, const spm_complex64_t beta, spm_complex64_t *y, int incy );
+
+int spm_zgemm_csc( const spm_side_t side, const spm_trans_t transA, const spm_trans_t transB, const spm_int_t K,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zgemm_csr( const spm_side_t side, const spm_trans_t transA, const spm_trans_t transB, const spm_int_t K,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zgemm_ijv( const spm_side_t side, const spm_trans_t transA, const spm_trans_t transB, const spm_int_t K,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zgemm( const spm_side_t side, const spm_trans_t transA, const spm_trans_t transB, const spm_int_t K,
+	       const spm_complex64_t alpha, const spmatrix_t      *A,
+	                                    const spm_complex64_t *B, int ldb,
+	       const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+
+int spm_zhemm_csc( const spm_side_t side, const spm_int_t N,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zhemm_csr( const spm_side_t side, const spm_int_t N,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zhemm_ijv( const spm_side_t side, const spm_int_t N,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zhemm( const spm_side_t side, const spm_int_t N,
+	       const spm_complex64_t alpha, const spmatrix_t      *A,
+	                                    const spm_complex64_t *B, int ldb,
+	       const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+
+int spm_zsymm_csc( const spm_side_t side, const spm_int_t N,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zsymm_csr( const spm_side_t side, const spm_int_t N,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zsymm_ijv( const spm_side_t side, const spm_int_t N,
+		   const spm_complex64_t alpha, const spmatrix_t      *A,
+		                                const spm_complex64_t *B, int ldb,
+		   const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+int spm_zsymm( const spm_side_t side, const spm_int_t N,
+	       const spm_complex64_t alpha, const spmatrix_t      *A,
+	                                    const spm_complex64_t *B, int ldb,
+	       const spm_complex64_t beta,        spm_complex64_t *C, int ldc );
+
 int spm_zspmv( spm_trans_t trans, spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *x, spm_int_t incx, spm_complex64_t beta, spm_complex64_t *y, spm_int_t incy );
 int spm_zspmm( spm_side_t side, spm_trans_t transA, spm_trans_t transB, spm_int_t K, spm_complex64_t alpha, const spmatrix_t *A, const spm_complex64_t *B, spm_int_t ldb, spm_complex64_t beta, spm_complex64_t *C, spm_int_t ldc );