done with cleanup
authorNikhil Jain <nikhil@illinois.edu>
Tue, 25 Oct 2011 02:56:44 +0000 (21:56 -0500)
committerNikhil Jain <nikhil@illinois.edu>
Tue, 25 Oct 2011 02:56:44 +0000 (21:56 -0500)
defines.h
particleGenerator.cpp
plummer.cpp
util.cpp
util.h

index 7db4aa6..8b81ba0 100644 (file)
--- a/defines.h
+++ b/defines.h
@@ -28,10 +28,8 @@ const double opening_geometry_factor = 2 / sqrt(3.0);
 #define COSMO_CONST(x) (x)
 #define CONVERT_TO_COSMO_TYPE
 
-//#define RRDEBUG CkPrintf
 #define RRDEBUG /* empty */
 
-//#define TB_DEBUG CkPrintf
 #define TB_DEBUG /* empty */
 
 #define NUM_PRIORITY_BITS (sizeof(int)*CHAR_BIT)
index 5d662e7..402289c 100644 (file)
@@ -60,7 +60,6 @@ void generateRandom(ofstream &file, int np){
   file.write((char *)&tnow, sizeof(real));
 
   for(int i = 0; i < np; i++){
-    //printf("writing particle: ");
     // position
     for(int j = 0; j < 3; j++){
       vals[j] = ((1.0*(rand()%max))/max)-0.5;
@@ -68,12 +67,9 @@ void generateRandom(ofstream &file, int np){
     // velocity
     for(int j = 3; j < 6; j++){
       vals[j] = (maxVel-minVel)*((1.0*(rand()%max))/max)+minVel;
-      //vals[j] = ((maxVel-minVel)*(1.0*(rand()%max))/max)+minVel;
     }
     vals[6] = mass;
     vals[7] = soft;
-    //printf("%lf,", mass);
-    //printf("%lf\n", soft);
     printf("%f %f %f %f\n", vals[0], vals[1], vals[2], sqrt(vals[3]*vals[3]+vals[4]*vals[4]+vals[5]*vals[5]));
     file.write((char *)vals,size);
   }
@@ -91,7 +87,6 @@ int main(int argc, char **argv){
   ofstream file;
   file.open(fileName, ios::out|ios::binary);
   generateRandom(file,numParticles);
-  //generateSequential(file,numParticles);
   file.close();
   return 0;
 }
index 0d3d7a0..0b88ac2 100644 (file)
@@ -65,11 +65,8 @@ Particle *testdata(int nbody)
       pickshell(p->velocity, vsc * v);
       cmv += p->velocity;
       
-      //cout << "particle generated" << endl;
    }
 
-   //cout << "done first loop" << endl;
-
    offset = 4.0;
 
    for (p = bodytab + halfnbody; p < bodytab+nbody; p++) {
@@ -97,8 +94,6 @@ Particle *testdata(int nbody)
  */
 
 void pickshell(Vector3D<Real> &vec, Real rad)
-   //Real vec[];                     /* coordinate vector chosen */
-   //Real rad;                       /* radius of chosen point */
 {
    register int k;
    Real rsq, rsc;
@@ -108,7 +103,6 @@ void pickshell(Vector3D<Real> &vec, Real rad)
      vec.y = xrand(-1.0,1.0);
      vec.z = xrand(-1.0,1.0);
      rsq = vec.lengthSquared();
-     //cout << "rsq " << vec.x << "," << vec.y << "," << vec.z << "," << rsq << endl;
    } while (rsq > 1.0);
 
    rsc = rad / sqrt(rsq);
index 3b77034..4744b1d 100644 (file)
--- a/util.cpp
+++ b/util.cpp
@@ -52,12 +52,6 @@ bool CompareParticleToKey(void *a, Key k){
 
 void findSplitters(Particle *particles, int start, int end, int *splitters, Key childKey, int childDepth){
   int nRankBits = LOG_BRANCH_FACTOR;
-  /*
-  CkPrintf("findSplitters parentKey %lx nRankBits %d startChildKey %lx parentDepth %d childDepth %d\n", 
-      parentKey, 
-      nRankBits, childKey,
-      parentDepth, childDepth);
-  */
   // particles of first child always begin at index 0
   splitters[0] = start;
   // for all other children, set the beginning of each 
@@ -65,7 +59,6 @@ void findSplitters(Particle *particles, int start, int end, int *splitters, Key
   for(int i = 1; i < BRANCH_FACTOR; i++){
     childKey++;
     Key testKey = (childKey << (TREE_KEY_BITS-(nRankBits*childDepth+1))); 
-    //CkPrintf("(%d) findSplitters [%d - %d] start %lx end %lx test %lx\n", CkMyPe(), start, end-1, particles[start].key, particles[end-1].key, testKey);
     int firstGEIdx = binary_search_ge<Particle>(testKey, particles, start, end, CompareParticleToKey); 
     splitters[i] = firstGEIdx;
     start = firstGEIdx;
diff --git a/util.h b/util.h
index 1d4f2ad..a8b0b0e 100644 (file)
--- a/util.h
+++ b/util.h
@@ -30,9 +30,7 @@ int binary_search_ge(Key check, T *particles, int start, int end, BinarySearchGE
   return lo;
 }
 
-//template<typename T>
 void getMomentsFromParticles(Node<ForceData> *bucket);
-//template<typename T>
 void getMomentsFromChildren(Node<ForceData> *node);
 void getBoundingBoxFromParticles(Node<ForceData> *node);
 void getBoundingBoxFromChildren(Node<ForceData> *node);