mapmat.h

Go to the documentation of this file.

#define NONE 0
#define RING 1
#define BUTTERFLY 2
#define NONBLOCKING 3 
#define NOEMPTY 4 
#define ALLTOALLV 5 
#define ALLREDUCE 6
#define SEQ 0 
#define OMP 1
#define GPU 2
 


typedef struct { 
  int           flag;                   // flag for communication scheme (NONE, RING, BUTTERFLY ...)
  int           m;                      // number local rows 
  int           nnz;                    // number non-zero per rows 
  int           *indices;               // column indices tab; size = m * nnz; can be a global or local numbering
  double        *values;                // non-zero values tab; size = m * nnz
  //--------local shaping---------------
  int           lcount;
  int           *lindices;              // local indices tab (monotony with global numbering); size = lcount  
#ifdef W_MPI
  MPI_Comm      comm;                   // MPI communicator
  //--------com shaping-----------------
  int           *com_indices, com_count;// communicated indices tab, and size
  int           steps;                  // number of steps in the communication scheme
  int           *nS, *nR;               // number of indices (to send and to receive); size = steps
  int           **R, **S;               // sending or receiving indices tab
#endif
}Mat;


int MatInit(Mat *A, int m, int nnz, int *indices, double *values, int flag
#ifdef W_MPI
, MPI_Comm comm
#endif
);

void MatSetIndices(Mat *A, int m, int nnz, int *indices);

void MatSetValues(Mat *A, int m, int nnz, double *values);

void MatFree(Mat *A);                                           

int MatLocalShape(Mat *A, int sflag);

#if W_MPI
int MatComShape(Mat *A, int flag,  MPI_Comm comm);
#endif

int MatVecProd(Mat *A, double *x, double *y, int pflag);

int TrMatVecProd(Mat *A, double *y, double* x, int pflag);


#if W_MPI
int TrMatVecProd_Naive(Mat *A, double *y, double* x, int pflag);
#endif

int MatLoad(Mat *A, char *filename);

int MatSave(Mat *A, char* filename);

#if W_MPI
int MatInfo(Mat *A, int master, char* filename);
#endif

int greedyreduce(Mat *A, double* x);


// Doxygen definitions