diff --git a/vec.cc b/vec.cc
index 9713039..6425bbd 100644
--- a/vec.cc
+++ b/vec.cc
@@ -969,10 +969,15 @@ return m;\
 
 
 INSTANTIZE_NONCOMPLEX(char)
+INSTANTIZE_NONCOMPLEX(unsigned char)
 INSTANTIZE_NONCOMPLEX(short)
+INSTANTIZE_NONCOMPLEX(unsigned short)
 INSTANTIZE_NONCOMPLEX(int)
+INSTANTIZE_NONCOMPLEX(unsigned int)
 INSTANTIZE_NONCOMPLEX(long)
+INSTANTIZE_NONCOMPLEX(unsigned long)
 INSTANTIZE_NONCOMPLEX(long long)
+INSTANTIZE_NONCOMPLEX(unsigned long long)
 INSTANTIZE_NONCOMPLEX(float)
 INSTANTIZE_NONCOMPLEX(double)
 
diff --git a/vec.h b/vec.h
index 70e0651..9f61b93 100644
--- a/vec.h
+++ b/vec.h
@@ -170,7 +170,16 @@ public:
 	NRVec& operator|=(const NRVec &rhs);
 
 	//! relational operators 
-	const bool operator!=(const NRVec &rhs) const {if(nn!=rhs.nn) return 1; return LA_traits<T>::gencmp(v,rhs.v,nn);}
+	const bool operator!=(const NRVec &rhs) const 
+		{
+		if(nn!=rhs.nn) return 1; 
+		if(LA_traits<T>::is_plaindata()) return LA_traits<T>::gencmp(v,rhs.v,nn);
+		else
+			{
+			for(int i=0; i<nn; ++i) if(v[i]!=rhs.v[i]) return 1;
+			return 0;
+			}
+		}
 	const bool operator==(const NRVec &rhs) const {return !(*this != rhs);};
 	const bool operator>(const NRVec &rhs) const;
 	const bool operator<(const NRVec &rhs) const;