Changeset 17386 for issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateNodes.cpp

Timestamp:

03/05/14 13:27:33 (11 years ago)

Author:

seroussi

Message:

BUG: fixed problem with 2d taylor hood

File:

• issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateNodes.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateNodes.cpp

@@ -296 +296 @@
                         /*P2 velocity*/
                         EdgesPartitioning(&my_edges,iomodel);
-                        FacesPartitioning(&my_faces,iomodel);
                         for(i=0;i<iomodel->numberofvertices;i++){
                                 if(iomodel->my_vertices[i]){
@@ -308 +307 @@
                         }
                         id0 = id0+iomodel->numberofvertices+iomodel->numberofedges;
-                        for(i=0;i<iomodel->numberoffaces;i++){
-                                if(iomodel->faces[i*6+5]==2){/*Vertical quads*/
-                                        if(my_faces[i]){
-                                                node = new Node(id0+i+1,iomodel->numberofvertices+iomodel->numberofedges+i,lid++,0,iomodel,analysis,FSvelocityEnum);
-                                                nodes->AddObject(node);
-                                        }
-                                }
-                                else if(iomodel->faces[i*6+5]==1){/*Triangular base/top*/
-                                        /*Nothing*/
-                                }
-                                else{
-                                        _error_("not supported");
+              if(iomodel->meshtype==Mesh3DEnum){
+                                FacesPartitioning(&my_faces,iomodel);
+                                for(i=0;i<iomodel->numberoffaces;i++){
+                                        if(iomodel->faces[i*6+5]==2){/*Vertical quads*/
+                                                if(my_faces[i]){
+                                                        node = new Node(id0+i+1,iomodel->numberofvertices+iomodel->numberofedges+i,lid++,0,iomodel,analysis,FSvelocityEnum);
+                                                        nodes->AddObject(node);
+                                                }
+                                        }
+                                        else if(iomodel->faces[i*6+5]==1){/*Triangular base/top*/
+                                                /*Nothing*/
+                                        }
+                                        else{
+                                                _error_("not supported");
+                                        }
                                 }
                         }