[cig-commits] [commit] knepley/fix-faults-parallel: Updated documentation (including README and release notes) about buried fault edges. (528d0bf)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Tue May 20 20:52:46 PDT 2014
Repository : https://github.com/geodynamics/pylith
On branch : knepley/fix-faults-parallel
Link : https://github.com/geodynamics/pylith/compare/8174f906ec3de3aed714ad6c753af1116600da7b...528d0bfc0b629761b6ba7419e807c4bdbe47ef9c
>---------------------------------------------------------------
commit 528d0bfc0b629761b6ba7419e807c4bdbe47ef9c
Author: Brad Aagaard <baagaard at usgs.gov>
Date: Tue May 20 18:50:36 2014 -0700
Updated documentation (including README and release notes) about buried fault edges.
Also switched to figure similar to JGR paper to illustrate insertion of cohesive cells.
>---------------------------------------------------------------
528d0bfc0b629761b6ba7419e807c4bdbe47ef9c
README | 14 +
doc/releasenotes/announce_v2.0.0.txt | 14 +
.../boundaryconditions/boundaryconditions.lyx | 50 ++-
.../boundaryconditions}/figs/cohesivecell.pdf | Bin 61338 -> 58493 bytes
.../boundaryconditions/figs/cohesivecell.tex | 361 +++++++++++++++++++++
.../boundaryconditions}/figs/figcolors.tex | 0
.../boundaryconditions/figs/quad4cohesivekin.fig | 122 -------
.../boundaryconditions/figs/quad4cohesivekin.pdf | Bin 5081 -> 0 bytes
8 files changed, 425 insertions(+), 136 deletions(-)
diff --git a/README b/README
index e463fca..abe5f05 100644
--- a/README
+++ b/README
@@ -77,6 +77,11 @@ name of the ParMetis/Metis partitioner was changed from "parmetis" to
[pylithapp.mesh_generator]
distributor.partitioner = metis
+Buried edges of faults are handled differently in v2.0. A separate
+nodeset/pset should be created and contain the vertices on the buried
+edges of the fault. See the Section 6.4.2 of the PyLith manual for
+more information.
+
----------------------------------------------------------------------
Version 2.0.0
----------------------------------------------------------------------
@@ -91,6 +96,15 @@ Version 2.0.0
results in a more efficient memory layout, resulting in better
performance.
+* Improved treatment of buried fault edges, so that the slip naturally
+ tapers to zero along the buried edges.
+
+ An additional nodeset/pset is used to designate the buried edges of
+ a fault. This allows the cohesive cells to be inserted up to the
+ edge of the fault without splitting the mesh at the fault edge. The
+ slip will naturally taper to zero at along the buried edges as a
+ result of how the cohesive cells are created.
+
* Switched from using Subversion to Git for version control.
The source code repository changed from a CIG maintained Subversion
diff --git a/doc/releasenotes/announce_v2.0.0.txt b/doc/releasenotes/announce_v2.0.0.txt
index 2e31d56..38085e8 100644
--- a/doc/releasenotes/announce_v2.0.0.txt
+++ b/doc/releasenotes/announce_v2.0.0.txt
@@ -54,6 +54,11 @@ name of the ParMetis/Metis partitioner was changed from "parmetis" to
[pylithapp.mesh_generator]
distributor.partitioner = metis
+Buried edges of faults are handled differently in v2.0. A separate
+nodeset/pset should be created and contain the vertices on the buried
+edges of the fault. See the Section 6.4.2 of the PyLith manual for
+more information.
+
======================================================================
RELEASE NOTES - PyLith v2.0.0
======================================================================
@@ -68,6 +73,15 @@ RELEASE NOTES - PyLith v2.0.0
results in a more efficient memory layout, resulting in better
performance.
+* Improved treatment of buried fault edges, so that the slip naturally
+ tapers to zero along the buried edges.
+
+ An additional nodeset/pset is used to designate the buried edges of
+ a fault. This allows the cohesive cells to be inserted up to the
+ edge of the fault without splitting the mesh at the fault edge. The
+ slip will naturally taper to zero at along the buried edges as a
+ result of how the cohesive cells are created.
+
* Switched from using Subversion to Git for version control.
The source code repository changed from a CIG maintained Subversion
diff --git a/doc/userguide/boundaryconditions/boundaryconditions.lyx b/doc/userguide/boundaryconditions/boundaryconditions.lyx
index 05dc580..5c50937 100644
--- a/doc/userguide/boundaryconditions/boundaryconditions.lyx
+++ b/doc/userguide/boundaryconditions/boundaryconditions.lyx
@@ -2248,12 +2248,15 @@ t mesh, additional degrees of freedom are added along with adjustment of
Figure
\begin_inset CommandInset ref
LatexCommand ref
-reference "fig:fault:cohesive:quad4"
+reference "fig:fault:cohesive:cells"
\end_inset
- illustrates the results of inserting a cohesive cell in a mesh consisting
- of two quadrilateral cells.
+ illustrates the results of inserting cohesive cells in a mesh consisting
+ of triangular cells.
+ This example also shows the distinction between how buried fault edges
+ are handled differently than fault edges that reach the edge of the domain,
+ such as the ground surface.
\end_layout
\begin_layout Standard
@@ -2268,7 +2271,7 @@ status open
\begin_layout Plain Layout
\align center
\begin_inset Graphics
- filename figs/quad4cohesivekin.pdf
+ filename figs/cohesivecell.pdf
\end_inset
@@ -2279,13 +2282,12 @@ status open
\begin_inset Caption
\begin_layout Plain Layout
-Example of cohesive cell inserted into a mesh of two quadrilateral cells.
- The zero thickness cohesive cell (shown with dashed lines) controls slip
- on the fault via the relative motion between vertices 3 and 7 and 2 and
- 6.
+Example of cohesive cells inserted into a mesh of triangular cells.
+ The zero thickness cohesive cells control slip on the fault via the relative
+ motion between the vertices on the positive and negative sides of the fault.
\begin_inset CommandInset label
LatexCommand label
-name "fig:fault:cohesive:quad4"
+name "fig:fault:cohesive:cells"
\end_inset
@@ -2302,6 +2304,19 @@ name "fig:fault:cohesive:quad4"
\end_layout
+\begin_layout Standard
+For faults that have buried edges, splitting the mesh apart and inserting
+ the cohesive cells becomes complex at the buried edges due to the ambiguity
+ of defining where the fault ends and how to insert the cohesive cell.
+ In PyLith v2.0.0 we have changed how the buried edges of the fault are managed.
+ An additional group of fault nodes is specified (e.g., via a nodeset from
+ CUBIT) that marks the buried edges of the fault.
+ This allows the cohesive cell insertion algorithm to djust the topology
+ so that cohseive cells are inserted up to the buried edge of the fault
+ but no additional degrees of freedom are added on the fault edge.
+ This naturally forces slip to zero along the buried edges.
+\end_layout
+
\begin_layout Subsection
Fault Parameters
\end_layout
@@ -2322,9 +2337,12 @@ material-id
\end_layout
\begin_layout Description
-name Name or label for the fault.
- This associates sets of vertices with the fault surface in addition to
- being used in error and diagnostic reports.
+label Name of group of vertices associated with the fault surface.
+ This label is also used in error and diagnostic reports.
+\end_layout
+
+\begin_layout Description
+edge Name of group of vertices marking the buried edges of the fault.
\end_layout
\begin_layout Description
@@ -2409,7 +2427,11 @@ fault = pylith.faults.FaultCohesiveKin ; default
\begin_inset Newline newline
\end_inset
-label = fault A ; Group of vertices defining the fault surface
+label = fault_A ; Group of vertices defining the fault surface
+\end_layout
+
+\begin_layout LyX-Code
+edge = fault_edge ; Group of vertices defining the buried edges
\end_layout
\begin_layout LyX-Code
@@ -2474,7 +2496,7 @@ In conventional implementations the additional degrees of freedom associated
as shown in Figure
\begin_inset CommandInset ref
LatexCommand ref
-reference "fig:fault:cohesive:quad4"
+reference "fig:fault:cohesive:cells"
\end_inset
diff --git a/doc/presentations/pylith1.7/figs/cohesivecell.pdf b/doc/userguide/boundaryconditions/figs/cohesivecell.pdf
similarity index 79%
copy from doc/presentations/pylith1.7/figs/cohesivecell.pdf
copy to doc/userguide/boundaryconditions/figs/cohesivecell.pdf
index eb50712..0d51b05 100644
Binary files a/doc/presentations/pylith1.7/figs/cohesivecell.pdf and b/doc/userguide/boundaryconditions/figs/cohesivecell.pdf differ
diff --git a/doc/userguide/boundaryconditions/figs/cohesivecell.tex b/doc/userguide/boundaryconditions/figs/cohesivecell.tex
new file mode 100644
index 0000000..9f65a1d
--- /dev/null
+++ b/doc/userguide/boundaryconditions/figs/cohesivecell.tex
@@ -0,0 +1,361 @@
+\documentclass{standalone}
+\usepackage[none]{hyphenat}
+\usepackage{helvet}
+\renewcommand{\familydefault}{phv}
+
+\makeatletter
+\renewcommand{\large}{\@setfontsize\large{10}{9}}
+\renewcommand{\normalsize}{\@setfontsize\normalsize{8}{9.5}}
+\makeatother
+
+\usepackage{tikz}
+
+\begin{document}
+\pagestyle{empty}
+
+\input{figcolors}
+\pgfdeclarelayer{background}
+\pgfsetlayers{background,main}
+
+
+\begin{tikzpicture}[>=latex,line width=0.75pt]
+
+\usetikzlibrary{shapes,calc}%
+
+\colorlet{colpositive}{red}
+\colorlet{colnegative}{blue}
+\colorlet{collagrange}{orange}
+\colorlet{coledge}{ltgreen}
+\colorlet{colother}{dkslate}
+
+\tikzstyle{title}=[text width=37mm,font=\bfseries\large]
+
+\tikzstyle{vertex}=[circle,color=colother,fill=colother,inner sep=1.5pt]
+\tikzstyle{vertexe}=[circle,color=coledge,fill=coledge,inner sep=1.5pt]
+\tikzstyle{vertexp}=[circle,color=colpositive,fill=colpositive,inner sep=1.5pt]
+\tikzstyle{vertexn}=[circle,color=colnegative,fill=colnegative,inner sep=1.5pt]
+\tikzstyle{vertexl}=[circle,color=collagrange,fill=collagrange,inner sep=1.5pt]
+\tikzstyle{edgel}=[color=collagrange,line width=1.0pt]
+
+\tikzstyle{fsurfp}=[color=colpositive,line width=1.5pt]
+\tikzstyle{fsurfn}=[color=colnegative,line width=1.5pt]
+
+\tikzstyle{cell}=[color=black,line width=0.5pt]
+\tikzstyle{cellp}=[fill=colpositive!30]
+\tikzstyle{celln}=[fill=colnegative!30]
+
+\tikzstyle{classify}=[font=\bfseries\large]
+
+% Reference points
+\coordinate (o1) at (-60mm,0mm);
+\coordinate (o2) at (-20mm,0mm);
+\coordinate (o3) at (+20mm,0mm);
+\coordinate (o4) at (+60mm,0mm);
+
+
+% ----------------------------------------------------------------------
+% Original mesh
+
+\node[title] at ($ (o1)+(0,40mm) $) {(a) Original mesh};
+
+% Vertices
+\node[vertexn] (v1) at ($ (o1)+(0mm,+30mm) $) {};
+\node[vertexn] (v2) at ($ (o1)+(0mm,+10mm) $) {};
+\node[vertexn] (v3) at ($ (o1)+(0mm,-10mm) $) {};
+\node[vertexe] (v4) at ($ (o1)+(0mm,-30mm) $) {};
+
+\node[vertex] (v5) at ($ (o1)+(-11.25mm,0.0mm) $) {};
+\node[vertex] (v6) at ($ (o1)+(+11.25mm,0.0mm) $) {};
+
+\node[vertex] (v7) at ($ (o1)+(-11.25mm,+20.0mm) $) {};
+\node[vertex] (v8) at ($ (o1)+(+11.25mm,+20.0mm) $) {};
+\node[vertex] (v9) at ($ (o1)+(-11.25mm,-20.0mm) $) {};
+\node[vertex] (v10) at ($ (o1)+(+11.25mm,-20.0mm) $) {};
+
+\node[vertex] (v11) at ($ (o1)+(-11.25mm,+30.0mm) $) {};
+\node[vertex] (v12) at ($ (o1)+(+11.25mm,+30.0mm) $) {};
+\node[vertex] (v13) at ($ (o1)+(-7.5mm,-40.0mm) $) {};
+\node[vertex] (v14) at ($ (o1)+(+7.5mm,-40.0mm) $) {};
+
+% Cells
+\draw[cell] (v1) -- (v11) -- (v7) -- (v1);
+\draw[cell] (v1) -- (v7) -- (v2) -- (v1);
+\draw[cell] (v2) -- (v7) -- (v5) -- (v2);
+\draw[cell] (v2) -- (v5) -- (v3) -- (v2);
+\draw[cell] (v3) -- (v5) -- (v9) -- (v3);
+\draw[cell] (v3) -- (v9) -- (v4) -- (v3);
+\draw[cell] (v4) -- (v9) -- (v13) -- (v4);
+\draw[cell] (v4) -- (v13) -- (v14) -- (v4);
+\draw[cell] (v4) -- (v14) -- (v10) -- (v4);
+\draw[cell] (v3) -- (v4) -- (v10) -- (v3);
+\draw[cell] (v3) -- (v10) -- (v6) -- (v3);
+\draw[cell] (v3) -- (v6) -- (v2) -- (v3);
+\draw[cell] (v2) -- (v6) -- (v8) -- (v2);
+\draw[cell] (v1) -- (v2) -- (v8) -- (v1);
+\draw[cell] (v1) -- (v8) -- (v12) -- (v1);
+
+\draw[fsurfn] (v1) -- (v2) -- (v3) -- (v4);
+
+\node(tlabel) [color=colnegative] at ($ (o1)+(-19mm,+15mm) $) {fault vertex};
+\path[color=colnegative,->] (tlabel) edge [bend right] (v2);
+
+\node(tlabel) [color=coledge] at ($ (o1)+(-10mm,-45mm) $) {fault edge vertex};
+\path[color=coledge,->] (tlabel) edge [bend right] (v4);
+
+\coordinate[color=black,label=right:{$\mathbf{n}$}] (nlabel) at ($ (o1)+(5mm,0mm) $);
+\path[color=black,->] ($ (o1) $) edge (nlabel);
+
+% ----------------------------------------------------------------------
+% Add vertices, create cohesive cell
+
+\node[title] at ($ (o2)+(0,40mm) $) {(b) Add colocated vertices};
+
+% Vertices
+\node[vertexe] (v4) at ($ (o2)+(0mm,-30mm) $) {};
+
+\node[vertexn] (v1) at ($ (o2)+(-1.5mm,+30mm) $) {};
+\node[vertexn] (v2) at ($ (o2)+(-1.5mm,+10mm) $) {};
+\node[vertexn] (v3) at ($ (o2)+(-1.5mm,-10mm) $) {};
+
+\node[vertexp] (v5) at ($ (o2)+(+1.5mm,+30mm) $) {};
+\node[vertexp] (v6) at ($ (o2)+(+1.5mm,+10mm) $) {};
+\node[vertexp] (v7) at ($ (o2)+(+1.5mm,-10mm) $) {};
+
+
+% surface and edges
+\draw[fsurfn] (v1) -- (v2) -- (v3) -- (v4);
+\draw[fsurfp] (v5) -- (v6) -- (v7) -- (v4);
+\draw[edgel] (v1) -- (v5);
+\draw[edgel] (v2) -- (v6);
+\draw[edgel] (v3) -- (v7);
+
+% Fault surface labels
+\node(tlabel) [color=colnegative] at ($ (o2)+(-10mm,+20mm) $) {$S_{f-}$};
+\path [color=colnegative,->] (tlabel.south) edge [bend right] ($ (o2)+(-2mm,+7.5mm) $);
+
+\node(tlabel) [color=colpositive] at ($ (o2)+(+10mm,+20mm) $) {$S_{f+}$};
+\path [color=colpositive,->] (tlabel.south) edge [bend left] ($ (o2)+(+2mm,+7.5mm) $);
+
+
+% Fault vertex labels
+\node(tlabel) [text width=25mm,align=right,color=colnegative] at ($ (o2)+(-17mm,-30mm) $) {Original fault vertex (negative side)};
+\path [color=colnegative,->] (tlabel.north east) edge [bend right] (v3);
+
+\node(tlabel) [text width=25mm,align=right,color=collagrange] at ($ (o2)+(-12mm,-40mm) $) {Add Lagrange multiplier edge};
+\path [color=collagrange,->] (tlabel.north east) edge [bend right] ($ (o2)+(0mm,-10mm) $);
+
+\node(tlabel) [text width=25mm,align=right,color=colpositive] at ($(o2)+(-7mm,-50mm) $) {Add vertex on positive side};
+\path [color=colpositive,->] (tlabel.north east) edge [bend right] (v7);
+
+
+
+% ----------------------------------------------------------------------
+% Update cells with faces on fault
+
+\node[title] at ($ (o3)+(0,40mm) $) {(c) Update cells with fault faces};
+
+\coordinate (p1) at ($ (o3)+(0mm,+30mm)-(1.5mm,0mm) $);
+\coordinate (p2) at ($ (o3)+(0mm,+30mm)+(1.5mm,0mm) $);
+
+\coordinate (p3) at ($ (o3)+(0mm,+10mm)-(1.5mm,0mm) $);
+\coordinate (p4) at ($ (o3)+(0mm,+10mm)+(1.5mm,0mm) $);
+
+\coordinate (p5) at ($ (o3)+(0mm,-10mm)-(1.5mm,0mm) $);
+\coordinate (p6) at ($ (o3)+(0mm,-10mm)+(1.5mm,0mm) $);
+
+\coordinate (p7) at ($ (o3)+(0mm,-30mm) $);
+
+\coordinate (p8) at ($ (o3)+(-11.25mm,+20.0mm) $);
+\coordinate (p9) at ($ (o3)+(+11.25mm,+20.0mm) $);
+\coordinate (p10) at ($ (o3)+(-11.25mm, +0.0mm) $);
+\coordinate (p11) at ($ (o3)+(+11.25mm, +0.0mm) $);
+\coordinate (p12) at ($ (o3)+(-11.25mm,-20.0mm) $);
+\coordinate (p13) at ($ (o3)+(+11.25mm,-20.0mm) $);
+
+
+% Vertices
+\node[vertexe] (v7) at (p7) {};
+
+\node[vertexn] (v1) at (p1) {};
+\node[vertexn] (v3) at (p3) {};
+\node[vertexn] (v5) at (p5) {};
+
+\node[vertexp] (v2) at (p2) {};
+\node[vertexp] (v4) at (p4) {};
+\node[vertexp] (v6) at (p6) {};
+
+\node[vertex] (v8) at (p8) {};
+\node[vertex] (v9) at (p9) {};
+\node[vertex] (v10) at (p10) {};
+\node[vertex] (v11) at (p11) {};
+\node[vertex] (v12) at (p12) {};
+\node[vertex] (v13) at (p13) {};
+
+% surface and edges
+\draw[fsurfp] (v2) -- (v4) -- (v6) -- (v7);
+\draw[fsurfn] (v1) -- (v3) -- (v5) -- (v7);
+\draw[edgel] (v1) -- (v2);
+\draw[edgel] (v3) -- (v4);
+\draw[edgel] (v5) -- (v6);
+
+
+\begin{pgfonlayer}{background}
+% cells
+\fill[celln] (p1) -- (p8) -- (p3) -- cycle;
+\fill[celln] (p3) -- (p10) -- (p5) -- cycle;
+\fill[celln] (p5) -- (p12) -- (p7) -- cycle;
+
+\fill[cellp] (p2) -- (p4) -- (p9) -- cycle;
+\fill[cellp] (p4) -- (p6) -- (p11) -- cycle;
+\fill[cellp] (p6) -- (p7) -- (p13) -- cycle;
+
+\draw[cell] (v1) -- (v8) -- (v3) -- (v1);
+\draw[cell] (v3) -- (v10) -- (v5) -- (v3);
+\draw[cell] (v5) -- (v12) -- (v7) -- (v5);
+
+\draw[cell] (v2) -- (v4) -- (v9) -- (v2);
+\draw[cell] (v4) -- (v6) -- (v11) -- (v4);
+\draw[cell] (v6) -- (v7) -- (v13) -- (v6);
+
+\end{pgfonlayer}
+
+\node[classify] at ($ 1/3*(p1)+1/3*(p8)+1/3*(p3) $) {-};
+\node[classify] at ($ 1/3*(p3)+1/3*(p10)+1/3*(p5) $) {-};
+\node[classify] at ($ 1/3*(p5)+1/3*(p12)+1/3*(p7) $) {-};
+
+\node[classify] at ($ 1/3*(p2)+1/3*(p4)+1/3*(p9) $) {+};
+\node[classify] at ($ 1/3*(p4)+1/3*(p6)+1/3*(p11) $) {+};
+\node[classify] at ($ 1/3*(p6)+1/3*(p7)+1/3*(p13) $) {+};
+
+% Fault vertex labels
+\node(tlabel) [text width=12mm,color=colnegative] at ($ (o3)+(-12mm,-40mm) $) {Cell on negative side};
+\path [color=colnegative,->] (tlabel.north) edge [bend left] ($ (v7)+(-5mm,+5mm) $);
+
+\node(tlabel) [text width=12mm,color=colpositive] at ($ (o3)+(+12mm,-40mm) $) {Cell on positive side};
+\path [color=colpositive,->] (tlabel.north) edge [bend right] ($ (v7)+(+5mm,+5mm) $);
+
+
+% ----------------------------------------------------------------------
+% Classify cells and update other cells
+
+\node[title] at ($ (o4)+(0,40mm) $) {(d) Classify cells and update remaining cells};
+
+\coordinate (p1) at ($ (o4)+(0mm,+30mm)-(1.5mm,0mm) $);
+\coordinate (p2) at ($ (o4)+(0mm,+30mm)+(1.5mm,0mm) $);
+
+\coordinate (p3) at ($ (o4)+(0mm,+10mm)-(1.5mm,0mm) $);
+\coordinate (p4) at ($ (o4)+(0mm,+10mm)+(1.5mm,0mm) $);
+
+\coordinate (p5) at ($ (o4)+(0mm,-10mm)-(1.5mm,0mm) $);
+\coordinate (p6) at ($ (o4)+(0mm,-10mm)+(1.5mm,0mm) $);
+
+\coordinate (p7) at ($ (o4)+(0mm,-30mm) $);
+
+\coordinate (p8) at ($ (o4)+(-11.25mm,+20.0mm) $);
+\coordinate (p9) at ($ (o4)+(+11.25mm,+20.0mm) $);
+\coordinate (p10) at ($ (o4)+(-11.25mm, +0.0mm) $);
+\coordinate (p11) at ($ (o4)+(+11.25mm, +0.0mm) $);
+\coordinate (p12) at ($ (o4)+(-11.25mm,-20.0mm) $);
+\coordinate (p13) at ($ (o4)+(+11.25mm,-20.0mm) $);
+
+\coordinate (p14) at ($ (o4)+(-11.25mm,+30.0mm) $);
+\coordinate (p15) at ($ (o4)+(+11.25mm,+30.0mm) $);
+
+\coordinate (p16) at ($ (o4)+(-7.5mm,-40.0mm) $);
+\coordinate (p17) at ($ (o4)+(+7.5mm,-40.0mm) $);
+
+
+% Vertices
+\node[vertexe] (v7) at (p7) {};
+
+\node[vertexn] (v1) at (p1) {};
+\node[vertexn] (v3) at (p3) {};
+\node[vertexn] (v5) at (p5) {};
+
+\node[vertexp] (v2) at (p2) {};
+\node[vertexp] (v4) at (p4) {};
+\node[vertexp] (v6) at (p6) {};
+
+\node[vertex] (v8) at (p8) {};
+\node[vertex] (v9) at (p9) {};
+\node[vertex] (v10) at (p10) {};
+\node[vertex] (v11) at (p11) {};
+\node[vertex] (v12) at (p12) {};
+\node[vertex] (v13) at (p13) {};
+
+\node[vertex] (v14) at (p14) {};
+\node[vertex] (v15) at (p15) {};
+\node[vertex] (v16) at (p16) {};
+\node[vertex] (v17) at (p17) {};
+
+% surface and edges
+\draw[fsurfp] (v2) -- (v4) -- (v6) -- (v7);
+\draw[fsurfn] (v1) -- (v3) -- (v5) -- (v7);
+\draw[edgel] (v1) -- (v2);
+\draw[edgel] (v3) -- (v4);
+\draw[edgel] (v5) -- (v6);
+
+
+\begin{pgfonlayer}{background}
+% cells
+\fill[celln] (p1) -- (p8) -- (p3) -- cycle;
+\fill[celln] (p3) -- (p10) -- (p5) -- cycle;
+\fill[celln] (p5) -- (p12) -- (p7) -- cycle;
+
+\fill[celln] (p1) -- (p14) -- (p8) -- cycle;
+\fill[celln] (p3) -- (p8) -- (p10) -- cycle;
+\fill[celln] (p5) -- (p10) -- (p12) -- cycle;
+
+\fill[cellp] (p2) -- (p4) -- (p9) -- cycle;
+\fill[cellp] (p4) -- (p6) -- (p11) -- cycle;
+\fill[cellp] (p6) -- (p7) -- (p13) -- cycle;
+
+\fill[cellp] (p2) -- (p9) -- (p15) -- cycle;
+\fill[cellp] (p4) -- (p11) -- (p9) -- cycle;
+\fill[cellp] (p6) -- (p13) -- (p11) -- cycle;
+
+% negative
+\draw[cell] (v1) -- (v8) -- (v3) -- (v1);
+\draw[cell] (v3) -- (v10) -- (v5) -- (v3);
+\draw[cell] (v5) -- (v12) -- (v7) -- (v5);
+
+\draw[cell] (v1) -- (v14) -- (v8) -- (v1);
+\draw[cell] (v3) -- (v8) -- (v10) -- (v3);
+\draw[cell] (v5) -- (v10) -- (v12) -- (v5);
+
+% positive
+\draw[cell] (p7) -- (v12) -- (v16) -- (v7);
+\draw[cell] (v2) -- (v4) -- (v9) -- (v2);
+\draw[cell] (v4) -- (v6) -- (v11) -- (v4);
+
+\draw[cell] (v2) -- (v9) -- (v15) -- (v2);
+\draw[cell] (v4) -- (v11) -- (v9) -- (v4);
+\draw[cell] (v6) -- (v13) -- (v11) -- (v6);
+
+% unclassified
+\draw[cell] (v6) -- (v7) -- (v13) -- (v6);
+\draw[cell] (v7) -- (v17) -- (v13) -- (v7);
+\draw[cell] (v7) -- (v16) -- (v17) -- (v7);
+
+\end{pgfonlayer}
+
+\node[classify] at ($ 1/3*(p1)+1/3*(p8)+1/3*(p3) $) {-};
+\node[classify] at ($ 1/3*(p3)+1/3*(p10)+1/3*(p5) $) {-};
+\node[classify] at ($ 1/3*(p5)+1/3*(p12)+1/3*(p7) $) {-};
+
+\node[classify] at ($ 1/3*(p1)+1/3*(p14)+1/3*(p8) $) {-};
+\node[classify] at ($ 1/3*(p3)+1/3*(p8)+1/3*(p10) $) {-};
+\node[classify] at ($ 1/3*(p5)+1/3*(p10)+1/3*(p12) $) {-};
+
+\node[classify] at ($ 1/3*(p2)+1/3*(p4)+1/3*(p9) $) {+};
+\node[classify] at ($ 1/3*(p4)+1/3*(p6)+1/3*(p11) $) {+};
+\node[classify] at ($ 1/3*(p6)+1/3*(p7)+1/3*(p13) $) {+};
+
+\node[classify] at ($ 1/3*(p2)+1/3*(p9)+1/3*(p15) $) {+};
+\node[classify] at ($ 1/3*(p4)+1/3*(p11)+1/3*(p9) $) {+};
+\node[classify] at ($ 1/3*(p6)+1/3*(p13)+1/3*(p11) $) {+};
+
+
+
+\end{tikzpicture}
+\end{document}
diff --git a/doc/presentations/pylith1.7/figs/figcolors.tex b/doc/userguide/boundaryconditions/figs/figcolors.tex
similarity index 100%
copy from doc/presentations/pylith1.7/figs/figcolors.tex
copy to doc/userguide/boundaryconditions/figs/figcolors.tex
diff --git a/doc/userguide/boundaryconditions/figs/quad4cohesivekin.fig b/doc/userguide/boundaryconditions/figs/quad4cohesivekin.fig
deleted file mode 100644
index a98a7d1..0000000
--- a/doc/userguide/boundaryconditions/figs/quad4cohesivekin.fig
+++ /dev/null
@@ -1,122 +0,0 @@
-#FIG 3.2
-Landscape
-Center
-Inches
-Letter
-100.00
-Single
--2
-1200 2
-0 32 #010101
-0 33 #f9f9f9
-0 34 #2f3548
-0 35 #727fa8
-0 36 #fa4040
-0 37 #c40003
-0 38 #fabb69
-0 39 #f17f00
-0 40 #fafa73
-0 41 #e1e100
-0 42 #60c70e
-0 43 #3b7d08
-0 44 #33bbfa
-0 45 #1e6e91
-0 46 #c98cfa
-0 47 #6101af
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 300 1125 37 37 300 1125 337 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 300 2325 37 37 300 2325 337 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1500 2325 37 37 1500 2325 1537 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1500 1125 37 37 1500 1125 1537 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2700 1125 37 37 2700 1125 2737 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2700 2325 37 37 2700 2325 2737 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4200 1125 37 37 4200 1125 4237 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4200 2325 37 37 4200 2325 4237 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6600 1125 37 37 6600 1125 6637 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6600 2325 37 37 6600 2325 6637 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1650 3300 37 37 1650 3300 1687 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1350 3300 37 37 1350 3300 1387 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1350 4200 37 37 1350 4200 1387 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2550 3300 37 37 2550 3300 2587 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2550 4200 37 37 2550 4200 2587 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1650 4200 37 37 1650 4200 1687 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 450 3300 37 37 450 3300 487 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 450 4200 37 37 450 4200 487 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4050 3300 37 37 4050 3300 4087 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4050 4200 37 37 4050 4200 4087 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4950 4200 37 37 4950 4200 4987 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4950 3300 37 37 4950 3300 4987 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5850 3300 37 37 5850 3300 5887 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5850 4200 37 37 5850 4200 5887 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6750 4200 37 37 6750 4200 6787 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6750 3300 37 37 6750 3300 6787 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 1125 37 37 5400 1125 5437 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 2325 37 37 5400 2325 5437 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5362 2325 37 37 5362 2325 5399 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5438 2325 37 37 5438 2325 5475 2325
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5362 1125 37 37 5362 1125 5399 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5438 1125 37 37 5438 1125 5475 1125
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5250 3300 37 37 5250 3300 5287 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5550 3300 37 37 5550 3300 5587 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5250 4200 37 37 5250 4200 5287 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5550 4200 37 37 5550 4200 5587 4200
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 3300 37 37 5400 3300 5437 3300
-1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 4200 37 37 5400 4200 5437 4200
-2 3 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 5
- 300 1125 300 2325 2700 2325 2700 1125 300 1125
-2 1 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 2
- 1500 1125 1500 2325
-2 3 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 5
- 4200 1125 4200 2325 6600 2325 6600 1125 4200 1125
-2 1 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 2
- 5400 1125 5400 2325
-2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
- 450 3300 1350 3300 1350 4200 450 4200 450 3300
-2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
- 1650 3300 2550 3300 2550 4200 1650 4200 1650 3300
-2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
- 4050 3300 4950 3300 4950 4200 4050 4200 4050 3300
-2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
- 5850 3300 6750 3300 6750 4200 5850 4200 5850 3300
-2 2 1 1 45 7 50 -1 -1 4.000 0 0 -1 0 0 5
- 5250 3300 5550 3300 5550 4200 5250 4200 5250 3300
-4 1 0 40 -1 18 12 0.0000 4 195 1245 1500 600 Original Mesh\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 4050 2550 0\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 6675 2550 4\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 6675 1050 5\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 4050 1050 1\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 2775 1050 5\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 150 1050 1\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 150 2550 0\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 2775 2550 4\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 300 3225 1\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 1350 4425 2\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 1650 4425 2\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 300 4425 0\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 1350 3225 3\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 1650 3225 3\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 2625 3225 5\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 2625 4425 4\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 6825 3225 5\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 3975 4425 0\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 3900 3225 1\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 6825 4425 4\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 1500 2550 2\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 1500 1050 3\001
-4 1 0 40 -1 18 12 0.0000 4 180 2160 3375 4800 Exploded view of meshes\001
-4 1 0 40 -1 18 12 0.0000 4 195 4560 5400 600 Mesh with Cohesive Cell using Lagrange Multipliers\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 5250 1275 3\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 5475 1275 7\001
-4 0 0 40 -1 18 10 0.0000 4 105 90 5250 2250 2\001
-4 0 0 40 -1 18 10 0.0000 4 105 75 5475 2250 6\001
-4 1 0 40 -1 18 10 0.0000 4 105 90 5400 4425 8\001
-4 1 0 40 -1 18 10 0.0000 4 105 75 5400 3225 9\001
-4 1 0 40 -1 18 10 0.0000 4 105 75 5625 3225 7\001
-4 1 0 40 -1 18 10 0.0000 4 105 90 5175 3225 3\001
-4 1 0 40 -1 18 10 0.0000 4 105 75 5925 3225 7\001
-4 1 0 40 -1 18 10 0.0000 4 105 90 4875 3225 3\001
-4 1 0 40 -1 18 10 0.0000 4 105 75 5925 4425 6\001
-4 1 0 40 -1 18 10 0.0000 4 105 75 5625 4425 6\001
-4 1 0 40 -1 18 10 0.0000 4 105 90 5175 4425 2\001
-4 1 0 40 -1 18 10 0.0000 4 105 90 4875 4425 2\001
-4 1 0 40 -1 18 10 0.0000 4 105 90 5400 2550 8\001
-4 1 0 40 -1 18 10 0.0000 4 105 75 5400 1050 9\001
diff --git a/doc/userguide/boundaryconditions/figs/quad4cohesivekin.pdf b/doc/userguide/boundaryconditions/figs/quad4cohesivekin.pdf
deleted file mode 100644
index 6a434d6..0000000
Binary files a/doc/userguide/boundaryconditions/figs/quad4cohesivekin.pdf and /dev/null differ
More information about the CIG-COMMITS
mailing list