Changeset 5412

Timestamp:

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

Author:

Mathieu Morlighem

Message:

Fixed misfit 3

Location:

issm/trunk/src/m/solutions

Files:

• adjointbalancedthickness_core.m (modified) (1 diff)
• adjointdiagnostic_core.m (modified) (2 diffs)

issm/trunk/src/m/solutions/adjointbalancedthickness_core.m ¶

@@ -15 +15 @@
         femmodel=solver_linear(femmodel);
 
+        %Call SurfaceAreax, because some it might be needed by PVector
+        femmodel.elements=SurfaceArea(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
+
+        %compute adjoint
         displaystring('\n%s',['      computing adjoint']);
         femmodel=SetCurrentConfiguration(femmodel,BalancedthicknessAnalysisEnum,AdjointBalancedthicknessAnalysisEnum);
         femmodel=solver_adjoint_linear(femmodel);
 
+        %Save results
+        femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,AdjointEnum);
         if(solution_type==AdjointSolutionEnum & ~control_analysis)
                 displaystring(verbose,'\n%s',['      saving results...']);

issm/trunk/src/m/solutions/adjointdiagnostic_core.m ¶

@@ -19 +19 @@
         femmodel=solver_diagnostic_nonlinear(femmodel,conserve_loads);
 
+        %Call SurfaceAreax, because some it might be needed by PVector
+        femmodel.elements=SurfaceArea(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
+
+        %compute adjoint
         displaystring('\n%s',['      computing adjoint']);
         if(isstokes), femmodel=SetCurrentConfiguration(femmodel,DiagnosticStokesAnalysisEnum,AdjointStokesAnalysisEnum);
@@ -24 +28 @@
         femmodel=solver_adjoint_linear(femmodel);
 
+        %save results
         if(solution_type==AdjointSolutionEnum & ~control_analysis)
                 displaystring(verbose,'\n%s',['      saving results...']);