    DevMemF dev_randTable(neededRands);
    // End Device Memory Pointers
 
+#ifdef DEBUG_TRACE
+   std::cout << "Finished main storage" << std::endl;
+#endif
+
    int percentComplete = 0; // Used to display progress to the user
    int percentSize = 0;
+
+#ifdef DEBUG_TRACE
+   std::cout << "Initializing random number generator" << std::endl;
+#endif
 
    // Set up the random number generator
    curandGenerator_t randGenerator;

      //cudaThreadSynchronize();
      //potent2Timer.stop();

      //fieldTimer.start();
      // calculate E field at Grid points
      field(dev_ex,dev_ey,dev_phi);

      //logger.flush();
      // END DEBUG

#ifdef DEBUG_TRACE
      std::cout << "Dens" << std::endl;
#endif
      //densTimer.start();
      // determine the charge density at the grid points
      dens(dev_rho, dev_rhoe,dev_rhoi, 
