Changeset 11289

Timestamp:

02/01/12 13:34:04 (13 years ago)

Author:

Mathieu Morlighem

Message:

Better way to set TAO parameters

File:

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

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

@@ -35 +35 @@
         if(ierr) _error_("Could not initialize Tao");
 
-        /*Line search options*/
-        /*TAO options: http://www.mcs.anl.gov/research/projects/tao/docs/manualpages/solver/TaoSetFromOptions.html*/
         info = PetscOptionsSetValue("-tao_monitor",""); if(info) _error_("STOP");
-        info = PetscOptionsSetValue("-tao_max_its","20"); if(info) _error_("STOP");
-        info = PetscOptionsSetValue("-tao_max_funcs","50"); if(info) _error_("STOP");
 
         /*Initialize argument*/
@@ -48 +44 @@
         info = TaoSetType(tao,"tao_blmvm"); if(info) _error_("STOP");
         //info = TaoSetType(tao,"tao_cg"); if(info) _error_("STOP");
+        TaoSetMaximumFunctionEvaluations(tao,10);
+        TaoSetMaximumIterations(tao,3);
 
         int size=femmodel->vertices->NumberOfVertices();
@@ -97 +95 @@
         VecFree(&gradient);
         CostFunctionx(fcn, femmodel->elements,femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials, femmodel->parameters);
+        _printf_(true,"f(x)=%g\n",*fcn);
         return 0;
 }