Some more timing output

bioem.cpp

@@ -95,6 +95,8 @@ int bioem::configure(int ac, char* av[])
 	// ****** And Precalculating necessary grids, map crosscorrelations and kernels  ********
 	// *************************************************************************************
 
+	HighResTimer timer;
+
 	if (mpi_rank == 0)
 	{
 		// *** Inizialzing default variables ***
@@ -207,11 +209,15 @@ int bioem::configure(int ac, char* av[])
 			cout << e.what() << "\n";
 			return 1;
 		}
-			//check for consitency in multiple MRCs
-			if(  RefMap.readMultMRC && not(RefMap.readMRC) ){
-			 cout << "For Multiple MRCs command --ReadMRC is necesary too";
-			 exit(1);
-			}
+
+		//check for consitency in multiple MRCs
+		if(RefMap.readMultMRC && not(RefMap.readMRC))
+		{
+			cout << "For Multiple MRCs command --ReadMRC is necesary too";
+			exit(1);
+		}
+
+		if (DebugOutput >= 2 && mpi_rank == 0) timer.ResetStart();
 		// ********************* Reading Parameter Input ***************************
 		param.readParameters(infile.c_str());
 
@@ -220,9 +226,9 @@ int bioem::configure(int ac, char* av[])
 
 		// ********************* Reading Particle Maps Input **********************
 		RefMap.readRefMaps(param, mapfile.c_str());
+		if (DebugOutput >= 2 && mpi_rank == 0) printf("Reading Input Data %f\n", timer.GetCurrentElapsedTime());
 	}
 
-	HighResTimer timer;
 #ifdef WITH_MPI
 	if (mpi_size > 1)
 	{
@@ -237,7 +243,7 @@ int bioem::configure(int ac, char* av[])
 		MPI_Bcast(&RefMap, sizeof(RefMap), MPI_BYTE, 0, MPI_COMM_WORLD);
 		if (mpi_rank != 0) RefMap.maps = (myfloat_t*) mallocchk(RefMap.refMapSize * sizeof(myfloat_t) * RefMap.ntotRefMap);
 		MPI_Bcast(RefMap.maps, RefMap.refMapSize * sizeof(myfloat_t) * RefMap.ntotRefMap, MPI_BYTE, 0, MPI_COMM_WORLD);
-		if (DebugOutput >= 2 && mpi_rank == 0) printf("MPI Broadcast Data %f\n", timer.GetCurrentElapsedTime());
+		if (DebugOutput >= 2 && mpi_rank == 0) printf("MPI Broadcast of Input Data %f\n", timer.GetCurrentElapsedTime());
 	}
 	#endif