[cig-commits] commit: Set up an adaptive sinker test problem

[Mercurial]

changeset:   191:ba42d6e531fb
user:        Walter Landry <wlandry at caltech.edu>
date:        Fri Apr 22 19:46:41 2011 -0700
files:       input/sinker.input src/FACStokes/FACStokes.C src/FACStokes/initializeLevelData.C src/main.C
description:
Set up an adaptive sinker test problem


diff -r 2a7cedcc2455 -r ba42d6e531fb input/sinker.input
--- a/input/sinker.input	Fri Apr 22 15:09:04 2011 -0700
+++ b/input/sinker.input	Fri Apr 22 19:46:41 2011 -0700
@@ -32,15 +32,15 @@ FACStokes {
         // This is the input for the cell-centered Stokes FAC solver
         // class in the SAMRAI library.
         enable_logging = TRUE   // Bool flag to switch logging on/off
-        max_cycles = 3         // Max number of FAC cycles to use
+        max_cycles = 5         // Max number of FAC cycles to use
         residual_tol = 1e-3     // Residual tolerance to solve for
-        num_pre_sweeps = 3      // Number of presmoothing sweeps to use
-        num_post_sweeps = 3     // Number of postsmoothing sweeps to use
+        num_pre_sweeps = 5      // Number of presmoothing sweeps to use
+        num_post_sweeps = 5     // Number of postsmoothing sweeps to use
         smoothing_choice = "Tackley"
         coarse_solver_choice = "Tackley"
         // smoothing_choice = "Gerya"
         // coarse_solver_choice = "Gerya"
-        coarse_solver_max_iterations = 3
+        coarse_solver_max_iterations = 25
         coarse_solver_tolerance = 1e-5
         // p_prolongation_method = "P_MDPI_REFINE"
         p_prolongation_method = "P_REFINE"
@@ -118,7 +118,7 @@ PatchHierarchy {
     //              [level 0 entry]
     //   etc....                       
     // }
-    max_levels = 3
+    max_levels = 5
     proper_nesting_buffer = 2, 2, 2, 2, 2, 2
     ratio_to_coarser {
         level_1            = 2, 2
diff -r 2a7cedcc2455 -r ba42d6e531fb src/FACStokes/FACStokes.C
--- a/src/FACStokes/FACStokes.C	Fri Apr 22 15:09:04 2011 -0700
+++ b/src/FACStokes/FACStokes.C	Fri Apr 22 19:46:41 2011 -0700
@@ -139,8 +139,8 @@ namespace SAMRAI {
                                                   operator */);
 
     d_adaptation_threshold=database->getDoubleWithDefault("adaption_threshold",
-                                                          1.0e-15);
-                                                          // 1e1);
+                                                          // 1.0e-15);
+                                                          2e-3);
 
     /*
      * Specify an implementation of solv::RobinBcCoefStrategy for the
diff -r 2a7cedcc2455 -r ba42d6e531fb src/FACStokes/initializeLevelData.C
--- a/src/FACStokes/initializeLevelData.C	Fri Apr 22 15:09:04 2011 -0700
+++ b/src/FACStokes/initializeLevelData.C	Fri Apr 22 19:46:41 2011 -0700
@@ -46,11 +46,11 @@ void SAMRAI::FACStokes::initializeLevelD
   }
 
   // const double inclusion_radius=0.5;
-  const double inclusion_viscosity=1e2;
+  const double inclusion_viscosity=2;
   const double background_viscosity=1;
 
-  const double background_density(1), block_density(1);
-  // const double background_density(1), block_density(1.03);
+  // const double background_density(1), block_density(1);
+  const double background_density(1), block_density(1.03);
 
   static double viscosity[129][129];
   static int max_x(0), max_y(0);
diff -r 2a7cedcc2455 -r ba42d6e531fb src/main.C
--- a/src/main.C	Fri Apr 22 15:09:04 2011 -0700
+++ b/src/main.C	Fri Apr 22 19:46:41 2011 -0700
@@ -302,7 +302,7 @@ int main(
       {
             tbox::Array<int> tag_buffer(patch_hierarchy->getMaxNumberOfLevels());
             for (int ln = 0; ln < tag_buffer.getSize(); ++ln) {
-               tag_buffer[ln] = 2;
+               tag_buffer[ln] = 1;
             }
             gridding_algorithm->regridAllFinerLevels(
                0,