source: issm/oecreview/Archive/26740-27031/ISSM-26942-26943.diff

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

@@ -94 +94 @@
         GetVectorFromControlInputsx(&X, femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,"value");
         GetVectorFromControlInputsx(&XL,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,"lowerbound");
         GetVectorFromControlInputsx(&XU,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,"upperbound");
+        #if PETSC_VERSION_LT(3,17,0)
         TaoSetInitialVector(tao,X->pvector->vector);
+        #else
+        //TaoSetSolution(tao,X->pvector->vector);
+        _error_("not implemented yet");
+        #endif
         TaoSetVariableBounds(tao,XL->pvector->vector,XU->pvector->vector);
         delete XL;
         delete XU;
@@ -109 +114 @@
         TaoView(tao,PETSC_VIEWER_STDOUT_WORLD);
 
         /*Save results*/
+        #if PETSC_VERSION_LT(3,17,0)
         TaoGetSolutionVector(tao,&X->pvector->vector);
+        #else
+        TaoGetSolution(tao,&X->pvector->vector);
+        #endif
         G=new Vector<IssmDouble>(0); VecFree(&G->pvector->vector);
+        #if PETSC_VERSION_LT(3,17,0)
         TaoGetGradientVector(tao,&G->pvector->vector);
+        #else
+        //TaoGetGradient(tao,&G->pvector->vector);
+        _error_("not implemented yet");
+        #endif
         SetControlInputsFromVectorx(femmodel,X);
         ControlInputSetGradientx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,G);
         femmodel->OutputControlsx(&femmodel->results);