Changeset 3453

Timestamp:

04/08/10 08:08:32 (15 years ago)

Author:

seroussi

Message:

added CreateConstraints for Stokes

File:

• issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateConstraintsDiagnosticStokes.cpp (modified) (3 diffs)

issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateConstraintsDiagnosticStokes.cpp

@@ -10 +10 @@
 #include "../IoModel.h"
 
-
 void    CreateConstraintsDiagnosticStokes(DataSet** pconstraints, IoModel* iomodel,ConstDataHandle iomodel_handle){
 
@@ -18 +17 @@
 
         /*spc intermediary data: */
-        int spc_sid;
-        int spc_node;
-        int spc_dof;
-        double spc_value;
-        int count;
-        
         double* gridonstokes=NULL;
-
 
         /*Create constraints: */
@@ -36 +28 @@
         IoModelFetchData(&gridonstokes,NULL,NULL,iomodel_handle,"gridonstokes");
 
-        count=0;
+        count=1; //matlab indexing
         /*Create spcs from x,y,z, as well as the spc values on those spcs: */
         for (i=0;i<iomodel->numberofvertices;i++){
-        #ifdef _PARALLEL_
-        /*keep only this partition's nodes:*/
-        if((iomodel->my_vertices[i]==1)){
-        #endif
+                /*keep only this partition's nodes:*/
+                if((iomodel->my_vertices[i]==1)){
 
-                if ((int)!gridonstokes[i]){
-        
-                        /*This grid will see its vx,vy and vz dofs spc'd to pattyn velocities: */
-                        spc_sid=count;
-                        spc_node=i+1;
-                        spc_dof=1; //we enforce first x translation degree of freedom
-                        spc_value=0;
+                        if ((int)!gridonstokes[i]){
 
-                        spc = new Spc(spc_sid,spc_node,spc_dof,spc_value);
-                        constraints->AddObject(spc);
-                        count++;
+                                /*This grid will see its vx,vy and vz dofs spc'd to pattyn velocities: */
+                                constraints->AddObject(new Spc(count,i+1,1,0));
+                                count++;
 
-                        spc_sid=count;
-                        spc_node=i+1;
-                        spc_dof=2; //we enforce first y translation degree of freedom
-                        spc_value=0;
-                        
-                        spc = new Spc(spc_sid,spc_node,spc_dof,spc_value);
-                        constraints->AddObject(spc);
-                        count++;
-                        
-                        spc_sid=count;
-                        spc_node=i+1;
-                        spc_dof=3; //we enforce first y translation degree of freedom
-                        spc_value=0;
-                        
-                        spc = new Spc(spc_sid,spc_node,spc_dof,spc_value);
-                        constraints->AddObject(spc);
-                        count++;
+                                constraints->AddObject(new Spc(count,i+1,2,0));
+                                count++;
 
-                }
-
-        #ifdef _PARALLEL_
-        } //if((my_vertices[i]==1))
-        #endif
+                                constraints->AddObject(new Spc(count,i+1,3,0));
+                                count++;
+                        }
+                } //if((my_vertices[i]==1))
         }