Changeset 14002

Timestamp:

11/26/12 13:29:44 (12 years ago)

Author:

Mathieu Morlighem

Message:

BUG: allow head_e and head_l to be -1, needed for ISMIP1105

Location:

issm/trunk-jpl/src/c

Files:

• classes/FemModel.cpp (modified) (1 diff)
• modules/ModelProcessorx/DiagnosticVert/CreateNodesDiagnosticVert.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/FemModel.cpp

@@ -502 +502 @@
 
                         /*Loop over elements that hold node number i*/
-                        _assert_(head_e[node->Sid()]!=-1 || head_l[node->Sid()]!=-1);
+                        //if(head_e[node->Sid()]==-1 && head_l[node->Sid()]==-1){
+                        //      printf("[%i] vertex %i\n",IssmComm::GetRank(),node->Sid()+1);
+                        //}
                         for(j=head_e[node->Sid()];j!=-1;j=next_e[j]){
                                 offset=count2offset_e[j];

issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateNodesDiagnosticVert.cpp

@@ -44 +44 @@
 
         for (i=0;i<numberofvertices;i++){
-                if(iomodel->my_vertices[i]==1){
+                if(iomodel->my_vertices[i]){
 
                         /*Add node to nodes dataset: */