Added an IsSorted() method for debugging purposes

SpMat.h

@@ -110,6 +110,7 @@ public:
   void Print() const {Print(8);}
   void WarningsOn() {_pw=true;}
   void WarningsOff() {_pw=false;}
+  bool IsSorted() const { return(is_sorted()); }  // Returns true if all _ri arrays are sorted. For debugging.
 
 
   T Peek(unsigned int r, unsigned int c) const;
@@ -220,6 +221,7 @@ private:
   bool same_sparsity(const SpMat<T>& M) const;
   SpMat<T>& add_same_sparsity_mat_to_me(const SpMat<T>& M, double s);
   SpMat<T>& add_diff_sparsity_mat_to_me(const SpMat<T>& M, double s);
+  bool is_sorted() const;
 };
 
 //@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@@ -1222,6 +1224,27 @@ SpMat<T>& SpMat<T>::add_diff_sparsity_mat_to_me(const SpMat<T>& M, double s)
   return(*this);
 }
 
+/////////////////////////////////////////////////////////////////////
+//
+// Returns true if all _ri arrays are sorted in ascending order
+//
+/////////////////////////////////////////////////////////////////////
+
+template<class T>
+bool SpMat<T>::is_sorted() const
+{
+  for (unsigned int ci=0; ci<_ri.size(); ci++) {
+    if (_ri[ci].size() > 1) {
+      unsigned int cri = _ri[ci][0];
+      for (unsigned int ii=1; ii<_ri[ci].size(); ii++) {
+	if (_ri[ci][ii] < cri) return(false);
+	else cri = _ri[ci][ii];
+      }
+    }
+  }
+  return(true);
+}
+
 /*
 template<class T>
 SpMat<T>& SpMat<T>::add_diff_sparsity_mat_to_me(const SpMat<T>& M, double s)