Changeset 8602
- Timestamp:
- 06/10/11 13:44:57 (14 years ago)
- Location:
- issm/trunk/src/m
- Files:
-
- 3 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/control_core.m
r6323 r8602 7 7 8 8 %recover parameters common to all solutions 9 nsteps=femmodel.parameters.NSteps;10 9 num_controls=femmodel.parameters.NumControls; 10 num_responses=femmodel.parameters.NumResponses; 11 responses=femmodel.parameters.CmResponses; 11 12 control_type=femmodel.parameters.ControlType; 12 13 solution_type=femmodel.parameters.SolutionType; 13 responses=femmodel.parameters.CmResponses;14 nsteps=femmodel.parameters.NSteps; 14 15 maxiter=femmodel.parameters.MaxIter; 15 16 cm_jump=femmodel.parameters.CmJump; … … 41 42 42 43 issmprintf(VerboseControl(),['control method step ' num2str(n) '/' num2str(femmodel.parameters.NSteps)]); 43 [femmodel.elements,femmodel.loads]=InputUpdateFromConstant(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,responses(n),CmResponseEnum);44 femmodel.parameters.StepResponses=responses(n,:); 44 45 45 46 %In case we are running a steady state control method, compute new temperature field using new parameter distribution: -
issm/trunk/src/m/solutions/objectivefunctionC.m
r8397 r8602 1 function J =objectivefunctionC(search_scalar,femmodel ,n);1 function J =objectivefunctionC(search_scalar,femmodel); 2 2 %OBJECTIVEFUNCTIONC - objective function that return a parameter for a certain function 3 3 4 4 conserve_loads=true; 5 J=0; 6 5 7 %recover some parameters 6 response=femmodel.parameters.CmResponses(n); 7 analysis_type=femmodel.parameters.AnalysisType; 8 solution_type=femmodel.parameters.SolutionType; 9 isstokes=femmodel.parameters.IsStokes; 8 num_responses = femmodel.parameters.NumResponses; 9 responses = femmodel.parameters.StepResponses; 10 analysis_type = femmodel.parameters.AnalysisType; 11 solution_type = femmodel.parameters.SolutionType; 12 isstokes = femmodel.parameters.IsStokes; 10 13 11 14 %set current configuration … … 33 36 34 37 %Compute misfit for this velocity field 35 J=CostFunction(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials, femmodel.parameters,response); 38 for i=1:num_responses 39 J=J+CostFunction(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials, femmodel.parameters,responses(i)); 40 end
Note:
See TracChangeset
for help on using the changeset viewer.