[cig-commits] r14489 - doc/cigma/manual

luis at geodynamics.org luis at geodynamics.org
Fri Mar 27 12:53:17 PDT 2009

Author: luis
Date: 2009-03-27 12:53:17 -0700 (Fri, 27 Mar 2009)
New Revision: 14489

Modified:
doc/cigma/manual/main2.lyx
Log:
Recover changes from stress

Modified: doc/cigma/manual/main2.lyx
===================================================================
--- doc/cigma/manual/main2.lyx	2009-03-27 18:59:10 UTC (rev 14488)
+++ doc/cigma/manual/main2.lyx	2009-03-27 19:53:17 UTC (rev 14489)
@@ -3762,6 +3762,459 @@
\end_layout

\begin_layout Section
+Laplace Problem
+\end_layout
+
+\begin_layout Standard
+Here we obtain a sequence of solutions over five refinement levels by solving
+ the Laplace problem
+\begin_inset Formula $\nabla^{2}\phi(x,y)=4x^{4}+4y^{4}$
+\end_inset
+
+ inside
+\begin_inset Formula $\Omega=[-1,1]^{2}$
+\end_inset
+
+, subject to
+\begin_inset Formula $\phi(x_{0},y_{0})=x_{0}^{2}+y_{0}^{2}$
+\end_inset
+
+ for points
+\begin_inset Formula $(x_{0},y_{0})$
+\end_inset
+
+ in the boundary
+\begin_inset Formula $\partial\Omega$
+\end_inset
+
+.
+ This problem can easily be solved using the Deal.II library, and is available
+ under Step 4 in the documentation.
+\end_layout
+
+\begin_layout LyX-Code
+$levels= +\begin_inset Quotes erd +\end_inset + +2 3 4 5 +\begin_inset Quotes erd +\end_inset + + +\end_layout + +\begin_layout LyX-Code +$ for i in ${levels}; do +\end_layout + +\begin_layout LyX-Code + cigma compare +\backslash + +\end_layout + +\begin_layout LyX-Code + square6.vtk square${i}.vtk
+\backslash
+
+\end_layout
+
+\begin_layout LyX-Code
+      -o square.h5:/error_6_${i} +\end_layout + +\begin_layout LyX-Code + vtk-residuals --output-log-values +\backslash + +\end_layout + +\begin_layout LyX-Code + -m square6.vtk +\backslash + +\end_layout + +\begin_layout LyX-Code + -i square.h5:/error_6_${i}
+\backslash
+
+\end_layout
+
+\begin_layout LyX-Code
+      -o log_error_square_6_{i}.vtk:log_error +\end_layout + +\begin_layout LyX-Code + done +\end_layout + +\begin_layout Standard +From output of the above commands we can collect the following table +\end_layout + +\begin_layout Standard +\begin_inset Tabular +<lyxtabular version="3" rows="5" columns="3"> +<features> +<column alignment="center" valignment="top" width="0"> +<column alignment="center" valignment="top" width="0"> +<column alignment="center" valignment="top" width="0"> +<row> +<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +Case +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +\begin_inset Formulah+\end_inset + + +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +\begin_inset FormulaL_{2}+\end_inset + + +\end_layout + +\end_inset +</cell> +</row> +<row> +<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +err_6_2 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.70717 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.08280 +\end_layout + +\end_inset +</cell> +</row> +<row> +<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +err_6_3 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.35355 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.02462 +\end_layout + +\end_inset +</cell> +</row> +<row> +<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +err_6_4 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.17677 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.00631 +\end_layout + +\end_inset +</cell> +</row> +<row> +<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +err_6_5 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.08838 +\end_layout + +\end_inset +</cell> +<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none"> +\begin_inset Text + +\begin_layout Plain Layout +0.00137 +\end_layout + +\end_inset +</cell> +</row> +</lyxtabular> + +\end_inset + + +\end_layout + +\begin_layout Standard +\begin_inset Float figure +placement H +wide false +sideways false +status collapsed + +\begin_layout Plain Layout +\begin_inset Graphics + filename figures2/laplace_square6.png + lyxscale 80 + scale 25 + +\end_inset + + +\begin_inset Caption + +\begin_layout Plain Layout +Caption for 5.4 here +\end_layout + +\end_inset + + +\end_layout + +\end_inset + + +\end_layout + +\begin_layout Standard +\begin_inset Float figure +placement H +wide false +sideways false +status collapsed + +\begin_layout Plain Layout +\begin_inset Graphics + filename figures2/log_res_square_6_2.png + lyxscale 40 + scale 18 + +\end_inset + + +\begin_inset Caption + +\begin_layout Plain Layout +Caption for 5.4 here +\end_layout + +\end_inset + + +\end_layout + +\end_inset + + +\end_layout + +\begin_layout Standard +\begin_inset Float figure +placement H +wide false +sideways false +status collapsed + +\begin_layout Plain Layout +\begin_inset Graphics + filename figures2/log_res_square_6_3.png + lyxscale 40 + scale 18 + +\end_inset + + +\begin_inset Caption + +\begin_layout Plain Layout +Caption for 5.4 here +\end_layout + +\end_inset + + +\end_layout + +\end_inset + + +\end_layout + +\begin_layout Standard +\begin_inset Float figure +placement H +wide false +sideways false +status collapsed + +\begin_layout Plain Layout +\begin_inset Graphics + filename figures2/log_res_square_6_4.png + lyxscale 40 + scale 18 + +\end_inset + + +\begin_inset Caption + +\begin_layout Plain Layout +Caption for 5.4 here +\end_layout + +\end_inset + + +\end_layout + +\end_inset + + +\end_layout + +\begin_layout Standard +\begin_inset Float figure +placement H +wide false +sideways false +status collapsed + +\begin_layout Plain Layout +\begin_inset Graphics + filename figures2/log_res_square_6_5.png + lyxscale 40 + scale 18 + +\end_inset + + +\begin_inset Caption + +\begin_layout Plain Layout +Caption for 5.4 here +\end_layout + +\end_inset + + +\end_layout + +\end_inset + + +\end_layout + +\begin_layout Standard +\begin_inset space ~ +\end_inset + + +\begin_inset Float figure +placement H +wide false +sideways false +status open + +\begin_layout Plain Layout +\begin_inset Graphics + filename figures2/alpha_square.png + scale 30 + +\end_inset + + +\begin_inset Caption + +\begin_layout Plain Layout +Caption for 5.4 here +\end_layout + +\end_inset + + +\end_layout + +\end_inset + + +\end_layout + +\begin_layout Section Mantle Convection \end_layout @@ -4250,11 +4703,15 @@ cigma compare ${bm}.pressure p64.vts:${p}  -o circ_inc.h5:/pressure_064
\end_layout

