Changeset 6032


Ignore:
Timestamp:
09/24/10 15:17:52 (15 years ago)
Author:
seroussi
Message:

fixed fset for tiling

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/objects/Node.cpp

    r5843 r6032  
    332332        int i;
    333333        int count=0;
     334        int count2=0;
    334335               
    335336        if(approximation_enum==NoneApproximationEnum){
    … …  
    349350                                        }
    350351                                }
    351                                 ISSMASSERT(count);
     352                                ISSMASSERT(count); //at least one dof should be the approximation requested
    352353                        }
    353354                        else for(i=0;i<this->indexing.gsize;i++) outdoflist[i]=indexing.gdoflist[i];
    … …  
    356357                        if(indexing.doftype){
    357358                                count=0;
    358                                 for(i=0;i<this->indexing.gsize;i++){
    359                                         if((indexing.doftype[i]==approximation_enum) && (indexing.f_set[i])){
    360                                                 outdoflist[count]=indexing.fdoflist[i];
    361                                                 count++;
     359                                count2=0;
     360                                for(i=0;i<this->indexing.gsize;i++){
     361                                        if(indexing.f_set[i]){
     362                                                if(indexing.doftype[i]==approximation_enum){
     363                                                        outdoflist[count]=indexing.fdoflist[count2];
     364                                                        count++;
     365                                                }
     366                                                count2++;
    362367                                        }
    363368                                }
    364                                 ISSMASSERT(count);
    365369                        }
    366370                        else for(i=0;i<this->indexing.fsize;i++) outdoflist[i]=indexing.fdoflist[i];
    … …  
    369373                        if(indexing.doftype){
    370374                                count=0;
    371                                 for(i=0;i<this->indexing.gsize;i++){
    372                                         if((indexing.doftype[i]==approximation_enum) && (indexing.s_set[i])){
    373                                                 outdoflist[count]=indexing.sdoflist[i];
    374                                                 count++;
     375                                count2=0;
     376                                for(i=0;i<this->indexing.gsize;i++){
     377                                        if(indexing.s_set[i]){
     378                                                if(indexing.doftype[i]==approximation_enum){
     379                                                        outdoflist[count]=indexing.sdoflist[count2];
     380                                                        count++;
     381                                                }
     382                                                count2++;
    375383                                        }
    376384                                }
    377                                 ISSMASSERT(count);
    378385                        }
    379386                        else for(i=0;i<this->indexing.ssize;i++) outdoflist[i]=indexing.sdoflist[i];
    … …  
    432439                                count2=0;
    433440                                for(i=0;i<this->indexing.gsize;i++){
    434                                         if((indexing.doftype[i]==approximation_enum) && (indexing.f_set[i])){
    435                                                 outdoflist[count]=count2;
    436                                                 count++;
     441                                        if(indexing.doftype[i]==approximation_enum){
     442                                                if(indexing.f_set[i]){
     443                                                        outdoflist[count]=count2;
     444                                                        count++;
     445                                                }
     446                                                count2++;
    437447                                        }
    438                                         if(indexing.doftype[i]==approximation_enum)count2++;
    439                                 }
    440                                 ISSMASSERT(count);
     448                                }
     449                                ISSMASSERT(count2);
    441450                        }
    442451                        else{
    … …  
    456465                                count2=0;
    457466                                for(i=0;i<this->indexing.gsize;i++){
    458                                         if((indexing.doftype[i]==approximation_enum) && (indexing.s_set[i])){
    459                                                 outdoflist[count]=count2;
    460                                                 count++;
     467                                        if(indexing.doftype[i]==approximation_enum){
     468                                                if(indexing.s_set[i]){
     469                                                        outdoflist[count]=count2;
     470                                                        count++;
     471                                                }
     472                                                count2++;
    461473                                        }
    462                                         if(indexing.doftype[i]==approximation_enum)count2++;
    463                                 }
    464                                 ISSMASSERT(count);
     474                                }
     475                                ISSMASSERT(count2);
    465476                        }
    466477                        else{
Note: See TracChangeset for help on using the changeset viewer.