diff --git a/conjgrad.h b/conjgrad.h
index 27c61ad..5775b20 100644
--- a/conjgrad.h
+++ b/conjgrad.h
@@ -46,7 +46,7 @@ double t,tt,bscal,ascal;
 
 NRVec<T> p,rr, *r;
 NRVec<T> q(m),s(m);
-if(issquare) r=&s; else r = new NRVec<T>(m);
+if(issquare) r=&s; else r = new NRVec<T>(n); //A^T s
 
 if(doguess)
 	{