source: issm/oecreview/Archive/16133-16554/ISSM-16483-16484.diff@ 16556

../trunk-jpl/src/c/analyses/controltao_core.cpp

@@ -81 +81 @@
         if(VerboseControl()) _printf0_("   Starting optimization\n");
         TaoSolve(tao);
         TaoView(tao,PETSC_VIEWER_STDOUT_WORLD);
+
+        /*Save results*/
         TaoGetSolutionVector(tao,&X->pvector->vector);
         G=new Vector<IssmDouble>(0); VecFree(&G->pvector->vector);
         TaoGetGradientVector(tao,&G->pvector->vector);
         SetControlInputsFromVectorx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,X);
         ControlInputSetGradientx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,G);
-        for(int i=0;i<num_controls;i++){
-                InputToResultx(femmodel,control_list[i]);
-        }
+        femmodel->OutputControlsx(&femmodel->results);
         femmodel->results->AddObject(new GenericExternalResult<double*>(femmodel->results->Size()+1,JEnum,user.J,maxiter+3,1,1,0));
 
         /*Finalize*/