Add != Zeros overload for Matrix, Vector.

Makefile.in

@@ -44,7 +44,7 @@ docs:
 	doxygen 
 
 
-TESTS=lu slice vector_resize gauss_jordan eigen-sqrt determinant chol_toon chol_lapack simplex sym_eigen fill so3 complex qr gr_svd diagonal_matrix gaussian_elimination
+TESTS=lu slice vector_resize gauss_jordan eigen-sqrt determinant chol_toon chol_lapack simplex sym_eigen fill so3 complex qr gr_svd diagonal_matrix gaussian_elimination zeros
 
 
 TEST_RESULT=$(TESTS:%=regressions/%.result)

internal/matrix.hh

@@ -283,6 +283,13 @@ public:
 		      return 1;
 	    return 0;
 	}
+
+	template<class Op>
+	bool operator!=(const Operator<Op>& op)
+	{
+		return op.notequal(*this);
+	}
+
 
 	///@}
 

internal/objects.h

@@ -145,6 +145,25 @@ template<> struct Operator<Internal::Zero> {
 	}
 	///@}
 
+	template<int R, int C, class P, class B>
+	bool notequal(Matrix<R,C,P,B>& m) const {
+		for(int r=0; r<m.num_rows(); r++)
+			for(int c=0; c<m.num_cols(); c++)
+				if(m[r][c] != 0)
+					return 1;
+
+		return 0;
+	}
+
+
+	template<int S, class P, class B>
+	bool notequal(Vector<S,P,B>& v) const {
+		for(int i=0; i<v.size(); i++)
+			if(v[i] != 0)
+				return 1;
+		return 0;
+	}
+
 	///Generate a sized Zero object for constructing dynamic vectors.
 	Operator<Internal::SizedZero> operator()(int s);
 	///Generate a sized Zero object for constructing dynamic matrices.

internal/vector.hh

@@ -357,7 +357,14 @@ public:
 		    return 1;
 		return 0;
 	}
-
+
+
+	template<class Op>
+	bool operator!=(const Operator<Op>& op)
+	{
+		return op.notequal(*this);
+	}		
+
 	/// @}
 
 	/// @name Misc

regressions/zeros.cc

@@ -0,0 +1,28 @@
+#include "regressions/regression.h"
+
+int main()
+{
+	Vector<5> v1 = Zeros;
+	cout << v1 << endl;
+
+	Vector<> v2 = Zeros(3);
+	cout << v2 << endl;
+
+	Matrix<3,6> m1 = Zeros;
+	cout << m1 << endl;
+
+	Matrix<> m2 = Zeros(2,3);
+	cout << m2 << endl;
+
+
+	Vector<4> v3 = Zeros;
+	cout << (v3 != Zeros) << endl;
+
+	v3 = Ones;
+	cout << (v3 != Zeros) << endl;
+
+	Matrix<4> m3 = Zeros;
+	cout << (m3 != Zeros) << endl;
+	m3 = Ones;
+	cout << (m3 != Zeros) << endl;
+}

regressions/zeros.txt

@@ -0,0 +1,15 @@
+0 0 0 0 0
+
+0 0 0
+
+0 0 0 0 0 0
+0 0 0 0 0 0
+0 0 0 0 0 0
+
+0 0 0
+0 0 0
+
+0
+1
+0
+1