[cig-commits] r20234 - in short/3D/PyLith/branches/v1.7-trunk: . doc/userguide/runpylith

brad at geodynamics.org brad at geodynamics.org
Wed May 30 15:14:46 PDT 2012


Author: brad
Date: 2012-05-30 15:14:45 -0700 (Wed, 30 May 2012)
New Revision: 20234

Modified:
   short/3D/PyLith/branches/v1.7-trunk/TODO
   short/3D/PyLith/branches/v1.7-trunk/doc/userguide/runpylith/runpylith.lyx
Log:
Added section on finite-element integration parameters.

Modified: short/3D/PyLith/branches/v1.7-trunk/TODO
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/TODO	2012-05-30 21:31:29 UTC (rev 20233)
+++ short/3D/PyLith/branches/v1.7-trunk/TODO	2012-05-30 22:14:45 UTC (rev 20234)
@@ -25,13 +25,11 @@
   - PETSc w/CUDA (Add to Section 5.1)
       + Add link to PETSc CUDA http://www.mcs.anl.gov/petsc/documentation/installation.html#CUDA
       + Add table with PETSc options (see solvertest/pc_fieldsplit_gpu.cfg).
-  - FIATSimplex, FIATLagrange
-     cell.shape -> cell.dimension
   - Examples/Tutorials
-    + Table with concepts/features
     + examples/3d/hex8/step20 (dike opening)
     + examples/2d/greensfns
     + examples/2d/subduction/step04.cfg [stress driven afterslip via friction]
+    + Table with concepts/features
 
   - Update Components appendix [DONE?, need to check]
 
@@ -58,7 +56,11 @@
     + fit to yield surface
     + allow tensile yield
     + Derivation of fit to yield surface (insert immediately before table)
+  - FIATSimplex, FIATLagrange
+     cell.shape -> cell.dimension
 
+  Fixed bug in using friction in quasi-static simulation in parallel.
+
 * configure
   
   + Check compatibility of PyLith options with PETSc

Modified: short/3D/PyLith/branches/v1.7-trunk/doc/userguide/runpylith/runpylith.lyx
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/doc/userguide/runpylith/runpylith.lyx	2012-05-30 21:31:29 UTC (rev 20233)
+++ short/3D/PyLith/branches/v1.7-trunk/doc/userguide/runpylith/runpylith.lyx	2012-05-30 22:14:45 UTC (rev 20234)
@@ -1591,6 +1591,68 @@
 \end_layout
 
 \begin_layout Subsection
+Finite-Element Integration Settings
+\end_layout
+
+\begin_layout Standard
+PyLith uses numerical quadrature to evaluate the finite-element integrals
+ for the residual and system Jacobian (see Chapter 
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "cha:Governing-Equations"
+
+\end_inset
+
+).
+ PyLith employs FIAT (finite element automatic tabulator) to compute the
+ basis functions and their derivatives at the quadrature points for various
+ quadrature schemes and cell shapes.
+ The parameters for Lagrange cells (lines, quadrilaterals, hexahedra) are
+ specified using the FIATLagrange object, whereas the parameters for Simplex
+ cells (lines, triangles, tetrahedra) are specified using the FIATSimplex
+ object.
+ Both objects use the same set of parameters and PyLith will setup the basis
+ functions and quadrature scheme appropriately for the two families of cells.
+ The quadrature scheme and basis functions must be set for each material
+ and boundary condition involving finite-element integrations (Dirichlet
+ boundary conditions are constraints and do not involve integrations).
+ Furthermore, the integration schemes can be set independently.
+ The current version of PyLith supports basis functions with linear variations
+ in the field (P1); support for higher order cells will be added in the
+ future.
+ The properties for the FIATLagrange and FIATSimplex objects are
+\end_layout
+
+\begin_layout Description
+dimension Dimension of the cell (0,1,2,3; default is 3).
+\end_layout
+
+\begin_layout Description
+degree Degree of the finite-element cell (default is 1).
+\end_layout
+
+\begin_layout Description
+order Order of quadrature rule (default is degree+1); hardwired to be equal
+ to degree for faults.
+\end_layout
+
+\begin_layout Description
+collocate_quad Collocate quadrature points with vertices (default is False);
+ hardwired to True for faults.
+\end_layout
+
+\begin_layout Standard
+See Section 
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "sec:material:parameters"
+
+\end_inset
+
+ for an example of setting these properties for a material.
+\end_layout
+
+\begin_layout Subsection
 \begin_inset CommandInset label
 LatexCommand label
 name "sec:petsc:options"



More information about the CIG-COMMITS mailing list