Ignore:
Timestamp:
09/15/10 13:09:01 (15 years ago)
Author:
Mathieu Morlighem
Message:

some renaming

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/objects/Numerics/ElementMatrix.cpp

    r5821 r5827  
    3030
    3131        this->row_fsize=0;
    32         this->row_finternaldoflist=NULL;
    33         this->row_fexternaldoflist=NULL;
     32        this->row_flocaldoflist=NULL;
     33        this->row_fglobaldoflist=NULL;
    3434        this->row_ssize=0;
    35         this->row_sinternaldoflist=NULL;
    36         this->row_sexternaldoflist=NULL;
     35        this->row_slocaldoflist=NULL;
     36        this->row_sglobaldoflist=NULL;
    3737
    3838        this->col_fsize=0;
    39         this->col_finternaldoflist=NULL;
    40         this->col_fexternaldoflist=NULL;
     39        this->col_flocaldoflist=NULL;
     40        this->col_fglobaldoflist=NULL;
    4141        this->col_ssize=0;
    42         this->col_sinternaldoflist=NULL;
    43         this->col_sexternaldoflist=NULL;
    44 
    45 }
    46 /*}}}*/
    47 /*FUNCTION ElementMatrix::ElementMatrix(int gsize,bool square,int* in_gexternaldoflist){{{1*/
    48 ElementMatrix::ElementMatrix(int gsize,bool square,int* in_gexternaldoflist){
     42        this->col_slocaldoflist=NULL;
     43        this->col_sglobaldoflist=NULL;
     44
     45}
     46/*}}}*/
     47/*FUNCTION ElementMatrix::ElementMatrix(int gsize,bool square,int* in_gglobaldoflist){{{1*/
     48ElementMatrix::ElementMatrix(int gsize,bool square,int* in_gglobaldoflist){
    4949
    5050        if(square=false)ISSMERROR(" calling square constructor with false square flag!");
    … …  
    5959       
    6060        if(this->nrows){
    61                 this->gexternaldoflist=(int*)xmalloc(this->nrows*sizeof(int));
    62                 memcpy(this->gexternaldoflist,in_gexternaldoflist,this->nrows*sizeof(int));
    63         }
    64         else{
    65                 this->gexternaldoflist=NULL;
     61                this->gglobaldoflist=(int*)xmalloc(this->nrows*sizeof(int));
     62                memcpy(this->gglobaldoflist,in_gglobaldoflist,this->nrows*sizeof(int));
     63        }
     64        else{
     65                this->gglobaldoflist=NULL;
    6666        }
    6767
    6868        /*don't do rows and cols, because kff is false:*/
    6969        this->row_fsize=0;
    70         this->row_finternaldoflist=NULL;
    71         this->row_fexternaldoflist=NULL;
     70        this->row_flocaldoflist=NULL;
     71        this->row_fglobaldoflist=NULL;
    7272        this->row_ssize=0;
    73         this->row_sinternaldoflist=NULL;
    74         this->row_sexternaldoflist=NULL;
     73        this->row_slocaldoflist=NULL;
     74        this->row_sglobaldoflist=NULL;
    7575       
    7676        this->col_fsize=0;
    77         this->col_finternaldoflist=NULL;
    78         this->col_fexternaldoflist=NULL;
     77        this->col_flocaldoflist=NULL;
     78        this->col_fglobaldoflist=NULL;
    7979        this->col_ssize=0;
    80         this->col_sinternaldoflist=NULL;
    81         this->col_sexternaldoflist=NULL;
    82 
    83 }
    84 /*}}}*/
    85 /*FUNCTION ElementMatrix::ElementMatrix(int gsize,bool square,int* finternaldoflist,int* fexternaldoflist,int fsize,int* sinternaldoflist,int* sexternaldoflist,int ssize){{{1*/
    86 ElementMatrix::ElementMatrix(int gsize,bool square,int* finternaldoflist,int* fexternaldoflist,int fsize,int* sinternaldoflist,int* sexternaldoflist,int ssize){
     80        this->col_slocaldoflist=NULL;
     81        this->col_sglobaldoflist=NULL;
     82
     83}
     84/*}}}*/
     85/*FUNCTION ElementMatrix::ElementMatrix(int gsize,bool square,int* flocaldoflist,int* fglobaldoflist,int fsize,int* slocaldoflist,int* sglobaldoflist,int ssize){{{1*/
     86ElementMatrix::ElementMatrix(int gsize,bool square,int* flocaldoflist,int* fglobaldoflist,int fsize,int* slocaldoflist,int* sglobaldoflist,int ssize){
    8787
    8888        if(square=false)ISSMERROR(" calling square constructor with false square flag!");
    … …  
    9999        this->row_fsize=fsize;
    100100        if(fsize){
    101                 this->row_finternaldoflist=(int*)xmalloc(fsize*sizeof(int));
    102                 this->row_fexternaldoflist=(int*)xmalloc(fsize*sizeof(int));
    103                 memcpy(this->row_finternaldoflist,finternaldoflist,fsize*sizeof(int));
    104                 memcpy(this->row_fexternaldoflist,fexternaldoflist,fsize*sizeof(int));
    105         }
    106         else{
    107                 this->row_finternaldoflist=NULL;
    108                 this->row_fexternaldoflist=NULL;
     101                this->row_flocaldoflist=(int*)xmalloc(fsize*sizeof(int));
     102                this->row_fglobaldoflist=(int*)xmalloc(fsize*sizeof(int));
     103                memcpy(this->row_flocaldoflist,flocaldoflist,fsize*sizeof(int));
     104                memcpy(this->row_fglobaldoflist,fglobaldoflist,fsize*sizeof(int));
     105        }
     106        else{
     107                this->row_flocaldoflist=NULL;
     108                this->row_fglobaldoflist=NULL;
    109109        }
    110110
    111111        this->row_ssize=ssize;
    112112        if(ssize){
    113                 this->row_sinternaldoflist=(int*)xmalloc(ssize*sizeof(int));
    114                 this->row_sexternaldoflist=(int*)xmalloc(ssize*sizeof(int));
    115                 memcpy(this->row_sinternaldoflist,sinternaldoflist,ssize*sizeof(int));
    116                 memcpy(this->row_sexternaldoflist,sexternaldoflist,ssize*sizeof(int));
    117         }
    118         else{
    119                 this->row_sinternaldoflist=NULL;
    120                 this->row_sexternaldoflist=NULL;
     113                this->row_slocaldoflist=(int*)xmalloc(ssize*sizeof(int));
     114                this->row_sglobaldoflist=(int*)xmalloc(ssize*sizeof(int));
     115                memcpy(this->row_slocaldoflist,slocaldoflist,ssize*sizeof(int));
     116                memcpy(this->row_sglobaldoflist,sglobaldoflist,ssize*sizeof(int));
     117        }
     118        else{
     119                this->row_slocaldoflist=NULL;
     120                this->row_sglobaldoflist=NULL;
    121121        }
    122122
    123123        /*don't do cols, we can't pick them up from the rows: */
    124124        this->col_fsize=0;
    125         this->col_finternaldoflist=NULL;
    126         this->col_fexternaldoflist=NULL;
     125        this->col_flocaldoflist=NULL;
     126        this->col_fglobaldoflist=NULL;
    127127        this->col_ssize=0;
    128         this->col_sinternaldoflist=NULL;
    129         this->col_sexternaldoflist=NULL;
     128        this->col_slocaldoflist=NULL;
     129        this->col_sglobaldoflist=NULL;
    130130
    131131        /*don't do g-set: */
    132         this->gexternaldoflist=NULL;
     132        this->gglobaldoflist=NULL;
    133133
    134134}
    … …  
    138138       
    139139        xfree((void**)&this->values);
    140         xfree((void**)&this->gexternaldoflist);
    141         xfree((void**)&this->row_finternaldoflist);
    142         xfree((void**)&this->row_fexternaldoflist);
    143         xfree((void**)&this->row_sinternaldoflist);
    144         xfree((void**)&this->row_sexternaldoflist);
    145         xfree((void**)&this->col_finternaldoflist);
    146         xfree((void**)&this->col_fexternaldoflist);
    147         xfree((void**)&this->col_sinternaldoflist);
    148         xfree((void**)&this->col_sexternaldoflist);
     140        xfree((void**)&this->gglobaldoflist);
     141        xfree((void**)&this->row_flocaldoflist);
     142        xfree((void**)&this->row_fglobaldoflist);
     143        xfree((void**)&this->row_slocaldoflist);
     144        xfree((void**)&this->row_sglobaldoflist);
     145        xfree((void**)&this->col_flocaldoflist);
     146        xfree((void**)&this->col_fglobaldoflist);
     147        xfree((void**)&this->col_slocaldoflist);
     148        xfree((void**)&this->col_sglobaldoflist);
    149149}
    150150/*}}}*/
    … …  
    167167
    168168        int i,j;
    169         double* internalvalues=NULL;
     169        double* localvalues=NULL;
    170170
    171171        if(this->square){
    … …  
    173173               
    174174                if(!this->kff){
    175                         /*add internal values into global  matrix, using the fexternaldoflist: */
    176                         MatSetValues(Kgg,this->nrows,this->gexternaldoflist,this->nrows,this->gexternaldoflist,(const double*)values,ADD_VALUES);
     175                        /*add local values into global  matrix, using the fglobaldoflist: */
     176                        MatSetValues(Kgg,this->nrows,this->gglobaldoflist,this->nrows,this->gglobaldoflist,(const double*)values,ADD_VALUES);
    177177                }
    178178                else{
    179179                        if(this->row_fsize){
    180180                                /*first, retrieve values that are in the f-set from the g-set values matrix: */
    181                                 internalvalues=(double*)xmalloc(this->row_fsize*this->row_fsize*sizeof(double));
     181                                localvalues=(double*)xmalloc(this->row_fsize*this->row_fsize*sizeof(double));
    182182                                for(i=0;i<this->row_fsize;i++){
    183183                                        for(j=0;j<this->row_fsize;j++){
    184                                                 *(internalvalues+this->row_fsize*i+j)=*(this->values+this->ncols*this->row_finternaldoflist[i]+this->row_finternaldoflist[j]);
     184                                                *(localvalues+this->row_fsize*i+j)=*(this->values+this->ncols*this->row_flocaldoflist[i]+this->row_flocaldoflist[j]);
    185185                                        }
    186186                                }
    187                                 /*add internal values into global  matrix, using the fexternaldoflist: */
    188                                 MatSetValues(Kff,this->row_fsize,this->row_fexternaldoflist,this->row_fsize,this->row_fexternaldoflist,(const double*)internalvalues,ADD_VALUES);
     187                                /*add local values into global  matrix, using the fglobaldoflist: */
     188                                MatSetValues(Kff,this->row_fsize,this->row_fglobaldoflist,this->row_fsize,this->row_fglobaldoflist,(const double*)localvalues,ADD_VALUES);
    189189
    190190                                /*Free ressources:*/
    191                                 xfree((void**)&internalvalues);
     191                                xfree((void**)&localvalues);
    192192                        }
    193193
    … …  
    195195                        if((this->row_ssize!=0) && (this->row_fsize!=0)){
    196196                                /*first, retrieve values that are in the f and s-set from the g-set values matrix: */
    197                                 internalvalues=(double*)xmalloc(this->row_fsize*this->row_ssize*sizeof(double));
     197                                localvalues=(double*)xmalloc(this->row_fsize*this->row_ssize*sizeof(double));
    198198                                for(i=0;i<this->row_fsize;i++){
    199199                                        for(j=0;j<this->row_ssize;j++){
    200                                                 *(internalvalues+this->row_ssize*i+j)=*(this->values+this->ncols*this->row_finternaldoflist[i]+this->row_sinternaldoflist[j]);
     200                                                *(localvalues+this->row_ssize*i+j)=*(this->values+this->ncols*this->row_flocaldoflist[i]+this->row_slocaldoflist[j]);
    201201                                        }
    202202                                }
    203                                 /*add internal values into global  matrix, using the fexternaldoflist: */
    204                                 MatSetValues(Kfs,this->row_fsize,this->row_fexternaldoflist,this->row_ssize,this->row_sexternaldoflist,(const double*)internalvalues,ADD_VALUES);
     203                                /*add local values into global  matrix, using the fglobaldoflist: */
     204                                MatSetValues(Kfs,this->row_fsize,this->row_fglobaldoflist,this->row_ssize,this->row_sglobaldoflist,(const double*)localvalues,ADD_VALUES);
    205205
    206206                                /*Free ressources:*/
    207                                 xfree((void**)&internalvalues);
     207                                xfree((void**)&localvalues);
    208208                        }
    209209                }
    … …  
    232232        }
    233233
    234         printf("   gexternaldoflist (%p): ",gexternaldoflist);
    235     if(gexternaldoflist) for(i=0;i<nrows;i++)printf("%i ",gexternaldoflist[i]); printf("\n");
     234        printf("   gglobaldoflist (%p): ",gglobaldoflist);
     235    if(gglobaldoflist) for(i=0;i<nrows;i++)printf("%i ",gglobaldoflist[i]); printf("\n");
    236236
    237237        printf("   row_fsize: %i\n",row_fsize);
    238         printf("   row_finternaldoflist (%p): ",row_finternaldoflist);
    239     if(row_finternaldoflist) for(i=0;i<row_fsize;i++)printf("%i ",row_finternaldoflist[i]); printf("\n");
    240         printf("   row_fexternaldoflist (%p): ",row_fexternaldoflist);
    241         if(row_fexternaldoflist)for(i=0;i<row_fsize;i++)printf("%i ",row_fexternaldoflist[i]); printf("\n");
     238        printf("   row_flocaldoflist (%p): ",row_flocaldoflist);
     239    if(row_flocaldoflist) for(i=0;i<row_fsize;i++)printf("%i ",row_flocaldoflist[i]); printf("\n");
     240        printf("   row_fglobaldoflist (%p): ",row_fglobaldoflist);
     241        if(row_fglobaldoflist)for(i=0;i<row_fsize;i++)printf("%i ",row_fglobaldoflist[i]); printf("\n");
    242242
    243243        printf("   row_ssize: %i\n",row_ssize);
    244         printf("   row_sinternaldoflist (%p): ",row_sinternaldoflist);
    245     if(row_sinternaldoflist)for(i=0;i<row_ssize;i++)printf("%i ",row_sinternaldoflist[i]); printf("\n");
    246         printf("   row_sexternaldoflist (%p): ",row_sexternaldoflist);
    247         if(row_sexternaldoflist)for(i=0;i<row_ssize;i++)printf("%i ",row_sexternaldoflist[i]); printf("\n");
     244        printf("   row_slocaldoflist (%p): ",row_slocaldoflist);
     245    if(row_slocaldoflist)for(i=0;i<row_ssize;i++)printf("%i ",row_slocaldoflist[i]); printf("\n");
     246        printf("   row_sglobaldoflist (%p): ",row_sglobaldoflist);
     247        if(row_sglobaldoflist)for(i=0;i<row_ssize;i++)printf("%i ",row_sglobaldoflist[i]); printf("\n");
    248248
    249249        if(!square){
    250250                printf("   col_fsize: %i\n",col_fsize);
    251                 printf("   col_finternaldoflist (%p): ",col_finternaldoflist);
    252                 if(col_finternaldoflist)for(i=0;i<col_fsize;i++)printf("%i ",col_finternaldoflist[i]); printf("\n");
    253                 printf("   col_fexternaldoflist (%p): ",col_fexternaldoflist);
    254                 if(col_fexternaldoflist)for(i=0;i<col_fsize;i++)printf("%i ",col_fexternaldoflist[i]); printf("\n");
     251                printf("   col_flocaldoflist (%p): ",col_flocaldoflist);
     252                if(col_flocaldoflist)for(i=0;i<col_fsize;i++)printf("%i ",col_flocaldoflist[i]); printf("\n");
     253                printf("   col_fglobaldoflist (%p): ",col_fglobaldoflist);
     254                if(col_fglobaldoflist)for(i=0;i<col_fsize;i++)printf("%i ",col_fglobaldoflist[i]); printf("\n");
    255255
    256256                printf("   col_ssize: %i\n",col_ssize);
    257                 printf("   col_sinternaldoflist (%p): ",col_sinternaldoflist);
    258                 if(col_sinternaldoflist)for(i=0;i<col_ssize;i++)printf("%i ",col_sinternaldoflist[i]); printf("\n");
    259                 printf("   col_sexternaldoflist (%p): ",col_sexternaldoflist);
    260                 if(col_sexternaldoflist)for(i=0;i<col_ssize;i++)printf("%i ",col_sexternaldoflist[i]); printf("\n");
    261         }
    262 }
    263 /*}}}*/
     257                printf("   col_slocaldoflist (%p): ",col_slocaldoflist);
     258                if(col_slocaldoflist)for(i=0;i<col_ssize;i++)printf("%i ",col_slocaldoflist[i]); printf("\n");
     259                printf("   col_sglobaldoflist (%p): ",col_sglobaldoflist);
     260                if(col_sglobaldoflist)for(i=0;i<col_ssize;i++)printf("%i ",col_sglobaldoflist[i]); printf("\n");
     261        }
     262}
     263/*}}}*/
Note: See TracChangeset for help on using the changeset viewer.