use debug code only when environment variable set

bioem.cpp

@@ -37,12 +37,12 @@ ostream& operator<<(ostream& os, const vector<T>& v)
 
 bioem::bioem()
 {
-	
+
 }
 
 bioem::~bioem()
 {
-	
+
 }
 
 int bioem::configure(int ac, char* av[])
@@ -119,18 +119,18 @@ int bioem::configure(int ac, char* av[])
             cout << "Reading model file in PDB format.\n";
             Model.readPDB=true;
         }
-        
+
         if (vm.count("DumpMaps"))
         {
             cout << "Dumping Maps after reading from file.\n";
             RefMap.dumpMap = true;
         }
-        
+
         if (vm.count("LoadMapDump"))
         {
             cout << "Loading Map dump.\n";
             RefMap.loadMap = true;
-        }        
+        }
 
         if (vm.count("Particlesfile"))
         {
@@ -163,11 +163,13 @@ int bioem::configure(int ac, char* av[])
 
     /****************** Precalculating Necessary Stuff *********************/
     precalculate();
-	
-	param.nTotGridAngles = 10;
-	param.nTotCTFs = 10;
-	//param.param_device.maxDisplaceCenter = 0;
-	
+
+	if (getenv("BIOEM_DEBUG_BREAK"))
+	{
+		param.nTotGridAngles = atoi(getenv("BIOEM_DEBUG_BREAK"));
+		param.nTotCTFs = atoi(getenv("BIOEM_DEBUG_BREAK"));
+	}
+
 	deviceInit();
 
     return(0);
@@ -231,7 +233,7 @@ int bioem::run()
 	deviceStartRun();
 
     /******************************** MAIN CYCLE ******************************************/
-	
+
     /*** Declaring Private variables for each thread *****/
     mycomplex_t* proj_mapFFT;
     bioem_map conv_map;
@@ -280,7 +282,7 @@ int bioem::run()
     }
     //deallocating fftw_complex vector
     fftw_free(proj_mapFFT);
-	
+
 	deviceFinishRun();
 
     /************* Writing Out Probabilities ***************/