Changeset 2019
- Timestamp:
- 08/27/09 14:10:01 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/parallel/ControlInitialization.cpp
r2016 r2019 102 102 VecDuplicatePatch(&ug_horiz,ug); FieldExtrudex( ug_horiz,fem_dh->elements,fem_dh->nodes, fem_dh->loads,fem_dh-> materials,"velocity",1); 103 103 104 //pressure105 if(debug)_printf_("%s\n"," computing pressure according to MacAyeal...");106 ComputePressurex( &pg,fem_dh->elements, fem_dh->nodes, fem_dh->loads, fem_dh->materials, numberofnodes);107 108 104 //vertical velocity 109 105 if(debug)_printf_("%s\n"," computing vertical velocities..."); … … 141 137 VecFree(&ug); 142 138 diagnostic_core_nonlinear(&ug,NULL,NULL,NULL,fem_ds,inputs,DiagnosticAnalysisEnum(),StokesAnalysisEnum()); 143 VecFree(&pg);144 139 145 140 //Add in inputs 146 141 inputs->Add("velocity",ug,numberofdofspernode_ds,numberofnodes); 147 VecFree(&ug);148 142 149 143 /*Assign output*/ … … 152 146 /*Free ressources:*/ 153 147 xfree((void**)&dofset); 148 VecFree(&ug); 149 VecFree(&ug_horiz); 150 VecFree(&ug_vert); 151 VecFree(&ug_stokes); 152 VecFree(&pg); 153 VecFree(&slopex); 154 VecFree(&slopey); 154 155 155 156 }
Note:
See TracChangeset
for help on using the changeset viewer.