Changeset 21540 for issm/trunk-jpl/src/c/classes/AdaptiveMeshRefinement.cpp

Timestamp:

02/09/17 17:55:28 (8 years ago)

Author:

tsantos

Message:

CHG: changes in AMR methods, minor. (NOTE: AMR not working yet)

File:

• issm/trunk-jpl/src/c/classes/AdaptiveMeshRefinement.cpp (modified) (6 diffs)

issm/trunk-jpl/src/c/classes/AdaptiveMeshRefinement.cpp

@@ -180 +180 @@
     /*Verify the new geometry*/
     this->CheckMesh(nvertices,nelements,nsegments,this->elementswidth,px,py,pz,pelements,psegments);
-     
+
+         _printf_("\trefinement process done!\n\n");
+
     delete nohangingnodesmesh;
 }
@@ -187 +189 @@
     
     /*Refine mesh hmax times*/
-    for(int hlevel=1;hlevel<=this->hmax;hlevel++){
+   _printf_("\n\trefinement process (hmax = " << this->hmax << ")\n");
+        _printf_("\tprogress:  ");
+        for(int hlevel=1;hlevel<=this->hmax;hlevel++){
         
         /*Set elements to be refined using some criteria*/
@@ -195 +199 @@
         /*Refine the mesh*/
         this->RefineMesh(gmesh, ElemVec);
-    }
-    
+                  
+                  _printf_("*  ");
+    }
+    _printf_("\n");
 }
 /*}}}*/
@@ -235 +241 @@
 /*}}}*/
 void AdaptiveMeshRefinement::RefineMeshToAvoidHangingNodes(TPZGeoMesh *gmesh){/*{{{*/
-    
+   
+         _printf_("\trefine to avoid hanging nodes...\n");
     /*Refine elements to avoid hanging nodes: non-uniform refinement*/
          const int NElem = gmesh->NElements();
@@ -389 +396 @@
     
     // itapopo inserir modo de encontrar criterio TESTING!!!! Come to false!
-         if(true) this->TagAllElements(gmesh,ElemVec); //uniform, refine all elements!
+         if(false) this->TagAllElements(gmesh,ElemVec); //uniform, refine all elements!
 
     /* Adaptive refinement. This refines some elements following some criteria*/
-    //this->TagElementsNearGroundingLine(gmesh, GLvec, hlevel, ElemVec);
+    this->TagElementsNearGroundingLine(gmesh, GLvec, hlevel, ElemVec);
 
 }
@@ -410 +417 @@
     
     /* Tag elements near grounding line */
-    double MaxRegion = 20000.; //itapopo
-    double alpha = 1.0;         //itapopo
+    double MaxRegion = 40000.; //itapopo
+    double alpha = log(1.5);         //itapopo
     double MaxDistance = MaxRegion / std::exp(alpha*(hlevel-1));