spmv_vbr.h

/*
 *  file :  spmv_vbr.h
 *  desc :  Implementation of variable block row format.
 *          C-wrappers for SPARSKIT's Fortran routines.
 */

#if !defined(INC_SPMV_VBR_H)
#define INC_SPMV_VBR_H

#define ZERO_BASED_INDEXING

#if !defined(INC_SPMV_H)
typedef struct tagSparseMatrixCSR SparseMatrixCSR;
#endif

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

  int mb, nb;           /* # of block rows, cols */

  int *ia;          /* length nr+1 */
  int *ja;          /* length ia[nr]-1 */
  int *ka;          /* length ia[nr]-1 (?) */
  int *kvstr;           /* length nr+1 */
  int *kvstc;           /* length nc+1 */
  double *a;            /* length nnz */
} SparseMatrixVBR;

extern void spmv_alloc_vbr (SparseMatrixVBR * B, int m, int n, int nnz,
                int mb, int nb, int num_blocks);

extern void spmv_convert_csr2vbr (const SparseMatrixCSR * A,
                  SparseMatrixVBR * B);

extern void spmv_convert_csr2vbr_autoblock (int threshold,
                        const SparseMatrixCSR * A,
                        SparseMatrixVBR * B);

extern void spmv_destroy_vbr (SparseMatrixVBR * B);

extern void spmv_vbr (const SparseMatrixVBR * B, const double *x, double *y);

extern void spmv_convert_csr2vbr_conformal (const SparseMatrixCSR * A,
                        double threshold,
                        SparseMatrixVBR * B);

/*
 *  SPARSKIT Fortran routines
 */
#if defined(F77_EXTERN_UNDERSCORES)
#define csrvbr csrvbr_
#define vbrmv vbrmv_
#endif

extern void csrvbr (const int *n,
            const int *ptr, const int *ind, const double *val,
            int *nr, int *nc, int *kvstr, int *kvstc,
            int *ib, int *jb, int *kb, double *b,
            int *job, int *iwk, const int *nkmax, const int *nzmax,
            int *ierr);

extern void vbrmv (const int *nr, const int *nc,
           const int *ia, const int *ja, const int *ka,
           const double *a, const int *kvstr, const int *kvstc,
           const double *x, double *b);


#endif

/* eof */

---