00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00038 #ifndef MM_IO_H
00039 #define MM_IO_H
00040
00041 #define MM_MAX_LINE_LENGTH 1025
00042 #define MatrixMarketBanner "%%MatrixMarket"
00043 #define MM_MAX_TOKEN_LENGTH 64
00044
00045 #ifdef __cplusplus
00046 #define EXTERN_C extern "C"
00047 #else
00048 #define EXTERN_C
00049 #endif
00050
00051
00052 typedef char MM_typecode[4];
00053
00054 EXTERN_C char *mm_typecode_to_str(MM_typecode matcode);
00055
00056 EXTERN_C int mm_read_banner(FILE *f, MM_typecode *matcode);
00057 EXTERN_C int mm_read_mtx_crd_size(FILE *f, int *M, int *N, int *nz);
00058 EXTERN_C int mm_read_mtx_array_size(FILE *f, int *M, int *N);
00059
00060 EXTERN_C int mm_write_banner(FILE *f, MM_typecode matcode);
00061 EXTERN_C int mm_write_mtx_crd_size(FILE *f, int M, int N, int nz);
00062 EXTERN_C int mm_write_mtx_array_size(FILE *f, int M, int N);
00063
00064
00065
00066
00067 #define mm_is_matrix(typecode) ((typecode)[0]=='M')
00068
00069 #define mm_is_sparse(typecode) ((typecode)[1]=='C')
00070 #define mm_is_coordinate(typecode)((typecode)[1]=='C')
00071 #define mm_is_dense(typecode) ((typecode)[1]=='A')
00072 #define mm_is_array(typecode) ((typecode)[1]=='A')
00073
00074 #define mm_is_complex(typecode) ((typecode)[2]=='C')
00075 #define mm_is_real(typecode) ((typecode)[2]=='R')
00076 #define mm_is_pattern(typecode) ((typecode)[2]=='P')
00077 #define mm_is_integer(typecode) ((typecode)[2]=='I')
00078
00079 #define mm_is_symmetric(typecode)((typecode)[3]=='S')
00080 #define mm_is_general(typecode) ((typecode)[3]=='G')
00081 #define mm_is_skew(typecode) ((typecode)[3]=='K')
00082 #define mm_is_hermitian(typecode)((typecode)[3]=='H')
00083
00084 int mm_is_valid(MM_typecode matcode);
00085
00086
00087
00088
00089 #define mm_set_matrix(typecode) ((*typecode)[0]='M')
00090 #define mm_set_coordinate(typecode) ((*typecode)[1]='C')
00091 #define mm_set_array(typecode) ((*typecode)[1]='A')
00092 #define mm_set_dense(typecode) mm_set_array(typecode)
00093 #define mm_set_sparse(typecode) mm_set_coordinate(typecode)
00094
00095 #define mm_set_complex(typecode)((*typecode)[2]='C')
00096 #define mm_set_real(typecode) ((*typecode)[2]='R')
00097 #define mm_set_pattern(typecode)((*typecode)[2]='P')
00098 #define mm_set_integer(typecode)((*typecode)[2]='I')
00099
00100
00101 #define mm_set_symmetric(typecode)((*typecode)[3]='S')
00102 #define mm_set_general(typecode)((*typecode)[3]='G')
00103 #define mm_set_skew(typecode) ((*typecode)[3]='K')
00104 #define mm_set_hermitian(typecode)((*typecode)[3]='H')
00105
00106 #define mm_clear_typecode(typecode) ((*typecode)[0]=(*typecode)[1]= \
00107 (*typecode)[2]=' ',(*typecode)[3]='G')
00108
00109 #define mm_initialize_typecode(typecode) mm_clear_typecode(typecode)
00110
00111
00112
00113
00114
00115 #define MM_COULD_NOT_READ_FILE 11
00116 #define MM_PREMATURE_EOF 12
00117 #define MM_NOT_MTX 13
00118 #define MM_NO_HEADER 14
00119 #define MM_UNSUPPORTED_TYPE 15
00120 #define MM_LINE_TOO_LONG 16
00121 #define MM_COULD_NOT_WRITE_FILE 17
00122
00123
00124
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140 #define MM_MTX_STR "matrix"
00141 #define MM_ARRAY_STR "array"
00142 #define MM_DENSE_STR "array"
00143 #define MM_COORDINATE_STR "coordinate"
00144 #define MM_SPARSE_STR "coordinate"
00145 #define MM_COMPLEX_STR "complex"
00146 #define MM_REAL_STR "real"
00147 #define MM_INT_STR "integer"
00148 #define MM_GENERAL_STR "general"
00149 #define MM_SYMM_STR "symmetric"
00150 #define MM_HERM_STR "hermitian"
00151 #define MM_SKEW_STR "skew-symmetric"
00152 #define MM_PATTERN_STR "pattern"
00153
00154
00155
00156
00157 int mm_write_mtx_crd(char fname[], int M, int N, int nz, int I[], int J[],
00158 double val[], MM_typecode matcode);
00159 int mm_read_mtx_crd_data(FILE *f, int M, int N, int nz, int I[], int J[],
00160 double val[], MM_typecode matcode);
00161 int mm_read_mtx_crd_entry(FILE *f, int *I, int *J, double *real, double *img,
00162 MM_typecode matcode);
00163
00164 int mm_read_unsymmetric_sparse(const char *fname, int *M_, int *N_, int *nz_,
00165 double **val_, int **I_, int **J_);
00166
00167
00168
00169 #endif