diff --git a/smat.h b/smat.h
index d642936..d6c0b18 100644
--- a/smat.h
+++ b/smat.h
@@ -92,6 +92,7 @@ inline NRSMat<T>::NRSMat(const T& a, const int n) : nn(n),
 {
 	*count =1;
 	if(a != (T)0) for(int i=0; i<NN2; i++) v[i] = a;
+	else memset(v, 0, NN2*sizeof(T));
 }
 
 template <typename T>