[cig-commits] commit: Modified input to use multigrid. Seems to work well for viscosity

Mercurial hg at geodynamics.org
Sat Mar 5 20:39:06 PST 2011


changeset:   105:935aa6a028a9
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Sat Mar 05 14:12:56 2011 -0800
files:       input/shear_corner.input
description:
Modified input to use multigrid.  Seems to work well for viscosity
difference of 10.


diff -r 5ee9bf364e5b -r 935aa6a028a9 input/shear_corner.input
--- a/input/shear_corner.input	Sat Mar 05 11:30:34 2011 -0800
+++ b/input/shear_corner.input	Sat Mar 05 14:12:56 2011 -0800
@@ -32,7 +32,7 @@ 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 = 1         // Max number of FAC cycles to use
+    max_cycles = 100         // Max number of FAC cycles to use
     residual_tol = 1e-8     // Residual tolerance to solve for
     num_pre_sweeps = 5      // Number of presmoothing sweeps to use
     num_post_sweeps = 5     // Number of postsmoothing sweeps to use
@@ -40,8 +40,8 @@ FACStokes {
     coarse_solver_choice = "Tackley"
     // smoothing_choice = "Gerya"
     // coarse_solver_choice = "Gerya"
-    coarse_solver_max_iterations = 25000
-    coarse_solver_tolerance = 1e-3
+    coarse_solver_max_iterations = 25
+    coarse_solver_tolerance = 1e-10
     p_prolongation_method = "P_REFINE" // Type of refinement
       					  // used in prolongation.
                                           // Suggested values are
@@ -90,8 +90,8 @@ CartesianGridGeometry {
   //                        domain is periodic.  Zero indicates not
   //                        periodic, non-zero value indicates periodicity.
   //                        [0]
-  domain_boxes = [(0,0), (63,63)]
-  // domain_boxes = [(0,0), (3,3)]
+  // domain_boxes = [(0,0), (63,63)]
+  domain_boxes = [(0,0), (3,3)]
   x_lo         = 0, 0
   x_up         = 1, 1
 }
@@ -99,10 +99,13 @@ StandardTagAndInitialize {
 StandardTagAndInitialize {
   tagging_method = "REFINE_BOXES"
   RefineBoxes {
-    level_0 = [(0,0),(63,63)]
-    // level_0 = [(0,0),(3,3)]
-    // level_1 = [(0,0),(7,7)]
-    // level_2 = [(0,0),(15,15)]
+    // level_0 = [(0,0),(63,63)]
+    level_0 = [(0,0),(3,3)]
+    level_1 = [(0,0),(7,7)]
+    level_2 = [(0,0),(15,15)]
+    level_3 = [(0,0),(31,31)]
+    level_4 = [(0,0),(63,63)]
+    level_5 = [(0,0),(127,127)]
     // level_1 = [(2,2),(3,3)]
     // level_2 = [(4,4),(5,5)]
     // level_3 = [(8,8),(20,20)]
@@ -136,7 +139,7 @@ PatchHierarchy {
    //              [level 0 entry]
    //   etc....                       
    // }
-   max_levels = 1
+   max_levels = 3
    ratio_to_coarser {
       level_1            = 2, 2
       level_2            = 2, 2



More information about the CIG-COMMITS mailing list