Changeset 1212
- Timestamp:
- 07/02/09 10:41:07 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/parallel/control_core.cpp
r1209 r1212 100 100 101 101 /*take care of Stokes: compute slope and get spc once for all*/ 102 if (isstokes){ 103 ControlPrepareStokes(fems,inputs); 104 } 102 if (isstokes) ControlPrepareStokes(fems,inputs); 105 103 106 104 /*Start looping: */ … … 112 110 113 111 /*Update parameters: */ 114 if (isstokes){ 115 UpdateFromInputsx(fem_ds->elements,fem_ds->nodes,fem_ds->loads, fem_ds->materials,inputs); 116 } 117 else{ 118 UpdateFromInputsx(fem_dh->elements,fem_dh->nodes,fem_dh->loads, fem_dh->materials,inputs); 119 } 112 if (isstokes) UpdateFromInputsx(fem_ds->elements,fem_ds->nodes,fem_ds->loads, fem_ds->materials,inputs); 113 else UpdateFromInputsx(fem_dh->elements,fem_dh->nodes,fem_dh->loads, fem_dh->materials,inputs); 120 114 121 115 _printf_("%s\n"," computing gradJ..."); 122 if (isstokes){ 123 grad_g=GradJCompute(inputs,fem_ds); 124 } 125 else{ 126 grad_g=GradJCompute(inputs,fem_dh); 127 } 116 if (isstokes) grad_g=GradJCompute(inputs,fem_ds); 117 else grad_g=GradJCompute(inputs,fem_dh); 128 118 _printf_("%s\n"," done."); 129 119 … … 136 126 137 127 _printf_("%s\n"," optimizing along gradient direction..."); 138 if (isstokes){ 139 optargs.femmodel=fem_ds; 140 } 141 else{ 142 optargs.femmodel=fem_dh; 143 } 128 if (isstokes) optargs.femmodel=fem_ds; 129 else optargs.femmodel=fem_dh; 144 130 optargs.param_g=param_g; optargs.grad_g=grad_g_double; optargs.inputs=inputs;optargs.n=n; 145 131 optpars.xmin=0; optpars.xmax=1; optpars.tolerance=tolx; optpars.maxiter=(int)maxiter[n];
Note:
See TracChangeset
for help on using the changeset viewer.