Changeset 3442
- Timestamp:
- 04/07/10 17:04:33 (15 years ago)
- Location:
- issm/trunk/src/c
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp
r3439 r3442 16 16 int i; 17 17 int count=0; 18 19 18 DataSet* constraints = NULL; 20 Spc* spc = NULL;21 int node1,node2;22 19 23 20 /*Create constraints: */ … … 30 27 /*Create spcs from x,y,z, as well as the spc values on those spcs: */ 31 28 for (i=0;i<iomodel->numberofvertices;i++){ 32 #ifdef _PARALLEL_ 33 /*keep only this partition's nodes:*/ 34 if((iomodel->my_grids[i]==1)){ 35 #endif 29 if(iomodel->my_vertices[i]){ 30 31 if ((int)iomodel->spcthickness[2*i]){ 32 33 constraints->AddObject(new Spc(count,i+1,1,*(iomodel->spcthickness+2*i+1)));//we enforce first translation degree of freedom, for temperature 34 count++; 35 } 36 } 37 } 36 38 37 if ((int)spcthickness[2*i]){ 38 39 /*This grid needs to be spc'd: */ 40 41 spc_sid=count; 42 spc_node=i+1; 43 spc_dof=1; //we enforce first translation degree of freedom, for temperature 44 spc_value=*(spcthickness+2*i+1); 45 46 spc = new Spc(spc_sid,spc_node,spc_dof,spc_value); 47 constraints->AddObject(spc); 48 count++; 49 } 50 51 #ifdef _PARALLEL_ 52 } //if((my_grids[i]==1)) 53 #endif 54 } 39 /*Free data: */ 40 xfree((void**)&iomodel->spcthickness); 55 41 56 42 /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these … … 58 44 constraints->Presort(); 59 45 60 /*Free data: */61 xfree((void**)&spcthickness);62 63 46 cleanup_and_return: 64 47 -
issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp
r3439 r3442 14 14 15 15 16 int i; 17 16 18 DataSet* loads = NULL; 17 19 Icefront* icefront = NULL; -
issm/trunk/src/c/objects/Node.cpp
r3441 r3442 135 135 /*On a 3d mesh, we may have collapsed elements, hence dead grids. Freeze them out: */ 136 136 if (!iomodel->gridonbed[i]){ 137 for(k=1;k<=n ode_numdofs;k++){138 node->FreezeDof(k);137 for(k=1;k<=numdofs;k++){ 138 this->FreezeDof(k); 139 139 } 140 140 }
Note:
See TracChangeset
for help on using the changeset viewer.