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

Timestamp:

08/02/13 11:47:37 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: implemented FS condensed mini

File:

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

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

@@ -130 +130 @@
                 case MINIcondensedEnum:
                         _assert_(approximation==FSApproximationEnum);
-                        /*P1 velocity (bubble statically condensed)*/
-                        for(i=0;i<iomodel->numberofvertices;i++){
-                                if(iomodel->my_vertices[i]){
-                                        nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i,iomodel,analysis,FSvelocityEnum));
-                                }
-                        }
-                        /*P1 pressure*/
-                        for(i=0;i<iomodel->numberofvertices;i++){
-                                if(iomodel->my_vertices[i]){
-                                        nodes->AddObject(new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,i,iomodel,analysis,FSpressureEnum));
+                        /*P1+ velocity (bubble statically condensed)*/
+                        for(i=0;i<iomodel->numberofvertices;i++){
+                                if(iomodel->my_vertices[i]){
+                                        nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i,iomodel,analysis,FSvelocityEnum));
+                                }
+                        }
+                        for(i=0;i<iomodel->numberofelements;i++){
+                                if(iomodel->my_elements[i]){
+                                        node = new Node(iomodel->nodecounter+iomodel->numberofvertices+i+1,iomodel->numberofvertices+i,0,iomodel,analysis,FSvelocityEnum);
+                                        node->Deactivate();
+                                        nodes->AddObject(node);
+                                }
+                        }
+                        /*P1 pressure*/
+                        for(i=0;i<iomodel->numberofvertices;i++){
+                                if(iomodel->my_vertices[i]){
+                                        nodes->AddObject(new Node(iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofelements+i+1,iomodel->numberofvertices+iomodel->numberofelements+i,i,iomodel,analysis,FSpressureEnum));
                                 }
                         }