bsr.h

#pragma once

#ifdef __cplusplus
extern "C" {
#endif

typedef
struct bsr_matrix
{
    // number or rows
    int nrows;
    // number or columns
    int ncols;
    // number or non-zero blocks
    int nnz;
    // block-size
    int b;

    // pointer ot row offsets
    int *rowptr;
    // pointer to column indices
    int *colidx;
    // pointer to double-precision values
    double *dbl;
    // pointer to single-precision values
    float *flt;

} bsr_matrix;

bsr_matrix* bsr_alloc();

void bsr_free(bsr_matrix *A);

void bsr_init(bsr_matrix *A, int nrows, int nnz, int b);

void bsr_vmspmv3(bsr_matrix *A, const double *x, double *y);

void bsr_vdspmv3(bsr_matrix *A, const double *x, double *y);

void bsr_downcast(bsr_matrix *M);

void bsr_info(bsr_matrix *A);

void bsr_sparsity(const bsr_matrix *A, const char *name);

void bsr_nonzeros(bsr_matrix *A, const char *name);

#ifdef __cplusplus
}
#endif