ret_val.elem().elem().elem().elem() = lsum; return ret_val;
return lsum;
norm2(const QDPType<TVec ,OLattice< TVec > >& s1){#ifdef DEBUG_BLAS  QDPIO::cout << "Using BJ sumsq all" << endl;#endif  int n_3vec = (all.end() - all.start() + 1)*Ns;  const REAL *s1ptr =  &(s1.elem(all.start()).elem(0).elem(0).real());      UnaryReturn< OLattice< TVec >, FnNorm2>::Type_t  ret_val;  DOUBLE lsum = 0;   local_sumsq(&lsum, (REAL *)s1ptr, n_3vec);   Internal::globalSum(lsum);  ret_val.elem().elem().elem().elem() = lsum;  return ret_val;}