Changeset 560


Ignore:
Timestamp:
05/21/09 10:47:14 (16 years ago)
Author:
seroussi
Message:

no penalties in hutter

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp

    r323 r560  
    8383        xfree((void**)&model->gridonhutter);
    8484
    85         //deal with mpcs for 2d-3d mesh transitions
    86        
    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 dof
    97                                         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 dof
    107                                         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 
    12485        /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these
    12586         * datasets, it will not be redone: */
Note: See TracChangeset for help on using the changeset viewer.