spmv_gencb.h

#if !defined(INCLUDE_SPMV_GCB_H)
#define INCLUDE_SPMV_GCB_H

typedef struct tagSparseMatrixGeneralCSR SparseMatrixGeneralCSR;

typedef enum
{ gcsrTypeBlockRows, gcsrTypeSparseRows } gcsr_type;

struct tagSparseMatrixGeneralCSR
{
  gcsr_type type;       /* type of block this is */

  int row;          /* first row in this block */
  int col_min;          /* min column index for this block */
  int m;            /* # of rows in this block */
  int n;            /* maximum width of this block */
  int nnz;          /* total # of non-zeros in this block */

  int *row_ptr;         /* row pointers to columns of each row (m+1) */
  int *col_ind;         /* column indices within each row (nnz) */
  double *val;          /* data values (nnz) */

  /* blocks are linked in a list */
  SparseMatrixGeneralCSR *next;

  /* variable fields, depending on type */
  int *row_ind;
};

typedef struct tagSparseMatrixCacheBlockedGCSR
{
  int m, n;         /* total matrix dimensions */
  int nnz;          /* total # of non-zeros */

  int num_blocks;
  SparseMatrixGeneralCSR *mats; /* matrix sub-blocks (num_blocks) */
} SparseMatrixCacheBlockedGCSR;

extern void convert_csr2cbgcsr (int R, int C,
                int m, int n, const int *ptr, const int *ind,
                const double *val,
                SparseMatrixCacheBlockedGCSR * B);
    /*
     *  convert_csr2cbgcsr( R, C, m, n, ptr, ind, val, B );
     *
     *  Given an m x n CSR sparse matrix (ptr, ind, val),
     *  converts it to R x C cache blocked matrix B.
     */

extern void destroy_cbgcsr (SparseMatrixCacheBlockedGCSR * B);

extern void spmv_cbgcsr (const SparseMatrixCacheBlockedGCSR * B,
             const double *x, double *y);
    /* computes y += B*x */

#endif

/*
 *  $Log: spmv_gencb.h,v $
 *  Revision 1.1.1.1  2005/06/29 22:07:51  rvuduc
 *  Initial import
 *
 *
 *  eof
 */

---