Changeset 3359 for issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateConstraintsPrognostic2.cpp

Timestamp:

03/31/10 16:43:20 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added discontinuous galerkin elements and prognostic2 solution, to be continued...

File:

• issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateConstraintsPrognostic2.cpp (modified) (2 diffs)

issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateConstraintsPrognostic2.cpp

@@ -2 +2 @@
  * CreateConstraintsPrognostic2.c:
  */
-
 
 #include "../../DataSet/DataSet.h"
@@ -11 +10 @@
 #include "../IoModel.h"
 
-
 void    CreateConstraintsPrognostic2(DataSet** pconstraints, IoModel* iomodel,ConstDataHandle iomodel_handle){
-
-
-        int i;
-        int count;
         
         DataSet* constraints = NULL;
-        Spc*    spc  = NULL;
-
-        /*spc intermediary data: */
-        int spc_sid;
-        int spc_node;
-        int spc_dof;
-        double spc_value;
-        
-        double* spcthickness=NULL;
         
         /*Create constraints: */
         constraints = new DataSet(ConstraintsEnum());
-
-        /*Fetch data: */
-        IoModelFetchData(&spcthickness,NULL,NULL,iomodel_handle,"spcthickness");
-
-        count=0;
-
-        /*Create spcs from x,y,z, as well as the spc values on those spcs: */
-        for (i=0;i<iomodel->numberofnodes;i++){
-        #ifdef _PARALLEL_
-        /*keep only this partition's nodes:*/
-        if((iomodel->my_grids[i]==1)){
-        #endif
-
-                if ((int)spcthickness[2*i]){
-        
-                        /*This grid needs to be spc'd: */
-
-                        spc_sid=count;
-                        spc_node=i+1;
-                        spc_dof=1; //we enforce first translation degree of freedom, for temperature
-                        spc_value=*(spcthickness+2*i+1);
-
-                        spc = new Spc(spc_sid,spc_node,spc_dof,spc_value);
-                        constraints->AddObject(spc);
-                        count++;
-                }
-
-        #ifdef _PARALLEL_
-        } //if((my_grids[i]==1))
-        #endif
-        }
-
-        /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these 
-         * datasets, it will not be redone: */
-        constraints->Presort();
-
-        /*Free data: */
-        xfree((void**)&spcthickness);
-        
-        cleanup_and_return:
         
         /*Assign output pointer: */