if ( s.hasOrderedRep() ) {
if ( all.hasOrderedRep() ) {
norm2(const QDPType<TVec ,OLattice< TVec > >& s1, const Subset& s){#ifdef DEBUG_BLAS  QDPIO::cout << "Using BJ sumsq" << endl;#endif  if ( s.hasOrderedRep() ) {#ifdef DEBUG_BLAS    QDPIO::cout << "BJ sumsq " << endl;#endif    int n_3vec = (s.end() - s.start() + 1)*Ns;    const REAL *s1ptr =  &(s1.elem(s.start()).elem(0).elem(0).real().elem());        // I am relying on this being a Double here     UnaryReturn< OLattice< TVec >, FnNorm2>::Type_t  lsum;    zero_rep(lsum.elem());        local_sumsq((DOUBLE *)&lsum.elem().elem().elem().elem().elem(),		(REAL *)s1ptr, 	       n_3vec);         Internal::globalSum(lsum.elem());        return lsum;  }  else {   return sum(localNorm2(s1),s);  }}