- Timestamp:
- 09/15/14 16:17:28 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/SystemMatricesx/SystemMatricesx.cpp
r17492 r18521 40 40 /*Get complete stiffness matrix without penalties*/ 41 41 for (i=0;i<femmodel->elements->Size();i++){ 42 element= dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));42 element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i)); 43 43 ElementMatrix* Ke = analysis->CreateKMatrix(element); 44 44 ElementVector* pe = analysis->CreatePVector(element); … … 50 50 51 51 for (i=0;i<femmodel->loads->Size();i++){ 52 load= dynamic_cast<Load*>(femmodel->loads->GetObjectByOffset(i));52 load=xDynamicCast<Load*>(femmodel->loads->GetObjectByOffset(i)); 53 53 if(load->InAnalysis(configuration_type)) load->CreateKMatrix(Kff_temp,NULL); 54 54 } … … 71 71 /*Fill stiffness matrix and load vector from elements*/ 72 72 for (i=0;i<femmodel->elements->Size();i++){ 73 element= dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));73 element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i)); 74 74 ElementMatrix* Ke = analysis->CreateKMatrix(element); 75 75 ElementVector* pe = analysis->CreatePVector(element); … … 85 85 /*Fill stiffness matrix and load vector from loads*/ 86 86 for(i=0;i<femmodel->loads->Size();i++){ 87 load= dynamic_cast<Load*>(femmodel->loads->GetObjectByOffset(i));87 load=xDynamicCast<Load*>(femmodel->loads->GetObjectByOffset(i)); 88 88 if(load->InAnalysis(configuration_type)){ 89 89 load->CreateKMatrix(Kff,Kfs); … … 95 95 if(femmodel->loads->IsPenalty(configuration_type)){ 96 96 for (i=0;i<femmodel->loads->Size();i++){ 97 load= dynamic_cast<Load*>(femmodel->loads->GetObjectByOffset(i));97 load=xDynamicCast<Load*>(femmodel->loads->GetObjectByOffset(i)); 98 98 if(load->InAnalysis(configuration_type)){ 99 99 load->PenaltyCreateKMatrix(Kff,Kfs,kmax); … … 106 106 if(pdf){ 107 107 for(i=0;i<femmodel->elements->Size();i++){ 108 element= dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));108 element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i)); 109 109 ElementVector* de=analysis->CreateDVector(element); 110 110 if(de) de->InsertIntoGlobal(df);
Note:
See TracChangeset
for help on using the changeset viewer.