semisparsemat ostream output

sparsemat.h

@@ -326,10 +326,9 @@ return *this;
 template <typename T>
 class SemiSparseMat {
 	friend class NRMat<T>;
-protected:
+public:
 	NRVec<T> diagonal;
 	SparseMat<T> offdiagonal;
-public:
         SemiSparseMat() {};
         SemiSparseMat(const SPMatindex n, const SPMatindex m) : offdiagonal(n,m) {if(n==m) diagonal.resize(n);};
         SPMatindex nrows() const {return offdiagonal.nrows();}
@@ -361,6 +360,14 @@ public:
 	void resize(const SPMatindex n, const SPMatindex m) {offdiagonal.resize(n,m); diagonal.resize(n==m?n:0);}
 };
 
+template <class T>
+std::ostream& operator<<(std::ostream &s, const SemiSparseMat<T> &x)
+                {
+		s << "diagonal= "<<x.diagonal;
+		s << "offdiagonal= "<<x.offdiagonal;
+                return s;
+                }
+
 
 }//namespace
 #endif

t.cc

@@ -4639,6 +4639,7 @@ davidson(aa,r,(NRVec<double> *)NULL,"eivecs",m,1,1e-5,0,300,300);
 cout <<r;
 if(n<=20)
 	{
+	cout<<aa;
 	NRMat a(aa);
 	NRVec<double> b(n);
 	cout <<a;