Merge branch 'master' of github.com:lattice/quda

lattice · Oct 13, 2011 · 3de6e8f · 3de6e8f
2 parents 24e2e03 + 7be4de4
commit 3de6e8f
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 3 deletions.
diff --git a/lib/inv_gcr_quda.cpp b/lib/inv_gcr_quda.cpp
@@ -168,7 +168,8 @@ void GCR::operator()(cudaColorSpinorField &x, cudaColorSpinorField &b)
 
   // create sloppy fields used for orthogonalization
   param.precision = invParam.cuda_prec_sloppy;
-  cudaColorSpinorField *p[Nkrylov], *Ap[Nkrylov];
+  cudaColorSpinorField **p = new cudaColorSpinorField*[Nkrylov];
+  cudaColorSpinorField **Ap = new cudaColorSpinorField*[Nkrylov];
   for (int i=0; i<Nkrylov; i++) {
     p[i] = new cudaColorSpinorField(x, param);
     Ap[i] = new cudaColorSpinorField(x, param);
@@ -367,6 +368,8 @@ void GCR::operator()(cudaColorSpinorField &x, cudaColorSpinorField &b)
     delete p[i];
     delete Ap[i];
   }
+  delete[] p;
+  delete[] Ap;
 
   delete alpha;
   for (int i=0; i<Nkrylov; i++) delete []beta[i];

diff --git a/tests/test_util.cpp b/tests/test_util.cpp
@@ -1053,7 +1053,7 @@ createSiteLinkCPU(void** link,  QudaPrecision precision, int phase)
 
 template <typename Float>
 int compareLink(Float **linkA, Float **linkB, int len) {
-  int fail_check = 16;
+  const int fail_check = 16;
   int fail[fail_check];
   for (int f=0; f<fail_check; f++) fail[f] = 0;
 
@@ -1207,7 +1207,7 @@ createHwCPU(void* hw,  QudaPrecision precision)
 
 template <typename Float>
 void compare_mom(Float *momA, Float *momB, int len) {
-  int fail_check = 16;
+  const int fail_check = 16;
   int fail[fail_check];
   for (int f=0; f<fail_check; f++) fail[f] = 0;