Changeset 17106


Ignore:
Timestamp:
01/13/14 19:45:49 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: use id0 instead of iomodel->nodecounter

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

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

    r17103 r17106  
    25962596
    25972597        /* Start  looping on the number of gaussian points: */
    2598         Gauss* gauss = element->NewGauss(5);
     2598        Gauss* gauss = element->NewGauss(10);
    25992599        for(int ig=gauss->begin();ig<gauss->end();ig++){
    26002600                gauss->GaussPoint(ig);
    … …  
    26802680                element->GetGroundedPart(&point1,&fraction1,&fraction2,&mainlyfloating);
    26812681                //gauss = element->NewGauss(point1,fraction1,fraction2,mainlyfloating,2);
    2682                 gauss=element->NewGaussBase(3);
     2682                gauss=element->NewGaussBase(10);
    26832683        }
    26842684        else{
    2685                 gauss=element->NewGaussBase(3);
     2685                gauss=element->NewGaussBase(10);
    26862686        }
    26872687
    … …  
    27712771
    27722772        /* Start  looping on the number of gaussian points: */
    2773         Gauss* gauss=element->NewGauss(5);
     2773        Gauss* gauss=element->NewGauss(10);
    27742774        for(int ig=gauss->begin();ig<gauss->end();ig++){
    27752775                gauss->GaussPoint(ig);
    … …  
    28442844
    28452845        /* Start  looping on the number of gaussian points: */
    2846         Gauss* gauss=element->NewGaussBase(5);
     2846        Gauss* gauss=element->NewGaussBase(10);
    28472847        for(int ig=gauss->begin();ig<gauss->end();ig++){
    28482848                gauss->GaussPoint(ig);

  • issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateNodes.cpp

    r17104 r17106  
    1313        /*Intermediaries*/
    1414        int   i,j,counter,vnodes,lid=0;
     15        int   id0 = iomodel->nodecounter;
    1516        bool *my_edges = NULL;
    1617        bool *my_nodes = NULL;
    … …  
    2122                        for(i=0;i<iomodel->numberofvertices;i++){
    2223                                if(iomodel->my_vertices[i]){
    23                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     24                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    2425                                }
    2526                        }
    … …  
    3132                                for(j=0;j<3;j++){
    3233                                        if(my_nodes[3*i+j]){
    33                                                 nodes->AddObject(new Node(iomodel->nodecounter+3*i+j+1,iomodel->nodecounter+3*i+j,lid++,iomodel->elements[+3*i+j]-1,iomodel,analysis,approximation));
     34                                                nodes->AddObject(new Node(id0+3*i+j+1,id0+3*i+j,lid++,iomodel->elements[+3*i+j]-1,iomodel,analysis,approximation));
    3435
    3536                                        }
    … …  
    4142                        for(i=0;i<iomodel->numberofvertices;i++){
    4243                                if(iomodel->my_vertices[i]){
    43                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     44                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    4445                                }
    4546                        }
    4647                        for(i=0;i<iomodel->numberofelements;i++){
    4748                                if(iomodel->my_elements[i]){
    48                                         nodes->AddObject(new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,approximation));
     49                                        nodes->AddObject(new Node(id0+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,approximation));
    4950                                }
    5051                        }
    … …  
    5455                        for(i=0;i<iomodel->numberofvertices;i++){
    5556                                if(iomodel->my_vertices[i]){
    56                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     57                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    5758                                }
    5859                        }
    5960                        for(i=0;i<iomodel->numberofelements;i++){
    6061                                if(iomodel->my_elements[i]){
    61                                         node = new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,approximation);
     62                                        node = new Node(id0+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,approximation);
    6263                                        node->Deactivate();
    6364                                        nodes->AddObject(node);
    … …  
    7071                        for(i=0;i<iomodel->numberofvertices;i++){
    7172                                if(iomodel->my_vertices[i]){
    72                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     73                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    7374                                }
    7475                        }
    … …  
    7879                                if(iomodel->edges[i*3+2]==2){
    7980                                        if(my_edges[i]){
    80                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,counter+1,lid++,0,iomodel,analysis,approximation);
     81                                                node = new Node(id0+iomodel->numberofvertices+i+1,counter+1,lid++,0,iomodel,analysis,approximation);
    8182                                                nodes->AddObject(node);
    8283                                        }
    … …  
    9091                        for(i=0;i<iomodel->numberofvertices;i++){
    9192                                if(iomodel->my_vertices[i]){
    92                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     93                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    9394                                }
    9495                        }
    … …  
    9899                                if(iomodel->edges[i*3+2]==2){
    99100                                        if(my_edges[i]){
    100                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+2*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
    101                                                 nodes->AddObject(node);
    102                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+2*i+2,counter+2,lid++,0,iomodel,analysis,approximation);
     101                                                node = new Node(id0+iomodel->numberofvertices+2*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
     102                                                nodes->AddObject(node);
     103                                                node = new Node(id0+iomodel->numberofvertices+2*i+2,counter+2,lid++,0,iomodel,analysis,approximation);
    103104                                                nodes->AddObject(node);
    104105                                        }
    … …  
    112113                        for(i=0;i<iomodel->numberofvertices;i++){
    113114                                if(iomodel->my_vertices[i]){
    114                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     115                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    115116                                }
    116117                        }
    … …  
    120121                                if(iomodel->edges[i*3+2]!=2){
    121122                                        if(my_edges[i]){
    122                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,counter+1,lid++,0,iomodel,analysis,approximation);
     123                                                node = new Node(id0+iomodel->numberofvertices+i+1,counter+1,lid++,0,iomodel,analysis,approximation);
    123124                                                nodes->AddObject(node);
    124125                                        }
    … …  
    132133                        for(i=0;i<iomodel->numberofvertices;i++){
    133134                                if(iomodel->my_vertices[i]){
    134                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     135                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    135136                                }
    136137                        }
    137138                        for(i=0;i<iomodel->numberofedges;i++){
    138139                                if(my_edges[i]){
    139                                         nodes->AddObject(new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,approximation));
     140                                        nodes->AddObject(new Node(id0+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,approximation));
    140141                                }
    141142                        }
    … …  
    145146                        for(i=0;i<iomodel->numberofvertices;i++){
    146147                                if(iomodel->my_vertices[i]){
    147                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,approximation));
     148                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,approximation));
    148149                                }
    149150                        }
    … …  
    152153                                if(iomodel->edges[i*3+2]==2){/*Vertical edges*/
    153154                                        if(my_edges[i]){
    154                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
    155                                                 nodes->AddObject(node);
    156                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,approximation);
    157                                                 nodes->AddObject(node);
    158                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,approximation);
     155                                                node = new Node(id0+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
     156                                                nodes->AddObject(node);
     157                                                node = new Node(id0+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,approximation);
     158                                                nodes->AddObject(node);
     159                                                node = new Node(id0+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,approximation);
    159160                                                nodes->AddObject(node);
    160161                                        }
    … …  
    163164                                else if(iomodel->edges[i*3+2]==1){/*Basal edges*/
    164165                                        if(my_edges[i]){
    165                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
    166                                                 nodes->AddObject(node);
    167                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,approximation);
    168                                                 nodes->AddObject(node);
    169                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,approximation);
    170                                                 nodes->AddObject(node);
    171                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+4,counter+4,lid++,0,iomodel,analysis,approximation);
     166                                                node = new Node(id0+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
     167                                                nodes->AddObject(node);
     168                                                node = new Node(id0+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,approximation);
     169                                                nodes->AddObject(node);
     170                                                node = new Node(id0+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,approximation);
     171                                                nodes->AddObject(node);
     172                                                node = new Node(id0+iomodel->numberofvertices+4*i+4,counter+4,lid++,0,iomodel,analysis,approximation);
    172173                                                nodes->AddObject(node);
    173174                                        }
    … …  
    176177                                else if(iomodel->edges[i*3+2]==3){/*Top edges*/
    177178                                        if(my_edges[i]){
    178                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
     179                                                node = new Node(id0+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,approximation);
    179180                                                nodes->AddObject(node);
    180181                                        }
    … …  
    193194                        for(i=0;i<iomodel->numberofvertices;i++){
    194195                                if(iomodel->my_vertices[i]){
    195                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
    196                                 }
    197                         }
    198                         /*P1 pressure*/
    199                         vnodes = iomodel->nodecounter+iomodel->numberofvertices;
     196                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
     197                                }
     198                        }
     199                        /*P1 pressure*/
     200                        vnodes = id0+iomodel->numberofvertices;
    200201                        for(i=0;i<iomodel->numberofvertices;i++){
    201202                                if(iomodel->my_vertices[i]){
    … …  
    209210                        for(i=0;i<iomodel->numberofvertices;i++){
    210211                                if(iomodel->my_vertices[i]){
    211                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
    212                                 }
    213                         }
    214                         /*P1 pressure*/
    215                         vnodes = iomodel->nodecounter+iomodel->numberofvertices;
     212                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
     213                                }
     214                        }
     215                        /*P1 pressure*/
     216                        vnodes = id0+iomodel->numberofvertices;
    216217                        for(i=0;i<iomodel->numberofvertices;i++){
    217218                                if(iomodel->my_vertices[i]){
    … …  
    225226                        for(i=0;i<iomodel->numberofvertices;i++){
    226227                                if(iomodel->my_vertices[i]){
    227                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
     228                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
    228229                                }
    229230                        }
    230231                        for(i=0;i<iomodel->numberofelements;i++){
    231232                                if(iomodel->my_elements[i]){
    232                                         node = new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,FSvelocityEnum);
     233                                        node = new Node(id0+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,FSvelocityEnum);
    233234                                        node->Deactivate();
    234235                                        nodes->AddObject(node);
    … …  
    236237                        }
    237238                        /*P1 pressure*/
    238                         vnodes = iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofelements;
     239                        vnodes = id0+iomodel->numberofvertices+iomodel->numberofelements;
    239240                        for(i=0;i<iomodel->numberofvertices;i++){
    240241                                if(iomodel->my_vertices[i]){
    … …  
    248249                        for(i=0;i<iomodel->numberofvertices;i++){
    249250                                if(iomodel->my_vertices[i]){
    250                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
     251                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
    251252                                }
    252253                        }
    253254                        for(i=0;i<iomodel->numberofelements;i++){
    254255                                if(iomodel->my_elements[i]){
    255                                         nodes->AddObject(new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,FSvelocityEnum));
    256                                 }
    257                         }
    258                         /*P1 pressure*/
    259                         vnodes = iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofelements;
     256                                        nodes->AddObject(new Node(id0+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,FSvelocityEnum));
     257                                }
     258                        }
     259                        /*P1 pressure*/
     260                        vnodes = id0+iomodel->numberofvertices+iomodel->numberofelements;
    260261                        for(i=0;i<iomodel->numberofvertices;i++){
    261262                                if(iomodel->my_vertices[i]){
    … …  
    270271                        for(i=0;i<iomodel->numberofvertices;i++){
    271272                                if(iomodel->my_vertices[i]){
    272                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
     273                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
    273274                                }
    274275                        }
    275276                        for(i=0;i<iomodel->numberofedges;i++){
    276277                                if(my_edges[i]){
    277                                         nodes->AddObject(new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,FSvelocityEnum));
    278                                 }
    279                         }
    280 
    281                         /*P1 pressure*/
    282                         vnodes = iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofedges;
     278                                        nodes->AddObject(new Node(id0+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,lid++,0,iomodel,analysis,FSvelocityEnum));
     279                                }
     280                        }
     281
     282                        /*P1 pressure*/
     283                        vnodes = id0+iomodel->numberofvertices+iomodel->numberofedges;
    283284                        for(i=0;i<iomodel->numberofvertices;i++){
    284285                                if(iomodel->my_vertices[i]){
    … …  
    293294                        for(i=0;i<iomodel->numberofvertices;i++){
    294295                                if(iomodel->my_vertices[i]){
    295                                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
     296                                        nodes->AddObject(new Node(id0+i+1,i,lid++,i,iomodel,analysis,FSvelocityEnum));
    296297                                }
    297298                        }
    … …  
    300301                                if(iomodel->edges[i*3+2]==2){/*Vertical edges*/
    301302                                        if(my_edges[i]){
    302                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,FSvelocityEnum);
    303                                                 nodes->AddObject(node);
    304                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,FSvelocityEnum);
    305                                                 nodes->AddObject(node);
    306                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,FSvelocityEnum);
     303                                                node = new Node(id0+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,FSvelocityEnum);
     304                                                nodes->AddObject(node);
     305                                                node = new Node(id0+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,FSvelocityEnum);
     306                                                nodes->AddObject(node);
     307                                                node = new Node(id0+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,FSvelocityEnum);
    307308                                                nodes->AddObject(node);
    308309                                        }
    … …  
    311312                                else if(iomodel->edges[i*3+2]==1){/*Basal edges*/
    312313                                        if(my_edges[i]){
    313                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,FSvelocityEnum);
    314                                                 nodes->AddObject(node);
    315                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,FSvelocityEnum);
    316                                                 nodes->AddObject(node);
    317                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,FSvelocityEnum);
    318                                                 nodes->AddObject(node);
    319                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+4,counter+4,lid++,0,iomodel,analysis,FSvelocityEnum);
     314                                                node = new Node(id0+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,FSvelocityEnum);
     315                                                nodes->AddObject(node);
     316                                                node = new Node(id0+iomodel->numberofvertices+4*i+2,counter+2,lid++,0,iomodel,analysis,FSvelocityEnum);
     317                                                nodes->AddObject(node);
     318                                                node = new Node(id0+iomodel->numberofvertices+4*i+3,counter+3,lid++,0,iomodel,analysis,FSvelocityEnum);
     319                                                nodes->AddObject(node);
     320                                                node = new Node(id0+iomodel->numberofvertices+4*i+4,counter+4,lid++,0,iomodel,analysis,FSvelocityEnum);
    320321                                                nodes->AddObject(node);
    321322                                        }
    … …  
    324325                                else if(iomodel->edges[i*3+2]==3){/*Top edges*/
    325326                                        if(my_edges[i]){
    326                                                 node = new Node(iomodel->nodecounter+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,FSvelocityEnum);
     327                                                node = new Node(id0+iomodel->numberofvertices+4*i+1,counter+1,lid++,0,iomodel,analysis,FSvelocityEnum);
    327328                                                nodes->AddObject(node);
    328329                                        }
    … …  
    334335                        }
    335336                        /*P1xP3 pressure*/
    336                         vnodes = iomodel->nodecounter+iomodel->numberofvertices+4*iomodel->numberofedges;
     337                        vnodes = id0+iomodel->numberofvertices+4*iomodel->numberofedges;
    337338                        for(i=0;i<iomodel->numberofvertices;i++){
    338339                                if(iomodel->my_vertices[i]){

  • issm/trunk-jpl/src/c/shared/Exceptions/Exceptions.cpp

    r16220 r17106  
    4545
    4646        file_line= what_line;
     47        this->Report();
    4748
    4849}/*}}}*/
Note: See TracChangeset for help on using the changeset viewer.