Changeset 2489

Timestamp:

10/22/09 11:53:24 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added ellapsed time in diagnostic parallel

File:

• issm/trunk/src/c/parallel/diagnostic.cpp (modified) (8 diffs)

issm/trunk/src/c/parallel/diagnostic.cpp

@@ -45 +45 @@
         DataSet* parameters=NULL;
 
+        /*time*/
+        double   start, finish;
+        double   start_core, finish_core;
+        double   start_init, finish_init;
+
         MODULEBOOT();
 
@@ -51 +56 @@
         #endif
 
+        /*Initialize Petsc and get start time*/
         PetscInitialize(&argc,&argv,(char *)0,"");  
+        MPI_Barrier(MPI_COMM_WORLD); start=MPI_Wtime();
 
         /*Size and rank: */
@@ -63 +70 @@
 
         /*Initialize model structure: */
+        MPI_Barrier(MPI_COMM_WORLD); start_init=MPI_Wtime();
         model=new Model();
 
@@ -106 +114 @@
         results=new DataSet(ResultsEnum());
 
+        MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
+
         /*are we running the solution sequence, or a qmu wrapper around it? : */
         if(!qmu_analysis){
@@ -111 +121 @@
                         /*run diagnostic analysis: */
                         _printf_("call computational core:\n");
+                        MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
                         diagnostic_core(results,model,inputs);
+                        MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
 
                         /*Add analysis_type to results: */
@@ -123 +135 @@
                         /*run control analysis: */
                         _printf_("call computational core:\n");
+                        MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
                         control_core(results,model,inputs);
+                        MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
 
                         /*Add analysis_type and control_type to results: */
@@ -145 +159 @@
 
                 #ifdef _HAVE_DAKOTA_ 
+                MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
                 Qmux(model,inputs,DiagnosticAnalysisEnum(),NoneAnalysisEnum());
+                MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
                 #else
                 throw ErrorException(__FUNCT__," Dakota not present, cannot do qmu!");
@@ -164 +180 @@
         delete processed_results;
 
+        /*Get finish time and close*/
+        MPI_Barrier(MPI_COMM_WORLD); finish = MPI_Wtime( );
+        _printf_("\n   %-34s %f seconds  \n","Model initialization elapsed time:",finish_init-start_init);
+        _printf_("   %-34s %f seconds  \n","Core solution elapsed time:",finish_core-start_core);
+        _printf_("   %-34s %f seconds\n\n","Total elapsed time:",finish-start);
         _printf_("closing MPI and Petsc\n");
         PetscFinalize();