liblinbox-dev 1.4.2-3 / usr / include / linbox / util / matrix-stream.inl

/* Copyright (C) LinBox
 *
 *
 *
 * ========LICENCE========
 * This file is part of the library LinBox.
 *
  * LinBox is free software: you can redistribute it and/or modify
 * it under the terms of the  GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.	 See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 * ========LICENCE========
 */

#ifndef __LINBOX_util_matrix_stream_INL
#define __LINBOX_util_matrix_stream_INL

#include "linbox/util/formats/matrix-stream-readers.h"

namespace LinBox
{

	template<class Field>
	bool MatrixStream<Field>::readWhiteSpace()
	{
		char c;
		while(in.get(c)) {
			if (isspace(c))
				switch(c) {
				case '\n': ++lineNumber; break;
				case '\r': if( in.peek() != '\n' ) ++lineNumber; break;
				}
			else {
				in.putback(c);
				return true;
			}

		}
		return false; // because eof or read error.
	}

	template<class Field>
	void MatrixStreamReader<Field>::saveTriple(size_t m, size_t n, const Element& v )
	{
		static std::pair<std::pair<size_t,size_t>,Element> temp;
		temp.first.first = m;
		temp.first.second = n;
		temp.second = v;
		savedTriples.push(temp);
	}

	template<class Field>
	MatrixStreamError MatrixStreamReader<Field>::init
	(const char* firstLine, std::istream* i, MatrixStream<Field>* m )
	{
		if( !i || !m || !firstLine ) throw "Bad istream or MatrixStream";
		sin = i;
		ms = m;
		return initImpl(firstLine);
	}

	template<class Field>
	MatrixStreamError MatrixStreamReader<Field>::nextTriple
	(size_t& m, size_t& n, Element& v)
	{
		if( savedTriples.size() == 0 ) {
			if( atEnd ) {
				if( lastError <= GOOD ) lastError = END_OF_MATRIX;
				return lastError;
			}
			if( lastError > GOOD ) return lastError;
			lastError =  nextTripleImpl(m,n,v);
			return lastError;
		}
		m = savedTriples.front().first.first;
		n = savedTriples.front().first.second;
		v = savedTriples.front().second;
		savedTriples.pop();
		return GOOD;
	}

	template<class Field>
	MatrixStreamError MatrixStreamReader<Field>::getArray
	(std::vector<Element> &array)
	{
		MatrixStreamError mse = GOOD;
		size_t c = 0,i,j;
		Element v;

		while( true ) {
			mse = nextTriple(i,j,v);
			if( mse > GOOD ) break;
			if( i > 0 ) {
				mse = getColumns(c);
				if( mse > GOOD ) break;
			}
			size_t loc = i*c+j;
			if( loc >= array.size() )
				array.resize(c ? (i+1)*c : loc+1);
			array[loc] = v;
		}
		if( mse > END_OF_MATRIX ) return mse;
		mse = getRows(j);
		if( mse > END_OF_MATRIX ) return mse;
		if( array.size() < j*c ) array.resize(j*c);
		return GOOD;
	}

	template<class Field>
	MatrixStreamError MatrixStreamReader<Field>::saveNext()
	{
		if( lastError > GOOD ) return lastError;
		if( atEnd ) {
			lastError = END_OF_MATRIX;
			return lastError;
		}
		size_t m, n;
		Element v;
		lastError = nextTripleImpl(m,n,v);
		if( lastError <= GOOD ) saveTriple(m,n,v);
		return lastError;
	}

	template<class Field>
	MatrixStreamError MatrixStreamReader<Field>::getRows(size_t& m)
	{
		MatrixStreamError toRet = GOOD;
		while( !knowM ) {
			if( atEnd ) return END_OF_MATRIX;
			toRet = saveNext();
			m = _m; // when no entries, it is good
			if( toRet > GOOD ) return toRet;
		}
		m = _m;
		return toRet;
	}

	template<class Field>
	MatrixStreamError MatrixStreamReader<Field>::getColumns(size_t& n)
	{
		MatrixStreamError toRet = GOOD;
		while( !knowN ) {
			if( atEnd ) return END_OF_MATRIX;
			toRet = saveNext();
			n = _n; // when no entries, it is good
			if( toRet > GOOD ) return toRet;
		}
		n = _n;
		return toRet;
	}

