Changeset 560
- Timestamp:
- 05/21/09 10:47:14 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp
r323 r560 83 83 xfree((void**)&model->gridonhutter); 84 84 85 //deal with mpcs for 2d-3d mesh transitions86 87 /*Fetch data: */88 ModelFetchData((void**)&model->penalties,&model->numpenalties,NULL,model_handle,"penalties","Matrix","Mat");89 90 if (strcmp(model->meshtype,"3d")==0){91 92 if (model->numpenalties){93 for (i=0;i<model->numpenalties;i++){94 for (j=0;j<model->numlayers-1;j++){95 96 //constrain first dof97 rgb_id=count;98 rgb_dof=1;99 rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j);100 rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j+1);101 102 rgb = new Rgb(rgb_id,rgb_nodeid1,rgb_nodeid2,rgb_dof);103 constraints->AddObject(rgb);104 count++;105 106 //constrain second dof107 rgb_id=count;108 rgb_dof=2;109 rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j);110 rgb_nodeid1=(int)*(model->penalties+(model->numlayers-1)*i+j+1);111 112 rgb = new Rgb(rgb_id,rgb_nodeid1,rgb_nodeid2,rgb_dof);113 constraints->AddObject(rgb);114 count++;115 116 }117 }118 }119 }120 121 /*Free data: */122 xfree((void**)&model->penalties);123 124 85 /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these 125 86 * datasets, it will not be redone: */
Note:
See TracChangeset
for help on using the changeset viewer.