[cig-commits] [commit] master: inclusion benchmark (8807ba9)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/485a82a477d4e359644ce65c4155a55cdde612d7...7b60f0be040eacf5c8619a8e155331ef57bec890

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

commit 8807ba93c97450203b4054ad35d849d0ed0b1faa
Author: Timo Heister <timo.heister at gmail.com>
Date:   Thu May 8 13:38:55 2014 -0400

    inclusion benchmark


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

8807ba93c97450203b4054ad35d849d0ed0b1faa
 .../inclusion/adaptive.prm                           | 19 +++++++++++++++++--
 .../inclusion.prm => benchmark/inclusion/base.prm    |  2 +-
 benchmark/inclusion/notes.txt                        |  2 ++
 benchmark/inclusion/run.sh                           | 20 ++++++++++++++++++++
 4 files changed, 40 insertions(+), 3 deletions(-)

diff --git a/cookbooks/future/inclusion.prm b/benchmark/inclusion/adaptive.prm
similarity index 77%
copy from cookbooks/future/inclusion.prm
copy to benchmark/inclusion/adaptive.prm
index 81ededc..ad95fb2 100644
--- a/cookbooks/future/inclusion.prm
+++ b/benchmark/inclusion/adaptive.prm
@@ -65,8 +65,15 @@ end
 
 
 subsection Mesh refinement
-  set Initial adaptive refinement              = 0
-  set Initial global refinement                = 6
+  set Initial adaptive refinement              = 13
+  set Initial global refinement                = 3
+
+  set Strategy = velocity#viscosity
+  set Run postprocessors on initial refinement = true
+  set Coarsening fraction                      = 0.0
+  set Refinement fraction                      = 0.5
+
+
 end
 
 
@@ -74,4 +81,12 @@ end
 
 subsection Postprocess
   set List of postprocessors = DuretzEtAl error, visualization
+
+
+
+  subsection Visualization
+    set List of output variables = viscosity
+    set Time between graphical output = 0
+  end
+
 end
diff --git a/cookbooks/future/inclusion.prm b/benchmark/inclusion/base.prm
similarity index 96%
copy from cookbooks/future/inclusion.prm
copy to benchmark/inclusion/base.prm
index 81ededc..cdeeb62 100644
--- a/cookbooks/future/inclusion.prm
+++ b/benchmark/inclusion/base.prm
@@ -66,7 +66,7 @@ end
 
 subsection Mesh refinement
   set Initial adaptive refinement              = 0
-  set Initial global refinement                = 6
+  set Initial global refinement                = 0
 end
 
 
diff --git a/benchmark/inclusion/notes.txt b/benchmark/inclusion/notes.txt
new file mode 100644
index 0000000..9dd6f56
--- /dev/null
+++ b/benchmark/inclusion/notes.txt
@@ -0,0 +1,2 @@
+This directory contains the "pure shear/inclusion" benchmark defined in Duretz et al 2011
+The results were published in Kronbichler,Heister,Bangerth,2012
diff --git a/benchmark/inclusion/run.sh b/benchmark/inclusion/run.sh
new file mode 100755
index 0000000..f007b62
--- /dev/null
+++ b/benchmark/inclusion/run.sh
@@ -0,0 +1,20 @@
+#!/bin/bash
+
+mkdir output
+
+# global ref:
+for r in "3" "4" "5" "6" "7" "8"
+do
+echo "ref $r:"
+cp base.prm temp.prm
+echo "subsection Mesh refinement" >>temp.prm
+echo "set Initial global refinement = $r" >> temp.prm
+echo "end" >> temp.prm
+../../aspect temp.prm | grep Error
+rm -f temp.prm
+done
+
+
+
+# adaptive ref:
+../../aspect adaptive.prm | egrep "freedom|Error"