matrix_coo.cc
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include "matrix_coo.hh"
	extern "C" {
	#include "mmio.h"
	}

	void MatrixCOO::read(const std::string & fn) {
	int nz;

	int ret_code;
	MM_typecode matcode;
	FILE * f;

	if ((f = fopen(fn.c_str(), "r")) == NULL) {
	printf("Could not open matrix");
	exit(1);
	}

	if (mm_read_banner(f, &matcode) != 0) {
	printf("Could not process Matrix Market banner.\n");
	exit(1);
	}

	// Matrix is sparse
	if (not(mm_is_matrix(matcode) and mm_is_coordinate(matcode))) {
	printf("Sorry, this application does not support ");
	printf("Market Market type: [%s]\n", mm_typecode_to_str(matcode));
	exit(1);
	}

	if ((ret_code = mm_read_mtx_crd_size(f, &m_m, &m_n, &nz)) != 0) {
	exit(1);
	}

	/* reserve memory for matrices */
	irn.resize(nz);
	jcn.resize(nz);
	a.resize(nz);

	/* NOTE: when reading in doubles, ANSI C requires the use of the "l" */
	/* specifier as in "%lg", "%lf", "%le", otherwise errors will occur */
	/* (ANSI C X3.159-1989, Sec. 4.9.6.2, p. 136 lines 13-15) */
	m_is_sym = mm_is_symmetric(matcode);
	for (int i = 0; i < nz; i++) {
	int I, J;
	double a_;

	fscanf(f, "%d %d %lg\n", &I, &J, &a_);
	I--; /* adjust from 1-based to 0-based */
	J--;

	irn[i] = I;
	jcn[i] = J;
	a[i] = a_;
	}

	if (f != stdin) {
	fclose(f);
	}
	}