#ifndef _PrintOp
#define _PrintOp



// Overload the cout operator for vectors
template <class T>
ostream & operator << (ostream &os, const Vector<T> &v)
{
  	int n = v.Size();

   	#ifndef _VectorTabDelim
   		if(n > 0)
		{
			/* Print Vector as comma-delimited */
			os << "[";
			for(int i=0; i<n-1; i++) os << v[i] << ",";
			os << v[n-1] << "]";
		}
		else os << "[]";
	#else
		if(n > 0)
		{
			/* Print Vector as tab-delimited */
			os << "[";
			for(int i=0; i<n-1; i++) os << v[i] << "\t";
			os << v[n-1] << "]";
		}
		else os << "[]";
	#endif

   	return os;
}


/* Overload the cout operator for Matrix objects */
template <class T>
ostream & operator << (ostream &os, const Matrix<T> &A)
{
	int m = A.Rows();
	int n = A.Cols();
	if(m > 0 && n > 0)
	{
		os << "\n";
		for(int i=0; i<m-1; i++)
		{
			os << "[ ";
			for(int j=0; j<n-1; j++) os << A(i,j) << "\t";
			os << A(i,n-1) << "]\n";
		}
		os << "[ ";
		for(int j=0; j<n-1; j++) os << A(m-1,j) << "\t";
		os << A(m-1,n-1) << "]";
	}
	else os << "[]";
	return os;
}

#endif