Changeset 16242

Timestamp:

09/24/13 13:25:58 (11 years ago)

Author:

Eric.Larour

Message:

CHG fixing check + verbosing tape statistics

File:

• issm/trunk-jpl/src/c/analyses/ad_core.cpp (modified) (8 diffs)

issm/trunk-jpl/src/c/analyses/ad_core.cpp

@@ -41 +41 @@
 
                 #ifdef _HAVE_ADOLC_
+                        if(VerboseAutodiff())_printf0_("   start ad core\n");
 
                         /*First, stop tracing: */
                         trace_off();
 
-                        /*preliminary checks: */
+                        /*retrieve num_dependents and num_independents only on rank 0: */
                         if (my_rank==0) {
                                 femmodel->parameters->FindParam(&num_dependents,AutodiffNumDependentsEnum);
                                 femmodel->parameters->FindParam(&num_independents,AutodiffNumIndependentsEnum);
-                                if(!(num_dependents*num_independents)) return;
                         }
-
-                        if(VerboseAutodiff())_printf0_("   start ad core\n");
+                                
+                        /*if we don't have dependents or independents, no point in running ad_core: */
+                        if(!(num_dependents*num_independents)) return;
+
+                        /*Print tape statistics so that user can kill this run if something is off already:*/
+                        tapestats(my_rank,tape_stats); //reading of tape statistics
+                        if(VerboseAutodiff()){
+                                _printf0_("   ADOLC statistics: \n");
+                                _printf0_("   "<<setw(45)<<left<<"Number of independents: " <<tape_stats[0] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<"Number of dependents: " <<tape_stats[1] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<"Maximum number of live active variables: " <<tape_stats[2] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<"Number of values in the taylor (value) stack: " <<tape_stats[3] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" Number of operations per buffer == OBUFSIZE (usrparms.h) " <<tape_stats[4] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" overall number of operations " <<tape_stats[5] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" operations file written or not " <<tape_stats[6] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" overall number of locations " <<tape_stats[7] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" locations file written or not " <<tape_stats[8] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" overall number of values " <<tape_stats[9] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" values file written or not " <<tape_stats[10] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" Number of locations per buffer == LBUFSIZE (usrparms.h) " <<tape_stats[11] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" Number of values per buffer == CBUFSIZE (usrparms.h) " <<tape_stats[12] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" Number of taylors per buffer <= TBUFSIZE (usrparms.h) " <<tape_stats[13] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" Number of eq_*_prod for sparsity pattern " <<tape_stats[14] << "\n");
+                                _printf0_("   "<<setw(45)<<left<<" represents the size of the stats vector " <<tape_stats[15] << "\n");
+                        }
+
 
                         /*retrieve state variable: */
@@ -73 +97 @@
                         anEDF_for_solverx_p->dp_y=xNew<double>(anEDF_for_solverx_p->max_m);
 
-                        if (strcmp(driver,"fos_forward")==0){
+                        if (strcmp(driver,"fos_forward")==0){ /*{{{*/
 
                                 int     anIndepIndex;
@@ -110 +134 @@
                                 xDelete(jacTimesTangentDir);
                                 xDelete(tangentDir);
-                        }
-                        else if ((strcmp(driver,"fov_forward")==0) || (strcmp(driver,"fov_forward_all")==0)){
+                        } /*}}}*/
+                        else if ((strcmp(driver,"fov_forward")==0) || (strcmp(driver,"fov_forward_all")==0)){ /*{{{*/
 
                                 int      tangentDirNum;
@@ -181 +205 @@
                                 xDelete(jacTimesSeed);
                                 xDelete(indepIndices);
-                        }
-                        else if (strcmp(driver,"fos_reverse")==0) {
+                        } /*}}}*/
+                        else if (strcmp(driver,"fos_reverse")==0) { /*{{{*/
 
                                 int     aDepIndex=0;
@@ -219 +243 @@
                                 xDelete(weightVectorTimesJac);
                                 xDelete(aWeightVector);
-                        }
-                        else if ((strcmp(driver,"fov_reverse")==0) || (strcmp(driver,"fov_reverse_all")==0)){
+                        } /*}}}*/
+                        else if ((strcmp(driver,"fov_reverse")==0) || (strcmp(driver,"fov_reverse_all")==0)){ /*{{{*/
 
                                 int* depIndices=NULL;
@@ -246 +270 @@
 
                                 /*set the forward method function pointers: */
-#ifdef _HAVE_GSL_
+                                #ifdef _HAVE_GSL_
                                 anEDF_for_solverx_p->fov_reverse=EDF_fov_reverse_for_solverx;
-#endif
+                                #endif
 
                                 /*allocate the space for the parameters to invoke the EDF fos_reverse :*/
@@ -282 +306 @@
                                 xDelete(weightsTimesJac);
                                 xDelete(depIndices);
-                        }
+                        } /*}}}*/
                         else _error_("driver: " << driver << " not yet supported!");
 
@@ -297 +321 @@
                         xDelete(anEDF_for_solverx_p->dpp_Z);
 
-                        /*Print statistics:*/
-                        tapestats(my_rank,tape_stats); //reading of tape statistics
-                        if(VerboseAutodiff()){
-                                _printf0_("   ADOLC statistics: \n");
-                                _printf0_("   "<<setw(45)<<left<<"Number of independents: " <<tape_stats[0] << "\n");
-                                _printf0_("   "<<setw(45)<<left<<"Number of dependents: " <<tape_stats[1] << "\n");
-                                _printf0_("   "<<setw(45)<<left<<"Maximal number of live active variables: " <<tape_stats[2] << "\n");
-                                _printf0_("   "<<setw(45)<<left<<"Size of value stack (number of overwrites): " <<tape_stats[3] << "\n");
-                                _printf0_("   "<<setw(45)<<left<<"Buffer size (a multiple of eight): " <<tape_stats[4] << "\n");
-                                _printf0_("   "<<setw(45)<<left<<"Total number of operations recorded: " <<tape_stats[5] << "\n");
-                        }
                         if(VerboseAutodiff())_printf0_("   end AD core\n");