Changeset 5414 for issm/trunk/src/c/solutions/adjointbalancedthickness_core.cpp

Timestamp:

08/19/10 11:22:37 (15 years ago)

Author:

Mathieu Morlighem

Message:

Fixed misfit 3

File:

• issm/trunk/src/c/solutions/adjointbalancedthickness_core.cpp (modified) (1 diff)

issm/trunk/src/c/solutions/adjointbalancedthickness_core.cpp

@@ -24 +24 @@
         femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
-        /*set analysis type to compute velocity: */
+        /*compute thickness */
         _printf_("%s\n","      computing thickness");
         femmodel->SetCurrentConfiguration(BalancedthicknessAnalysisEnum);
         solver_linear(femmodel);
 
-        /*Update inputs using adjoint solution, and same type of setup as diagnostic solution: */
+        /*Call SurfaceAreax, because some it might be needed by PVector*/
+        SurfaceAreax(NULL,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
+
+        /*compute adjoint*/
         _printf_("%s\n","      computing adjoint");
         femmodel->SetCurrentConfiguration(BalancedthicknessAnalysisEnum,AdjointBalancedthicknessAnalysisEnum);
         solver_adjoint_linear(femmodel);
         
+        /*Save results*/
         if(solution_type==AdjointSolutionEnum && !control_analysis){
                 if(verbose)_printf_("saving results:\n");