[cig-commits] r11489 - short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes

[brad at geodynamics.org]

Author: brad
Date: 2008-03-20 11:55:48 -0700 (Thu, 20 Mar 2008)
New Revision: 11489

Removed:
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/classify_pts.lagrit
Modified:
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit
Log:
Decided improving quality of cells was not important.

Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/classify_pts.lagrit
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/classify_pts.lagrit	2008-03-20 06:20:56 UTC (rev 11488)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/classify_pts.lagrit	2008-03-20 18:55:48 UTC (rev 11489)
@@ -1,40 +0,0 @@
-* -*- LaGriT -*-
-*
-* ----------------------------------------------------------------------
-*
-*                           Brad T. Aagaard
-*                        U.S. Geological Survey
-*
-* <LicenseText>
-*
-* ----------------------------------------------------------------------
-*
-
-* Classify points according to types.
-*
-* Point sets:
-*   pts_interface = points on interface
-*   pts_notinterface = point not on interface
-*   pts_minaspect = all points associated with elements < aspect ratio
-*   pts_intminaspect = all points in interior associated 
-*       with elements < aspect ratio
-* Element sets:
-*   el_minaspect = elements with aspect ratio < value
-
-pset / pts_interface / surface / -interface-
-pset / pts_parents / attribute / itp1 / 1,0,0 / 41 / eq
-pset / pts_all / seq / 1,0,0
-pset / pts_notinterface / not / pts_all,pts_interface,pts_parents
-
-quality / aspect / y
-eltset / el_minaspect / aratio / lt / 0.6
-pset / pts_minaspect / eltset / el_minaspect
-pset / pts_interior / attribute / itp1 / 1,0,0 / 0 / eq
-pset / pts_intminaspect / inter / pts_minaspect,pts_interior
-pset / pts_notinterfaceaspect / union / pts_notinterface,pts_minaspect
-
-* ----------------------------------------------------------------------
-finish
-
-
-* End of file

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit	2008-03-20 06:20:56 UTC (rev 11488)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit	2008-03-20 18:55:48 UTC (rev 11489)
@@ -199,46 +199,6 @@
 resetpts/cell_color
 
 * ----------------------------------------------------------------------
-* Improve element size and quality
-* ----------------------------------------------------------------------
-*
-define / MO_NAME / box
-mode / recon / geom
-
-massage / len_max / len_min / 1.0e-4 / 1.0e-4 / 1,0,0 / nosmooth
-crush_thin_tets / MO_NAME / 0.4
-rmpoint / compress
-resetpts / itp
-
-*smooth / position / laplace
-input / classify_pts.lagrit
-smooth / position / geometry / pset,get,pts_intminaspect
-
-massage / len_max / len_min / 1.0e-4 / 1.0e-4 / 1,0,0 / &
-  nosmooth
-crush_thin_tets / MO_NAME / 0.4
-rmpoint / compress
-resetpts / itp
-
-** Plump most distorted elements
-input / classify_pts.lagrit
-smooth / position / geometry / pset,get,pts_intminaspect
-
-* Cleanup --------------------------------------------------
-eltset / el_minaspect / delete
-pset / pts_minaspect / delete
-pset / pts_interior / delete
-pset / pts_interface / delete
-pset / pts_parents / delete
-pset / pts_all / delete
-pset / pts_notinterface / delete
-
-filter / 1,0,0 / 10.0
-rmpoint / compress
-
-quality / aspect
-*
-* ----------------------------------------------------------------------
 * Remove parent/child hierarchy
 * ----------------------------------------------------------------------
 resetpts
@@ -249,6 +209,12 @@
 * ----------------------------------------------------------------------
 *
 * Identify nodes on fault surface
+*
+* Nodes on fault surface are defined as being the intersection of the
+* points on the hanging wall and the points on the footwall. This is a
+* more robust definition than using the points on the "fault surface"
+* because smoothing, etc of the points can cause points to move
+* epsilon off the "fault surface".
 eltset / el_hangelast / itetclr / eq / 3
 eltset / el_hangvisco / itetclr / eq / 5
 eltset / el_footelast / itetclr / eq / 4
@@ -260,7 +226,7 @@
 pset / tmp / inter / pts_hang,pts_foot
 pset / fault / geom / xyz / pset,get,tmp / &
    -0.001e+3   -0.001e+3 -16.001e+3 / &
-   16.001e+3   16.001e+3  +0.001e+3 /
+   20.001e+3   16.001e+3  +0.001e+3 /
 pset / tmp / delete
 pset / pts_hang / delete
 pset / pts_foot / delete