Checking in debug version temporarily.
authorTerry L. Wilmarth <wilmarth@uiuc.edu>
Tue, 2 Aug 2005 15:59:04 +0000 (15:59 +0000)
committerTerry L. Wilmarth <wilmarth@uiuc.edu>
Tue, 2 Aug 2005 15:59:04 +0000 (15:59 +0000)
src/libs/ck-libs/TMRC2D/femrefine.C
src/libs/ck-libs/TMRC2D/tri.C

index cb6d0d192ea994762e8749765ee48906a4c212e0..1948e030259cfe7befd8b2d0585ac4f44caf1c17 100644 (file)
@@ -267,6 +267,7 @@ void FEM_Modify_IDXL(FEM_Refine_Operation_Data *data,refineData &op){
 };
 
 void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
+  CmiMemoryCheck();
   int meshID = data->meshID;
   int nodeID = data->nodeID;
   int sparseID = data->sparseID;
@@ -279,7 +280,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
   double *coord = data->coord;
   AllocTable2d<int> *sparseConnTable, *sparseBoundaryTable;
   CkVec<FEM_Attribute *> *sparseattrs = data->sparseattrs;
-  
+  CmiMemoryCheck();  
   /*
   FEM_DataAttribute *boundaryAttr = (FEM_DataAttribute *)data->node->lookup(FEM_BOUNDARY,"split");
   if(boundaryAttr != NULL){
@@ -296,6 +297,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
   // current number of nodes in the mesh
   int *connData = connTable->getData();
   int flags=op.flag;
+  CmiMemoryCheck();
     
   if((flags & 0x1) || (flags & 0x2)){
     //new node 
@@ -309,6 +311,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
       data->node->setLength(D+1);
       data->cur_nodes = D+1;
     }  
+  CmiMemoryCheck();
     for(int i=0;i<attrs->size();i++){
       FEM_Attribute *a = (FEM_Attribute *)(*attrs)[i];
       if(a->getAttr()<FEM_ATTRIB_TAG_MAX){
@@ -334,8 +337,9 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
          }
        }
       }        
+      CmiMemoryCheck();
     }
-
+    CmiMemoryCheck();
     int AandB[2];
     AandB[0]=A;
     AandB[1]=B;
@@ -349,6 +353,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
     /* add the new sparse element <D,B> and modify the connectivity of the old
        one from <A,B> to <A,D> and change the hashtable to reflect that change
     */
+  CmiMemoryCheck();
     if(sparseID != -1){
       int oldsidx = nodes2sparse->get(intdual(A,B))-1;
       int newsidx = data->sparse->size();
@@ -376,6 +381,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
       nodes2sparse->put(intdual(D,B)) = newsidx+1;
     }
   }
+  CmiMemoryCheck();
   //add a new triangle
   /*TODO: replace  FEM_ELEM with parameter*/
   int newTri = op._new;
@@ -385,6 +391,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
     data->elem->setLength(newTri+1);
   }    
   D = newnodes->get(intdual(A,B));
+  CmiMemoryCheck();
   for(int j=0;j<elemattrs->size();j++){
     if((*elemattrs)[j]->getAttr() == FEM_CONN){
       DEBUGINT(CkPrintf("elem attr conn code %d \n",(*elemattrs)[j]->getAttr()));
@@ -419,6 +426,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
       }
     }
   }
+  CmiMemoryCheck();
   if(sparseID != -1){ /* add sparse element (edge between C and D) */
     int cdidx = data->sparse->size();
     data->sparse->setLength(cdidx+1);
@@ -440,6 +448,7 @@ void FEM_Refine_Operation(FEM_Refine_Operation_Data *data,refineData &op){
     }
     nodes2sparse->put(intdual(C,D)) = cdidx+1;
   }
+  CmiMemoryCheck();
 }
 
 
index 50bf0420f2c08f935a922ca6da89a467c38a3a95..31e68c3d2f53b7d9f501647837a94eaf4e178417 100644 (file)
@@ -826,6 +826,7 @@ void chunk::multipleRefine(double *desiredArea, refineClient *client)
 #ifdef TDEBUG1
   sanityCheck(); // quietly make sure mesh is in shape
 #endif
+  CmiMemoryCheck();
   theClient = client; // initialize refine client associated with this chunk
   
   if (refineStack) delete [] refineStack;
@@ -855,12 +856,14 @@ void chunk::multipleRefine(double *desiredArea, refineClient *client)
   CkPrintf("TMRC2D: [%d] multipleRefine DONE.\n", cid);
 #endif
   CkWaitQD();
+  CmiMemoryCheck();
 }
 
 void chunk::multipleCoarsen(double *desiredArea, refineClient *client)
 {
   int i;
   double precThrshld, area, qFactor;
+  CmiMemoryCheck();
 #ifdef TDEBUG2
   CkPrintf("TMRC2D: [%d] multipleCoarsen....\n", cid);
 #endif
@@ -908,6 +911,7 @@ void chunk::multipleCoarsen(double *desiredArea, refineClient *client)
     }
   }
 */
+  CmiMemoryCheck();
 }
 
 void chunk::newMesh(int meshID_,int nEl, int nGhost, const int *conn_,