Changeset 13621 for issm/trunk-jpl/src/c/solutions/control_core.cpp
- Timestamp:
- 10/11/12 11:20:39 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/solutions/control_core.cpp
r13540 r13621 16 16 17 17 int i,n; 18 18 19 19 /*parameters: */ 20 20 int num_controls,num_responses; … … 32 32 IssmDouble* maxiter=NULL; 33 33 IssmDouble* cm_jump=NULL; 34 34 35 35 /*intermediary: */ 36 36 IssmDouble search_scalar=1; … … 72 72 J=xNew<IssmDouble>(nsteps); 73 73 step_responses=xNew<int>(num_responses); 74 74 75 75 /*Initialize some of the BrentSearch arguments: */ 76 76 optargs.femmodel=femmodel; 77 77 optpars.xmin=0; optpars.xmax=1; 78 78 79 79 /*Start looping: */ 80 80 for(n=0;n<nsteps;n++){ … … 84 84 for(i=0;i<num_responses;i++) step_responses[i]=reCast<int,IssmDouble>(responses[n*num_responses+i]); 85 85 femmodel->parameters->SetParam(step_responses,1,num_responses,StepResponsesEnum); 86 86 87 87 /*In steady state inversion, compute new temperature field now*/ 88 88 if(solution_type==SteadystateSolutionEnum) solutioncore(femmodel); … … 104 104 if(VerboseControl()) _pprintLine_(" updating parameter using optimized search scalar"); //true means update save controls 105 105 InputControlUpdatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,search_scalar,true); 106 106 107 107 if(controlconvergence(J[n],tol_cm)) break; 108 108 }
Note:
See TracChangeset
for help on using the changeset viewer.