+\begin_layout Plain Layout
+\begin_inset Note Note
+status open
+
\begin_layout Section
Strikeslip Benchmark Convergence
\end_layout

-\begin_layout Standard
+\begin_layout Plain Layout
This benchmark problem computes the viscoelastic (Maxwell) relaxation of
stresses from a single, finite strike-slip earthquake in 3D without gravity.
In order to obtain several data points, we use the CUBIT mesh generator
@@ -4353,460 +4810,11 @@
done
\end_layout

-\begin_layout Section
-Laplace Problem
-\end_layout
-
-\begin_layout Standard
-Here we obtain a sequence of solutions over five refinement levels by solving
- the Laplace problem
-\begin_inset Formula $\nabla^{2}\phi(x,y)=4x^{4}+4y^{4}$
\end_inset

- inside
-\begin_inset Formula $\Omega=[-1,1]^{2}$
-\end_inset

-, subject to
-\begin_inset Formula $\phi(x_{0},y_{0})=x_{0}^{2}+y_{0}^{2}$
-\end_inset
-
- for points
-\begin_inset Formula $(x_{0},y_{0})$
-\end_inset
-
- in the boundary
-\begin_inset Formula $\partial\Omega$
-\end_inset
-
-.
- This problem can easily be solved using the Deal.II library, and is available
- under Step 4 in the documentation.
\end_layout

