Changeset 25972
- Timestamp:
- 02/03/21 19:49:27 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r25958 r25972 3015 3015 IssmDouble* values = xNew<IssmDouble>(numnodes); 3016 3016 3017 if(this->Id()==104 || this->Id()==266 || this->Id()==264 || this->Id()==247 || this->Id()==263 || this->Id()==249 || this->Id()==250){3018 if(enum_type==MaskIceLevelsetEnum) printf("-----------------\n");3019 }3020 3021 3017 /*Use the dof list to index into the solution vector: */ 3022 3018 for(int i=0;i<numnodes;i++){ 3023 3019 values[i]=solution[doflist[i]]; 3024 if(this->Id()==104 || this->Id()==266 || this->Id()==264 || this->Id()==247 || this->Id()==263 || this->Id()==249 || this->Id()==250){3025 if(enum_type==MaskIceLevelsetEnum) printf("%g\n",values[i]);3026 if(enum_type==MaskIceLevelsetEnum && values[i]>0.5 && values[i]<1.5 )this->nodes[i]->DeepEcho();3027 }3028 3020 if(xIsNan<IssmDouble>(values[i])) _error_("NaN found in solution vector"); 3029 3021 if(xIsInf<IssmDouble>(values[i])) _error_("Inf found in solution vector, SID = " << this->Sid()); 3030 3022 } 3031 3032 3023 3033 3024 /*Add input to the element: */
Note:
See TracChangeset
for help on using the changeset viewer.