Changeset 3442


Ignore:
Timestamp:
04/07/10 17:04:33 (15 years ago)
Author:
Eric.Larour
Message:

debug

Location:
issm/trunk/src/c
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp

    r3439 r3442  
    1616        int i;
    1717        int count=0;
    18        
    1918        DataSet* constraints = NULL;
    20         Spc*    spc  = NULL;
    21         int     node1,node2;
    2219
    2320        /*Create constraints: */
     
    3027        /*Create spcs from x,y,z, as well as the spc values on those spcs: */
    3128        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        }
    3638
    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);
    5541
    5642        /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these
     
    5844        constraints->Presort();
    5945
    60         /*Free data: */
    61         xfree((void**)&spcthickness);
    62        
    6346        cleanup_and_return:
    6447       
  • issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp

    r3439 r3442  
    1414
    1515
     16        int i;
     17       
    1618        DataSet*    loads    = NULL;
    1719        Icefront*   icefront = NULL;
  • issm/trunk/src/c/objects/Node.cpp

    r3441 r3442  
    135135                /*On a 3d mesh, we may have collapsed elements, hence dead grids. Freeze them out: */
    136136                if (!iomodel->gridonbed[i]){
    137                         for(k=1;k<=node_numdofs;k++){
    138                                 node->FreezeDof(k);
     137                        for(k=1;k<=numdofs;k++){
     138                                this->FreezeDof(k);
    139139                        }
    140140                }
Note: See TracChangeset for help on using the changeset viewer.