-\begin_layout LyX-Code
-$levels= -\begin_inset Quotes erd -\end_inset - -2 3 4 5 -\begin_inset Quotes erd -\end_inset - - -\end_layout - -\begin_layout LyX-Code -$ for i in ${levels}; do -\end_layout - -\begin_layout LyX-Code - cigma compare -\backslash - -\end_layout - -\begin_layout LyX-Code - square6.vtk square${i}.vtk
-\backslash
-
-\end_layout
-
-\begin_layout LyX-Code
-      -o square.h5:/error_6_${i} -\end_layout - -\begin_layout LyX-Code - vtk-residuals --output-log-values -\backslash - -\end_layout - -\begin_layout LyX-Code - -m square6.vtk -\backslash - -\end_layout - -\begin_layout LyX-Code - -i square.h5:/error_6_${i}
-\backslash
-
-\end_layout
-
-\begin_layout LyX-Code
-      -o log_error_square_6_{i}.vtk:log_error -\end_layout - -\begin_layout LyX-Code - done -\end_layout - -\begin_layout Standard -From output of the above commands we can collect the following table -\end_layout - -\begin_layout Standard -\begin_inset Tabular -<lyxtabular version="3" rows="5" columns="3"> -<features> -<column alignment="center" valignment="top" width="0"> -<column alignment="center" valignment="top" width="0"> -<column alignment="center" valignment="top" width="0"> -<row> -<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none"> -\begin_inset Text - -\begin_layout Plain Layout -Case -\end_layout - -\end_inset -</cell> -<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none"> -\begin_inset Text - -\begin_layout Plain Layout -\begin_inset Formulah-\end_inset - - -\end_layout - -\end_inset -</cell> -<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none"> -\begin_inset Text - -\begin_layout Plain Layout - -\begin_inset FormulaL_{2}\$
-\end_inset
-
-
-\end_layout
-
-\end_inset
-</cell>
-</row>
-<row>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-err_6_2
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.70717
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.08280
-\end_layout
-
-\end_inset
-</cell>
-</row>
-<row>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-err_6_3
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.35355
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.02462
-\end_layout
-
-\end_inset
-</cell>
-</row>
-<row>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-err_6_4
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.17677
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.00631
-\end_layout
-
-\end_inset
-</cell>
-</row>
-<row>
-<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-err_6_5
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.08838
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-0.00137
-\end_layout
-
-\end_inset
-</cell>
-</row>
-</lyxtabular>
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status collapsed
-
-\begin_layout Plain Layout
-\begin_inset Graphics
-	filename figures2/laplace_square6.png
-	lyxscale 80
-	scale 25
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-Caption for 5.4 here
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status collapsed
-
-\begin_layout Plain Layout
-\begin_inset Graphics
-	filename figures2/log_res_square_6_2.png
-	lyxscale 40
-	scale 18
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-Caption for 5.4 here
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status collapsed
-
-\begin_layout Plain Layout
-\begin_inset Graphics
-	filename figures2/log_res_square_6_3.png
-	lyxscale 40
-	scale 18
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-Caption for 5.4 here
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status collapsed
-
-\begin_layout Plain Layout
-\begin_inset Graphics
-	filename figures2/log_res_square_6_4.png
-	lyxscale 40
-	scale 18
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-Caption for 5.4 here
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status collapsed
-
-\begin_layout Plain Layout
-\begin_inset Graphics
-	filename figures2/log_res_square_6_5.png
-	lyxscale 40
-	scale 18
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-Caption for 5.4 here
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset space ~
-\end_inset
-
-
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\begin_inset Graphics
-	filename figures2/alpha_square.png
-	scale 30
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-Caption for 5.4 here
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
\begin_layout Chapter
\start_of_appendix
Input and Output