[cig-commits] commit: Make it use just a single grid for now. Seems to work for viscosity

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


changeset:   104:5ee9bf364e5b
user:        Walter Landry <wlandry at caltech.edu>
date:        Sat Mar 05 11:30:34 2011 -0800
files:       input/shear_corner.input
description:
Make it use just a single grid for now.  Seems to work for viscosity
variations up to 10^10.


diff -r c82543f37f62 -r 5ee9bf364e5b input/shear_corner.input
--- a/input/shear_corner.input	Sat Mar 05 11:29:15 2011 -0800
+++ b/input/shear_corner.input	Sat Mar 05 11:30:34 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 = 100         // Max number of FAC cycles to use
+    max_cycles = 1         // 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 = 25
-    coarse_solver_tolerance = 1e-10
+    coarse_solver_max_iterations = 25000
+    coarse_solver_tolerance = 1e-3
     p_prolongation_method = "P_REFINE" // Type of refinement
       					  // used in prolongation.
                                           // Suggested values are
@@ -90,7 +90,8 @@ CartesianGridGeometry {
   //                        domain is periodic.  Zero indicates not
   //                        periodic, non-zero value indicates periodicity.
   //                        [0]
-  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
 }
@@ -98,11 +99,12 @@ StandardTagAndInitialize {
 StandardTagAndInitialize {
   tagging_method = "REFINE_BOXES"
   RefineBoxes {
-    level_0 = [(0,0),(3,3)]
+    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_1 = [(2,2),(3,3)]
-    level_2 = [(4,4),(5,5)]
+    // level_1 = [(2,2),(3,3)]
+    // level_2 = [(4,4),(5,5)]
     // level_3 = [(8,8),(20,20)]
     //etc.
   }
@@ -134,7 +136,7 @@ PatchHierarchy {
    //              [level 0 entry]
    //   etc....                       
    // }
-   max_levels = 3
+   max_levels = 1
    ratio_to_coarser {
       level_1            = 2, 2
       level_2            = 2, 2
@@ -142,7 +144,7 @@ PatchHierarchy {
       level_4            = 2, 2
    }
    largest_patch_size {
-      level_0 = 32, 32
+      level_0 = 64, 64
       // level_0 = 4, 4
       // all finer levels will use same values as level_0...
    }



More information about the CIG-COMMITS mailing list