Changeset 11513

Timestamp:

02/23/12 14:29:41 (13 years ago)

Author:

Mathieu Morlighem

Message:

Fixed TAO Monitor for parallel runs

File:

• issm/trunk-jpl/src/c/solutions/controltao_core.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/solutions/controltao_core.cpp

@@ -17 +17 @@
 
 /*Local prototype*/
-int FormFunctionGradient(TaoSolver tao,Vec,double*,Vec,void*);
-int IssmMonitor(TaoSolver,void *);
+int FormFunctionGradient(TaoSolver,Vec,double*,Vec,void*);
+int IssmMonitor(TaoSolver,void*);
 typedef struct {
         FemModel* femmodel;
@@ -60 +60 @@
 
         /*Prepare all TAO parameters*/
-        TaoSetMonitor(tao,IssmMonitor,NULL,NULL);
+        TaoSetMonitor(tao,IssmMonitor,&user,NULL);
         TaoSetMaximumFunctionEvaluations(tao,maxiter);
         TaoSetMaximumIterations(tao,nsteps);
@@ -146 +146 @@
 int IssmMonitor(TaoSolver tao, void *userCtx){
 
-        AppCtx *user = (AppCtx *)userCtx;
-        int     its;
-        double  f,gnorm,cnorm,xdiff;
+        int       i,its,num_responses;
+        double    f,gnorm,cnorm,xdiff;
+        AppCtx   *user      = (AppCtx *)userCtx;
+        FemModel *femmodel  = user->femmodel;
+        Element  *element   = NULL;
+        int      *responses = NULL;
+
+        femmodel->parameters->FindParam(&num_responses,InversionNumCostFunctionsEnum);
+        femmodel->parameters->FindParam(&responses,NULL,NULL,StepResponsesEnum);
 
         TaoGetSolutionStatus(tao, &its, &f, &gnorm, &cnorm, &xdiff, NULL);
-        if(its==0)printf("Iter       Function      Residual\n");
-        printf("%4i   %12.7g  %12.7g\n",its,f,gnorm);
-        printf("iter: %i\n",tao->niter);
+        if(its==0) _printf_(true,"Iter       Function      Residual  |  List of contributions\n");
+        if(its==0) _printf_(true,"-----------------------------------+-----------------------\n");
+        _printf_(true,"%4i   %12.7g  %12.7g  | ",its,f,gnorm);
+
+        /*Retrieve objective functions independently*/
+        for(i=0;i<num_responses;i++){
+                Responsex(&f,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,EnumToStringx(responses[i]),false,i);
+                _printf_(true," %12.7g ","",f);
+        }
+        _printf_(true,"\n");
+
+        /*Clean-up and return*/
+        xfree((void**)&responses);
         return 0;
 }