Changeset 3824

Timestamp:

05/18/10 15:55:12 (15 years ago)

Author:

seroussi

Message:

nan in lower and upper elements

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -58 +58 @@
         penta_matice_id=index+1; //refers to the corresponding ice material object
         penta_matpar_id=iomodel->numberofelements+1; //refers to the constant material parameters object
-        penta_elements_ids[0]=(int)(iomodel->upperelements[index]);
-        penta_elements_ids[1]=(int)(iomodel->lowerelements[index]);
+
+        if isnan(iomodel->upperelements[index]){
+                penta_elements_ids[0]=this->id; //upper penta is the same penta
+        }
+        else{
+                penta_elements_ids[0]=(int)(iomodel->upperelements[index]);
+        }
+        if isnan(iomodel->lowerelements[index]){
+                penta_elements_ids[1]=this->id; //lower penta is the same penta
+        }
+        else{
+                penta_elements_ids[1]=(int)(iomodel->lowerelements[index]);
+        }
 
         this->InitHookNodes(penta_node_ids); this->nodes=NULL;