Changeset 6021
- Timestamp:
- 09/24/10 10:42:23 (15 years ago)
- Location:
- issm/trunk/src/c/objects/Numerics
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Numerics/ElementMatrix.cpp
r5989 r6021 26 26 this->ncols=0; 27 27 this->values=NULL; 28 this-> square=false;28 this->dofsymmetrical=false; 29 29 this->kff=false; 30 30 … … 68 68 69 69 /*General Case: Ke1 and Ke2 are not empty*/ 70 if(!Ke1-> square || !Ke2->square) ISSMERROR("merging 2 non squarematrices not implemented yet");70 if(!Ke1->dofsymmetrical || !Ke2->dofsymmetrical) ISSMERROR("merging 2 non dofsymmetrical matrices not implemented yet"); 71 71 72 72 /*Initialize itransformation matrix Ke[P[i]] = Ke2[i]*/ … … 90 90 this->nrows=gsize; 91 91 this->ncols=gsize; 92 this-> square=true;92 this->dofsymmetrical=true; 93 93 this->kff=Ke1->kff; 94 94 … … 195 195 196 196 /*get Matrix size and properties*/ 197 this-> square=true;197 this->dofsymmetrical=true; 198 198 this->nrows=GetNumberOfDofs(nodes,numnodes,GsetEnum,approximation); 199 199 this->ncols=this->nrows; … … 223 223 } 224 224 225 /*Because this matrix is " square" don't do cols, we can pick them up from the rows: */225 /*Because this matrix is "dofsymmetrical" don't do cols, we can pick them up from the rows: */ 226 226 this->col_fsize=0; 227 227 this->col_flocaldoflist=NULL; … … 268 268 double* localvalues=NULL; 269 269 270 if(this-> square){270 if(this->dofsymmetrical){ 271 271 /*only use row dofs to add values into global matrices: */ 272 272 … … 309 309 } 310 310 else{ 311 ISSMERROR(" non squarematrix AddToGlobal routine not support yet!");311 ISSMERROR(" non dofsymmetrical matrix AddToGlobal routine not support yet!"); 312 312 } 313 313 … … 334 334 printf(" nrows: %i\n",nrows); 335 335 printf(" ncols: %i\n",ncols); 336 printf(" square: %s\n",square?"true":"false");336 printf(" dofsymmetrical: %s\n",dofsymmetrical?"true":"false"); 337 337 printf(" kff: %s\n",kff?"true":"false"); 338 338 … … 359 359 if(row_sglobaldoflist)for(i=0;i<row_ssize;i++)printf("%i ",row_sglobaldoflist[i]); printf("\n"); 360 360 361 if(! square){361 if(!dofsymmetrical){ 362 362 printf(" col_fsize: %i\n",col_fsize); 363 363 printf(" col_flocaldoflist (%p): ",col_flocaldoflist); … … 381 381 this->nrows =Ke->nrows; 382 382 this->ncols =Ke->ncols; 383 this-> square=Ke->square;383 this->dofsymmetrical=Ke->dofsymmetrical; 384 384 this->kff =Ke->kff; 385 385 -
issm/trunk/src/c/objects/Numerics/ElementMatrix.h
r5989 r6021 23 23 int nrows; 24 24 int ncols; 25 bool square;25 bool dofsymmetrical; 26 26 bool kff; 27 27 double* values;
Note:
See TracChangeset
for help on using the changeset viewer.