I have made the Cshift work successfully with open mp threading in

every routine. Collapse(2) is now working under clang-omp++.

benchmarks/Grid_wilson.cc

@@ -77,7 +77,7 @@ int main (int argc, char ** argv)
   WilsonMatrix Dw(Umu,mass);
   
   std::cout << "Calling Dw"<<std::endl;
-  int ncall=10000;
+  int ncall=1000;
   double t0=usecond();
   for(int i=0;i<ncall;i++){
     Dw.multiply(src,result);