[cig-commits] [commit] master: add test (311a53d)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Mon May 19 20:14:36 PDT 2014


Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/76d4275352ef2cae5de9a073acd1c03a92c2670c...4f3d06fd1f3754419813db37ec9ef7f0f6f3cb15

>---------------------------------------------------------------

commit 311a53d02af603c087776988520da5b4d63239c4
Author: Timo Heister <timo.heister at gmail.com>
Date:   Mon May 19 20:41:39 2014 -0400

    add test


>---------------------------------------------------------------

311a53d02af603c087776988520da5b4d63239c4
 .../latent-heat.prm => tests/direct_solver_1.prm   |  9 ++-
 tests/direct_solver_1/screen-output                | 88 ++++++++++++++++++++++
 tests/direct_solver_1/statistics                   | 18 +++++
 3 files changed, 112 insertions(+), 3 deletions(-)

diff --git a/cookbooks/latent-heat.prm b/tests/direct_solver_1.prm
similarity index 95%
copy from cookbooks/latent-heat.prm
copy to tests/direct_solver_1.prm
index 8288722..9174f37 100644
--- a/cookbooks/latent-heat.prm
+++ b/tests/direct_solver_1.prm
@@ -5,11 +5,14 @@
 
 set Dimension = 2
 
+set Use direct solver for Stokes system = true
+
 set Start time                             = 0
-set End time                               = 5e17
+set End time                               = 6e15
 set Use years in output instead of seconds = false
 
 set Output directory                       = output
+set Pressure normalization                 = surface
 
 
 subsection Geometry model
@@ -120,7 +123,7 @@ end
 
 subsection Mesh refinement
   set Initial adaptive refinement        = 0 
-  set Initial global refinement          = 7 
+  set Initial global refinement          = 5
   set Time steps between mesh refinement = 0
 
 end
@@ -155,7 +158,7 @@ subsection Postprocess
     # a steady state). For following the development of the system or checking
     # if the solution already reached steady state, this parameter can be set 
     # to a smaller value.
-    set Time between graphical output = 5e17
+    set Time between graphical output = 0
     set List of output variables      = density
   end
 end
diff --git a/tests/direct_solver_1/screen-output b/tests/direct_solver_1/screen-output
new file mode 100644
index 0000000..c97de0b
--- /dev/null
+++ b/tests/direct_solver_1/screen-output
@@ -0,0 +1,88 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
+--     . running in DEBUG mode
+--     . running with 1 MPI process
+--     . using Trilinos
+-----------------------------------------------------------------------------
+
+Number of active cells: 1,024 (on 6 levels)
+Number of degrees of freedom: 13,764 (9,539+4,225)
+
+*** Timestep 0:  t=0 seconds
+   Solving temperature system... 0 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00000
+
+*** Timestep 1:  t=7.26222e+14 seconds
+   Solving temperature system... 177 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00001
+
+*** Timestep 2:  t=1.4246e+15 seconds
+   Solving temperature system... 20 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00002
+
+*** Timestep 3:  t=2.11579e+15 seconds
+   Solving temperature system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00003
+
+*** Timestep 4:  t=2.79859e+15 seconds
+   Solving temperature system... 19 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00004
+
+*** Timestep 5:  t=3.48103e+15 seconds
+   Solving temperature system... 17 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00005
+
+*** Timestep 6:  t=4.16434e+15 seconds
+   Solving temperature system... 18 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00006
+
+*** Timestep 7:  t=4.85169e+15 seconds
+   Solving temperature system... 16 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00007
+
+*** Timestep 8:  t=5.54591e+15 seconds
+   Solving temperature system... 17 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00008
+
+*** Timestep 9:  t=6e+15 seconds
+   Solving temperature system... 13 iterations.
+
+   Postprocessing:
+     Writing graphical output: output-direct_solver_1/solution-00009
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start    |      5.58s |            |
+|                                             |            |            |
+| Section                         | no. calls |  wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system          |        10 |      1.64s |        29% |
+| Assemble temperature system     |        10 |       1.2s |        22% |
+| Build temperature preconditioner|        10 |    0.0827s |       1.5% |
+| Solve Stokes system             |        10 |      1.39s |        25% |
+| Solve temperature system        |        10 |     0.143s |       2.6% |
+| Initialization                  |         2 |    0.0294s |      0.53% |
+| Postprocessing                  |        10 |       0.7s |        13% |
+| Setup dof systems               |         1 |     0.087s |       1.6% |
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/direct_solver_1/statistics b/tests/direct_solver_1/statistics
new file mode 100644
index 0000000..d0f3311
--- /dev/null
+++ b/tests/direct_solver_1/statistics
@@ -0,0 +1,18 @@
+# 1: Time step number
+# 2: Time (seconds)
+# 3: Number of mesh cells
+# 4: Number of Stokes degrees of freedom
+# 5: Number of temperature degrees of freedom
+# 6: Iterations for temperature solver
+# 7: Time step size (seconds)
+# 8: Visualization file name
+0 0.0000e+00 1024 9539 4225   0 7.2622e+14 output-direct_solver_1/solution-00000 
+1 7.2622e+14 1024 9539 4225 177 6.9838e+14 output-direct_solver_1/solution-00001 
+2 1.4246e+15 1024 9539 4225  20 6.9119e+14 output-direct_solver_1/solution-00002 
+3 2.1158e+15 1024 9539 4225  18 6.8280e+14 output-direct_solver_1/solution-00003 
+4 2.7986e+15 1024 9539 4225  19 6.8244e+14 output-direct_solver_1/solution-00004 
+5 3.4810e+15 1024 9539 4225  17 6.8331e+14 output-direct_solver_1/solution-00005 
+6 4.1643e+15 1024 9539 4225  18 6.8734e+14 output-direct_solver_1/solution-00006 
+7 4.8517e+15 1024 9539 4225  16 6.9422e+14 output-direct_solver_1/solution-00007 
+8 5.5459e+15 1024 9539 4225  17 4.5409e+14 output-direct_solver_1/solution-00008 
+9 6.0000e+15 1024 9539 4225  13 7.1111e+14 output-direct_solver_1/solution-00009 



More information about the CIG-COMMITS mailing list