1 /* 2 * Matrix Market I/O library for ANSI C 3 * 4 * See http://math.nist.gov/MatrixMarket for details. 5 * 6 * 7 */ 8 9 #ifndef MM_IO_H 10 #define MM_IO_H 11 12 #define MM_MAX_LINE_LENGTH 1025 13 #define MatrixMarketBanner "%%MatrixMarket" 14 #define MM_MAX_TOKEN_LENGTH 64 15 16 #include <sys/types.h> 17 18 typedef char MM_typecode[4]; 19 20 char *mm_typecode_to_str(MM_typecode matcode); 21 22 int mm_read_banner(FILE *f, MM_typecode *matcode); 23 int mm_read_mtx_crd_size(FILE *f, int64_t *M, int64_t *N, int64_t *nz, int64_t *numlinesread); // ABAB: Modified 24 int mm_read_mtx_array_size(FILE *f, int *M, int *N); 25 26 int mm_write_banner(FILE *f, MM_typecode matcode); 27 int mm_write_mtx_crd_size(FILE *f, int M, int N, int nz); 28 int mm_write_mtx_array_size(FILE *f, int M, int N); 29 30 31 /********************* MM_typecode query fucntions ***************************/ 32 33 #define mm_is_matrix(typecode) ((typecode)[0]=='M') 34 35 #define mm_is_sparse(typecode) ((typecode)[1]=='C') 36 #define mm_is_coordinate(typecode)((typecode)[1]=='C') 37 #define mm_is_dense(typecode) ((typecode)[1]=='A') 38 #define mm_is_array(typecode) ((typecode)[1]=='A') 39 40 #define mm_is_complex(typecode) ((typecode)[2]=='C') 41 #define mm_is_real(typecode) ((typecode)[2]=='R') 42 #define mm_is_pattern(typecode) ((typecode)[2]=='P') 43 #define mm_is_integer(typecode) ((typecode)[2]=='I') 44 45 #define mm_is_symmetric(typecode)((typecode)[3]=='S') 46 #define mm_is_general(typecode) ((typecode)[3]=='G') 47 #define mm_is_skew(typecode) ((typecode)[3]=='K') 48 #define mm_is_hermitian(typecode)((typecode)[3]=='H') 49 50 int mm_is_valid(MM_typecode matcode); /* too complex for a macro */ 51 52 53 /********************* MM_typecode modify fucntions ***************************/ 54 55 #define mm_set_matrix(typecode) ((*typecode)[0]='M') 56 #define mm_set_coordinate(typecode) ((*typecode)[1]='C') 57 #define mm_set_array(typecode) ((*typecode)[1]='A') 58 #define mm_set_dense(typecode) mm_set_array(typecode) 59 #define mm_set_sparse(typecode) mm_set_coordinate(typecode) 60 61 #define mm_set_complex(typecode)((*typecode)[2]='C') 62 #define mm_set_real(typecode) ((*typecode)[2]='R') 63 #define mm_set_pattern(typecode)((*typecode)[2]='P') 64 #define mm_set_integer(typecode)((*typecode)[2]='I') 65 66 67 #define mm_set_symmetric(typecode)((*typecode)[3]='S') 68 #define mm_set_general(typecode)((*typecode)[3]='G') 69 #define mm_set_skew(typecode) ((*typecode)[3]='K') 70 #define mm_set_hermitian(typecode)((*typecode)[3]='H') 71 72 #define mm_clear_typecode(typecode) ((*typecode)[0]=(*typecode)[1]= \ 73 (*typecode)[2]=' ',(*typecode)[3]='G') 74 75 #define mm_initialize_typecode(typecode) mm_clear_typecode(typecode) 76 77 78 /********************* Matrix Market error codes ***************************/ 79 80 81 #define MM_COULD_NOT_READ_FILE 11 82 #define MM_PREMATURE_EOF 12 83 #define MM_NOT_MTX 13 84 #define MM_NO_HEADER 14 85 #define MM_UNSUPPORTED_TYPE 15 86 #define MM_LINE_TOO_LONG 16 87 #define MM_COULD_NOT_WRITE_FILE 17 88 89 90 /******************** Matrix Market internal definitions ******************** 91 92 MM_matrix_typecode: 4-character sequence 93 94 ojbect sparse/ data storage 95 dense type scheme 96 97 string position: [0] [1] [2] [3] 98 99 Matrix typecode: M(atrix) C(oord) R(eal) G(eneral) 100 A(array) C(omplex) H(ermitian) 101 P(attern) S(ymmetric) 102 I(nteger) K(kew) 103 104 ***********************************************************************/ 105 106 #define MM_MTX_STR "matrix" 107 #define MM_ARRAY_STR "array" 108 #define MM_DENSE_STR "array" 109 #define MM_COORDINATE_STR "coordinate" 110 #define MM_SPARSE_STR "coordinate" 111 #define MM_COMPLEX_STR "complex" 112 #define MM_REAL_STR "real" 113 #define MM_INT_STR "integer" 114 #define MM_GENERAL_STR "general" 115 #define MM_SYMM_STR "symmetric" 116 #define MM_HERM_STR "hermitian" 117 #define MM_SKEW_STR "skew-symmetric" 118 #define MM_PATTERN_STR "pattern" 119 120 121 /* high level routines */ 122 123 int mm_write_mtx_crd(char fname[], int M, int N, int nz, int I[], int J[], 124 double val[], MM_typecode matcode); 125 int mm_read_mtx_crd_data(FILE *f, int M, int N, int nz, int I[], int J[], 126 double val[], MM_typecode matcode); 127 int mm_read_mtx_crd_entry(FILE *f, int *I, int *J, double *real, double *img, 128 MM_typecode matcode); 129 130 131 #endif 132