- Timestamp:
- 09/07/20 03:12:14 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/SystemMatricesx/SystemMatricesx.cpp
r24679 r25539 41 41 42 42 /*Get complete stiffness matrix without penalties*/ 43 for (i=0;i<femmodel->elements->Size();i++){44 element =xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i));43 for(Object* & object : femmodel->elements->objects){ 44 element = xDynamicCast<Element*>(object); 45 45 if(!element->AnyFSet() && analysisenum!=StressbalanceAnalysisEnum) continue; 46 46 ElementMatrix* Ke = analysis->CreateKMatrix(element); … … 52 52 } 53 53 54 for (i=0;i<femmodel->loads->Size();i++){55 load =xDynamicCast<Load*>(femmodel->loads->GetObjectByOffset(i));54 for(Object* & object : femmodel->loads->objects){ 55 load = xDynamicCast<Load*>(object); 56 56 load->CreateKMatrix(Kff_temp,NULL); 57 57 } … … 83 83 84 84 /*Fill stiffness matrix and load vector from elements*/ 85 for (i=0;i<femmodel->elements->Size();i++){86 element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i));85 for(Object* & object : femmodel->elements->objects){ 86 element = xDynamicCast<Element*>(object); 87 87 if(!element->AnyFSet() && analysisenum!=StressbalanceAnalysisEnum) continue; 88 88 ElementMatrix* Ke = analysis->CreateKMatrix(element); … … 98 98 99 99 /*Fill stiffness matrix and load vector from loads*/ 100 for( i=0;i<femmodel->loads->Size();i++){101 load=xDynamicCast<Load*>(femmodel->loads->GetObjectByOffset(i));100 for(Object* & object : femmodel->loads->objects){ 101 load = xDynamicCast<Load*>(object); 102 102 load->CreateKMatrix(Kff,Kfs); 103 103 load->CreatePVector(pf); … … 106 106 /*Now deal with penalties (only in loads)*/ 107 107 if(ispenalty){ 108 for (i=0;i<femmodel->loads->Size();i++){109 load =xDynamicCast<Load*>(femmodel->loads->GetObjectByOffset(i));108 for(Object* & object : femmodel->loads->objects){ 109 load = xDynamicCast<Load*>(object); 110 110 load->PenaltyCreateKMatrix(Kff,Kfs,kmax); 111 111 load->PenaltyCreatePVector(pf,kmax); … … 115 115 /*Create dof vector for stiffness matrix preconditioning*/ 116 116 if(pdf){ 117 for(i=0;i<femmodel->elements->Size();i++){118 element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i));119 120 121 117 for(Object* & object : femmodel->elements->objects){ 118 element = xDynamicCast<Element*>(object); 119 ElementVector* de=analysis->CreateDVector(element); 120 if(de) de->InsertIntoGlobal(df); 121 delete de; 122 122 } 123 123 }
Note:
See TracChangeset
for help on using the changeset viewer.