Bug fix for non-mpi runs

This fixes a bug that was introduced at some point for non-mpi VR. There was the incorrect addressing of an array using a index i rather than the correct Part[i].GetID() index. This has now been fixed. MPI did not have this error as i=Part[i].GetID() for the mpi runs at this point.

Conflicts:
	src/search.cxx

src/search.cxx

@@ -248,7 +248,7 @@ Int_t* SearchFullSet(Options &opt, const Int_t nbodies, vector<Particle> &Part,
 // #else
             numingroup=BuildNumInGroup(nbodies, numgroups, pfof);
             for (i=0;i<nbodies;i++) {
-                Part[i].SetType((numingroup[pfof[i]]>=MINSUBSIZE));
+                Part[i].SetType((numingroup[pfof[Part[i].GetID()]]>=MINSUBSIZE));
                 numlocalden += (Part[i].GetType()>0);
             }
             if (opt.fofbgtype>FOF6D) {
@@ -270,6 +270,7 @@ Int_t* SearchFullSet(Options &opt, const Int_t nbodies, vector<Particle> &Part,
         }
         for (i=0;i<nbodies;i++) {numinstrucs+=(pfof[i]>0);}
         LOG(debug) << "Number of particles in large subhalo searchable structures " << numinstrucs;
+        LOG(debug) << "Number of particles that will have local velocity densitites calculated "<< numlocalden;
         if (numlocalden>0) GetVelocityDensity(opt, nbodies, Part.data(), tree);
         for (i=0;i<nbodies;i++) Part[i].SetType(storetype[i]);
         delete[] storetype;