	template<class Field>
	void MatrixStream<Field>::init()
	{
		lineNumber = 1;

		//Skip comments
		readWhiteSpace();
		while( !in.eof() && in.peek() == '#' ) {
			char c;
			while( in.get(c) ) {
				if( c == '\n' ) break;
				if( c == '\r' ) {
					if( in.peek() == '\n' )
						in.get();
					break;
				}
			}
			++lineNumber;
			readWhiteSpace();
		}

		//Get first line
		firstLine = new char[FIRST_LINE_LIMIT];
		in.getline(firstLine,FIRST_LINE_LIMIT);
		firstLine[in.gcount()] = '\0';

		//Initialize readers
		currentError = NO_FORMAT;
		__MATRIX_STREAM_READERDEFS;
		delete[] firstLine;

		if( !reader ) return;
		else if( currentError > GOOD )
			errorLineNumber = lineNumber;
	}

	template<class Field>
	void MatrixStream<Field>::addReader( MatrixStreamReader<Field>* r )
	{
		if( currentError == GOOD ) {
			delete r;
			return;
		}

		MatrixStreamError mse = r->init( firstLine, &in, this );
		if( mse < currentError ) {
			if( reader ) delete reader;
			reader = r;
			currentError = mse;
		}
		else delete r ;
		return ;
	}

	template<class Field>
	MatrixStream<Field>::MatrixStream(const Field& fld, std::istream& i ) :
		reader(NULL),in(i),readAnythingYet(false),f(fld)
	{
		init();
		if( currentError > GOOD ){
#ifndef NDEBUG
			reportError(__func__,__LINE__);
#endif
			throw currentError;
		}
	}

	template<class Field>
	void MatrixStream<Field>::newmatrix()
	{
		readAnythingYet = false;
		init();
		if( currentError > GOOD ) throw currentError;
	}

	template<class Field>
	bool MatrixStream<Field>::nextTriple(size_t& m, size_t& n, Element& v)
	{
		if( currentError > GOOD ) return false;

		do {
			currentError = reader->nextTriple(m,n,v);
		} while( f.isZero(v) && currentError == GOOD );

		if( currentError != GOOD ) {
			errorLineNumber = lineNumber;
			return false;
		}

		readAnythingYet = true;
		return true;
	}

	template<class Field>
	bool MatrixStream<Field>::getArray(std::vector<Element> &array)
	{
		if( currentError > GOOD || readAnythingYet ) return false;
		currentError = reader->getArray(array);

		if( currentError != GOOD ) {
			errorLineNumber = lineNumber;
			return false;
		}

		readAnythingYet = true;
		return true;
	}

	template<class Field>
	bool MatrixStream<Field>::getRows(size_t& m)
	{
		MatrixStreamError mse = reader->getRows(m);

		if( currentError > GOOD )
			return (mse == GOOD);
		else if( mse > GOOD ) {
			currentError = mse;
			errorLineNumber = lineNumber;
			return false;
		}
		else return true;
	}

	template<class Field>
	bool MatrixStream<Field>::getColumns(size_t& n)
	{
		MatrixStreamError mse = reader->getColumns(n);

		if( currentError > GOOD )
			return (mse == GOOD);
		else if( mse > GOOD ) {
			currentError = mse;
			errorLineNumber = lineNumber;
			return false;
		}
		else return true;
	}

	template<class Field>
	bool MatrixStream<Field>::getDimensions( size_t& m, size_t& n )
	{
		bool r = getRows(m);

		bool c = getColumns(n); // need getColumns to be called

		return( r && c );
	}

	template<class Field>
	MatrixStreamError MatrixStream<Field>::reportError
	( const char* func, int line ) const
	{
		std::cerr << std::endl
		<< "ERROR (" << func << ":" << line << "): "
			<< "Problem reading matrix:" << std::endl;
		switch( getError() ) {
		case END_OF_MATRIX:
			std::cerr << "There is no more data in the matrix file.";
			break;
		case END_OF_FILE:
			std::cerr << "An EOF was encountered unexpectedly in reading the data.";
			break;
		case BAD_FORMAT:
			std::cerr << "There is a formatting error in the matrix.";
			break;
		case NO_FORMAT:
			std::cerr << "The matrix format is not recognized or supported.";
			break;
		case GOOD: break;
		default: break;
		}
		std::cerr << std::endl << "At line number: " << lineNumber << std::endl
		<< "Matrix format is " << getFormat() << std::endl;
		return currentError;
	}

	template<class Field>
	int MatrixStream<Field>::getLineNumber() const {
		if( currentError > GOOD ) return errorLineNumber;
		else return lineNumber;
	}

} // end of namespace LinBox

#endif // __LINBOX_util_matrix_stream_INL

// Local Variables:
// mode: C++
// tab-width: 8
// indent-tabs-mode: nil
// c-basic-offset: 8
// End:
// vim:sts=8:sw=8:ts=8:noet:sr:cino=>s,f0,{0,g0,(0,\:0,t0,+0,=s