[cig-commits] r17169 - in short/3D/PyLith/branches/v1.5-stable: libsrc/feassemble playpen/quadratic playpen/quadratic/twotet10 tests_auto/3dnew tests_auto/3dnew/tet10 tests_auto/3dnew/tet4
brad at geodynamics.org
brad at geodynamics.org
Thu Sep 2 17:33:01 PDT 2010
Author: brad
Date: 2010-09-02 17:33:01 -0700 (Thu, 02 Sep 2010)
New Revision: 17169
Added:
short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/
short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation.cfg
short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation.jpg
short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_slip.spatialdb
short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_sliptime.spatialdb
short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/pylithapp.cfg
short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/twotet10.mesh
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axial_disp.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp.cfg
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_gendb.py
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_soln.py
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation.cfg
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_slipLL.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_sliptime.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/geometry.jou
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/matprops.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/mesh.exo
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/mesh.jou
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/testpylith.py
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axial_disp.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp.cfg
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_gendb.py
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_soln.py
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/geometry.jou
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/matprops.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/mesh.exo
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/mesh.jou
short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/testpylith.py
Modified:
short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1D.cc
short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din2D.cc
short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din3D.cc
short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2D.cc
short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2Din3D.cc
short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature3D.cc
Log:
Temporarily reverted to using isoparameteric property of cells for computing quadrature rather than cell geometry due to corners in Sieve BFS ordering being at the end.
Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1D.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1D.cc 2010-09-03 00:31:48 UTC (rev 17168)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1D.cc 2010-09-03 00:33:01 UTC (rev 17169)
@@ -27,7 +27,7 @@
#include <cassert> // USES assert()
-//#define ISOPARAMETRIC
+#define ISOPARAMETRIC
// ----------------------------------------------------------------------
// Constructor
Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din2D.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din2D.cc 2010-09-03 00:31:48 UTC (rev 17168)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din2D.cc 2010-09-03 00:33:01 UTC (rev 17169)
@@ -27,7 +27,7 @@
#include <cassert> // USES assert()
-//#define ISOPARAMETRIC
+#define ISOPARAMETRIC
// ----------------------------------------------------------------------
// Constructor
Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din3D.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din3D.cc 2010-09-03 00:31:48 UTC (rev 17168)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature1Din3D.cc 2010-09-03 00:33:01 UTC (rev 17169)
@@ -27,7 +27,7 @@
#include <cassert> // USES assert()
-//#define ISOPARAMETRIC
+#define ISOPARAMETRIC
// ----------------------------------------------------------------------
// Constructor
Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2D.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2D.cc 2010-09-03 00:31:48 UTC (rev 17168)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2D.cc 2010-09-03 00:33:01 UTC (rev 17169)
@@ -27,7 +27,7 @@
#include <cassert> // USES assert()
-//#define ISOPARAMETRIC
+#define ISOPARAMETRIC
// ----------------------------------------------------------------------
// Constructor
Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2Din3D.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2Din3D.cc 2010-09-03 00:31:48 UTC (rev 17168)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature2Din3D.cc 2010-09-03 00:33:01 UTC (rev 17169)
@@ -29,7 +29,7 @@
#include <cassert> // USES assert()
#include <stdexcept> // USES std::runtime_error()
-//#define ISOPARAMETRIC
+#define ISOPARAMETRIC
// ----------------------------------------------------------------------
// Constructor
Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature3D.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature3D.cc 2010-09-03 00:31:48 UTC (rev 17168)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/Quadrature3D.cc 2010-09-03 00:33:01 UTC (rev 17169)
@@ -27,7 +27,7 @@
#include <cassert> // USES assert()
-//#define ISOPARAMETRIC
+#define ISOPARAMETRIC
// ----------------------------------------------------------------------
// Constructor
Added: short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation.cfg 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,107 @@
+# -*- Python -*-
+
+# The settings in this file (dislocation.cfg) will be read if it is
+# specified on the command line:
+# 'pylith dislocation.cfg'
+
+# If this file is given, the settings in this file will override all
+# others except for parameters specified directly on the command line.
+
+# These settings define an implicit problem using a fault with
+# kinematically-specified slip.
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Journal settings in addition to those given in 'pylithapp.cfg'
+[pylithapp.journal.info]
+quadrature2d = 1
+faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+# Specify the problem settings.
+# This is a time-dependent problem so we use that facility.
+[pylithapp.timedependent]
+
+# We want an implicit formulation.
+formulation = pylith.problems.Implicit
+
+# Set bc to an array of 1 boundary condition: 'bc'.
+bc = [bc]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# Provide information on the boundary conditions.
+[pylithapp.timedependent.bc.bc]
+
+# We use the default spatial database for the Dirichlet BC which sets
+# the displacements to zero for all time.
+
+# We are fixing the 0 (x), 1 (y), and 2 (z) degrees of freedom.
+bc_dof = [0, 1, 2]
+
+# The nodes associated with this boundary condition have the name
+# 'end points' in the mesh file.
+label = end points
+
+# We are assigning the label 'Dirichlet BC' to the database.
+db_initial.label = Dirichlet BC
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# Provide information on the fault (interface).
+[pylithapp.timedependent.interfaces]
+
+# Define fault properties.
+[pylithapp.timedependent.interfaces.fault]
+
+# The nodes associated with this fault have the name 'fault' in the mesh file.
+label = fault
+
+# NOTE: It is possible to assign an ID number to a fault (e.g.,
+# 'id = 10'). Care must be taken when doing this, however, because the
+# assigned ID will become the material ID for the cohesive element.
+# This ID must not conflict with any of the material ID numbers for
+# volume elements. The default ID for a fault is 100. If you have a
+# fault in your mesh you must:
+# 1. If you create your own fault ID, make sure it does not conflict
+# with any of you material ID's.
+# 2. If you use the default fault ID, make sure that none of your
+# material ID's are equal to 100.
+
+# The quadrature for a 3D fault is 2D.
+quadrature.cell.shape = triangle
+quadrature.cell.degree = 2
+
+# Give the spatial databases specifying the fault slip.
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+
+# Database specifying fault slip (0.01 m of left-lateral slip).
+slip.iohandler.filename = dislocation_slip.spatialdb
+
+# Database specifying time at which fault slips (0.0 s).
+slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output.writer]
+filename = dislocation.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = dislocation-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.material.output.writer]
+filename = dislocation-statevars.vtk
Added: short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation.jpg
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation.jpg
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_slip.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_slip.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_slip.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,42 @@
+#SPATIAL.ascii 1
+
+// This database is used to specify the fault slip for the
+// dislocation (kinematic fault) example.
+SimpleDB {
+
+ // There are 3 values specified in the database, corresponding to the
+ // slip values for the left-lateral, reverse-slip, and fault-opening
+ // components.
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+
+ // The fault slip values have units of meters.
+ value-units = m m m
+
+ // The values are specified at one spatial location.
+ num-locs = 1
+
+ // The dimension of the spatial distribution is 0, since the same data
+ // is specified for all points in the set.
+ data-dim = 0
+
+ // The spatial dimension of the database is 3.
+ space-dim = 3
+
+ // We are specifying the data in a Cartesian coordinate system.
+ cs-data = cartesian {
+
+ // Our units are already in meters, so we can just multiply by one.
+ to-meters = 1.0
+
+ // We are using a 3D Cartesian coordinate system.
+ space-dim = 3
+ }
+}
+
+// This is where the data is specified.
+// As described in dislocation.cfg, this database will be used to specify
+// uniform data for a fault (pylith.faults.SingleFault).
+// The entries are:
+// X-Coord, Y-Coord, Z-Coord, left-lateral slip, reverse-slip, fault-opening slip
+0.0 0.0 0.0 0.01 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_sliptime.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_sliptime.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/dislocation_sliptime.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,41 @@
+#SPATIAL.ascii 1
+
+// This database is used to specify the fault slip time for the
+// dislocation (kinematic fault) example.
+SimpleDB {
+
+ // There is one value specified in the database, corresponding to the
+ // time at which fault slip begins.
+ num-values = 1
+ value-names = slip-time
+
+ // The fault slip time has units of meters.
+ value-units = s
+
+ // The value is specified at one spatial location.
+ num-locs = 1
+
+ // The dimension of the spatial distribution is 0, since the same data
+ // is specified for all points in the set.
+ data-dim = 0
+
+ // The spatial dimension of the database is 3.
+ space-dim = 3
+
+ // We are specifying the data in a Cartesian coordinate system.
+ cs-data = cartesian {
+
+ // Our units are already in meters, so we can just multiply by one.
+ to-meters = 1.0
+
+ // We are using a 3D Cartesian coordinate system.
+ space-dim = 3
+ }
+}
+
+// This is where the data is specified.
+// As described in dislocation.cfg, this database will be used to specify
+// uniform data for a fault (pylith.faults.SingleFault).
+// The entries are:
+// X-Coord, Y-Coord, Z-Coord, slip time
+0.0 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/pylithapp.cfg 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,100 @@
+# -*- Python -*-
+
+# The settings in this file (pylithapp.cfg) will be read automatically
+# by pylith, as long as the file is placed in the run directory.
+
+# The settings in this file will override any settings in:
+# PREFIX/etc/pylithapp.cfg
+# $HOME/.pyre/pylithapp/pylithapp.cfg
+
+# The settings in this file will be overridden by any .cfg file given
+# on the command line or by any command line settings.
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# The settings below turn on journal info for the specified components.
+# If you want less output to stdout, you can turn these off.
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshioascii = 1
+homogeneous = 1
+implicitelasticity = 1
+quadrature2d = 1
+fiatsimplex = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+# The settings below control the mesh generation (importing mesh info).
+# Turn on debugging output for mesh generation.
+[pylithapp.mesh_generator]
+debug = 1
+
+# This component specification means we are using PyLith ASCII format,
+# and we then specify the filename and number of space dimensions for
+# the mesh.
+[pylithapp.mesh_generator.reader]
+filename = twotet10.mesh
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+# Specify the problem settings.
+# This is a time-dependent problem, so we select this as our problem type.
+# We select a total time of 0 sec, and a time step size of 1 sec, so we
+# are performing a single time step.
+# The spatial dimension for this problem is 2.
+# For an implicit formulation (using implicit.cfg), we will perform 1
+# implicit time step from t = -1.0 to t = 0.0 (elastic solution step).
+[pylithapp.timedependent]
+dimension = 3
+normalizer.length_scale = 1.0*m
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+dt = 1.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+# Specify the material information for the problem.
+# The material type is 3D isotropic elastic.
+[pylithapp.timedependent.materials]
+material = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.material]
+
+# We give a label of 'Elastic material' to this material.
+label = Elastic material
+
+# The cells associated with this material are given a material ID of 1
+# in the mesh file.
+id = 1
+
+# The properties for this material are given in the spatial database file
+# 'matprops.spatialdb'.
+db_properties.iohandler.filename = matprops.spatialdb
+
+# We are doing 3D quadrature for a tetrahedron.
+quadrature.cell.shape = tetrahedron
+quadrature.cell.degree = 2
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+# We are using all of the default settings for PETSc except for
+# specifying the block Jacobi preconditioner. Additional PETSc
+# command-line arguments may be found in the PETSc documentation.
+[pylithapp.petsc]
+pc_type = jacobi
+
+# start_in_debugger = true
+# debugger_timeout = 100
Added: short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/twotet10.mesh
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/twotet10.mesh (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/playpen/quadratic/twotet10/twotet10.mesh 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,62 @@
+// Global mesh object.
+// This defines a mesh composed of two tetrahedral elements.
+mesh = {
+
+ dimension = 3
+ use-index-zero = true
+ vertices = {
+ dimension = 3
+ count = 14
+ coordinates = {
+ 0 -1.0 0.0 0.0
+ 1 -0.5 -0.5 0.0
+ 2 -0.5 0.5 0.0
+ 3 -0.5 0.0 0.5
+ 4 0.0 -1.0 0.0
+ 5 0.0 0.0 0.0
+ 6 0.0 1.0 0.0
+ 7 0.0 -0.5 0.5
+ 8 0.0 0.5 0.5
+ 9 0.0 0.0 1.0
+ 10 0.5 -0.5 0.0
+ 11 0.5 0.5 0.0
+ 12 0.5 0.0 0.5
+ 13 1.0 0.0 0.0
+ }
+ }
+
+ cells = {
+
+ count = 2
+ num-corners = 10
+ simplices = {
+ 0 1 5 2 3 7 8 0 4 6 9
+ 1 10 11 5 7 12 8 4 13 6 9
+ }
+
+ material-ids = {
+ 0 1
+ 1 1
+ }
+ }
+
+ group = {
+ name = fault
+ type = vertices
+ count = 6
+ indices = {
+ 4 5 6 7 8 9
+ }
+ }
+
+ group = {
+ name = end points
+ type = vertices
+ count = 2
+ indices = {
+ 0
+ 13
+ }
+ }
+
+}
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/Makefile.am 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,70 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+TESTS = testpylith.py
+
+check_SCRIPTS = testpylith.py
+
+dist_noinst_PYTHON = \
+ TestHex8.py \
+ TestLgDeformRigidBody.py \
+ rigidbody_soln.py \
+ rigidbody_gendb.py \
+ TestLgDeformTraction.py \
+ lgdeformtraction_soln.py \
+ testpylith.py
+
+dist_noinst_DATA = \
+ geometry.jou \
+ mesh.jou \
+ mesh.exo \
+ matprops.spatialdb \
+ lgdeformrigidbody.cfg \
+ lgdeformtraction.cfg
+
+noinst_TMP = \
+ rigidbody_disp.spatialdb \
+ lgdeformrigidbody_t0000000.vtk \
+ lgdeformrigidbody-elastic_info.vtk \
+ lgdeformrigidbody-elastic_t0000000.vtk \
+ lgdeformrigidbody-viscoelastic_info.vtk \
+ lgdeformrigidbody-viscoelastic_t0000000.vtk \
+ lgdeformtraction_t0000000.vtk \
+ lgdeformtraction-elastic_info.vtk \
+ lgdeformtraction-elastic_t0000000.vtk \
+ lgdeformtraction-viscoelastic_info.vtk \
+ lgdeformtraction-viscoelastic_t0000000.vtk
+
+
+TESTS_ENVIRONMENT = $(PYTHON)
+
+
+# 'export' the input files by performing a mock install
+export_datadir = $(top_builddir)/tests/3dnew/hex8
+export-data: $(dist_noinst_PYTHON) $(dist_noinst_DATA)
+ for f in $(dist_noinst_PYTHON) $(dist_noinst_DATA); do $(install_sh_DATA) $(srcdir)/$$f $(export_datadir); done
+
+clean-data:
+ if [ "X$(top_srcdir)" != "X$(top_builddir)" ]; then for f in $(dist_noinst_PYTHON) $(dist_noinst_DATA) $(noinst_TMP); do $(RM) $(RM_FLAGS) $(export_datadir)/$$f; done; fi
+
+
+BUILT_SOURCES = export-data
+clean-local: clean-data
+CLEANFILES = *.pyc
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axial_disp.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axial_disp.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axial_disp.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,580 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = displacement-x displacement-y displacement-z
+ value-units = m m m
+ num-locs = 567
+ data-dim = 3
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1
+ space-dim = 3
+}
+}
+ -4.000000e+03 -4.000000e+03 -6.000000e+03 -7.037037e-01 1.851852e-01 2.777778e-01
+ -3.000000e+03 -4.000000e+03 -6.000000e+03 -5.277778e-01 1.851852e-01 2.777778e-01
+ -2.000000e+03 -4.000000e+03 -6.000000e+03 -3.518519e-01 1.851852e-01 2.777778e-01
+ -1.000000e+03 -4.000000e+03 -6.000000e+03 -1.759259e-01 1.851852e-01 2.777778e-01
+ 0.000000e+00 -4.000000e+03 -6.000000e+03 0.000000e+00 1.851852e-01 2.777778e-01
+ 1.000000e+03 -4.000000e+03 -6.000000e+03 1.759259e-01 1.851852e-01 2.777778e-01
+ 2.000000e+03 -4.000000e+03 -6.000000e+03 3.518519e-01 1.851852e-01 2.777778e-01
+ 3.000000e+03 -4.000000e+03 -6.000000e+03 5.277778e-01 1.851852e-01 2.777778e-01
+ 4.000000e+03 -4.000000e+03 -6.000000e+03 7.037037e-01 1.851852e-01 2.777778e-01
+ -4.000000e+03 -3.000000e+03 -6.000000e+03 -7.037037e-01 1.388889e-01 2.777778e-01
+ -3.000000e+03 -3.000000e+03 -6.000000e+03 -5.277778e-01 1.388889e-01 2.777778e-01
+ -2.000000e+03 -3.000000e+03 -6.000000e+03 -3.518519e-01 1.388889e-01 2.777778e-01
+ -1.000000e+03 -3.000000e+03 -6.000000e+03 -1.759259e-01 1.388889e-01 2.777778e-01
+ 0.000000e+00 -3.000000e+03 -6.000000e+03 0.000000e+00 1.388889e-01 2.777778e-01
+ 1.000000e+03 -3.000000e+03 -6.000000e+03 1.759259e-01 1.388889e-01 2.777778e-01
+ 2.000000e+03 -3.000000e+03 -6.000000e+03 3.518519e-01 1.388889e-01 2.777778e-01
+ 3.000000e+03 -3.000000e+03 -6.000000e+03 5.277778e-01 1.388889e-01 2.777778e-01
+ 4.000000e+03 -3.000000e+03 -6.000000e+03 7.037037e-01 1.388889e-01 2.777778e-01
+ -4.000000e+03 -2.000000e+03 -6.000000e+03 -7.037037e-01 9.259259e-02 2.777778e-01
+ -3.000000e+03 -2.000000e+03 -6.000000e+03 -5.277778e-01 9.259259e-02 2.777778e-01
+ -2.000000e+03 -2.000000e+03 -6.000000e+03 -3.518519e-01 9.259259e-02 2.777778e-01
+ -1.000000e+03 -2.000000e+03 -6.000000e+03 -1.759259e-01 9.259259e-02 2.777778e-01
+ 0.000000e+00 -2.000000e+03 -6.000000e+03 0.000000e+00 9.259259e-02 2.777778e-01
+ 1.000000e+03 -2.000000e+03 -6.000000e+03 1.759259e-01 9.259259e-02 2.777778e-01
+ 2.000000e+03 -2.000000e+03 -6.000000e+03 3.518519e-01 9.259259e-02 2.777778e-01
+ 3.000000e+03 -2.000000e+03 -6.000000e+03 5.277778e-01 9.259259e-02 2.777778e-01
+ 4.000000e+03 -2.000000e+03 -6.000000e+03 7.037037e-01 9.259259e-02 2.777778e-01
+ -4.000000e+03 -1.000000e+03 -6.000000e+03 -7.037037e-01 4.629630e-02 2.777778e-01
+ -3.000000e+03 -1.000000e+03 -6.000000e+03 -5.277778e-01 4.629630e-02 2.777778e-01
+ -2.000000e+03 -1.000000e+03 -6.000000e+03 -3.518519e-01 4.629630e-02 2.777778e-01
+ -1.000000e+03 -1.000000e+03 -6.000000e+03 -1.759259e-01 4.629630e-02 2.777778e-01
+ 0.000000e+00 -1.000000e+03 -6.000000e+03 0.000000e+00 4.629630e-02 2.777778e-01
+ 1.000000e+03 -1.000000e+03 -6.000000e+03 1.759259e-01 4.629630e-02 2.777778e-01
+ 2.000000e+03 -1.000000e+03 -6.000000e+03 3.518519e-01 4.629630e-02 2.777778e-01
+ 3.000000e+03 -1.000000e+03 -6.000000e+03 5.277778e-01 4.629630e-02 2.777778e-01
+ 4.000000e+03 -1.000000e+03 -6.000000e+03 7.037037e-01 4.629630e-02 2.777778e-01
+ -4.000000e+03 0.000000e+00 -6.000000e+03 -7.037037e-01 -0.000000e+00 2.777778e-01
+ -3.000000e+03 0.000000e+00 -6.000000e+03 -5.277778e-01 -0.000000e+00 2.777778e-01
+ -2.000000e+03 0.000000e+00 -6.000000e+03 -3.518519e-01 -0.000000e+00 2.777778e-01
+ -1.000000e+03 0.000000e+00 -6.000000e+03 -1.759259e-01 -0.000000e+00 2.777778e-01
+ 0.000000e+00 0.000000e+00 -6.000000e+03 0.000000e+00 0.000000e+00 2.777778e-01
+ 1.000000e+03 0.000000e+00 -6.000000e+03 1.759259e-01 0.000000e+00 2.777778e-01
+ 2.000000e+03 0.000000e+00 -6.000000e+03 3.518519e-01 0.000000e+00 2.777778e-01
+ 3.000000e+03 0.000000e+00 -6.000000e+03 5.277778e-01 0.000000e+00 2.777778e-01
+ 4.000000e+03 0.000000e+00 -6.000000e+03 7.037037e-01 0.000000e+00 2.777778e-01
+ -4.000000e+03 1.000000e+03 -6.000000e+03 -7.037037e-01 -4.629630e-02 2.777778e-01
+ -3.000000e+03 1.000000e+03 -6.000000e+03 -5.277778e-01 -4.629630e-02 2.777778e-01
+ -2.000000e+03 1.000000e+03 -6.000000e+03 -3.518519e-01 -4.629630e-02 2.777778e-01
+ -1.000000e+03 1.000000e+03 -6.000000e+03 -1.759259e-01 -4.629630e-02 2.777778e-01
+ 0.000000e+00 1.000000e+03 -6.000000e+03 0.000000e+00 -4.629630e-02 2.777778e-01
+ 1.000000e+03 1.000000e+03 -6.000000e+03 1.759259e-01 -4.629630e-02 2.777778e-01
+ 2.000000e+03 1.000000e+03 -6.000000e+03 3.518519e-01 -4.629630e-02 2.777778e-01
+ 3.000000e+03 1.000000e+03 -6.000000e+03 5.277778e-01 -4.629630e-02 2.777778e-01
+ 4.000000e+03 1.000000e+03 -6.000000e+03 7.037037e-01 -4.629630e-02 2.777778e-01
+ -4.000000e+03 2.000000e+03 -6.000000e+03 -7.037037e-01 -9.259259e-02 2.777778e-01
+ -3.000000e+03 2.000000e+03 -6.000000e+03 -5.277778e-01 -9.259259e-02 2.777778e-01
+ -2.000000e+03 2.000000e+03 -6.000000e+03 -3.518519e-01 -9.259259e-02 2.777778e-01
+ -1.000000e+03 2.000000e+03 -6.000000e+03 -1.759259e-01 -9.259259e-02 2.777778e-01
+ 0.000000e+00 2.000000e+03 -6.000000e+03 0.000000e+00 -9.259259e-02 2.777778e-01
+ 1.000000e+03 2.000000e+03 -6.000000e+03 1.759259e-01 -9.259259e-02 2.777778e-01
+ 2.000000e+03 2.000000e+03 -6.000000e+03 3.518519e-01 -9.259259e-02 2.777778e-01
+ 3.000000e+03 2.000000e+03 -6.000000e+03 5.277778e-01 -9.259259e-02 2.777778e-01
+ 4.000000e+03 2.000000e+03 -6.000000e+03 7.037037e-01 -9.259259e-02 2.777778e-01
+ -4.000000e+03 3.000000e+03 -6.000000e+03 -7.037037e-01 -1.388889e-01 2.777778e-01
+ -3.000000e+03 3.000000e+03 -6.000000e+03 -5.277778e-01 -1.388889e-01 2.777778e-01
+ -2.000000e+03 3.000000e+03 -6.000000e+03 -3.518519e-01 -1.388889e-01 2.777778e-01
+ -1.000000e+03 3.000000e+03 -6.000000e+03 -1.759259e-01 -1.388889e-01 2.777778e-01
+ 0.000000e+00 3.000000e+03 -6.000000e+03 0.000000e+00 -1.388889e-01 2.777778e-01
+ 1.000000e+03 3.000000e+03 -6.000000e+03 1.759259e-01 -1.388889e-01 2.777778e-01
+ 2.000000e+03 3.000000e+03 -6.000000e+03 3.518519e-01 -1.388889e-01 2.777778e-01
+ 3.000000e+03 3.000000e+03 -6.000000e+03 5.277778e-01 -1.388889e-01 2.777778e-01
+ 4.000000e+03 3.000000e+03 -6.000000e+03 7.037037e-01 -1.388889e-01 2.777778e-01
+ -4.000000e+03 4.000000e+03 -6.000000e+03 -7.037037e-01 -1.851852e-01 2.777778e-01
+ -3.000000e+03 4.000000e+03 -6.000000e+03 -5.277778e-01 -1.851852e-01 2.777778e-01
+ -2.000000e+03 4.000000e+03 -6.000000e+03 -3.518519e-01 -1.851852e-01 2.777778e-01
+ -1.000000e+03 4.000000e+03 -6.000000e+03 -1.759259e-01 -1.851852e-01 2.777778e-01
+ 0.000000e+00 4.000000e+03 -6.000000e+03 0.000000e+00 -1.851852e-01 2.777778e-01
+ 1.000000e+03 4.000000e+03 -6.000000e+03 1.759259e-01 -1.851852e-01 2.777778e-01
+ 2.000000e+03 4.000000e+03 -6.000000e+03 3.518519e-01 -1.851852e-01 2.777778e-01
+ 3.000000e+03 4.000000e+03 -6.000000e+03 5.277778e-01 -1.851852e-01 2.777778e-01
+ 4.000000e+03 4.000000e+03 -6.000000e+03 7.037037e-01 -1.851852e-01 2.777778e-01
+ -4.000000e+03 -4.000000e+03 -5.000000e+03 -7.037037e-01 1.851852e-01 2.314815e-01
+ -3.000000e+03 -4.000000e+03 -5.000000e+03 -5.277778e-01 1.851852e-01 2.314815e-01
+ -2.000000e+03 -4.000000e+03 -5.000000e+03 -3.518519e-01 1.851852e-01 2.314815e-01
+ -1.000000e+03 -4.000000e+03 -5.000000e+03 -1.759259e-01 1.851852e-01 2.314815e-01
+ 0.000000e+00 -4.000000e+03 -5.000000e+03 0.000000e+00 1.851852e-01 2.314815e-01
+ 1.000000e+03 -4.000000e+03 -5.000000e+03 1.759259e-01 1.851852e-01 2.314815e-01
+ 2.000000e+03 -4.000000e+03 -5.000000e+03 3.518519e-01 1.851852e-01 2.314815e-01
+ 3.000000e+03 -4.000000e+03 -5.000000e+03 5.277778e-01 1.851852e-01 2.314815e-01
+ 4.000000e+03 -4.000000e+03 -5.000000e+03 7.037037e-01 1.851852e-01 2.314815e-01
+ -4.000000e+03 -3.000000e+03 -5.000000e+03 -7.037037e-01 1.388889e-01 2.314815e-01
+ -3.000000e+03 -3.000000e+03 -5.000000e+03 -5.277778e-01 1.388889e-01 2.314815e-01
+ -2.000000e+03 -3.000000e+03 -5.000000e+03 -3.518519e-01 1.388889e-01 2.314815e-01
+ -1.000000e+03 -3.000000e+03 -5.000000e+03 -1.759259e-01 1.388889e-01 2.314815e-01
+ 0.000000e+00 -3.000000e+03 -5.000000e+03 0.000000e+00 1.388889e-01 2.314815e-01
+ 1.000000e+03 -3.000000e+03 -5.000000e+03 1.759259e-01 1.388889e-01 2.314815e-01
+ 2.000000e+03 -3.000000e+03 -5.000000e+03 3.518519e-01 1.388889e-01 2.314815e-01
+ 3.000000e+03 -3.000000e+03 -5.000000e+03 5.277778e-01 1.388889e-01 2.314815e-01
+ 4.000000e+03 -3.000000e+03 -5.000000e+03 7.037037e-01 1.388889e-01 2.314815e-01
+ -4.000000e+03 -2.000000e+03 -5.000000e+03 -7.037037e-01 9.259259e-02 2.314815e-01
+ -3.000000e+03 -2.000000e+03 -5.000000e+03 -5.277778e-01 9.259259e-02 2.314815e-01
+ -2.000000e+03 -2.000000e+03 -5.000000e+03 -3.518519e-01 9.259259e-02 2.314815e-01
+ -1.000000e+03 -2.000000e+03 -5.000000e+03 -1.759259e-01 9.259259e-02 2.314815e-01
+ 0.000000e+00 -2.000000e+03 -5.000000e+03 0.000000e+00 9.259259e-02 2.314815e-01
+ 1.000000e+03 -2.000000e+03 -5.000000e+03 1.759259e-01 9.259259e-02 2.314815e-01
+ 2.000000e+03 -2.000000e+03 -5.000000e+03 3.518519e-01 9.259259e-02 2.314815e-01
+ 3.000000e+03 -2.000000e+03 -5.000000e+03 5.277778e-01 9.259259e-02 2.314815e-01
+ 4.000000e+03 -2.000000e+03 -5.000000e+03 7.037037e-01 9.259259e-02 2.314815e-01
+ -4.000000e+03 -1.000000e+03 -5.000000e+03 -7.037037e-01 4.629630e-02 2.314815e-01
+ -3.000000e+03 -1.000000e+03 -5.000000e+03 -5.277778e-01 4.629630e-02 2.314815e-01
+ -2.000000e+03 -1.000000e+03 -5.000000e+03 -3.518519e-01 4.629630e-02 2.314815e-01
+ -1.000000e+03 -1.000000e+03 -5.000000e+03 -1.759259e-01 4.629630e-02 2.314815e-01
+ 0.000000e+00 -1.000000e+03 -5.000000e+03 0.000000e+00 4.629630e-02 2.314815e-01
+ 1.000000e+03 -1.000000e+03 -5.000000e+03 1.759259e-01 4.629630e-02 2.314815e-01
+ 2.000000e+03 -1.000000e+03 -5.000000e+03 3.518519e-01 4.629630e-02 2.314815e-01
+ 3.000000e+03 -1.000000e+03 -5.000000e+03 5.277778e-01 4.629630e-02 2.314815e-01
+ 4.000000e+03 -1.000000e+03 -5.000000e+03 7.037037e-01 4.629630e-02 2.314815e-01
+ -4.000000e+03 0.000000e+00 -5.000000e+03 -7.037037e-01 -0.000000e+00 2.314815e-01
+ -3.000000e+03 0.000000e+00 -5.000000e+03 -5.277778e-01 -0.000000e+00 2.314815e-01
+ -2.000000e+03 0.000000e+00 -5.000000e+03 -3.518519e-01 -0.000000e+00 2.314815e-01
+ -1.000000e+03 0.000000e+00 -5.000000e+03 -1.759259e-01 -0.000000e+00 2.314815e-01
+ 0.000000e+00 0.000000e+00 -5.000000e+03 0.000000e+00 0.000000e+00 2.314815e-01
+ 1.000000e+03 0.000000e+00 -5.000000e+03 1.759259e-01 0.000000e+00 2.314815e-01
+ 2.000000e+03 0.000000e+00 -5.000000e+03 3.518519e-01 0.000000e+00 2.314815e-01
+ 3.000000e+03 0.000000e+00 -5.000000e+03 5.277778e-01 0.000000e+00 2.314815e-01
+ 4.000000e+03 0.000000e+00 -5.000000e+03 7.037037e-01 0.000000e+00 2.314815e-01
+ -4.000000e+03 1.000000e+03 -5.000000e+03 -7.037037e-01 -4.629630e-02 2.314815e-01
+ -3.000000e+03 1.000000e+03 -5.000000e+03 -5.277778e-01 -4.629630e-02 2.314815e-01
+ -2.000000e+03 1.000000e+03 -5.000000e+03 -3.518519e-01 -4.629630e-02 2.314815e-01
+ -1.000000e+03 1.000000e+03 -5.000000e+03 -1.759259e-01 -4.629630e-02 2.314815e-01
+ 0.000000e+00 1.000000e+03 -5.000000e+03 0.000000e+00 -4.629630e-02 2.314815e-01
+ 1.000000e+03 1.000000e+03 -5.000000e+03 1.759259e-01 -4.629630e-02 2.314815e-01
+ 2.000000e+03 1.000000e+03 -5.000000e+03 3.518519e-01 -4.629630e-02 2.314815e-01
+ 3.000000e+03 1.000000e+03 -5.000000e+03 5.277778e-01 -4.629630e-02 2.314815e-01
+ 4.000000e+03 1.000000e+03 -5.000000e+03 7.037037e-01 -4.629630e-02 2.314815e-01
+ -4.000000e+03 2.000000e+03 -5.000000e+03 -7.037037e-01 -9.259259e-02 2.314815e-01
+ -3.000000e+03 2.000000e+03 -5.000000e+03 -5.277778e-01 -9.259259e-02 2.314815e-01
+ -2.000000e+03 2.000000e+03 -5.000000e+03 -3.518519e-01 -9.259259e-02 2.314815e-01
+ -1.000000e+03 2.000000e+03 -5.000000e+03 -1.759259e-01 -9.259259e-02 2.314815e-01
+ 0.000000e+00 2.000000e+03 -5.000000e+03 0.000000e+00 -9.259259e-02 2.314815e-01
+ 1.000000e+03 2.000000e+03 -5.000000e+03 1.759259e-01 -9.259259e-02 2.314815e-01
+ 2.000000e+03 2.000000e+03 -5.000000e+03 3.518519e-01 -9.259259e-02 2.314815e-01
+ 3.000000e+03 2.000000e+03 -5.000000e+03 5.277778e-01 -9.259259e-02 2.314815e-01
+ 4.000000e+03 2.000000e+03 -5.000000e+03 7.037037e-01 -9.259259e-02 2.314815e-01
+ -4.000000e+03 3.000000e+03 -5.000000e+03 -7.037037e-01 -1.388889e-01 2.314815e-01
+ -3.000000e+03 3.000000e+03 -5.000000e+03 -5.277778e-01 -1.388889e-01 2.314815e-01
+ -2.000000e+03 3.000000e+03 -5.000000e+03 -3.518519e-01 -1.388889e-01 2.314815e-01
+ -1.000000e+03 3.000000e+03 -5.000000e+03 -1.759259e-01 -1.388889e-01 2.314815e-01
+ 0.000000e+00 3.000000e+03 -5.000000e+03 0.000000e+00 -1.388889e-01 2.314815e-01
+ 1.000000e+03 3.000000e+03 -5.000000e+03 1.759259e-01 -1.388889e-01 2.314815e-01
+ 2.000000e+03 3.000000e+03 -5.000000e+03 3.518519e-01 -1.388889e-01 2.314815e-01
+ 3.000000e+03 3.000000e+03 -5.000000e+03 5.277778e-01 -1.388889e-01 2.314815e-01
+ 4.000000e+03 3.000000e+03 -5.000000e+03 7.037037e-01 -1.388889e-01 2.314815e-01
+ -4.000000e+03 4.000000e+03 -5.000000e+03 -7.037037e-01 -1.851852e-01 2.314815e-01
+ -3.000000e+03 4.000000e+03 -5.000000e+03 -5.277778e-01 -1.851852e-01 2.314815e-01
+ -2.000000e+03 4.000000e+03 -5.000000e+03 -3.518519e-01 -1.851852e-01 2.314815e-01
+ -1.000000e+03 4.000000e+03 -5.000000e+03 -1.759259e-01 -1.851852e-01 2.314815e-01
+ 0.000000e+00 4.000000e+03 -5.000000e+03 0.000000e+00 -1.851852e-01 2.314815e-01
+ 1.000000e+03 4.000000e+03 -5.000000e+03 1.759259e-01 -1.851852e-01 2.314815e-01
+ 2.000000e+03 4.000000e+03 -5.000000e+03 3.518519e-01 -1.851852e-01 2.314815e-01
+ 3.000000e+03 4.000000e+03 -5.000000e+03 5.277778e-01 -1.851852e-01 2.314815e-01
+ 4.000000e+03 4.000000e+03 -5.000000e+03 7.037037e-01 -1.851852e-01 2.314815e-01
+ -4.000000e+03 -4.000000e+03 -4.000000e+03 -7.037037e-01 1.851852e-01 1.851852e-01
+ -3.000000e+03 -4.000000e+03 -4.000000e+03 -5.277778e-01 1.851852e-01 1.851852e-01
+ -2.000000e+03 -4.000000e+03 -4.000000e+03 -3.518519e-01 1.851852e-01 1.851852e-01
+ -1.000000e+03 -4.000000e+03 -4.000000e+03 -1.759259e-01 1.851852e-01 1.851852e-01
+ 0.000000e+00 -4.000000e+03 -4.000000e+03 0.000000e+00 1.851852e-01 1.851852e-01
+ 1.000000e+03 -4.000000e+03 -4.000000e+03 1.759259e-01 1.851852e-01 1.851852e-01
+ 2.000000e+03 -4.000000e+03 -4.000000e+03 3.518519e-01 1.851852e-01 1.851852e-01
+ 3.000000e+03 -4.000000e+03 -4.000000e+03 5.277778e-01 1.851852e-01 1.851852e-01
+ 4.000000e+03 -4.000000e+03 -4.000000e+03 7.037037e-01 1.851852e-01 1.851852e-01
+ -4.000000e+03 -3.000000e+03 -4.000000e+03 -7.037037e-01 1.388889e-01 1.851852e-01
+ -3.000000e+03 -3.000000e+03 -4.000000e+03 -5.277778e-01 1.388889e-01 1.851852e-01
+ -2.000000e+03 -3.000000e+03 -4.000000e+03 -3.518519e-01 1.388889e-01 1.851852e-01
+ -1.000000e+03 -3.000000e+03 -4.000000e+03 -1.759259e-01 1.388889e-01 1.851852e-01
+ 0.000000e+00 -3.000000e+03 -4.000000e+03 0.000000e+00 1.388889e-01 1.851852e-01
+ 1.000000e+03 -3.000000e+03 -4.000000e+03 1.759259e-01 1.388889e-01 1.851852e-01
+ 2.000000e+03 -3.000000e+03 -4.000000e+03 3.518519e-01 1.388889e-01 1.851852e-01
+ 3.000000e+03 -3.000000e+03 -4.000000e+03 5.277778e-01 1.388889e-01 1.851852e-01
+ 4.000000e+03 -3.000000e+03 -4.000000e+03 7.037037e-01 1.388889e-01 1.851852e-01
+ -4.000000e+03 -2.000000e+03 -4.000000e+03 -7.037037e-01 9.259259e-02 1.851852e-01
+ -3.000000e+03 -2.000000e+03 -4.000000e+03 -5.277778e-01 9.259259e-02 1.851852e-01
+ -2.000000e+03 -2.000000e+03 -4.000000e+03 -3.518519e-01 9.259259e-02 1.851852e-01
+ -1.000000e+03 -2.000000e+03 -4.000000e+03 -1.759259e-01 9.259259e-02 1.851852e-01
+ 0.000000e+00 -2.000000e+03 -4.000000e+03 0.000000e+00 9.259259e-02 1.851852e-01
+ 1.000000e+03 -2.000000e+03 -4.000000e+03 1.759259e-01 9.259259e-02 1.851852e-01
+ 2.000000e+03 -2.000000e+03 -4.000000e+03 3.518519e-01 9.259259e-02 1.851852e-01
+ 3.000000e+03 -2.000000e+03 -4.000000e+03 5.277778e-01 9.259259e-02 1.851852e-01
+ 4.000000e+03 -2.000000e+03 -4.000000e+03 7.037037e-01 9.259259e-02 1.851852e-01
+ -4.000000e+03 -1.000000e+03 -4.000000e+03 -7.037037e-01 4.629630e-02 1.851852e-01
+ -3.000000e+03 -1.000000e+03 -4.000000e+03 -5.277778e-01 4.629630e-02 1.851852e-01
+ -2.000000e+03 -1.000000e+03 -4.000000e+03 -3.518519e-01 4.629630e-02 1.851852e-01
+ -1.000000e+03 -1.000000e+03 -4.000000e+03 -1.759259e-01 4.629630e-02 1.851852e-01
+ 0.000000e+00 -1.000000e+03 -4.000000e+03 0.000000e+00 4.629630e-02 1.851852e-01
+ 1.000000e+03 -1.000000e+03 -4.000000e+03 1.759259e-01 4.629630e-02 1.851852e-01
+ 2.000000e+03 -1.000000e+03 -4.000000e+03 3.518519e-01 4.629630e-02 1.851852e-01
+ 3.000000e+03 -1.000000e+03 -4.000000e+03 5.277778e-01 4.629630e-02 1.851852e-01
+ 4.000000e+03 -1.000000e+03 -4.000000e+03 7.037037e-01 4.629630e-02 1.851852e-01
+ -4.000000e+03 0.000000e+00 -4.000000e+03 -7.037037e-01 -0.000000e+00 1.851852e-01
+ -3.000000e+03 0.000000e+00 -4.000000e+03 -5.277778e-01 -0.000000e+00 1.851852e-01
+ -2.000000e+03 0.000000e+00 -4.000000e+03 -3.518519e-01 -0.000000e+00 1.851852e-01
+ -1.000000e+03 0.000000e+00 -4.000000e+03 -1.759259e-01 -0.000000e+00 1.851852e-01
+ 0.000000e+00 0.000000e+00 -4.000000e+03 0.000000e+00 0.000000e+00 1.851852e-01
+ 1.000000e+03 0.000000e+00 -4.000000e+03 1.759259e-01 0.000000e+00 1.851852e-01
+ 2.000000e+03 0.000000e+00 -4.000000e+03 3.518519e-01 0.000000e+00 1.851852e-01
+ 3.000000e+03 0.000000e+00 -4.000000e+03 5.277778e-01 0.000000e+00 1.851852e-01
+ 4.000000e+03 0.000000e+00 -4.000000e+03 7.037037e-01 0.000000e+00 1.851852e-01
+ -4.000000e+03 1.000000e+03 -4.000000e+03 -7.037037e-01 -4.629630e-02 1.851852e-01
+ -3.000000e+03 1.000000e+03 -4.000000e+03 -5.277778e-01 -4.629630e-02 1.851852e-01
+ -2.000000e+03 1.000000e+03 -4.000000e+03 -3.518519e-01 -4.629630e-02 1.851852e-01
+ -1.000000e+03 1.000000e+03 -4.000000e+03 -1.759259e-01 -4.629630e-02 1.851852e-01
+ 0.000000e+00 1.000000e+03 -4.000000e+03 0.000000e+00 -4.629630e-02 1.851852e-01
+ 1.000000e+03 1.000000e+03 -4.000000e+03 1.759259e-01 -4.629630e-02 1.851852e-01
+ 2.000000e+03 1.000000e+03 -4.000000e+03 3.518519e-01 -4.629630e-02 1.851852e-01
+ 3.000000e+03 1.000000e+03 -4.000000e+03 5.277778e-01 -4.629630e-02 1.851852e-01
+ 4.000000e+03 1.000000e+03 -4.000000e+03 7.037037e-01 -4.629630e-02 1.851852e-01
+ -4.000000e+03 2.000000e+03 -4.000000e+03 -7.037037e-01 -9.259259e-02 1.851852e-01
+ -3.000000e+03 2.000000e+03 -4.000000e+03 -5.277778e-01 -9.259259e-02 1.851852e-01
+ -2.000000e+03 2.000000e+03 -4.000000e+03 -3.518519e-01 -9.259259e-02 1.851852e-01
+ -1.000000e+03 2.000000e+03 -4.000000e+03 -1.759259e-01 -9.259259e-02 1.851852e-01
+ 0.000000e+00 2.000000e+03 -4.000000e+03 0.000000e+00 -9.259259e-02 1.851852e-01
+ 1.000000e+03 2.000000e+03 -4.000000e+03 1.759259e-01 -9.259259e-02 1.851852e-01
+ 2.000000e+03 2.000000e+03 -4.000000e+03 3.518519e-01 -9.259259e-02 1.851852e-01
+ 3.000000e+03 2.000000e+03 -4.000000e+03 5.277778e-01 -9.259259e-02 1.851852e-01
+ 4.000000e+03 2.000000e+03 -4.000000e+03 7.037037e-01 -9.259259e-02 1.851852e-01
+ -4.000000e+03 3.000000e+03 -4.000000e+03 -7.037037e-01 -1.388889e-01 1.851852e-01
+ -3.000000e+03 3.000000e+03 -4.000000e+03 -5.277778e-01 -1.388889e-01 1.851852e-01
+ -2.000000e+03 3.000000e+03 -4.000000e+03 -3.518519e-01 -1.388889e-01 1.851852e-01
+ -1.000000e+03 3.000000e+03 -4.000000e+03 -1.759259e-01 -1.388889e-01 1.851852e-01
+ 0.000000e+00 3.000000e+03 -4.000000e+03 0.000000e+00 -1.388889e-01 1.851852e-01
+ 1.000000e+03 3.000000e+03 -4.000000e+03 1.759259e-01 -1.388889e-01 1.851852e-01
+ 2.000000e+03 3.000000e+03 -4.000000e+03 3.518519e-01 -1.388889e-01 1.851852e-01
+ 3.000000e+03 3.000000e+03 -4.000000e+03 5.277778e-01 -1.388889e-01 1.851852e-01
+ 4.000000e+03 3.000000e+03 -4.000000e+03 7.037037e-01 -1.388889e-01 1.851852e-01
+ -4.000000e+03 4.000000e+03 -4.000000e+03 -7.037037e-01 -1.851852e-01 1.851852e-01
+ -3.000000e+03 4.000000e+03 -4.000000e+03 -5.277778e-01 -1.851852e-01 1.851852e-01
+ -2.000000e+03 4.000000e+03 -4.000000e+03 -3.518519e-01 -1.851852e-01 1.851852e-01
+ -1.000000e+03 4.000000e+03 -4.000000e+03 -1.759259e-01 -1.851852e-01 1.851852e-01
+ 0.000000e+00 4.000000e+03 -4.000000e+03 0.000000e+00 -1.851852e-01 1.851852e-01
+ 1.000000e+03 4.000000e+03 -4.000000e+03 1.759259e-01 -1.851852e-01 1.851852e-01
+ 2.000000e+03 4.000000e+03 -4.000000e+03 3.518519e-01 -1.851852e-01 1.851852e-01
+ 3.000000e+03 4.000000e+03 -4.000000e+03 5.277778e-01 -1.851852e-01 1.851852e-01
+ 4.000000e+03 4.000000e+03 -4.000000e+03 7.037037e-01 -1.851852e-01 1.851852e-01
+ -4.000000e+03 -4.000000e+03 -3.000000e+03 -7.037037e-01 1.851852e-01 1.388889e-01
+ -3.000000e+03 -4.000000e+03 -3.000000e+03 -5.277778e-01 1.851852e-01 1.388889e-01
+ -2.000000e+03 -4.000000e+03 -3.000000e+03 -3.518519e-01 1.851852e-01 1.388889e-01
+ -1.000000e+03 -4.000000e+03 -3.000000e+03 -1.759259e-01 1.851852e-01 1.388889e-01
+ 0.000000e+00 -4.000000e+03 -3.000000e+03 0.000000e+00 1.851852e-01 1.388889e-01
+ 1.000000e+03 -4.000000e+03 -3.000000e+03 1.759259e-01 1.851852e-01 1.388889e-01
+ 2.000000e+03 -4.000000e+03 -3.000000e+03 3.518519e-01 1.851852e-01 1.388889e-01
+ 3.000000e+03 -4.000000e+03 -3.000000e+03 5.277778e-01 1.851852e-01 1.388889e-01
+ 4.000000e+03 -4.000000e+03 -3.000000e+03 7.037037e-01 1.851852e-01 1.388889e-01
+ -4.000000e+03 -3.000000e+03 -3.000000e+03 -7.037037e-01 1.388889e-01 1.388889e-01
+ -3.000000e+03 -3.000000e+03 -3.000000e+03 -5.277778e-01 1.388889e-01 1.388889e-01
+ -2.000000e+03 -3.000000e+03 -3.000000e+03 -3.518519e-01 1.388889e-01 1.388889e-01
+ -1.000000e+03 -3.000000e+03 -3.000000e+03 -1.759259e-01 1.388889e-01 1.388889e-01
+ 0.000000e+00 -3.000000e+03 -3.000000e+03 0.000000e+00 1.388889e-01 1.388889e-01
+ 1.000000e+03 -3.000000e+03 -3.000000e+03 1.759259e-01 1.388889e-01 1.388889e-01
+ 2.000000e+03 -3.000000e+03 -3.000000e+03 3.518519e-01 1.388889e-01 1.388889e-01
+ 3.000000e+03 -3.000000e+03 -3.000000e+03 5.277778e-01 1.388889e-01 1.388889e-01
+ 4.000000e+03 -3.000000e+03 -3.000000e+03 7.037037e-01 1.388889e-01 1.388889e-01
+ -4.000000e+03 -2.000000e+03 -3.000000e+03 -7.037037e-01 9.259259e-02 1.388889e-01
+ -3.000000e+03 -2.000000e+03 -3.000000e+03 -5.277778e-01 9.259259e-02 1.388889e-01
+ -2.000000e+03 -2.000000e+03 -3.000000e+03 -3.518519e-01 9.259259e-02 1.388889e-01
+ -1.000000e+03 -2.000000e+03 -3.000000e+03 -1.759259e-01 9.259259e-02 1.388889e-01
+ 0.000000e+00 -2.000000e+03 -3.000000e+03 0.000000e+00 9.259259e-02 1.388889e-01
+ 1.000000e+03 -2.000000e+03 -3.000000e+03 1.759259e-01 9.259259e-02 1.388889e-01
+ 2.000000e+03 -2.000000e+03 -3.000000e+03 3.518519e-01 9.259259e-02 1.388889e-01
+ 3.000000e+03 -2.000000e+03 -3.000000e+03 5.277778e-01 9.259259e-02 1.388889e-01
+ 4.000000e+03 -2.000000e+03 -3.000000e+03 7.037037e-01 9.259259e-02 1.388889e-01
+ -4.000000e+03 -1.000000e+03 -3.000000e+03 -7.037037e-01 4.629630e-02 1.388889e-01
+ -3.000000e+03 -1.000000e+03 -3.000000e+03 -5.277778e-01 4.629630e-02 1.388889e-01
+ -2.000000e+03 -1.000000e+03 -3.000000e+03 -3.518519e-01 4.629630e-02 1.388889e-01
+ -1.000000e+03 -1.000000e+03 -3.000000e+03 -1.759259e-01 4.629630e-02 1.388889e-01
+ 0.000000e+00 -1.000000e+03 -3.000000e+03 0.000000e+00 4.629630e-02 1.388889e-01
+ 1.000000e+03 -1.000000e+03 -3.000000e+03 1.759259e-01 4.629630e-02 1.388889e-01
+ 2.000000e+03 -1.000000e+03 -3.000000e+03 3.518519e-01 4.629630e-02 1.388889e-01
+ 3.000000e+03 -1.000000e+03 -3.000000e+03 5.277778e-01 4.629630e-02 1.388889e-01
+ 4.000000e+03 -1.000000e+03 -3.000000e+03 7.037037e-01 4.629630e-02 1.388889e-01
+ -4.000000e+03 0.000000e+00 -3.000000e+03 -7.037037e-01 -0.000000e+00 1.388889e-01
+ -3.000000e+03 0.000000e+00 -3.000000e+03 -5.277778e-01 -0.000000e+00 1.388889e-01
+ -2.000000e+03 0.000000e+00 -3.000000e+03 -3.518519e-01 -0.000000e+00 1.388889e-01
+ -1.000000e+03 0.000000e+00 -3.000000e+03 -1.759259e-01 -0.000000e+00 1.388889e-01
+ 0.000000e+00 0.000000e+00 -3.000000e+03 0.000000e+00 0.000000e+00 1.388889e-01
+ 1.000000e+03 0.000000e+00 -3.000000e+03 1.759259e-01 0.000000e+00 1.388889e-01
+ 2.000000e+03 0.000000e+00 -3.000000e+03 3.518519e-01 0.000000e+00 1.388889e-01
+ 3.000000e+03 0.000000e+00 -3.000000e+03 5.277778e-01 0.000000e+00 1.388889e-01
+ 4.000000e+03 0.000000e+00 -3.000000e+03 7.037037e-01 0.000000e+00 1.388889e-01
+ -4.000000e+03 1.000000e+03 -3.000000e+03 -7.037037e-01 -4.629630e-02 1.388889e-01
+ -3.000000e+03 1.000000e+03 -3.000000e+03 -5.277778e-01 -4.629630e-02 1.388889e-01
+ -2.000000e+03 1.000000e+03 -3.000000e+03 -3.518519e-01 -4.629630e-02 1.388889e-01
+ -1.000000e+03 1.000000e+03 -3.000000e+03 -1.759259e-01 -4.629630e-02 1.388889e-01
+ 0.000000e+00 1.000000e+03 -3.000000e+03 0.000000e+00 -4.629630e-02 1.388889e-01
+ 1.000000e+03 1.000000e+03 -3.000000e+03 1.759259e-01 -4.629630e-02 1.388889e-01
+ 2.000000e+03 1.000000e+03 -3.000000e+03 3.518519e-01 -4.629630e-02 1.388889e-01
+ 3.000000e+03 1.000000e+03 -3.000000e+03 5.277778e-01 -4.629630e-02 1.388889e-01
+ 4.000000e+03 1.000000e+03 -3.000000e+03 7.037037e-01 -4.629630e-02 1.388889e-01
+ -4.000000e+03 2.000000e+03 -3.000000e+03 -7.037037e-01 -9.259259e-02 1.388889e-01
+ -3.000000e+03 2.000000e+03 -3.000000e+03 -5.277778e-01 -9.259259e-02 1.388889e-01
+ -2.000000e+03 2.000000e+03 -3.000000e+03 -3.518519e-01 -9.259259e-02 1.388889e-01
+ -1.000000e+03 2.000000e+03 -3.000000e+03 -1.759259e-01 -9.259259e-02 1.388889e-01
+ 0.000000e+00 2.000000e+03 -3.000000e+03 0.000000e+00 -9.259259e-02 1.388889e-01
+ 1.000000e+03 2.000000e+03 -3.000000e+03 1.759259e-01 -9.259259e-02 1.388889e-01
+ 2.000000e+03 2.000000e+03 -3.000000e+03 3.518519e-01 -9.259259e-02 1.388889e-01
+ 3.000000e+03 2.000000e+03 -3.000000e+03 5.277778e-01 -9.259259e-02 1.388889e-01
+ 4.000000e+03 2.000000e+03 -3.000000e+03 7.037037e-01 -9.259259e-02 1.388889e-01
+ -4.000000e+03 3.000000e+03 -3.000000e+03 -7.037037e-01 -1.388889e-01 1.388889e-01
+ -3.000000e+03 3.000000e+03 -3.000000e+03 -5.277778e-01 -1.388889e-01 1.388889e-01
+ -2.000000e+03 3.000000e+03 -3.000000e+03 -3.518519e-01 -1.388889e-01 1.388889e-01
+ -1.000000e+03 3.000000e+03 -3.000000e+03 -1.759259e-01 -1.388889e-01 1.388889e-01
+ 0.000000e+00 3.000000e+03 -3.000000e+03 0.000000e+00 -1.388889e-01 1.388889e-01
+ 1.000000e+03 3.000000e+03 -3.000000e+03 1.759259e-01 -1.388889e-01 1.388889e-01
+ 2.000000e+03 3.000000e+03 -3.000000e+03 3.518519e-01 -1.388889e-01 1.388889e-01
+ 3.000000e+03 3.000000e+03 -3.000000e+03 5.277778e-01 -1.388889e-01 1.388889e-01
+ 4.000000e+03 3.000000e+03 -3.000000e+03 7.037037e-01 -1.388889e-01 1.388889e-01
+ -4.000000e+03 4.000000e+03 -3.000000e+03 -7.037037e-01 -1.851852e-01 1.388889e-01
+ -3.000000e+03 4.000000e+03 -3.000000e+03 -5.277778e-01 -1.851852e-01 1.388889e-01
+ -2.000000e+03 4.000000e+03 -3.000000e+03 -3.518519e-01 -1.851852e-01 1.388889e-01
+ -1.000000e+03 4.000000e+03 -3.000000e+03 -1.759259e-01 -1.851852e-01 1.388889e-01
+ 0.000000e+00 4.000000e+03 -3.000000e+03 0.000000e+00 -1.851852e-01 1.388889e-01
+ 1.000000e+03 4.000000e+03 -3.000000e+03 1.759259e-01 -1.851852e-01 1.388889e-01
+ 2.000000e+03 4.000000e+03 -3.000000e+03 3.518519e-01 -1.851852e-01 1.388889e-01
+ 3.000000e+03 4.000000e+03 -3.000000e+03 5.277778e-01 -1.851852e-01 1.388889e-01
+ 4.000000e+03 4.000000e+03 -3.000000e+03 7.037037e-01 -1.851852e-01 1.388889e-01
+ -4.000000e+03 -4.000000e+03 -2.000000e+03 -7.037037e-01 1.851852e-01 9.259259e-02
+ -3.000000e+03 -4.000000e+03 -2.000000e+03 -5.277778e-01 1.851852e-01 9.259259e-02
+ -2.000000e+03 -4.000000e+03 -2.000000e+03 -3.518519e-01 1.851852e-01 9.259259e-02
+ -1.000000e+03 -4.000000e+03 -2.000000e+03 -1.759259e-01 1.851852e-01 9.259259e-02
+ 0.000000e+00 -4.000000e+03 -2.000000e+03 0.000000e+00 1.851852e-01 9.259259e-02
+ 1.000000e+03 -4.000000e+03 -2.000000e+03 1.759259e-01 1.851852e-01 9.259259e-02
+ 2.000000e+03 -4.000000e+03 -2.000000e+03 3.518519e-01 1.851852e-01 9.259259e-02
+ 3.000000e+03 -4.000000e+03 -2.000000e+03 5.277778e-01 1.851852e-01 9.259259e-02
+ 4.000000e+03 -4.000000e+03 -2.000000e+03 7.037037e-01 1.851852e-01 9.259259e-02
+ -4.000000e+03 -3.000000e+03 -2.000000e+03 -7.037037e-01 1.388889e-01 9.259259e-02
+ -3.000000e+03 -3.000000e+03 -2.000000e+03 -5.277778e-01 1.388889e-01 9.259259e-02
+ -2.000000e+03 -3.000000e+03 -2.000000e+03 -3.518519e-01 1.388889e-01 9.259259e-02
+ -1.000000e+03 -3.000000e+03 -2.000000e+03 -1.759259e-01 1.388889e-01 9.259259e-02
+ 0.000000e+00 -3.000000e+03 -2.000000e+03 0.000000e+00 1.388889e-01 9.259259e-02
+ 1.000000e+03 -3.000000e+03 -2.000000e+03 1.759259e-01 1.388889e-01 9.259259e-02
+ 2.000000e+03 -3.000000e+03 -2.000000e+03 3.518519e-01 1.388889e-01 9.259259e-02
+ 3.000000e+03 -3.000000e+03 -2.000000e+03 5.277778e-01 1.388889e-01 9.259259e-02
+ 4.000000e+03 -3.000000e+03 -2.000000e+03 7.037037e-01 1.388889e-01 9.259259e-02
+ -4.000000e+03 -2.000000e+03 -2.000000e+03 -7.037037e-01 9.259259e-02 9.259259e-02
+ -3.000000e+03 -2.000000e+03 -2.000000e+03 -5.277778e-01 9.259259e-02 9.259259e-02
+ -2.000000e+03 -2.000000e+03 -2.000000e+03 -3.518519e-01 9.259259e-02 9.259259e-02
+ -1.000000e+03 -2.000000e+03 -2.000000e+03 -1.759259e-01 9.259259e-02 9.259259e-02
+ 0.000000e+00 -2.000000e+03 -2.000000e+03 0.000000e+00 9.259259e-02 9.259259e-02
+ 1.000000e+03 -2.000000e+03 -2.000000e+03 1.759259e-01 9.259259e-02 9.259259e-02
+ 2.000000e+03 -2.000000e+03 -2.000000e+03 3.518519e-01 9.259259e-02 9.259259e-02
+ 3.000000e+03 -2.000000e+03 -2.000000e+03 5.277778e-01 9.259259e-02 9.259259e-02
+ 4.000000e+03 -2.000000e+03 -2.000000e+03 7.037037e-01 9.259259e-02 9.259259e-02
+ -4.000000e+03 -1.000000e+03 -2.000000e+03 -7.037037e-01 4.629630e-02 9.259259e-02
+ -3.000000e+03 -1.000000e+03 -2.000000e+03 -5.277778e-01 4.629630e-02 9.259259e-02
+ -2.000000e+03 -1.000000e+03 -2.000000e+03 -3.518519e-01 4.629630e-02 9.259259e-02
+ -1.000000e+03 -1.000000e+03 -2.000000e+03 -1.759259e-01 4.629630e-02 9.259259e-02
+ 0.000000e+00 -1.000000e+03 -2.000000e+03 0.000000e+00 4.629630e-02 9.259259e-02
+ 1.000000e+03 -1.000000e+03 -2.000000e+03 1.759259e-01 4.629630e-02 9.259259e-02
+ 2.000000e+03 -1.000000e+03 -2.000000e+03 3.518519e-01 4.629630e-02 9.259259e-02
+ 3.000000e+03 -1.000000e+03 -2.000000e+03 5.277778e-01 4.629630e-02 9.259259e-02
+ 4.000000e+03 -1.000000e+03 -2.000000e+03 7.037037e-01 4.629630e-02 9.259259e-02
+ -4.000000e+03 0.000000e+00 -2.000000e+03 -7.037037e-01 -0.000000e+00 9.259259e-02
+ -3.000000e+03 0.000000e+00 -2.000000e+03 -5.277778e-01 -0.000000e+00 9.259259e-02
+ -2.000000e+03 0.000000e+00 -2.000000e+03 -3.518519e-01 -0.000000e+00 9.259259e-02
+ -1.000000e+03 0.000000e+00 -2.000000e+03 -1.759259e-01 -0.000000e+00 9.259259e-02
+ 0.000000e+00 0.000000e+00 -2.000000e+03 0.000000e+00 0.000000e+00 9.259259e-02
+ 1.000000e+03 0.000000e+00 -2.000000e+03 1.759259e-01 0.000000e+00 9.259259e-02
+ 2.000000e+03 0.000000e+00 -2.000000e+03 3.518519e-01 0.000000e+00 9.259259e-02
+ 3.000000e+03 0.000000e+00 -2.000000e+03 5.277778e-01 0.000000e+00 9.259259e-02
+ 4.000000e+03 0.000000e+00 -2.000000e+03 7.037037e-01 0.000000e+00 9.259259e-02
+ -4.000000e+03 1.000000e+03 -2.000000e+03 -7.037037e-01 -4.629630e-02 9.259259e-02
+ -3.000000e+03 1.000000e+03 -2.000000e+03 -5.277778e-01 -4.629630e-02 9.259259e-02
+ -2.000000e+03 1.000000e+03 -2.000000e+03 -3.518519e-01 -4.629630e-02 9.259259e-02
+ -1.000000e+03 1.000000e+03 -2.000000e+03 -1.759259e-01 -4.629630e-02 9.259259e-02
+ 0.000000e+00 1.000000e+03 -2.000000e+03 0.000000e+00 -4.629630e-02 9.259259e-02
+ 1.000000e+03 1.000000e+03 -2.000000e+03 1.759259e-01 -4.629630e-02 9.259259e-02
+ 2.000000e+03 1.000000e+03 -2.000000e+03 3.518519e-01 -4.629630e-02 9.259259e-02
+ 3.000000e+03 1.000000e+03 -2.000000e+03 5.277778e-01 -4.629630e-02 9.259259e-02
+ 4.000000e+03 1.000000e+03 -2.000000e+03 7.037037e-01 -4.629630e-02 9.259259e-02
+ -4.000000e+03 2.000000e+03 -2.000000e+03 -7.037037e-01 -9.259259e-02 9.259259e-02
+ -3.000000e+03 2.000000e+03 -2.000000e+03 -5.277778e-01 -9.259259e-02 9.259259e-02
+ -2.000000e+03 2.000000e+03 -2.000000e+03 -3.518519e-01 -9.259259e-02 9.259259e-02
+ -1.000000e+03 2.000000e+03 -2.000000e+03 -1.759259e-01 -9.259259e-02 9.259259e-02
+ 0.000000e+00 2.000000e+03 -2.000000e+03 0.000000e+00 -9.259259e-02 9.259259e-02
+ 1.000000e+03 2.000000e+03 -2.000000e+03 1.759259e-01 -9.259259e-02 9.259259e-02
+ 2.000000e+03 2.000000e+03 -2.000000e+03 3.518519e-01 -9.259259e-02 9.259259e-02
+ 3.000000e+03 2.000000e+03 -2.000000e+03 5.277778e-01 -9.259259e-02 9.259259e-02
+ 4.000000e+03 2.000000e+03 -2.000000e+03 7.037037e-01 -9.259259e-02 9.259259e-02
+ -4.000000e+03 3.000000e+03 -2.000000e+03 -7.037037e-01 -1.388889e-01 9.259259e-02
+ -3.000000e+03 3.000000e+03 -2.000000e+03 -5.277778e-01 -1.388889e-01 9.259259e-02
+ -2.000000e+03 3.000000e+03 -2.000000e+03 -3.518519e-01 -1.388889e-01 9.259259e-02
+ -1.000000e+03 3.000000e+03 -2.000000e+03 -1.759259e-01 -1.388889e-01 9.259259e-02
+ 0.000000e+00 3.000000e+03 -2.000000e+03 0.000000e+00 -1.388889e-01 9.259259e-02
+ 1.000000e+03 3.000000e+03 -2.000000e+03 1.759259e-01 -1.388889e-01 9.259259e-02
+ 2.000000e+03 3.000000e+03 -2.000000e+03 3.518519e-01 -1.388889e-01 9.259259e-02
+ 3.000000e+03 3.000000e+03 -2.000000e+03 5.277778e-01 -1.388889e-01 9.259259e-02
+ 4.000000e+03 3.000000e+03 -2.000000e+03 7.037037e-01 -1.388889e-01 9.259259e-02
+ -4.000000e+03 4.000000e+03 -2.000000e+03 -7.037037e-01 -1.851852e-01 9.259259e-02
+ -3.000000e+03 4.000000e+03 -2.000000e+03 -5.277778e-01 -1.851852e-01 9.259259e-02
+ -2.000000e+03 4.000000e+03 -2.000000e+03 -3.518519e-01 -1.851852e-01 9.259259e-02
+ -1.000000e+03 4.000000e+03 -2.000000e+03 -1.759259e-01 -1.851852e-01 9.259259e-02
+ 0.000000e+00 4.000000e+03 -2.000000e+03 0.000000e+00 -1.851852e-01 9.259259e-02
+ 1.000000e+03 4.000000e+03 -2.000000e+03 1.759259e-01 -1.851852e-01 9.259259e-02
+ 2.000000e+03 4.000000e+03 -2.000000e+03 3.518519e-01 -1.851852e-01 9.259259e-02
+ 3.000000e+03 4.000000e+03 -2.000000e+03 5.277778e-01 -1.851852e-01 9.259259e-02
+ 4.000000e+03 4.000000e+03 -2.000000e+03 7.037037e-01 -1.851852e-01 9.259259e-02
+ -4.000000e+03 -4.000000e+03 -1.000000e+03 -7.037037e-01 1.851852e-01 4.629630e-02
+ -3.000000e+03 -4.000000e+03 -1.000000e+03 -5.277778e-01 1.851852e-01 4.629630e-02
+ -2.000000e+03 -4.000000e+03 -1.000000e+03 -3.518519e-01 1.851852e-01 4.629630e-02
+ -1.000000e+03 -4.000000e+03 -1.000000e+03 -1.759259e-01 1.851852e-01 4.629630e-02
+ 0.000000e+00 -4.000000e+03 -1.000000e+03 0.000000e+00 1.851852e-01 4.629630e-02
+ 1.000000e+03 -4.000000e+03 -1.000000e+03 1.759259e-01 1.851852e-01 4.629630e-02
+ 2.000000e+03 -4.000000e+03 -1.000000e+03 3.518519e-01 1.851852e-01 4.629630e-02
+ 3.000000e+03 -4.000000e+03 -1.000000e+03 5.277778e-01 1.851852e-01 4.629630e-02
+ 4.000000e+03 -4.000000e+03 -1.000000e+03 7.037037e-01 1.851852e-01 4.629630e-02
+ -4.000000e+03 -3.000000e+03 -1.000000e+03 -7.037037e-01 1.388889e-01 4.629630e-02
+ -3.000000e+03 -3.000000e+03 -1.000000e+03 -5.277778e-01 1.388889e-01 4.629630e-02
+ -2.000000e+03 -3.000000e+03 -1.000000e+03 -3.518519e-01 1.388889e-01 4.629630e-02
+ -1.000000e+03 -3.000000e+03 -1.000000e+03 -1.759259e-01 1.388889e-01 4.629630e-02
+ 0.000000e+00 -3.000000e+03 -1.000000e+03 0.000000e+00 1.388889e-01 4.629630e-02
+ 1.000000e+03 -3.000000e+03 -1.000000e+03 1.759259e-01 1.388889e-01 4.629630e-02
+ 2.000000e+03 -3.000000e+03 -1.000000e+03 3.518519e-01 1.388889e-01 4.629630e-02
+ 3.000000e+03 -3.000000e+03 -1.000000e+03 5.277778e-01 1.388889e-01 4.629630e-02
+ 4.000000e+03 -3.000000e+03 -1.000000e+03 7.037037e-01 1.388889e-01 4.629630e-02
+ -4.000000e+03 -2.000000e+03 -1.000000e+03 -7.037037e-01 9.259259e-02 4.629630e-02
+ -3.000000e+03 -2.000000e+03 -1.000000e+03 -5.277778e-01 9.259259e-02 4.629630e-02
+ -2.000000e+03 -2.000000e+03 -1.000000e+03 -3.518519e-01 9.259259e-02 4.629630e-02
+ -1.000000e+03 -2.000000e+03 -1.000000e+03 -1.759259e-01 9.259259e-02 4.629630e-02
+ 0.000000e+00 -2.000000e+03 -1.000000e+03 0.000000e+00 9.259259e-02 4.629630e-02
+ 1.000000e+03 -2.000000e+03 -1.000000e+03 1.759259e-01 9.259259e-02 4.629630e-02
+ 2.000000e+03 -2.000000e+03 -1.000000e+03 3.518519e-01 9.259259e-02 4.629630e-02
+ 3.000000e+03 -2.000000e+03 -1.000000e+03 5.277778e-01 9.259259e-02 4.629630e-02
+ 4.000000e+03 -2.000000e+03 -1.000000e+03 7.037037e-01 9.259259e-02 4.629630e-02
+ -4.000000e+03 -1.000000e+03 -1.000000e+03 -7.037037e-01 4.629630e-02 4.629630e-02
+ -3.000000e+03 -1.000000e+03 -1.000000e+03 -5.277778e-01 4.629630e-02 4.629630e-02
+ -2.000000e+03 -1.000000e+03 -1.000000e+03 -3.518519e-01 4.629630e-02 4.629630e-02
+ -1.000000e+03 -1.000000e+03 -1.000000e+03 -1.759259e-01 4.629630e-02 4.629630e-02
+ 0.000000e+00 -1.000000e+03 -1.000000e+03 0.000000e+00 4.629630e-02 4.629630e-02
+ 1.000000e+03 -1.000000e+03 -1.000000e+03 1.759259e-01 4.629630e-02 4.629630e-02
+ 2.000000e+03 -1.000000e+03 -1.000000e+03 3.518519e-01 4.629630e-02 4.629630e-02
+ 3.000000e+03 -1.000000e+03 -1.000000e+03 5.277778e-01 4.629630e-02 4.629630e-02
+ 4.000000e+03 -1.000000e+03 -1.000000e+03 7.037037e-01 4.629630e-02 4.629630e-02
+ -4.000000e+03 0.000000e+00 -1.000000e+03 -7.037037e-01 -0.000000e+00 4.629630e-02
+ -3.000000e+03 0.000000e+00 -1.000000e+03 -5.277778e-01 -0.000000e+00 4.629630e-02
+ -2.000000e+03 0.000000e+00 -1.000000e+03 -3.518519e-01 -0.000000e+00 4.629630e-02
+ -1.000000e+03 0.000000e+00 -1.000000e+03 -1.759259e-01 -0.000000e+00 4.629630e-02
+ 0.000000e+00 0.000000e+00 -1.000000e+03 0.000000e+00 0.000000e+00 4.629630e-02
+ 1.000000e+03 0.000000e+00 -1.000000e+03 1.759259e-01 0.000000e+00 4.629630e-02
+ 2.000000e+03 0.000000e+00 -1.000000e+03 3.518519e-01 0.000000e+00 4.629630e-02
+ 3.000000e+03 0.000000e+00 -1.000000e+03 5.277778e-01 0.000000e+00 4.629630e-02
+ 4.000000e+03 0.000000e+00 -1.000000e+03 7.037037e-01 0.000000e+00 4.629630e-02
+ -4.000000e+03 1.000000e+03 -1.000000e+03 -7.037037e-01 -4.629630e-02 4.629630e-02
+ -3.000000e+03 1.000000e+03 -1.000000e+03 -5.277778e-01 -4.629630e-02 4.629630e-02
+ -2.000000e+03 1.000000e+03 -1.000000e+03 -3.518519e-01 -4.629630e-02 4.629630e-02
+ -1.000000e+03 1.000000e+03 -1.000000e+03 -1.759259e-01 -4.629630e-02 4.629630e-02
+ 0.000000e+00 1.000000e+03 -1.000000e+03 0.000000e+00 -4.629630e-02 4.629630e-02
+ 1.000000e+03 1.000000e+03 -1.000000e+03 1.759259e-01 -4.629630e-02 4.629630e-02
+ 2.000000e+03 1.000000e+03 -1.000000e+03 3.518519e-01 -4.629630e-02 4.629630e-02
+ 3.000000e+03 1.000000e+03 -1.000000e+03 5.277778e-01 -4.629630e-02 4.629630e-02
+ 4.000000e+03 1.000000e+03 -1.000000e+03 7.037037e-01 -4.629630e-02 4.629630e-02
+ -4.000000e+03 2.000000e+03 -1.000000e+03 -7.037037e-01 -9.259259e-02 4.629630e-02
+ -3.000000e+03 2.000000e+03 -1.000000e+03 -5.277778e-01 -9.259259e-02 4.629630e-02
+ -2.000000e+03 2.000000e+03 -1.000000e+03 -3.518519e-01 -9.259259e-02 4.629630e-02
+ -1.000000e+03 2.000000e+03 -1.000000e+03 -1.759259e-01 -9.259259e-02 4.629630e-02
+ 0.000000e+00 2.000000e+03 -1.000000e+03 0.000000e+00 -9.259259e-02 4.629630e-02
+ 1.000000e+03 2.000000e+03 -1.000000e+03 1.759259e-01 -9.259259e-02 4.629630e-02
+ 2.000000e+03 2.000000e+03 -1.000000e+03 3.518519e-01 -9.259259e-02 4.629630e-02
+ 3.000000e+03 2.000000e+03 -1.000000e+03 5.277778e-01 -9.259259e-02 4.629630e-02
+ 4.000000e+03 2.000000e+03 -1.000000e+03 7.037037e-01 -9.259259e-02 4.629630e-02
+ -4.000000e+03 3.000000e+03 -1.000000e+03 -7.037037e-01 -1.388889e-01 4.629630e-02
+ -3.000000e+03 3.000000e+03 -1.000000e+03 -5.277778e-01 -1.388889e-01 4.629630e-02
+ -2.000000e+03 3.000000e+03 -1.000000e+03 -3.518519e-01 -1.388889e-01 4.629630e-02
+ -1.000000e+03 3.000000e+03 -1.000000e+03 -1.759259e-01 -1.388889e-01 4.629630e-02
+ 0.000000e+00 3.000000e+03 -1.000000e+03 0.000000e+00 -1.388889e-01 4.629630e-02
+ 1.000000e+03 3.000000e+03 -1.000000e+03 1.759259e-01 -1.388889e-01 4.629630e-02
+ 2.000000e+03 3.000000e+03 -1.000000e+03 3.518519e-01 -1.388889e-01 4.629630e-02
+ 3.000000e+03 3.000000e+03 -1.000000e+03 5.277778e-01 -1.388889e-01 4.629630e-02
+ 4.000000e+03 3.000000e+03 -1.000000e+03 7.037037e-01 -1.388889e-01 4.629630e-02
+ -4.000000e+03 4.000000e+03 -1.000000e+03 -7.037037e-01 -1.851852e-01 4.629630e-02
+ -3.000000e+03 4.000000e+03 -1.000000e+03 -5.277778e-01 -1.851852e-01 4.629630e-02
+ -2.000000e+03 4.000000e+03 -1.000000e+03 -3.518519e-01 -1.851852e-01 4.629630e-02
+ -1.000000e+03 4.000000e+03 -1.000000e+03 -1.759259e-01 -1.851852e-01 4.629630e-02
+ 0.000000e+00 4.000000e+03 -1.000000e+03 0.000000e+00 -1.851852e-01 4.629630e-02
+ 1.000000e+03 4.000000e+03 -1.000000e+03 1.759259e-01 -1.851852e-01 4.629630e-02
+ 2.000000e+03 4.000000e+03 -1.000000e+03 3.518519e-01 -1.851852e-01 4.629630e-02
+ 3.000000e+03 4.000000e+03 -1.000000e+03 5.277778e-01 -1.851852e-01 4.629630e-02
+ 4.000000e+03 4.000000e+03 -1.000000e+03 7.037037e-01 -1.851852e-01 4.629630e-02
+ -4.000000e+03 -4.000000e+03 0.000000e+00 -7.037037e-01 1.851852e-01 -0.000000e+00
+ -3.000000e+03 -4.000000e+03 0.000000e+00 -5.277778e-01 1.851852e-01 -0.000000e+00
+ -2.000000e+03 -4.000000e+03 0.000000e+00 -3.518519e-01 1.851852e-01 -0.000000e+00
+ -1.000000e+03 -4.000000e+03 0.000000e+00 -1.759259e-01 1.851852e-01 -0.000000e+00
+ 0.000000e+00 -4.000000e+03 0.000000e+00 0.000000e+00 1.851852e-01 0.000000e+00
+ 1.000000e+03 -4.000000e+03 0.000000e+00 1.759259e-01 1.851852e-01 0.000000e+00
+ 2.000000e+03 -4.000000e+03 0.000000e+00 3.518519e-01 1.851852e-01 0.000000e+00
+ 3.000000e+03 -4.000000e+03 0.000000e+00 5.277778e-01 1.851852e-01 0.000000e+00
+ 4.000000e+03 -4.000000e+03 0.000000e+00 7.037037e-01 1.851852e-01 0.000000e+00
+ -4.000000e+03 -3.000000e+03 0.000000e+00 -7.037037e-01 1.388889e-01 -0.000000e+00
+ -3.000000e+03 -3.000000e+03 0.000000e+00 -5.277778e-01 1.388889e-01 -0.000000e+00
+ -2.000000e+03 -3.000000e+03 0.000000e+00 -3.518519e-01 1.388889e-01 -0.000000e+00
+ -1.000000e+03 -3.000000e+03 0.000000e+00 -1.759259e-01 1.388889e-01 -0.000000e+00
+ 0.000000e+00 -3.000000e+03 0.000000e+00 0.000000e+00 1.388889e-01 0.000000e+00
+ 1.000000e+03 -3.000000e+03 0.000000e+00 1.759259e-01 1.388889e-01 0.000000e+00
+ 2.000000e+03 -3.000000e+03 0.000000e+00 3.518519e-01 1.388889e-01 0.000000e+00
+ 3.000000e+03 -3.000000e+03 0.000000e+00 5.277778e-01 1.388889e-01 0.000000e+00
+ 4.000000e+03 -3.000000e+03 0.000000e+00 7.037037e-01 1.388889e-01 0.000000e+00
+ -4.000000e+03 -2.000000e+03 0.000000e+00 -7.037037e-01 9.259259e-02 -0.000000e+00
+ -3.000000e+03 -2.000000e+03 0.000000e+00 -5.277778e-01 9.259259e-02 -0.000000e+00
+ -2.000000e+03 -2.000000e+03 0.000000e+00 -3.518519e-01 9.259259e-02 -0.000000e+00
+ -1.000000e+03 -2.000000e+03 0.000000e+00 -1.759259e-01 9.259259e-02 -0.000000e+00
+ 0.000000e+00 -2.000000e+03 0.000000e+00 0.000000e+00 9.259259e-02 0.000000e+00
+ 1.000000e+03 -2.000000e+03 0.000000e+00 1.759259e-01 9.259259e-02 0.000000e+00
+ 2.000000e+03 -2.000000e+03 0.000000e+00 3.518519e-01 9.259259e-02 0.000000e+00
+ 3.000000e+03 -2.000000e+03 0.000000e+00 5.277778e-01 9.259259e-02 0.000000e+00
+ 4.000000e+03 -2.000000e+03 0.000000e+00 7.037037e-01 9.259259e-02 0.000000e+00
+ -4.000000e+03 -1.000000e+03 0.000000e+00 -7.037037e-01 4.629630e-02 -0.000000e+00
+ -3.000000e+03 -1.000000e+03 0.000000e+00 -5.277778e-01 4.629630e-02 -0.000000e+00
+ -2.000000e+03 -1.000000e+03 0.000000e+00 -3.518519e-01 4.629630e-02 -0.000000e+00
+ -1.000000e+03 -1.000000e+03 0.000000e+00 -1.759259e-01 4.629630e-02 -0.000000e+00
+ 0.000000e+00 -1.000000e+03 0.000000e+00 0.000000e+00 4.629630e-02 0.000000e+00
+ 1.000000e+03 -1.000000e+03 0.000000e+00 1.759259e-01 4.629630e-02 0.000000e+00
+ 2.000000e+03 -1.000000e+03 0.000000e+00 3.518519e-01 4.629630e-02 0.000000e+00
+ 3.000000e+03 -1.000000e+03 0.000000e+00 5.277778e-01 4.629630e-02 0.000000e+00
+ 4.000000e+03 -1.000000e+03 0.000000e+00 7.037037e-01 4.629630e-02 0.000000e+00
+ -4.000000e+03 0.000000e+00 0.000000e+00 -7.037037e-01 0.000000e+00 0.000000e+00
+ -3.000000e+03 0.000000e+00 0.000000e+00 -5.277778e-01 0.000000e+00 0.000000e+00
+ -2.000000e+03 0.000000e+00 0.000000e+00 -3.518519e-01 0.000000e+00 0.000000e+00
+ -1.000000e+03 0.000000e+00 0.000000e+00 -1.759259e-01 0.000000e+00 0.000000e+00
+ 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
+ 1.000000e+03 0.000000e+00 0.000000e+00 1.759259e-01 0.000000e+00 0.000000e+00
+ 2.000000e+03 0.000000e+00 0.000000e+00 3.518519e-01 0.000000e+00 0.000000e+00
+ 3.000000e+03 0.000000e+00 0.000000e+00 5.277778e-01 0.000000e+00 0.000000e+00
+ 4.000000e+03 0.000000e+00 0.000000e+00 7.037037e-01 0.000000e+00 0.000000e+00
+ -4.000000e+03 1.000000e+03 0.000000e+00 -7.037037e-01 -4.629630e-02 0.000000e+00
+ -3.000000e+03 1.000000e+03 0.000000e+00 -5.277778e-01 -4.629630e-02 0.000000e+00
+ -2.000000e+03 1.000000e+03 0.000000e+00 -3.518519e-01 -4.629630e-02 0.000000e+00
+ -1.000000e+03 1.000000e+03 0.000000e+00 -1.759259e-01 -4.629630e-02 0.000000e+00
+ 0.000000e+00 1.000000e+03 0.000000e+00 0.000000e+00 -4.629630e-02 0.000000e+00
+ 1.000000e+03 1.000000e+03 0.000000e+00 1.759259e-01 -4.629630e-02 0.000000e+00
+ 2.000000e+03 1.000000e+03 0.000000e+00 3.518519e-01 -4.629630e-02 0.000000e+00
+ 3.000000e+03 1.000000e+03 0.000000e+00 5.277778e-01 -4.629630e-02 0.000000e+00
+ 4.000000e+03 1.000000e+03 0.000000e+00 7.037037e-01 -4.629630e-02 0.000000e+00
+ -4.000000e+03 2.000000e+03 0.000000e+00 -7.037037e-01 -9.259259e-02 0.000000e+00
+ -3.000000e+03 2.000000e+03 0.000000e+00 -5.277778e-01 -9.259259e-02 0.000000e+00
+ -2.000000e+03 2.000000e+03 0.000000e+00 -3.518519e-01 -9.259259e-02 0.000000e+00
+ -1.000000e+03 2.000000e+03 0.000000e+00 -1.759259e-01 -9.259259e-02 0.000000e+00
+ 0.000000e+00 2.000000e+03 0.000000e+00 0.000000e+00 -9.259259e-02 0.000000e+00
+ 1.000000e+03 2.000000e+03 0.000000e+00 1.759259e-01 -9.259259e-02 0.000000e+00
+ 2.000000e+03 2.000000e+03 0.000000e+00 3.518519e-01 -9.259259e-02 0.000000e+00
+ 3.000000e+03 2.000000e+03 0.000000e+00 5.277778e-01 -9.259259e-02 0.000000e+00
+ 4.000000e+03 2.000000e+03 0.000000e+00 7.037037e-01 -9.259259e-02 0.000000e+00
+ -4.000000e+03 3.000000e+03 0.000000e+00 -7.037037e-01 -1.388889e-01 0.000000e+00
+ -3.000000e+03 3.000000e+03 0.000000e+00 -5.277778e-01 -1.388889e-01 0.000000e+00
+ -2.000000e+03 3.000000e+03 0.000000e+00 -3.518519e-01 -1.388889e-01 0.000000e+00
+ -1.000000e+03 3.000000e+03 0.000000e+00 -1.759259e-01 -1.388889e-01 0.000000e+00
+ 0.000000e+00 3.000000e+03 0.000000e+00 0.000000e+00 -1.388889e-01 0.000000e+00
+ 1.000000e+03 3.000000e+03 0.000000e+00 1.759259e-01 -1.388889e-01 0.000000e+00
+ 2.000000e+03 3.000000e+03 0.000000e+00 3.518519e-01 -1.388889e-01 0.000000e+00
+ 3.000000e+03 3.000000e+03 0.000000e+00 5.277778e-01 -1.388889e-01 0.000000e+00
+ 4.000000e+03 3.000000e+03 0.000000e+00 7.037037e-01 -1.388889e-01 0.000000e+00
+ -4.000000e+03 4.000000e+03 0.000000e+00 -7.037037e-01 -1.851852e-01 0.000000e+00
+ -3.000000e+03 4.000000e+03 0.000000e+00 -5.277778e-01 -1.851852e-01 0.000000e+00
+ -2.000000e+03 4.000000e+03 0.000000e+00 -3.518519e-01 -1.851852e-01 0.000000e+00
+ -1.000000e+03 4.000000e+03 0.000000e+00 -1.759259e-01 -1.851852e-01 0.000000e+00
+ 0.000000e+00 4.000000e+03 0.000000e+00 0.000000e+00 -1.851852e-01 0.000000e+00
+ 1.000000e+03 4.000000e+03 0.000000e+00 1.759259e-01 -1.851852e-01 0.000000e+00
+ 2.000000e+03 4.000000e+03 0.000000e+00 3.518519e-01 -1.851852e-01 0.000000e+00
+ 3.000000e+03 4.000000e+03 0.000000e+00 5.277778e-01 -1.851852e-01 0.000000e+00
+ 4.000000e+03 4.000000e+03 0.000000e+00 7.037037e-01 -1.851852e-01 0.000000e+00
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp.cfg 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,136 @@
+# -*- Python -*-
+[pylithapp]
+
+[pylithapp.launcher] # WARNING: THIS IS NOT PORTABLE
+command = mpirun -np ${nodes}
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+#timedependent = 1
+#implicit = 1
+#petsc = 1
+#solverlinear = 1
+#meshiocubit = 1
+#implicitelasticity = 1
+#quadrature3d = 1
+#fiatlagrange = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = mesh.exo
+use_nodeset_names = True
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+bc = [x_neg,x_pos,y_neg,z_neg]
+
+[pylithapp.timedependent.formulation]
+output = [domain,subdomain]
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic,viscoelastic]
+materials.elastic = pylith.materials.ElasticIsotropic3D
+materials.viscoelastic = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.degree= 2
+
+[pylithapp.timedependent.materials.viscoelastic]
+label = Elastic material 2
+id = 2
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.degree= 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC +x eface
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC -x face
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [1]
+label = face_yneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC -y face
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC -z edge
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+pc_type = asm
+
+# Change the preconditioner settings.
+sub_pc_factor_shift_type = none
+
+ksp_rtol = 1.0e-8
+ksp_max_it = 100
+ksp_gmres_restart = 50
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+#log_summary = true
+# start_in_debugger = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain.writer]
+filename = axialdisp.vtk
+
+[pylithapp.problem.formulation.output.subdomain]
+label = face_zpos
+writer.filename = axialdisp-groundsurf.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = axialdisp-elastic.vtk
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = axialdisp-viscoelastic.vtk
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_gendb.py
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_gendb.py (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_gendb.py 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,95 @@
+#!/usr/bin/env python
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+## @file tests/3d/hex8/axialdisp_gendb.py
+##
+## @brief Python script to generate spatial database with displacement
+## boundary conditions for the axial displacement test.
+
+import numpy
+
+class GenerateDB(object):
+ """
+ Python object to generate spatial database with displacement
+ boundary conditions for the axial displacement test.
+ """
+
+ def __init__(self):
+ """
+ Constructor.
+ """
+ return
+
+
+ def run(self):
+ """
+ Generate the database.
+ """
+ # Domain
+ x = numpy.arange(-4000.0, 4000.1, 500.0)
+ y = numpy.arange(-4000.0, 4000.1, 500.0)
+ z = numpy.arange(-6000.0, 0.1, 500.0)
+ nxpts = x.shape[0]
+ nypts = y.shape[0]
+ nzpts = z.shape[0]
+
+ xx = numpy.tile(x, (1,nypts*nzpts))
+ yy = numpy.tile(y, (nxpts,nzpts)).transpose()
+ zz = numpy.tile(z, (nxpts*nypts,1)).transpose()
+
+ xyz = numpy.zeros( (nxpts*nypts*nzpts, 3), dtype=numpy.float64)
+ xyz[:,0] = numpy.ravel(xx)
+ xyz[:,1] = numpy.ravel(yy)
+ xyz[:,2] = numpy.ravel(zz)
+
+ from axialdisp_soln import AnalyticalSoln
+ soln = AnalyticalSoln()
+ disp = soln.displacement(xyz)
+
+ from spatialdata.geocoords.CSCart import CSCart
+ cs = CSCart()
+ cs.inventory.spaceDim = 3
+ cs._configure()
+ data = {'points': xyz,
+ 'coordsys': cs,
+ 'data_dim': 3,
+ 'values': [{'name': "displacement-x",
+ 'units': "m",
+ 'data': numpy.ravel(disp[:,0])},
+ {'name': "displacement-y",
+ 'units': "m",
+ 'data': numpy.ravel(disp[:,1])},
+ {'name': "displacement-z",
+ 'units': "m",
+ 'data': numpy.ravel(disp[:,2])}]}
+
+ from spatialdata.spatialdb.SimpleIOAscii import SimpleIOAscii
+ io = SimpleIOAscii()
+ io.inventory.filename = "axial_disp.spatialdb"
+ io._configure()
+ io.write(data)
+ return
+
+
+# ======================================================================
+if __name__ == "__main__":
+ app = GenerateDB()
+ app.run()
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_soln.py
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_soln.py (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/axialdisp_soln.py 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,107 @@
+#!/usr/bin/env python
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+## @file tests/2d/quad4/axialdisp_soln.py
+##
+## @brief Analytical solution to axial displacement problem.
+
+## 3-D axial extension test with linear hexahedral cells.
+
+import numpy
+
+# Physical properties
+p_density = 2500.0
+p_vs = 3000.0
+p_vp = 5291.502622129181
+
+p_mu = p_density*p_vs**2
+p_lambda = p_density*p_vp**2 - 2*p_mu
+
+# Uniform stress field (plane strain)
+sxx = 1.0e+7
+syy = 0.0
+szz = 0.0
+sxy = 0.0
+syz = 0.0
+sxz = 0.0
+
+# Uniform strain field
+exx = 1.0/(2*p_mu) * (sxx - p_lambda/(3*p_lambda+2*p_mu) * (sxx+syy+szz))
+eyy = 1.0/(2*p_mu) * (syy - p_lambda/(3*p_lambda+2*p_mu) * (sxx+syy+szz))
+ezz = 1.0/(2*p_mu) * (szz - p_lambda/(3*p_lambda+2*p_mu) * (sxx+syy+szz))
+
+exy = 1.0/(2*p_mu) * (sxy)
+eyz = 1.0/(2*p_mu) * (syz)
+exz = 1.0/(2*p_mu) * (sxz)
+
+#print exx,eyy,exy,ezz
+#print -exx*p_lambda/(p_lambda+2*p_mu)
+
+# ----------------------------------------------------------------------
+class AnalyticalSoln(object):
+ """
+ Analytical solution to axial/shear displacement problem.
+ """
+
+ def __init__(self):
+ return
+
+
+ def displacement(self, locs):
+ """
+ Compute displacement field at locations.
+ """
+ (npts, dim) = locs.shape
+ disp = numpy.zeros( (npts, 3), dtype=numpy.float64)
+ disp[:,0] = exx*locs[:,0] + exy*locs[:,1] + exz*locs[:,2]
+ disp[:,1] = eyy*locs[:,1] + exy*locs[:,0] + eyz*locs[:,2]
+ disp[:,2] = ezz*locs[:,2] + exz*locs[:,0] + eyz*locs[:,1]
+ return disp
+
+
+ def strain(self, locs):
+ """
+ Compute strain field at locations.
+ """
+ (npts, dim) = locs.shape
+ strain = numpy.zeros( (npts, 6), dtype=numpy.float64)
+ strain[:,0] = exx
+ strain[:,1] = eyy
+ strain[:,2] = ezz
+ strain[:,3] = exy
+ strain[:,4] = eyz
+ strain[:,5] = exz
+ return strain
+
+
+ def stress(self, locs):
+ """
+ Compute stress field at locations.
+ """
+ (npts, dim) = locs.shape
+ stress = numpy.zeros( (npts, 6), dtype=numpy.float64)
+ stress[:,0] = sxx
+ stress[:,1] = syy
+ stress[:,2] = szz
+ stress[:,3] = sxy
+ stress[:,4] = syz
+ stress[:,5] = sxz
+ return stress
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation.cfg 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,142 @@
+# -*- Python -*-
+[pylithapp]
+
+[pylithapp.launcher] # WARNING: THIS IS NOT PORTABLE
+command = mpirun -np ${nodes}
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+#timedependent = 1
+#implicit = 1
+#petsc = 1
+#solverlinear = 1
+#meshiocubit = 1
+#implicitelasticity = 1
+#quadrature3d = 1
+#fiatlagrange = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = mesh.exo
+use_nodeset_names = True
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+bc = [x_neg,x_pos]
+interfaces = [fault]
+
+[pylithapp.timedependent.formulation]
+output = [domain,subdomain]
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic,viscoelastic]
+materials.elastic = pylith.materials.ElasticIsotropic3D
+materials.viscoelastic = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.degree= 2
+
+[pylithapp.timedependent.materials.viscoelastic]
+label = Elastic material 2
+id = 2
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.degree = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0, 1, 2]
+label = face_xpos
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC +x edge
+db_initial.values = [displacement-x,displacement-y,displacement-z]
+db_initial.data = [0.0*m,+0.5*m,0.0*m]
+
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0, 1, 2]
+label = face_xneg
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Dirichlet BC -x edge
+db_initial.values = [displacement-x,displacement-y,displacement-z]
+db_initial.data = [0.0*m,-0.5*m,0.0*m]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault]
+id = 100
+label = fault
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.degree= 2
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = dislocation_slipLL.spatialdb
+slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+pc_type = asm
+
+# Change the preconditioner settings.
+sub_pc_factor_shift_type = none
+
+ksp_rtol = 1.0e-8
+ksp_max_it = 100
+ksp_gmres_restart = 50
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+#log_summary = true
+# start_in_debugger = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain.writer]
+filename = dislocation.vtk
+
+[pylithapp.problem.formulation.output.subdomain]
+label = face_zpos
+writer.filename = dislocation-groundsurf.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = dislocation-elastic.vtk
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = dislocation-viscoelastic.vtk
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_slipLL.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_slipLL.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_slipLL.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+ value-units = m m m
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 1.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_sliptime.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_sliptime.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/dislocation_sliptime.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 1
+ value-names = slip-time
+ value-units = s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 -1.0
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/geometry.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/geometry.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/geometry.jou 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,45 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 8 km x 8 km x 6 km
+# -4 km <= x <= 4 km
+# -4 km <= y <= 4 km
+# -6 km <= z <= 0 km
+reset
+brick x 8000 y 8000 z 6000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -3000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane xplane offset +2000
+surface 8 name "fault_surface_b"
+create planar surface with plane zplane offset -2000
+surface 9 name "material_interface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+webcut volume 1 with plane surface fault_surface_b
+webcut volume 1 with plane surface material_interface
+webcut volume 5 with plane surface material_interface
+webcut volume 6 with plane surface material_interface
+volume 1 name "elastic_xpos"
+volume 6 name "elastic_mid"
+volume 5 name "elastic_xneg"
+volume 7 name "visco_xpos"
+volume 9 name "visco_mid"
+volume 8 name "visco_xneg"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/matprops.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/matprops.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/matprops.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = density vs vp
+ value-units = kg/m**3 m/s m/s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2500.0 3000.0 5291.502622129181
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/mesh.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/mesh.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/mesh.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/mesh.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/mesh.jou 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,113 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 1000
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+volume all scheme tetmesh
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 5 6
+block 1 name "elastic"
+block 1 element type tetra10
+block 2 volume 7 9 8
+block 2 name "viscoelastic"
+block 2 element type tetra10
+
+# ----------------------------------------------------------------------
+# Create nodesets for faults
+# ----------------------------------------------------------------------
+group "fault" add node in surface fault_surface
+group "fault" add node in fault_surface at A
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+group "fault_b" add node in surface fault_surface_b
+group "fault_b" add node in surface fault_surface_b at A
+nodeset 11 group fault_b
+nodeset 11 name "fault_b"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 31
+group "face_xpos" add node in surface 39
+nodeset 20 group face_xpos
+nodeset 20 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 41
+group "face_xneg" add node in surface 49
+nodeset 21 group face_xneg
+nodeset 21 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 32
+group "face_ypos" add node in surface 44
+group "face_ypos" add node in surface 54
+group "face_ypos" add node in surface 38
+group "face_ypos" add node in surface 46
+group "face_ypos" add node in surface 56
+nodeset 22 group face_ypos
+nodeset 22 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 34
+group "face_yneg" add node in surface 42
+group "face_yneg" add node in surface 52
+group "face_yneg" add node in surface 36
+group "face_yneg" add node in surface 48
+group "face_yneg" add node in surface 58
+nodeset 23 group face_yneg
+nodeset 23 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 18
+group "face_zpos" add node in surface 27
+group "face_zpos" add node in surface 23
+nodeset 24 group face_zpos
+nodeset 24 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 17
+group "face_zneg" add node in surface 28
+group "face_zneg" add node in surface 21
+nodeset 25 group face_zneg
+nodeset 25 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 26 group face_zneg_nofault
+nodeset 26 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "mesh.exo" dimension 3 overwrite
+
+
+
+
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/testpylith.py
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/testpylith.py (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/testpylith.py 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,51 @@
+#!/usr/bin/env python
+#
+# ======================================================================
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ======================================================================
+#
+
+__requires__ = "PyLith"
+
+import unittest
+
+def suite():
+ """
+ Create test suite.
+ """
+ suite = unittest.TestSuite()
+
+ from TestLgDeformRigidBody import TestRigidBody
+ suite.addTest(unittest.makeSuite(TestRigidBody))
+
+ from TestLgDeformTraction import TestTraction
+ suite.addTest(unittest.makeSuite(TestTraction))
+
+ return suite
+
+
+def main():
+ """
+ Run test suite.
+ """
+ unittest.TextTestRunner(verbosity=2).run(suite())
+ return
+
+
+# ----------------------------------------------------------------------
+if __name__ == '__main__':
+ main()
+
+
+# End of file
Property changes on: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet10/testpylith.py
___________________________________________________________________
Name: svn:executable
+ *
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/Makefile.am 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,70 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+TESTS = testpylith.py
+
+check_SCRIPTS = testpylith.py
+
+dist_noinst_PYTHON = \
+ TestHex8.py \
+ TestLgDeformRigidBody.py \
+ rigidbody_soln.py \
+ rigidbody_gendb.py \
+ TestLgDeformTraction.py \
+ lgdeformtraction_soln.py \
+ testpylith.py
+
+dist_noinst_DATA = \
+ geometry.jou \
+ mesh.jou \
+ mesh.exo \
+ matprops.spatialdb \
+ lgdeformrigidbody.cfg \
+ lgdeformtraction.cfg
+
+noinst_TMP = \
+ rigidbody_disp.spatialdb \
+ lgdeformrigidbody_t0000000.vtk \
+ lgdeformrigidbody-elastic_info.vtk \
+ lgdeformrigidbody-elastic_t0000000.vtk \
+ lgdeformrigidbody-viscoelastic_info.vtk \
+ lgdeformrigidbody-viscoelastic_t0000000.vtk \
+ lgdeformtraction_t0000000.vtk \
+ lgdeformtraction-elastic_info.vtk \
+ lgdeformtraction-elastic_t0000000.vtk \
+ lgdeformtraction-viscoelastic_info.vtk \
+ lgdeformtraction-viscoelastic_t0000000.vtk
+
+
+TESTS_ENVIRONMENT = $(PYTHON)
+
+
+# 'export' the input files by performing a mock install
+export_datadir = $(top_builddir)/tests/3dnew/hex8
+export-data: $(dist_noinst_PYTHON) $(dist_noinst_DATA)
+ for f in $(dist_noinst_PYTHON) $(dist_noinst_DATA); do $(install_sh_DATA) $(srcdir)/$$f $(export_datadir); done
+
+clean-data:
+ if [ "X$(top_srcdir)" != "X$(top_builddir)" ]; then for f in $(dist_noinst_PYTHON) $(dist_noinst_DATA) $(noinst_TMP); do $(RM) $(RM_FLAGS) $(export_datadir)/$$f; done; fi
+
+
+BUILT_SOURCES = export-data
+clean-local: clean-data
+CLEANFILES = *.pyc
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axial_disp.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axial_disp.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axial_disp.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,580 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = displacement-x displacement-y displacement-z
+ value-units = m m m
+ num-locs = 567
+ data-dim = 3
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1
+ space-dim = 3
+}
+}
+ -4.000000e+03 -4.000000e+03 -6.000000e+03 -7.037037e-01 1.851852e-01 2.777778e-01
+ -3.000000e+03 -4.000000e+03 -6.000000e+03 -5.277778e-01 1.851852e-01 2.777778e-01
+ -2.000000e+03 -4.000000e+03 -6.000000e+03 -3.518519e-01 1.851852e-01 2.777778e-01
+ -1.000000e+03 -4.000000e+03 -6.000000e+03 -1.759259e-01 1.851852e-01 2.777778e-01
+ 0.000000e+00 -4.000000e+03 -6.000000e+03 0.000000e+00 1.851852e-01 2.777778e-01
+ 1.000000e+03 -4.000000e+03 -6.000000e+03 1.759259e-01 1.851852e-01 2.777778e-01
+ 2.000000e+03 -4.000000e+03 -6.000000e+03 3.518519e-01 1.851852e-01 2.777778e-01
+ 3.000000e+03 -4.000000e+03 -6.000000e+03 5.277778e-01 1.851852e-01 2.777778e-01
+ 4.000000e+03 -4.000000e+03 -6.000000e+03 7.037037e-01 1.851852e-01 2.777778e-01
+ -4.000000e+03 -3.000000e+03 -6.000000e+03 -7.037037e-01 1.388889e-01 2.777778e-01
+ -3.000000e+03 -3.000000e+03 -6.000000e+03 -5.277778e-01 1.388889e-01 2.777778e-01
+ -2.000000e+03 -3.000000e+03 -6.000000e+03 -3.518519e-01 1.388889e-01 2.777778e-01
+ -1.000000e+03 -3.000000e+03 -6.000000e+03 -1.759259e-01 1.388889e-01 2.777778e-01
+ 0.000000e+00 -3.000000e+03 -6.000000e+03 0.000000e+00 1.388889e-01 2.777778e-01
+ 1.000000e+03 -3.000000e+03 -6.000000e+03 1.759259e-01 1.388889e-01 2.777778e-01
+ 2.000000e+03 -3.000000e+03 -6.000000e+03 3.518519e-01 1.388889e-01 2.777778e-01
+ 3.000000e+03 -3.000000e+03 -6.000000e+03 5.277778e-01 1.388889e-01 2.777778e-01
+ 4.000000e+03 -3.000000e+03 -6.000000e+03 7.037037e-01 1.388889e-01 2.777778e-01
+ -4.000000e+03 -2.000000e+03 -6.000000e+03 -7.037037e-01 9.259259e-02 2.777778e-01
+ -3.000000e+03 -2.000000e+03 -6.000000e+03 -5.277778e-01 9.259259e-02 2.777778e-01
+ -2.000000e+03 -2.000000e+03 -6.000000e+03 -3.518519e-01 9.259259e-02 2.777778e-01
+ -1.000000e+03 -2.000000e+03 -6.000000e+03 -1.759259e-01 9.259259e-02 2.777778e-01
+ 0.000000e+00 -2.000000e+03 -6.000000e+03 0.000000e+00 9.259259e-02 2.777778e-01
+ 1.000000e+03 -2.000000e+03 -6.000000e+03 1.759259e-01 9.259259e-02 2.777778e-01
+ 2.000000e+03 -2.000000e+03 -6.000000e+03 3.518519e-01 9.259259e-02 2.777778e-01
+ 3.000000e+03 -2.000000e+03 -6.000000e+03 5.277778e-01 9.259259e-02 2.777778e-01
+ 4.000000e+03 -2.000000e+03 -6.000000e+03 7.037037e-01 9.259259e-02 2.777778e-01
+ -4.000000e+03 -1.000000e+03 -6.000000e+03 -7.037037e-01 4.629630e-02 2.777778e-01
+ -3.000000e+03 -1.000000e+03 -6.000000e+03 -5.277778e-01 4.629630e-02 2.777778e-01
+ -2.000000e+03 -1.000000e+03 -6.000000e+03 -3.518519e-01 4.629630e-02 2.777778e-01
+ -1.000000e+03 -1.000000e+03 -6.000000e+03 -1.759259e-01 4.629630e-02 2.777778e-01
+ 0.000000e+00 -1.000000e+03 -6.000000e+03 0.000000e+00 4.629630e-02 2.777778e-01
+ 1.000000e+03 -1.000000e+03 -6.000000e+03 1.759259e-01 4.629630e-02 2.777778e-01
+ 2.000000e+03 -1.000000e+03 -6.000000e+03 3.518519e-01 4.629630e-02 2.777778e-01
+ 3.000000e+03 -1.000000e+03 -6.000000e+03 5.277778e-01 4.629630e-02 2.777778e-01
+ 4.000000e+03 -1.000000e+03 -6.000000e+03 7.037037e-01 4.629630e-02 2.777778e-01
+ -4.000000e+03 0.000000e+00 -6.000000e+03 -7.037037e-01 -0.000000e+00 2.777778e-01
+ -3.000000e+03 0.000000e+00 -6.000000e+03 -5.277778e-01 -0.000000e+00 2.777778e-01
+ -2.000000e+03 0.000000e+00 -6.000000e+03 -3.518519e-01 -0.000000e+00 2.777778e-01
+ -1.000000e+03 0.000000e+00 -6.000000e+03 -1.759259e-01 -0.000000e+00 2.777778e-01
+ 0.000000e+00 0.000000e+00 -6.000000e+03 0.000000e+00 0.000000e+00 2.777778e-01
+ 1.000000e+03 0.000000e+00 -6.000000e+03 1.759259e-01 0.000000e+00 2.777778e-01
+ 2.000000e+03 0.000000e+00 -6.000000e+03 3.518519e-01 0.000000e+00 2.777778e-01
+ 3.000000e+03 0.000000e+00 -6.000000e+03 5.277778e-01 0.000000e+00 2.777778e-01
+ 4.000000e+03 0.000000e+00 -6.000000e+03 7.037037e-01 0.000000e+00 2.777778e-01
+ -4.000000e+03 1.000000e+03 -6.000000e+03 -7.037037e-01 -4.629630e-02 2.777778e-01
+ -3.000000e+03 1.000000e+03 -6.000000e+03 -5.277778e-01 -4.629630e-02 2.777778e-01
+ -2.000000e+03 1.000000e+03 -6.000000e+03 -3.518519e-01 -4.629630e-02 2.777778e-01
+ -1.000000e+03 1.000000e+03 -6.000000e+03 -1.759259e-01 -4.629630e-02 2.777778e-01
+ 0.000000e+00 1.000000e+03 -6.000000e+03 0.000000e+00 -4.629630e-02 2.777778e-01
+ 1.000000e+03 1.000000e+03 -6.000000e+03 1.759259e-01 -4.629630e-02 2.777778e-01
+ 2.000000e+03 1.000000e+03 -6.000000e+03 3.518519e-01 -4.629630e-02 2.777778e-01
+ 3.000000e+03 1.000000e+03 -6.000000e+03 5.277778e-01 -4.629630e-02 2.777778e-01
+ 4.000000e+03 1.000000e+03 -6.000000e+03 7.037037e-01 -4.629630e-02 2.777778e-01
+ -4.000000e+03 2.000000e+03 -6.000000e+03 -7.037037e-01 -9.259259e-02 2.777778e-01
+ -3.000000e+03 2.000000e+03 -6.000000e+03 -5.277778e-01 -9.259259e-02 2.777778e-01
+ -2.000000e+03 2.000000e+03 -6.000000e+03 -3.518519e-01 -9.259259e-02 2.777778e-01
+ -1.000000e+03 2.000000e+03 -6.000000e+03 -1.759259e-01 -9.259259e-02 2.777778e-01
+ 0.000000e+00 2.000000e+03 -6.000000e+03 0.000000e+00 -9.259259e-02 2.777778e-01
+ 1.000000e+03 2.000000e+03 -6.000000e+03 1.759259e-01 -9.259259e-02 2.777778e-01
+ 2.000000e+03 2.000000e+03 -6.000000e+03 3.518519e-01 -9.259259e-02 2.777778e-01
+ 3.000000e+03 2.000000e+03 -6.000000e+03 5.277778e-01 -9.259259e-02 2.777778e-01
+ 4.000000e+03 2.000000e+03 -6.000000e+03 7.037037e-01 -9.259259e-02 2.777778e-01
+ -4.000000e+03 3.000000e+03 -6.000000e+03 -7.037037e-01 -1.388889e-01 2.777778e-01
+ -3.000000e+03 3.000000e+03 -6.000000e+03 -5.277778e-01 -1.388889e-01 2.777778e-01
+ -2.000000e+03 3.000000e+03 -6.000000e+03 -3.518519e-01 -1.388889e-01 2.777778e-01
+ -1.000000e+03 3.000000e+03 -6.000000e+03 -1.759259e-01 -1.388889e-01 2.777778e-01
+ 0.000000e+00 3.000000e+03 -6.000000e+03 0.000000e+00 -1.388889e-01 2.777778e-01
+ 1.000000e+03 3.000000e+03 -6.000000e+03 1.759259e-01 -1.388889e-01 2.777778e-01
+ 2.000000e+03 3.000000e+03 -6.000000e+03 3.518519e-01 -1.388889e-01 2.777778e-01
+ 3.000000e+03 3.000000e+03 -6.000000e+03 5.277778e-01 -1.388889e-01 2.777778e-01
+ 4.000000e+03 3.000000e+03 -6.000000e+03 7.037037e-01 -1.388889e-01 2.777778e-01
+ -4.000000e+03 4.000000e+03 -6.000000e+03 -7.037037e-01 -1.851852e-01 2.777778e-01
+ -3.000000e+03 4.000000e+03 -6.000000e+03 -5.277778e-01 -1.851852e-01 2.777778e-01
+ -2.000000e+03 4.000000e+03 -6.000000e+03 -3.518519e-01 -1.851852e-01 2.777778e-01
+ -1.000000e+03 4.000000e+03 -6.000000e+03 -1.759259e-01 -1.851852e-01 2.777778e-01
+ 0.000000e+00 4.000000e+03 -6.000000e+03 0.000000e+00 -1.851852e-01 2.777778e-01
+ 1.000000e+03 4.000000e+03 -6.000000e+03 1.759259e-01 -1.851852e-01 2.777778e-01
+ 2.000000e+03 4.000000e+03 -6.000000e+03 3.518519e-01 -1.851852e-01 2.777778e-01
+ 3.000000e+03 4.000000e+03 -6.000000e+03 5.277778e-01 -1.851852e-01 2.777778e-01
+ 4.000000e+03 4.000000e+03 -6.000000e+03 7.037037e-01 -1.851852e-01 2.777778e-01
+ -4.000000e+03 -4.000000e+03 -5.000000e+03 -7.037037e-01 1.851852e-01 2.314815e-01
+ -3.000000e+03 -4.000000e+03 -5.000000e+03 -5.277778e-01 1.851852e-01 2.314815e-01
+ -2.000000e+03 -4.000000e+03 -5.000000e+03 -3.518519e-01 1.851852e-01 2.314815e-01
+ -1.000000e+03 -4.000000e+03 -5.000000e+03 -1.759259e-01 1.851852e-01 2.314815e-01
+ 0.000000e+00 -4.000000e+03 -5.000000e+03 0.000000e+00 1.851852e-01 2.314815e-01
+ 1.000000e+03 -4.000000e+03 -5.000000e+03 1.759259e-01 1.851852e-01 2.314815e-01
+ 2.000000e+03 -4.000000e+03 -5.000000e+03 3.518519e-01 1.851852e-01 2.314815e-01
+ 3.000000e+03 -4.000000e+03 -5.000000e+03 5.277778e-01 1.851852e-01 2.314815e-01
+ 4.000000e+03 -4.000000e+03 -5.000000e+03 7.037037e-01 1.851852e-01 2.314815e-01
+ -4.000000e+03 -3.000000e+03 -5.000000e+03 -7.037037e-01 1.388889e-01 2.314815e-01
+ -3.000000e+03 -3.000000e+03 -5.000000e+03 -5.277778e-01 1.388889e-01 2.314815e-01
+ -2.000000e+03 -3.000000e+03 -5.000000e+03 -3.518519e-01 1.388889e-01 2.314815e-01
+ -1.000000e+03 -3.000000e+03 -5.000000e+03 -1.759259e-01 1.388889e-01 2.314815e-01
+ 0.000000e+00 -3.000000e+03 -5.000000e+03 0.000000e+00 1.388889e-01 2.314815e-01
+ 1.000000e+03 -3.000000e+03 -5.000000e+03 1.759259e-01 1.388889e-01 2.314815e-01
+ 2.000000e+03 -3.000000e+03 -5.000000e+03 3.518519e-01 1.388889e-01 2.314815e-01
+ 3.000000e+03 -3.000000e+03 -5.000000e+03 5.277778e-01 1.388889e-01 2.314815e-01
+ 4.000000e+03 -3.000000e+03 -5.000000e+03 7.037037e-01 1.388889e-01 2.314815e-01
+ -4.000000e+03 -2.000000e+03 -5.000000e+03 -7.037037e-01 9.259259e-02 2.314815e-01
+ -3.000000e+03 -2.000000e+03 -5.000000e+03 -5.277778e-01 9.259259e-02 2.314815e-01
+ -2.000000e+03 -2.000000e+03 -5.000000e+03 -3.518519e-01 9.259259e-02 2.314815e-01
+ -1.000000e+03 -2.000000e+03 -5.000000e+03 -1.759259e-01 9.259259e-02 2.314815e-01
+ 0.000000e+00 -2.000000e+03 -5.000000e+03 0.000000e+00 9.259259e-02 2.314815e-01
+ 1.000000e+03 -2.000000e+03 -5.000000e+03 1.759259e-01 9.259259e-02 2.314815e-01
+ 2.000000e+03 -2.000000e+03 -5.000000e+03 3.518519e-01 9.259259e-02 2.314815e-01
+ 3.000000e+03 -2.000000e+03 -5.000000e+03 5.277778e-01 9.259259e-02 2.314815e-01
+ 4.000000e+03 -2.000000e+03 -5.000000e+03 7.037037e-01 9.259259e-02 2.314815e-01
+ -4.000000e+03 -1.000000e+03 -5.000000e+03 -7.037037e-01 4.629630e-02 2.314815e-01
+ -3.000000e+03 -1.000000e+03 -5.000000e+03 -5.277778e-01 4.629630e-02 2.314815e-01
+ -2.000000e+03 -1.000000e+03 -5.000000e+03 -3.518519e-01 4.629630e-02 2.314815e-01
+ -1.000000e+03 -1.000000e+03 -5.000000e+03 -1.759259e-01 4.629630e-02 2.314815e-01
+ 0.000000e+00 -1.000000e+03 -5.000000e+03 0.000000e+00 4.629630e-02 2.314815e-01
+ 1.000000e+03 -1.000000e+03 -5.000000e+03 1.759259e-01 4.629630e-02 2.314815e-01
+ 2.000000e+03 -1.000000e+03 -5.000000e+03 3.518519e-01 4.629630e-02 2.314815e-01
+ 3.000000e+03 -1.000000e+03 -5.000000e+03 5.277778e-01 4.629630e-02 2.314815e-01
+ 4.000000e+03 -1.000000e+03 -5.000000e+03 7.037037e-01 4.629630e-02 2.314815e-01
+ -4.000000e+03 0.000000e+00 -5.000000e+03 -7.037037e-01 -0.000000e+00 2.314815e-01
+ -3.000000e+03 0.000000e+00 -5.000000e+03 -5.277778e-01 -0.000000e+00 2.314815e-01
+ -2.000000e+03 0.000000e+00 -5.000000e+03 -3.518519e-01 -0.000000e+00 2.314815e-01
+ -1.000000e+03 0.000000e+00 -5.000000e+03 -1.759259e-01 -0.000000e+00 2.314815e-01
+ 0.000000e+00 0.000000e+00 -5.000000e+03 0.000000e+00 0.000000e+00 2.314815e-01
+ 1.000000e+03 0.000000e+00 -5.000000e+03 1.759259e-01 0.000000e+00 2.314815e-01
+ 2.000000e+03 0.000000e+00 -5.000000e+03 3.518519e-01 0.000000e+00 2.314815e-01
+ 3.000000e+03 0.000000e+00 -5.000000e+03 5.277778e-01 0.000000e+00 2.314815e-01
+ 4.000000e+03 0.000000e+00 -5.000000e+03 7.037037e-01 0.000000e+00 2.314815e-01
+ -4.000000e+03 1.000000e+03 -5.000000e+03 -7.037037e-01 -4.629630e-02 2.314815e-01
+ -3.000000e+03 1.000000e+03 -5.000000e+03 -5.277778e-01 -4.629630e-02 2.314815e-01
+ -2.000000e+03 1.000000e+03 -5.000000e+03 -3.518519e-01 -4.629630e-02 2.314815e-01
+ -1.000000e+03 1.000000e+03 -5.000000e+03 -1.759259e-01 -4.629630e-02 2.314815e-01
+ 0.000000e+00 1.000000e+03 -5.000000e+03 0.000000e+00 -4.629630e-02 2.314815e-01
+ 1.000000e+03 1.000000e+03 -5.000000e+03 1.759259e-01 -4.629630e-02 2.314815e-01
+ 2.000000e+03 1.000000e+03 -5.000000e+03 3.518519e-01 -4.629630e-02 2.314815e-01
+ 3.000000e+03 1.000000e+03 -5.000000e+03 5.277778e-01 -4.629630e-02 2.314815e-01
+ 4.000000e+03 1.000000e+03 -5.000000e+03 7.037037e-01 -4.629630e-02 2.314815e-01
+ -4.000000e+03 2.000000e+03 -5.000000e+03 -7.037037e-01 -9.259259e-02 2.314815e-01
+ -3.000000e+03 2.000000e+03 -5.000000e+03 -5.277778e-01 -9.259259e-02 2.314815e-01
+ -2.000000e+03 2.000000e+03 -5.000000e+03 -3.518519e-01 -9.259259e-02 2.314815e-01
+ -1.000000e+03 2.000000e+03 -5.000000e+03 -1.759259e-01 -9.259259e-02 2.314815e-01
+ 0.000000e+00 2.000000e+03 -5.000000e+03 0.000000e+00 -9.259259e-02 2.314815e-01
+ 1.000000e+03 2.000000e+03 -5.000000e+03 1.759259e-01 -9.259259e-02 2.314815e-01
+ 2.000000e+03 2.000000e+03 -5.000000e+03 3.518519e-01 -9.259259e-02 2.314815e-01
+ 3.000000e+03 2.000000e+03 -5.000000e+03 5.277778e-01 -9.259259e-02 2.314815e-01
+ 4.000000e+03 2.000000e+03 -5.000000e+03 7.037037e-01 -9.259259e-02 2.314815e-01
+ -4.000000e+03 3.000000e+03 -5.000000e+03 -7.037037e-01 -1.388889e-01 2.314815e-01
+ -3.000000e+03 3.000000e+03 -5.000000e+03 -5.277778e-01 -1.388889e-01 2.314815e-01
+ -2.000000e+03 3.000000e+03 -5.000000e+03 -3.518519e-01 -1.388889e-01 2.314815e-01
+ -1.000000e+03 3.000000e+03 -5.000000e+03 -1.759259e-01 -1.388889e-01 2.314815e-01
+ 0.000000e+00 3.000000e+03 -5.000000e+03 0.000000e+00 -1.388889e-01 2.314815e-01
+ 1.000000e+03 3.000000e+03 -5.000000e+03 1.759259e-01 -1.388889e-01 2.314815e-01
+ 2.000000e+03 3.000000e+03 -5.000000e+03 3.518519e-01 -1.388889e-01 2.314815e-01
+ 3.000000e+03 3.000000e+03 -5.000000e+03 5.277778e-01 -1.388889e-01 2.314815e-01
+ 4.000000e+03 3.000000e+03 -5.000000e+03 7.037037e-01 -1.388889e-01 2.314815e-01
+ -4.000000e+03 4.000000e+03 -5.000000e+03 -7.037037e-01 -1.851852e-01 2.314815e-01
+ -3.000000e+03 4.000000e+03 -5.000000e+03 -5.277778e-01 -1.851852e-01 2.314815e-01
+ -2.000000e+03 4.000000e+03 -5.000000e+03 -3.518519e-01 -1.851852e-01 2.314815e-01
+ -1.000000e+03 4.000000e+03 -5.000000e+03 -1.759259e-01 -1.851852e-01 2.314815e-01
+ 0.000000e+00 4.000000e+03 -5.000000e+03 0.000000e+00 -1.851852e-01 2.314815e-01
+ 1.000000e+03 4.000000e+03 -5.000000e+03 1.759259e-01 -1.851852e-01 2.314815e-01
+ 2.000000e+03 4.000000e+03 -5.000000e+03 3.518519e-01 -1.851852e-01 2.314815e-01
+ 3.000000e+03 4.000000e+03 -5.000000e+03 5.277778e-01 -1.851852e-01 2.314815e-01
+ 4.000000e+03 4.000000e+03 -5.000000e+03 7.037037e-01 -1.851852e-01 2.314815e-01
+ -4.000000e+03 -4.000000e+03 -4.000000e+03 -7.037037e-01 1.851852e-01 1.851852e-01
+ -3.000000e+03 -4.000000e+03 -4.000000e+03 -5.277778e-01 1.851852e-01 1.851852e-01
+ -2.000000e+03 -4.000000e+03 -4.000000e+03 -3.518519e-01 1.851852e-01 1.851852e-01
+ -1.000000e+03 -4.000000e+03 -4.000000e+03 -1.759259e-01 1.851852e-01 1.851852e-01
+ 0.000000e+00 -4.000000e+03 -4.000000e+03 0.000000e+00 1.851852e-01 1.851852e-01
+ 1.000000e+03 -4.000000e+03 -4.000000e+03 1.759259e-01 1.851852e-01 1.851852e-01
+ 2.000000e+03 -4.000000e+03 -4.000000e+03 3.518519e-01 1.851852e-01 1.851852e-01
+ 3.000000e+03 -4.000000e+03 -4.000000e+03 5.277778e-01 1.851852e-01 1.851852e-01
+ 4.000000e+03 -4.000000e+03 -4.000000e+03 7.037037e-01 1.851852e-01 1.851852e-01
+ -4.000000e+03 -3.000000e+03 -4.000000e+03 -7.037037e-01 1.388889e-01 1.851852e-01
+ -3.000000e+03 -3.000000e+03 -4.000000e+03 -5.277778e-01 1.388889e-01 1.851852e-01
+ -2.000000e+03 -3.000000e+03 -4.000000e+03 -3.518519e-01 1.388889e-01 1.851852e-01
+ -1.000000e+03 -3.000000e+03 -4.000000e+03 -1.759259e-01 1.388889e-01 1.851852e-01
+ 0.000000e+00 -3.000000e+03 -4.000000e+03 0.000000e+00 1.388889e-01 1.851852e-01
+ 1.000000e+03 -3.000000e+03 -4.000000e+03 1.759259e-01 1.388889e-01 1.851852e-01
+ 2.000000e+03 -3.000000e+03 -4.000000e+03 3.518519e-01 1.388889e-01 1.851852e-01
+ 3.000000e+03 -3.000000e+03 -4.000000e+03 5.277778e-01 1.388889e-01 1.851852e-01
+ 4.000000e+03 -3.000000e+03 -4.000000e+03 7.037037e-01 1.388889e-01 1.851852e-01
+ -4.000000e+03 -2.000000e+03 -4.000000e+03 -7.037037e-01 9.259259e-02 1.851852e-01
+ -3.000000e+03 -2.000000e+03 -4.000000e+03 -5.277778e-01 9.259259e-02 1.851852e-01
+ -2.000000e+03 -2.000000e+03 -4.000000e+03 -3.518519e-01 9.259259e-02 1.851852e-01
+ -1.000000e+03 -2.000000e+03 -4.000000e+03 -1.759259e-01 9.259259e-02 1.851852e-01
+ 0.000000e+00 -2.000000e+03 -4.000000e+03 0.000000e+00 9.259259e-02 1.851852e-01
+ 1.000000e+03 -2.000000e+03 -4.000000e+03 1.759259e-01 9.259259e-02 1.851852e-01
+ 2.000000e+03 -2.000000e+03 -4.000000e+03 3.518519e-01 9.259259e-02 1.851852e-01
+ 3.000000e+03 -2.000000e+03 -4.000000e+03 5.277778e-01 9.259259e-02 1.851852e-01
+ 4.000000e+03 -2.000000e+03 -4.000000e+03 7.037037e-01 9.259259e-02 1.851852e-01
+ -4.000000e+03 -1.000000e+03 -4.000000e+03 -7.037037e-01 4.629630e-02 1.851852e-01
+ -3.000000e+03 -1.000000e+03 -4.000000e+03 -5.277778e-01 4.629630e-02 1.851852e-01
+ -2.000000e+03 -1.000000e+03 -4.000000e+03 -3.518519e-01 4.629630e-02 1.851852e-01
+ -1.000000e+03 -1.000000e+03 -4.000000e+03 -1.759259e-01 4.629630e-02 1.851852e-01
+ 0.000000e+00 -1.000000e+03 -4.000000e+03 0.000000e+00 4.629630e-02 1.851852e-01
+ 1.000000e+03 -1.000000e+03 -4.000000e+03 1.759259e-01 4.629630e-02 1.851852e-01
+ 2.000000e+03 -1.000000e+03 -4.000000e+03 3.518519e-01 4.629630e-02 1.851852e-01
+ 3.000000e+03 -1.000000e+03 -4.000000e+03 5.277778e-01 4.629630e-02 1.851852e-01
+ 4.000000e+03 -1.000000e+03 -4.000000e+03 7.037037e-01 4.629630e-02 1.851852e-01
+ -4.000000e+03 0.000000e+00 -4.000000e+03 -7.037037e-01 -0.000000e+00 1.851852e-01
+ -3.000000e+03 0.000000e+00 -4.000000e+03 -5.277778e-01 -0.000000e+00 1.851852e-01
+ -2.000000e+03 0.000000e+00 -4.000000e+03 -3.518519e-01 -0.000000e+00 1.851852e-01
+ -1.000000e+03 0.000000e+00 -4.000000e+03 -1.759259e-01 -0.000000e+00 1.851852e-01
+ 0.000000e+00 0.000000e+00 -4.000000e+03 0.000000e+00 0.000000e+00 1.851852e-01
+ 1.000000e+03 0.000000e+00 -4.000000e+03 1.759259e-01 0.000000e+00 1.851852e-01
+ 2.000000e+03 0.000000e+00 -4.000000e+03 3.518519e-01 0.000000e+00 1.851852e-01
+ 3.000000e+03 0.000000e+00 -4.000000e+03 5.277778e-01 0.000000e+00 1.851852e-01
+ 4.000000e+03 0.000000e+00 -4.000000e+03 7.037037e-01 0.000000e+00 1.851852e-01
+ -4.000000e+03 1.000000e+03 -4.000000e+03 -7.037037e-01 -4.629630e-02 1.851852e-01
+ -3.000000e+03 1.000000e+03 -4.000000e+03 -5.277778e-01 -4.629630e-02 1.851852e-01
+ -2.000000e+03 1.000000e+03 -4.000000e+03 -3.518519e-01 -4.629630e-02 1.851852e-01
+ -1.000000e+03 1.000000e+03 -4.000000e+03 -1.759259e-01 -4.629630e-02 1.851852e-01
+ 0.000000e+00 1.000000e+03 -4.000000e+03 0.000000e+00 -4.629630e-02 1.851852e-01
+ 1.000000e+03 1.000000e+03 -4.000000e+03 1.759259e-01 -4.629630e-02 1.851852e-01
+ 2.000000e+03 1.000000e+03 -4.000000e+03 3.518519e-01 -4.629630e-02 1.851852e-01
+ 3.000000e+03 1.000000e+03 -4.000000e+03 5.277778e-01 -4.629630e-02 1.851852e-01
+ 4.000000e+03 1.000000e+03 -4.000000e+03 7.037037e-01 -4.629630e-02 1.851852e-01
+ -4.000000e+03 2.000000e+03 -4.000000e+03 -7.037037e-01 -9.259259e-02 1.851852e-01
+ -3.000000e+03 2.000000e+03 -4.000000e+03 -5.277778e-01 -9.259259e-02 1.851852e-01
+ -2.000000e+03 2.000000e+03 -4.000000e+03 -3.518519e-01 -9.259259e-02 1.851852e-01
+ -1.000000e+03 2.000000e+03 -4.000000e+03 -1.759259e-01 -9.259259e-02 1.851852e-01
+ 0.000000e+00 2.000000e+03 -4.000000e+03 0.000000e+00 -9.259259e-02 1.851852e-01
+ 1.000000e+03 2.000000e+03 -4.000000e+03 1.759259e-01 -9.259259e-02 1.851852e-01
+ 2.000000e+03 2.000000e+03 -4.000000e+03 3.518519e-01 -9.259259e-02 1.851852e-01
+ 3.000000e+03 2.000000e+03 -4.000000e+03 5.277778e-01 -9.259259e-02 1.851852e-01
+ 4.000000e+03 2.000000e+03 -4.000000e+03 7.037037e-01 -9.259259e-02 1.851852e-01
+ -4.000000e+03 3.000000e+03 -4.000000e+03 -7.037037e-01 -1.388889e-01 1.851852e-01
+ -3.000000e+03 3.000000e+03 -4.000000e+03 -5.277778e-01 -1.388889e-01 1.851852e-01
+ -2.000000e+03 3.000000e+03 -4.000000e+03 -3.518519e-01 -1.388889e-01 1.851852e-01
+ -1.000000e+03 3.000000e+03 -4.000000e+03 -1.759259e-01 -1.388889e-01 1.851852e-01
+ 0.000000e+00 3.000000e+03 -4.000000e+03 0.000000e+00 -1.388889e-01 1.851852e-01
+ 1.000000e+03 3.000000e+03 -4.000000e+03 1.759259e-01 -1.388889e-01 1.851852e-01
+ 2.000000e+03 3.000000e+03 -4.000000e+03 3.518519e-01 -1.388889e-01 1.851852e-01
+ 3.000000e+03 3.000000e+03 -4.000000e+03 5.277778e-01 -1.388889e-01 1.851852e-01
+ 4.000000e+03 3.000000e+03 -4.000000e+03 7.037037e-01 -1.388889e-01 1.851852e-01
+ -4.000000e+03 4.000000e+03 -4.000000e+03 -7.037037e-01 -1.851852e-01 1.851852e-01
+ -3.000000e+03 4.000000e+03 -4.000000e+03 -5.277778e-01 -1.851852e-01 1.851852e-01
+ -2.000000e+03 4.000000e+03 -4.000000e+03 -3.518519e-01 -1.851852e-01 1.851852e-01
+ -1.000000e+03 4.000000e+03 -4.000000e+03 -1.759259e-01 -1.851852e-01 1.851852e-01
+ 0.000000e+00 4.000000e+03 -4.000000e+03 0.000000e+00 -1.851852e-01 1.851852e-01
+ 1.000000e+03 4.000000e+03 -4.000000e+03 1.759259e-01 -1.851852e-01 1.851852e-01
+ 2.000000e+03 4.000000e+03 -4.000000e+03 3.518519e-01 -1.851852e-01 1.851852e-01
+ 3.000000e+03 4.000000e+03 -4.000000e+03 5.277778e-01 -1.851852e-01 1.851852e-01
+ 4.000000e+03 4.000000e+03 -4.000000e+03 7.037037e-01 -1.851852e-01 1.851852e-01
+ -4.000000e+03 -4.000000e+03 -3.000000e+03 -7.037037e-01 1.851852e-01 1.388889e-01
+ -3.000000e+03 -4.000000e+03 -3.000000e+03 -5.277778e-01 1.851852e-01 1.388889e-01
+ -2.000000e+03 -4.000000e+03 -3.000000e+03 -3.518519e-01 1.851852e-01 1.388889e-01
+ -1.000000e+03 -4.000000e+03 -3.000000e+03 -1.759259e-01 1.851852e-01 1.388889e-01
+ 0.000000e+00 -4.000000e+03 -3.000000e+03 0.000000e+00 1.851852e-01 1.388889e-01
+ 1.000000e+03 -4.000000e+03 -3.000000e+03 1.759259e-01 1.851852e-01 1.388889e-01
+ 2.000000e+03 -4.000000e+03 -3.000000e+03 3.518519e-01 1.851852e-01 1.388889e-01
+ 3.000000e+03 -4.000000e+03 -3.000000e+03 5.277778e-01 1.851852e-01 1.388889e-01
+ 4.000000e+03 -4.000000e+03 -3.000000e+03 7.037037e-01 1.851852e-01 1.388889e-01
+ -4.000000e+03 -3.000000e+03 -3.000000e+03 -7.037037e-01 1.388889e-01 1.388889e-01
+ -3.000000e+03 -3.000000e+03 -3.000000e+03 -5.277778e-01 1.388889e-01 1.388889e-01
+ -2.000000e+03 -3.000000e+03 -3.000000e+03 -3.518519e-01 1.388889e-01 1.388889e-01
+ -1.000000e+03 -3.000000e+03 -3.000000e+03 -1.759259e-01 1.388889e-01 1.388889e-01
+ 0.000000e+00 -3.000000e+03 -3.000000e+03 0.000000e+00 1.388889e-01 1.388889e-01
+ 1.000000e+03 -3.000000e+03 -3.000000e+03 1.759259e-01 1.388889e-01 1.388889e-01
+ 2.000000e+03 -3.000000e+03 -3.000000e+03 3.518519e-01 1.388889e-01 1.388889e-01
+ 3.000000e+03 -3.000000e+03 -3.000000e+03 5.277778e-01 1.388889e-01 1.388889e-01
+ 4.000000e+03 -3.000000e+03 -3.000000e+03 7.037037e-01 1.388889e-01 1.388889e-01
+ -4.000000e+03 -2.000000e+03 -3.000000e+03 -7.037037e-01 9.259259e-02 1.388889e-01
+ -3.000000e+03 -2.000000e+03 -3.000000e+03 -5.277778e-01 9.259259e-02 1.388889e-01
+ -2.000000e+03 -2.000000e+03 -3.000000e+03 -3.518519e-01 9.259259e-02 1.388889e-01
+ -1.000000e+03 -2.000000e+03 -3.000000e+03 -1.759259e-01 9.259259e-02 1.388889e-01
+ 0.000000e+00 -2.000000e+03 -3.000000e+03 0.000000e+00 9.259259e-02 1.388889e-01
+ 1.000000e+03 -2.000000e+03 -3.000000e+03 1.759259e-01 9.259259e-02 1.388889e-01
+ 2.000000e+03 -2.000000e+03 -3.000000e+03 3.518519e-01 9.259259e-02 1.388889e-01
+ 3.000000e+03 -2.000000e+03 -3.000000e+03 5.277778e-01 9.259259e-02 1.388889e-01
+ 4.000000e+03 -2.000000e+03 -3.000000e+03 7.037037e-01 9.259259e-02 1.388889e-01
+ -4.000000e+03 -1.000000e+03 -3.000000e+03 -7.037037e-01 4.629630e-02 1.388889e-01
+ -3.000000e+03 -1.000000e+03 -3.000000e+03 -5.277778e-01 4.629630e-02 1.388889e-01
+ -2.000000e+03 -1.000000e+03 -3.000000e+03 -3.518519e-01 4.629630e-02 1.388889e-01
+ -1.000000e+03 -1.000000e+03 -3.000000e+03 -1.759259e-01 4.629630e-02 1.388889e-01
+ 0.000000e+00 -1.000000e+03 -3.000000e+03 0.000000e+00 4.629630e-02 1.388889e-01
+ 1.000000e+03 -1.000000e+03 -3.000000e+03 1.759259e-01 4.629630e-02 1.388889e-01
+ 2.000000e+03 -1.000000e+03 -3.000000e+03 3.518519e-01 4.629630e-02 1.388889e-01
+ 3.000000e+03 -1.000000e+03 -3.000000e+03 5.277778e-01 4.629630e-02 1.388889e-01
+ 4.000000e+03 -1.000000e+03 -3.000000e+03 7.037037e-01 4.629630e-02 1.388889e-01
+ -4.000000e+03 0.000000e+00 -3.000000e+03 -7.037037e-01 -0.000000e+00 1.388889e-01
+ -3.000000e+03 0.000000e+00 -3.000000e+03 -5.277778e-01 -0.000000e+00 1.388889e-01
+ -2.000000e+03 0.000000e+00 -3.000000e+03 -3.518519e-01 -0.000000e+00 1.388889e-01
+ -1.000000e+03 0.000000e+00 -3.000000e+03 -1.759259e-01 -0.000000e+00 1.388889e-01
+ 0.000000e+00 0.000000e+00 -3.000000e+03 0.000000e+00 0.000000e+00 1.388889e-01
+ 1.000000e+03 0.000000e+00 -3.000000e+03 1.759259e-01 0.000000e+00 1.388889e-01
+ 2.000000e+03 0.000000e+00 -3.000000e+03 3.518519e-01 0.000000e+00 1.388889e-01
+ 3.000000e+03 0.000000e+00 -3.000000e+03 5.277778e-01 0.000000e+00 1.388889e-01
+ 4.000000e+03 0.000000e+00 -3.000000e+03 7.037037e-01 0.000000e+00 1.388889e-01
+ -4.000000e+03 1.000000e+03 -3.000000e+03 -7.037037e-01 -4.629630e-02 1.388889e-01
+ -3.000000e+03 1.000000e+03 -3.000000e+03 -5.277778e-01 -4.629630e-02 1.388889e-01
+ -2.000000e+03 1.000000e+03 -3.000000e+03 -3.518519e-01 -4.629630e-02 1.388889e-01
+ -1.000000e+03 1.000000e+03 -3.000000e+03 -1.759259e-01 -4.629630e-02 1.388889e-01
+ 0.000000e+00 1.000000e+03 -3.000000e+03 0.000000e+00 -4.629630e-02 1.388889e-01
+ 1.000000e+03 1.000000e+03 -3.000000e+03 1.759259e-01 -4.629630e-02 1.388889e-01
+ 2.000000e+03 1.000000e+03 -3.000000e+03 3.518519e-01 -4.629630e-02 1.388889e-01
+ 3.000000e+03 1.000000e+03 -3.000000e+03 5.277778e-01 -4.629630e-02 1.388889e-01
+ 4.000000e+03 1.000000e+03 -3.000000e+03 7.037037e-01 -4.629630e-02 1.388889e-01
+ -4.000000e+03 2.000000e+03 -3.000000e+03 -7.037037e-01 -9.259259e-02 1.388889e-01
+ -3.000000e+03 2.000000e+03 -3.000000e+03 -5.277778e-01 -9.259259e-02 1.388889e-01
+ -2.000000e+03 2.000000e+03 -3.000000e+03 -3.518519e-01 -9.259259e-02 1.388889e-01
+ -1.000000e+03 2.000000e+03 -3.000000e+03 -1.759259e-01 -9.259259e-02 1.388889e-01
+ 0.000000e+00 2.000000e+03 -3.000000e+03 0.000000e+00 -9.259259e-02 1.388889e-01
+ 1.000000e+03 2.000000e+03 -3.000000e+03 1.759259e-01 -9.259259e-02 1.388889e-01
+ 2.000000e+03 2.000000e+03 -3.000000e+03 3.518519e-01 -9.259259e-02 1.388889e-01
+ 3.000000e+03 2.000000e+03 -3.000000e+03 5.277778e-01 -9.259259e-02 1.388889e-01
+ 4.000000e+03 2.000000e+03 -3.000000e+03 7.037037e-01 -9.259259e-02 1.388889e-01
+ -4.000000e+03 3.000000e+03 -3.000000e+03 -7.037037e-01 -1.388889e-01 1.388889e-01
+ -3.000000e+03 3.000000e+03 -3.000000e+03 -5.277778e-01 -1.388889e-01 1.388889e-01
+ -2.000000e+03 3.000000e+03 -3.000000e+03 -3.518519e-01 -1.388889e-01 1.388889e-01
+ -1.000000e+03 3.000000e+03 -3.000000e+03 -1.759259e-01 -1.388889e-01 1.388889e-01
+ 0.000000e+00 3.000000e+03 -3.000000e+03 0.000000e+00 -1.388889e-01 1.388889e-01
+ 1.000000e+03 3.000000e+03 -3.000000e+03 1.759259e-01 -1.388889e-01 1.388889e-01
+ 2.000000e+03 3.000000e+03 -3.000000e+03 3.518519e-01 -1.388889e-01 1.388889e-01
+ 3.000000e+03 3.000000e+03 -3.000000e+03 5.277778e-01 -1.388889e-01 1.388889e-01
+ 4.000000e+03 3.000000e+03 -3.000000e+03 7.037037e-01 -1.388889e-01 1.388889e-01
+ -4.000000e+03 4.000000e+03 -3.000000e+03 -7.037037e-01 -1.851852e-01 1.388889e-01
+ -3.000000e+03 4.000000e+03 -3.000000e+03 -5.277778e-01 -1.851852e-01 1.388889e-01
+ -2.000000e+03 4.000000e+03 -3.000000e+03 -3.518519e-01 -1.851852e-01 1.388889e-01
+ -1.000000e+03 4.000000e+03 -3.000000e+03 -1.759259e-01 -1.851852e-01 1.388889e-01
+ 0.000000e+00 4.000000e+03 -3.000000e+03 0.000000e+00 -1.851852e-01 1.388889e-01
+ 1.000000e+03 4.000000e+03 -3.000000e+03 1.759259e-01 -1.851852e-01 1.388889e-01
+ 2.000000e+03 4.000000e+03 -3.000000e+03 3.518519e-01 -1.851852e-01 1.388889e-01
+ 3.000000e+03 4.000000e+03 -3.000000e+03 5.277778e-01 -1.851852e-01 1.388889e-01
+ 4.000000e+03 4.000000e+03 -3.000000e+03 7.037037e-01 -1.851852e-01 1.388889e-01
+ -4.000000e+03 -4.000000e+03 -2.000000e+03 -7.037037e-01 1.851852e-01 9.259259e-02
+ -3.000000e+03 -4.000000e+03 -2.000000e+03 -5.277778e-01 1.851852e-01 9.259259e-02
+ -2.000000e+03 -4.000000e+03 -2.000000e+03 -3.518519e-01 1.851852e-01 9.259259e-02
+ -1.000000e+03 -4.000000e+03 -2.000000e+03 -1.759259e-01 1.851852e-01 9.259259e-02
+ 0.000000e+00 -4.000000e+03 -2.000000e+03 0.000000e+00 1.851852e-01 9.259259e-02
+ 1.000000e+03 -4.000000e+03 -2.000000e+03 1.759259e-01 1.851852e-01 9.259259e-02
+ 2.000000e+03 -4.000000e+03 -2.000000e+03 3.518519e-01 1.851852e-01 9.259259e-02
+ 3.000000e+03 -4.000000e+03 -2.000000e+03 5.277778e-01 1.851852e-01 9.259259e-02
+ 4.000000e+03 -4.000000e+03 -2.000000e+03 7.037037e-01 1.851852e-01 9.259259e-02
+ -4.000000e+03 -3.000000e+03 -2.000000e+03 -7.037037e-01 1.388889e-01 9.259259e-02
+ -3.000000e+03 -3.000000e+03 -2.000000e+03 -5.277778e-01 1.388889e-01 9.259259e-02
+ -2.000000e+03 -3.000000e+03 -2.000000e+03 -3.518519e-01 1.388889e-01 9.259259e-02
+ -1.000000e+03 -3.000000e+03 -2.000000e+03 -1.759259e-01 1.388889e-01 9.259259e-02
+ 0.000000e+00 -3.000000e+03 -2.000000e+03 0.000000e+00 1.388889e-01 9.259259e-02
+ 1.000000e+03 -3.000000e+03 -2.000000e+03 1.759259e-01 1.388889e-01 9.259259e-02
+ 2.000000e+03 -3.000000e+03 -2.000000e+03 3.518519e-01 1.388889e-01 9.259259e-02
+ 3.000000e+03 -3.000000e+03 -2.000000e+03 5.277778e-01 1.388889e-01 9.259259e-02
+ 4.000000e+03 -3.000000e+03 -2.000000e+03 7.037037e-01 1.388889e-01 9.259259e-02
+ -4.000000e+03 -2.000000e+03 -2.000000e+03 -7.037037e-01 9.259259e-02 9.259259e-02
+ -3.000000e+03 -2.000000e+03 -2.000000e+03 -5.277778e-01 9.259259e-02 9.259259e-02
+ -2.000000e+03 -2.000000e+03 -2.000000e+03 -3.518519e-01 9.259259e-02 9.259259e-02
+ -1.000000e+03 -2.000000e+03 -2.000000e+03 -1.759259e-01 9.259259e-02 9.259259e-02
+ 0.000000e+00 -2.000000e+03 -2.000000e+03 0.000000e+00 9.259259e-02 9.259259e-02
+ 1.000000e+03 -2.000000e+03 -2.000000e+03 1.759259e-01 9.259259e-02 9.259259e-02
+ 2.000000e+03 -2.000000e+03 -2.000000e+03 3.518519e-01 9.259259e-02 9.259259e-02
+ 3.000000e+03 -2.000000e+03 -2.000000e+03 5.277778e-01 9.259259e-02 9.259259e-02
+ 4.000000e+03 -2.000000e+03 -2.000000e+03 7.037037e-01 9.259259e-02 9.259259e-02
+ -4.000000e+03 -1.000000e+03 -2.000000e+03 -7.037037e-01 4.629630e-02 9.259259e-02
+ -3.000000e+03 -1.000000e+03 -2.000000e+03 -5.277778e-01 4.629630e-02 9.259259e-02
+ -2.000000e+03 -1.000000e+03 -2.000000e+03 -3.518519e-01 4.629630e-02 9.259259e-02
+ -1.000000e+03 -1.000000e+03 -2.000000e+03 -1.759259e-01 4.629630e-02 9.259259e-02
+ 0.000000e+00 -1.000000e+03 -2.000000e+03 0.000000e+00 4.629630e-02 9.259259e-02
+ 1.000000e+03 -1.000000e+03 -2.000000e+03 1.759259e-01 4.629630e-02 9.259259e-02
+ 2.000000e+03 -1.000000e+03 -2.000000e+03 3.518519e-01 4.629630e-02 9.259259e-02
+ 3.000000e+03 -1.000000e+03 -2.000000e+03 5.277778e-01 4.629630e-02 9.259259e-02
+ 4.000000e+03 -1.000000e+03 -2.000000e+03 7.037037e-01 4.629630e-02 9.259259e-02
+ -4.000000e+03 0.000000e+00 -2.000000e+03 -7.037037e-01 -0.000000e+00 9.259259e-02
+ -3.000000e+03 0.000000e+00 -2.000000e+03 -5.277778e-01 -0.000000e+00 9.259259e-02
+ -2.000000e+03 0.000000e+00 -2.000000e+03 -3.518519e-01 -0.000000e+00 9.259259e-02
+ -1.000000e+03 0.000000e+00 -2.000000e+03 -1.759259e-01 -0.000000e+00 9.259259e-02
+ 0.000000e+00 0.000000e+00 -2.000000e+03 0.000000e+00 0.000000e+00 9.259259e-02
+ 1.000000e+03 0.000000e+00 -2.000000e+03 1.759259e-01 0.000000e+00 9.259259e-02
+ 2.000000e+03 0.000000e+00 -2.000000e+03 3.518519e-01 0.000000e+00 9.259259e-02
+ 3.000000e+03 0.000000e+00 -2.000000e+03 5.277778e-01 0.000000e+00 9.259259e-02
+ 4.000000e+03 0.000000e+00 -2.000000e+03 7.037037e-01 0.000000e+00 9.259259e-02
+ -4.000000e+03 1.000000e+03 -2.000000e+03 -7.037037e-01 -4.629630e-02 9.259259e-02
+ -3.000000e+03 1.000000e+03 -2.000000e+03 -5.277778e-01 -4.629630e-02 9.259259e-02
+ -2.000000e+03 1.000000e+03 -2.000000e+03 -3.518519e-01 -4.629630e-02 9.259259e-02
+ -1.000000e+03 1.000000e+03 -2.000000e+03 -1.759259e-01 -4.629630e-02 9.259259e-02
+ 0.000000e+00 1.000000e+03 -2.000000e+03 0.000000e+00 -4.629630e-02 9.259259e-02
+ 1.000000e+03 1.000000e+03 -2.000000e+03 1.759259e-01 -4.629630e-02 9.259259e-02
+ 2.000000e+03 1.000000e+03 -2.000000e+03 3.518519e-01 -4.629630e-02 9.259259e-02
+ 3.000000e+03 1.000000e+03 -2.000000e+03 5.277778e-01 -4.629630e-02 9.259259e-02
+ 4.000000e+03 1.000000e+03 -2.000000e+03 7.037037e-01 -4.629630e-02 9.259259e-02
+ -4.000000e+03 2.000000e+03 -2.000000e+03 -7.037037e-01 -9.259259e-02 9.259259e-02
+ -3.000000e+03 2.000000e+03 -2.000000e+03 -5.277778e-01 -9.259259e-02 9.259259e-02
+ -2.000000e+03 2.000000e+03 -2.000000e+03 -3.518519e-01 -9.259259e-02 9.259259e-02
+ -1.000000e+03 2.000000e+03 -2.000000e+03 -1.759259e-01 -9.259259e-02 9.259259e-02
+ 0.000000e+00 2.000000e+03 -2.000000e+03 0.000000e+00 -9.259259e-02 9.259259e-02
+ 1.000000e+03 2.000000e+03 -2.000000e+03 1.759259e-01 -9.259259e-02 9.259259e-02
+ 2.000000e+03 2.000000e+03 -2.000000e+03 3.518519e-01 -9.259259e-02 9.259259e-02
+ 3.000000e+03 2.000000e+03 -2.000000e+03 5.277778e-01 -9.259259e-02 9.259259e-02
+ 4.000000e+03 2.000000e+03 -2.000000e+03 7.037037e-01 -9.259259e-02 9.259259e-02
+ -4.000000e+03 3.000000e+03 -2.000000e+03 -7.037037e-01 -1.388889e-01 9.259259e-02
+ -3.000000e+03 3.000000e+03 -2.000000e+03 -5.277778e-01 -1.388889e-01 9.259259e-02
+ -2.000000e+03 3.000000e+03 -2.000000e+03 -3.518519e-01 -1.388889e-01 9.259259e-02
+ -1.000000e+03 3.000000e+03 -2.000000e+03 -1.759259e-01 -1.388889e-01 9.259259e-02
+ 0.000000e+00 3.000000e+03 -2.000000e+03 0.000000e+00 -1.388889e-01 9.259259e-02
+ 1.000000e+03 3.000000e+03 -2.000000e+03 1.759259e-01 -1.388889e-01 9.259259e-02
+ 2.000000e+03 3.000000e+03 -2.000000e+03 3.518519e-01 -1.388889e-01 9.259259e-02
+ 3.000000e+03 3.000000e+03 -2.000000e+03 5.277778e-01 -1.388889e-01 9.259259e-02
+ 4.000000e+03 3.000000e+03 -2.000000e+03 7.037037e-01 -1.388889e-01 9.259259e-02
+ -4.000000e+03 4.000000e+03 -2.000000e+03 -7.037037e-01 -1.851852e-01 9.259259e-02
+ -3.000000e+03 4.000000e+03 -2.000000e+03 -5.277778e-01 -1.851852e-01 9.259259e-02
+ -2.000000e+03 4.000000e+03 -2.000000e+03 -3.518519e-01 -1.851852e-01 9.259259e-02
+ -1.000000e+03 4.000000e+03 -2.000000e+03 -1.759259e-01 -1.851852e-01 9.259259e-02
+ 0.000000e+00 4.000000e+03 -2.000000e+03 0.000000e+00 -1.851852e-01 9.259259e-02
+ 1.000000e+03 4.000000e+03 -2.000000e+03 1.759259e-01 -1.851852e-01 9.259259e-02
+ 2.000000e+03 4.000000e+03 -2.000000e+03 3.518519e-01 -1.851852e-01 9.259259e-02
+ 3.000000e+03 4.000000e+03 -2.000000e+03 5.277778e-01 -1.851852e-01 9.259259e-02
+ 4.000000e+03 4.000000e+03 -2.000000e+03 7.037037e-01 -1.851852e-01 9.259259e-02
+ -4.000000e+03 -4.000000e+03 -1.000000e+03 -7.037037e-01 1.851852e-01 4.629630e-02
+ -3.000000e+03 -4.000000e+03 -1.000000e+03 -5.277778e-01 1.851852e-01 4.629630e-02
+ -2.000000e+03 -4.000000e+03 -1.000000e+03 -3.518519e-01 1.851852e-01 4.629630e-02
+ -1.000000e+03 -4.000000e+03 -1.000000e+03 -1.759259e-01 1.851852e-01 4.629630e-02
+ 0.000000e+00 -4.000000e+03 -1.000000e+03 0.000000e+00 1.851852e-01 4.629630e-02
+ 1.000000e+03 -4.000000e+03 -1.000000e+03 1.759259e-01 1.851852e-01 4.629630e-02
+ 2.000000e+03 -4.000000e+03 -1.000000e+03 3.518519e-01 1.851852e-01 4.629630e-02
+ 3.000000e+03 -4.000000e+03 -1.000000e+03 5.277778e-01 1.851852e-01 4.629630e-02
+ 4.000000e+03 -4.000000e+03 -1.000000e+03 7.037037e-01 1.851852e-01 4.629630e-02
+ -4.000000e+03 -3.000000e+03 -1.000000e+03 -7.037037e-01 1.388889e-01 4.629630e-02
+ -3.000000e+03 -3.000000e+03 -1.000000e+03 -5.277778e-01 1.388889e-01 4.629630e-02
+ -2.000000e+03 -3.000000e+03 -1.000000e+03 -3.518519e-01 1.388889e-01 4.629630e-02
+ -1.000000e+03 -3.000000e+03 -1.000000e+03 -1.759259e-01 1.388889e-01 4.629630e-02
+ 0.000000e+00 -3.000000e+03 -1.000000e+03 0.000000e+00 1.388889e-01 4.629630e-02
+ 1.000000e+03 -3.000000e+03 -1.000000e+03 1.759259e-01 1.388889e-01 4.629630e-02
+ 2.000000e+03 -3.000000e+03 -1.000000e+03 3.518519e-01 1.388889e-01 4.629630e-02
+ 3.000000e+03 -3.000000e+03 -1.000000e+03 5.277778e-01 1.388889e-01 4.629630e-02
+ 4.000000e+03 -3.000000e+03 -1.000000e+03 7.037037e-01 1.388889e-01 4.629630e-02
+ -4.000000e+03 -2.000000e+03 -1.000000e+03 -7.037037e-01 9.259259e-02 4.629630e-02
+ -3.000000e+03 -2.000000e+03 -1.000000e+03 -5.277778e-01 9.259259e-02 4.629630e-02
+ -2.000000e+03 -2.000000e+03 -1.000000e+03 -3.518519e-01 9.259259e-02 4.629630e-02
+ -1.000000e+03 -2.000000e+03 -1.000000e+03 -1.759259e-01 9.259259e-02 4.629630e-02
+ 0.000000e+00 -2.000000e+03 -1.000000e+03 0.000000e+00 9.259259e-02 4.629630e-02
+ 1.000000e+03 -2.000000e+03 -1.000000e+03 1.759259e-01 9.259259e-02 4.629630e-02
+ 2.000000e+03 -2.000000e+03 -1.000000e+03 3.518519e-01 9.259259e-02 4.629630e-02
+ 3.000000e+03 -2.000000e+03 -1.000000e+03 5.277778e-01 9.259259e-02 4.629630e-02
+ 4.000000e+03 -2.000000e+03 -1.000000e+03 7.037037e-01 9.259259e-02 4.629630e-02
+ -4.000000e+03 -1.000000e+03 -1.000000e+03 -7.037037e-01 4.629630e-02 4.629630e-02
+ -3.000000e+03 -1.000000e+03 -1.000000e+03 -5.277778e-01 4.629630e-02 4.629630e-02
+ -2.000000e+03 -1.000000e+03 -1.000000e+03 -3.518519e-01 4.629630e-02 4.629630e-02
+ -1.000000e+03 -1.000000e+03 -1.000000e+03 -1.759259e-01 4.629630e-02 4.629630e-02
+ 0.000000e+00 -1.000000e+03 -1.000000e+03 0.000000e+00 4.629630e-02 4.629630e-02
+ 1.000000e+03 -1.000000e+03 -1.000000e+03 1.759259e-01 4.629630e-02 4.629630e-02
+ 2.000000e+03 -1.000000e+03 -1.000000e+03 3.518519e-01 4.629630e-02 4.629630e-02
+ 3.000000e+03 -1.000000e+03 -1.000000e+03 5.277778e-01 4.629630e-02 4.629630e-02
+ 4.000000e+03 -1.000000e+03 -1.000000e+03 7.037037e-01 4.629630e-02 4.629630e-02
+ -4.000000e+03 0.000000e+00 -1.000000e+03 -7.037037e-01 -0.000000e+00 4.629630e-02
+ -3.000000e+03 0.000000e+00 -1.000000e+03 -5.277778e-01 -0.000000e+00 4.629630e-02
+ -2.000000e+03 0.000000e+00 -1.000000e+03 -3.518519e-01 -0.000000e+00 4.629630e-02
+ -1.000000e+03 0.000000e+00 -1.000000e+03 -1.759259e-01 -0.000000e+00 4.629630e-02
+ 0.000000e+00 0.000000e+00 -1.000000e+03 0.000000e+00 0.000000e+00 4.629630e-02
+ 1.000000e+03 0.000000e+00 -1.000000e+03 1.759259e-01 0.000000e+00 4.629630e-02
+ 2.000000e+03 0.000000e+00 -1.000000e+03 3.518519e-01 0.000000e+00 4.629630e-02
+ 3.000000e+03 0.000000e+00 -1.000000e+03 5.277778e-01 0.000000e+00 4.629630e-02
+ 4.000000e+03 0.000000e+00 -1.000000e+03 7.037037e-01 0.000000e+00 4.629630e-02
+ -4.000000e+03 1.000000e+03 -1.000000e+03 -7.037037e-01 -4.629630e-02 4.629630e-02
+ -3.000000e+03 1.000000e+03 -1.000000e+03 -5.277778e-01 -4.629630e-02 4.629630e-02
+ -2.000000e+03 1.000000e+03 -1.000000e+03 -3.518519e-01 -4.629630e-02 4.629630e-02
+ -1.000000e+03 1.000000e+03 -1.000000e+03 -1.759259e-01 -4.629630e-02 4.629630e-02
+ 0.000000e+00 1.000000e+03 -1.000000e+03 0.000000e+00 -4.629630e-02 4.629630e-02
+ 1.000000e+03 1.000000e+03 -1.000000e+03 1.759259e-01 -4.629630e-02 4.629630e-02
+ 2.000000e+03 1.000000e+03 -1.000000e+03 3.518519e-01 -4.629630e-02 4.629630e-02
+ 3.000000e+03 1.000000e+03 -1.000000e+03 5.277778e-01 -4.629630e-02 4.629630e-02
+ 4.000000e+03 1.000000e+03 -1.000000e+03 7.037037e-01 -4.629630e-02 4.629630e-02
+ -4.000000e+03 2.000000e+03 -1.000000e+03 -7.037037e-01 -9.259259e-02 4.629630e-02
+ -3.000000e+03 2.000000e+03 -1.000000e+03 -5.277778e-01 -9.259259e-02 4.629630e-02
+ -2.000000e+03 2.000000e+03 -1.000000e+03 -3.518519e-01 -9.259259e-02 4.629630e-02
+ -1.000000e+03 2.000000e+03 -1.000000e+03 -1.759259e-01 -9.259259e-02 4.629630e-02
+ 0.000000e+00 2.000000e+03 -1.000000e+03 0.000000e+00 -9.259259e-02 4.629630e-02
+ 1.000000e+03 2.000000e+03 -1.000000e+03 1.759259e-01 -9.259259e-02 4.629630e-02
+ 2.000000e+03 2.000000e+03 -1.000000e+03 3.518519e-01 -9.259259e-02 4.629630e-02
+ 3.000000e+03 2.000000e+03 -1.000000e+03 5.277778e-01 -9.259259e-02 4.629630e-02
+ 4.000000e+03 2.000000e+03 -1.000000e+03 7.037037e-01 -9.259259e-02 4.629630e-02
+ -4.000000e+03 3.000000e+03 -1.000000e+03 -7.037037e-01 -1.388889e-01 4.629630e-02
+ -3.000000e+03 3.000000e+03 -1.000000e+03 -5.277778e-01 -1.388889e-01 4.629630e-02
+ -2.000000e+03 3.000000e+03 -1.000000e+03 -3.518519e-01 -1.388889e-01 4.629630e-02
+ -1.000000e+03 3.000000e+03 -1.000000e+03 -1.759259e-01 -1.388889e-01 4.629630e-02
+ 0.000000e+00 3.000000e+03 -1.000000e+03 0.000000e+00 -1.388889e-01 4.629630e-02
+ 1.000000e+03 3.000000e+03 -1.000000e+03 1.759259e-01 -1.388889e-01 4.629630e-02
+ 2.000000e+03 3.000000e+03 -1.000000e+03 3.518519e-01 -1.388889e-01 4.629630e-02
+ 3.000000e+03 3.000000e+03 -1.000000e+03 5.277778e-01 -1.388889e-01 4.629630e-02
+ 4.000000e+03 3.000000e+03 -1.000000e+03 7.037037e-01 -1.388889e-01 4.629630e-02
+ -4.000000e+03 4.000000e+03 -1.000000e+03 -7.037037e-01 -1.851852e-01 4.629630e-02
+ -3.000000e+03 4.000000e+03 -1.000000e+03 -5.277778e-01 -1.851852e-01 4.629630e-02
+ -2.000000e+03 4.000000e+03 -1.000000e+03 -3.518519e-01 -1.851852e-01 4.629630e-02
+ -1.000000e+03 4.000000e+03 -1.000000e+03 -1.759259e-01 -1.851852e-01 4.629630e-02
+ 0.000000e+00 4.000000e+03 -1.000000e+03 0.000000e+00 -1.851852e-01 4.629630e-02
+ 1.000000e+03 4.000000e+03 -1.000000e+03 1.759259e-01 -1.851852e-01 4.629630e-02
+ 2.000000e+03 4.000000e+03 -1.000000e+03 3.518519e-01 -1.851852e-01 4.629630e-02
+ 3.000000e+03 4.000000e+03 -1.000000e+03 5.277778e-01 -1.851852e-01 4.629630e-02
+ 4.000000e+03 4.000000e+03 -1.000000e+03 7.037037e-01 -1.851852e-01 4.629630e-02
+ -4.000000e+03 -4.000000e+03 0.000000e+00 -7.037037e-01 1.851852e-01 -0.000000e+00
+ -3.000000e+03 -4.000000e+03 0.000000e+00 -5.277778e-01 1.851852e-01 -0.000000e+00
+ -2.000000e+03 -4.000000e+03 0.000000e+00 -3.518519e-01 1.851852e-01 -0.000000e+00
+ -1.000000e+03 -4.000000e+03 0.000000e+00 -1.759259e-01 1.851852e-01 -0.000000e+00
+ 0.000000e+00 -4.000000e+03 0.000000e+00 0.000000e+00 1.851852e-01 0.000000e+00
+ 1.000000e+03 -4.000000e+03 0.000000e+00 1.759259e-01 1.851852e-01 0.000000e+00
+ 2.000000e+03 -4.000000e+03 0.000000e+00 3.518519e-01 1.851852e-01 0.000000e+00
+ 3.000000e+03 -4.000000e+03 0.000000e+00 5.277778e-01 1.851852e-01 0.000000e+00
+ 4.000000e+03 -4.000000e+03 0.000000e+00 7.037037e-01 1.851852e-01 0.000000e+00
+ -4.000000e+03 -3.000000e+03 0.000000e+00 -7.037037e-01 1.388889e-01 -0.000000e+00
+ -3.000000e+03 -3.000000e+03 0.000000e+00 -5.277778e-01 1.388889e-01 -0.000000e+00
+ -2.000000e+03 -3.000000e+03 0.000000e+00 -3.518519e-01 1.388889e-01 -0.000000e+00
+ -1.000000e+03 -3.000000e+03 0.000000e+00 -1.759259e-01 1.388889e-01 -0.000000e+00
+ 0.000000e+00 -3.000000e+03 0.000000e+00 0.000000e+00 1.388889e-01 0.000000e+00
+ 1.000000e+03 -3.000000e+03 0.000000e+00 1.759259e-01 1.388889e-01 0.000000e+00
+ 2.000000e+03 -3.000000e+03 0.000000e+00 3.518519e-01 1.388889e-01 0.000000e+00
+ 3.000000e+03 -3.000000e+03 0.000000e+00 5.277778e-01 1.388889e-01 0.000000e+00
+ 4.000000e+03 -3.000000e+03 0.000000e+00 7.037037e-01 1.388889e-01 0.000000e+00
+ -4.000000e+03 -2.000000e+03 0.000000e+00 -7.037037e-01 9.259259e-02 -0.000000e+00
+ -3.000000e+03 -2.000000e+03 0.000000e+00 -5.277778e-01 9.259259e-02 -0.000000e+00
+ -2.000000e+03 -2.000000e+03 0.000000e+00 -3.518519e-01 9.259259e-02 -0.000000e+00
+ -1.000000e+03 -2.000000e+03 0.000000e+00 -1.759259e-01 9.259259e-02 -0.000000e+00
+ 0.000000e+00 -2.000000e+03 0.000000e+00 0.000000e+00 9.259259e-02 0.000000e+00
+ 1.000000e+03 -2.000000e+03 0.000000e+00 1.759259e-01 9.259259e-02 0.000000e+00
+ 2.000000e+03 -2.000000e+03 0.000000e+00 3.518519e-01 9.259259e-02 0.000000e+00
+ 3.000000e+03 -2.000000e+03 0.000000e+00 5.277778e-01 9.259259e-02 0.000000e+00
+ 4.000000e+03 -2.000000e+03 0.000000e+00 7.037037e-01 9.259259e-02 0.000000e+00
+ -4.000000e+03 -1.000000e+03 0.000000e+00 -7.037037e-01 4.629630e-02 -0.000000e+00
+ -3.000000e+03 -1.000000e+03 0.000000e+00 -5.277778e-01 4.629630e-02 -0.000000e+00
+ -2.000000e+03 -1.000000e+03 0.000000e+00 -3.518519e-01 4.629630e-02 -0.000000e+00
+ -1.000000e+03 -1.000000e+03 0.000000e+00 -1.759259e-01 4.629630e-02 -0.000000e+00
+ 0.000000e+00 -1.000000e+03 0.000000e+00 0.000000e+00 4.629630e-02 0.000000e+00
+ 1.000000e+03 -1.000000e+03 0.000000e+00 1.759259e-01 4.629630e-02 0.000000e+00
+ 2.000000e+03 -1.000000e+03 0.000000e+00 3.518519e-01 4.629630e-02 0.000000e+00
+ 3.000000e+03 -1.000000e+03 0.000000e+00 5.277778e-01 4.629630e-02 0.000000e+00
+ 4.000000e+03 -1.000000e+03 0.000000e+00 7.037037e-01 4.629630e-02 0.000000e+00
+ -4.000000e+03 0.000000e+00 0.000000e+00 -7.037037e-01 0.000000e+00 0.000000e+00
+ -3.000000e+03 0.000000e+00 0.000000e+00 -5.277778e-01 0.000000e+00 0.000000e+00
+ -2.000000e+03 0.000000e+00 0.000000e+00 -3.518519e-01 0.000000e+00 0.000000e+00
+ -1.000000e+03 0.000000e+00 0.000000e+00 -1.759259e-01 0.000000e+00 0.000000e+00
+ 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
+ 1.000000e+03 0.000000e+00 0.000000e+00 1.759259e-01 0.000000e+00 0.000000e+00
+ 2.000000e+03 0.000000e+00 0.000000e+00 3.518519e-01 0.000000e+00 0.000000e+00
+ 3.000000e+03 0.000000e+00 0.000000e+00 5.277778e-01 0.000000e+00 0.000000e+00
+ 4.000000e+03 0.000000e+00 0.000000e+00 7.037037e-01 0.000000e+00 0.000000e+00
+ -4.000000e+03 1.000000e+03 0.000000e+00 -7.037037e-01 -4.629630e-02 0.000000e+00
+ -3.000000e+03 1.000000e+03 0.000000e+00 -5.277778e-01 -4.629630e-02 0.000000e+00
+ -2.000000e+03 1.000000e+03 0.000000e+00 -3.518519e-01 -4.629630e-02 0.000000e+00
+ -1.000000e+03 1.000000e+03 0.000000e+00 -1.759259e-01 -4.629630e-02 0.000000e+00
+ 0.000000e+00 1.000000e+03 0.000000e+00 0.000000e+00 -4.629630e-02 0.000000e+00
+ 1.000000e+03 1.000000e+03 0.000000e+00 1.759259e-01 -4.629630e-02 0.000000e+00
+ 2.000000e+03 1.000000e+03 0.000000e+00 3.518519e-01 -4.629630e-02 0.000000e+00
+ 3.000000e+03 1.000000e+03 0.000000e+00 5.277778e-01 -4.629630e-02 0.000000e+00
+ 4.000000e+03 1.000000e+03 0.000000e+00 7.037037e-01 -4.629630e-02 0.000000e+00
+ -4.000000e+03 2.000000e+03 0.000000e+00 -7.037037e-01 -9.259259e-02 0.000000e+00
+ -3.000000e+03 2.000000e+03 0.000000e+00 -5.277778e-01 -9.259259e-02 0.000000e+00
+ -2.000000e+03 2.000000e+03 0.000000e+00 -3.518519e-01 -9.259259e-02 0.000000e+00
+ -1.000000e+03 2.000000e+03 0.000000e+00 -1.759259e-01 -9.259259e-02 0.000000e+00
+ 0.000000e+00 2.000000e+03 0.000000e+00 0.000000e+00 -9.259259e-02 0.000000e+00
+ 1.000000e+03 2.000000e+03 0.000000e+00 1.759259e-01 -9.259259e-02 0.000000e+00
+ 2.000000e+03 2.000000e+03 0.000000e+00 3.518519e-01 -9.259259e-02 0.000000e+00
+ 3.000000e+03 2.000000e+03 0.000000e+00 5.277778e-01 -9.259259e-02 0.000000e+00
+ 4.000000e+03 2.000000e+03 0.000000e+00 7.037037e-01 -9.259259e-02 0.000000e+00
+ -4.000000e+03 3.000000e+03 0.000000e+00 -7.037037e-01 -1.388889e-01 0.000000e+00
+ -3.000000e+03 3.000000e+03 0.000000e+00 -5.277778e-01 -1.388889e-01 0.000000e+00
+ -2.000000e+03 3.000000e+03 0.000000e+00 -3.518519e-01 -1.388889e-01 0.000000e+00
+ -1.000000e+03 3.000000e+03 0.000000e+00 -1.759259e-01 -1.388889e-01 0.000000e+00
+ 0.000000e+00 3.000000e+03 0.000000e+00 0.000000e+00 -1.388889e-01 0.000000e+00
+ 1.000000e+03 3.000000e+03 0.000000e+00 1.759259e-01 -1.388889e-01 0.000000e+00
+ 2.000000e+03 3.000000e+03 0.000000e+00 3.518519e-01 -1.388889e-01 0.000000e+00
+ 3.000000e+03 3.000000e+03 0.000000e+00 5.277778e-01 -1.388889e-01 0.000000e+00
+ 4.000000e+03 3.000000e+03 0.000000e+00 7.037037e-01 -1.388889e-01 0.000000e+00
+ -4.000000e+03 4.000000e+03 0.000000e+00 -7.037037e-01 -1.851852e-01 0.000000e+00
+ -3.000000e+03 4.000000e+03 0.000000e+00 -5.277778e-01 -1.851852e-01 0.000000e+00
+ -2.000000e+03 4.000000e+03 0.000000e+00 -3.518519e-01 -1.851852e-01 0.000000e+00
+ -1.000000e+03 4.000000e+03 0.000000e+00 -1.759259e-01 -1.851852e-01 0.000000e+00
+ 0.000000e+00 4.000000e+03 0.000000e+00 0.000000e+00 -1.851852e-01 0.000000e+00
+ 1.000000e+03 4.000000e+03 0.000000e+00 1.759259e-01 -1.851852e-01 0.000000e+00
+ 2.000000e+03 4.000000e+03 0.000000e+00 3.518519e-01 -1.851852e-01 0.000000e+00
+ 3.000000e+03 4.000000e+03 0.000000e+00 5.277778e-01 -1.851852e-01 0.000000e+00
+ 4.000000e+03 4.000000e+03 0.000000e+00 7.037037e-01 -1.851852e-01 0.000000e+00
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp.cfg 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,134 @@
+# -*- Python -*-
+[pylithapp]
+
+[pylithapp.launcher] # WARNING: THIS IS NOT PORTABLE
+command = mpirun -np ${nodes}
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+#timedependent = 1
+#implicit = 1
+#petsc = 1
+#solverlinear = 1
+#meshiocubit = 1
+#implicitelasticity = 1
+#quadrature3d = 1
+#fiatlagrange = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = mesh.exo
+use_nodeset_names = True
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+bc = [x_neg,x_pos,y_neg,z_neg]
+
+[pylithapp.timedependent.formulation]
+output = [domain,subdomain]
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic,viscoelastic]
+materials.elastic = pylith.materials.ElasticIsotropic3D
+materials.viscoelastic = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+
+[pylithapp.timedependent.materials.viscoelastic]
+label = Elastic material 2
+id = 2
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC +x eface
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC -x face
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [1]
+label = face_yneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC -y face
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Dirichlet BC -z edge
+db_initial.iohandler.filename = axial_disp.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+pc_type = asm
+
+# Change the preconditioner settings.
+sub_pc_factor_shift_type = none
+
+ksp_rtol = 1.0e-8
+ksp_max_it = 100
+ksp_gmres_restart = 50
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+#log_summary = true
+# start_in_debugger = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain.writer]
+filename = axialdisp.vtk
+
+[pylithapp.problem.formulation.output.subdomain]
+label = face_zpos
+writer.filename = axialdisp-groundsurf.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = axialdisp-elastic.vtk
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = axialdisp-viscoelastic.vtk
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_gendb.py
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_gendb.py (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_gendb.py 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,95 @@
+#!/usr/bin/env python
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+## @file tests/3d/hex8/axialdisp_gendb.py
+##
+## @brief Python script to generate spatial database with displacement
+## boundary conditions for the axial displacement test.
+
+import numpy
+
+class GenerateDB(object):
+ """
+ Python object to generate spatial database with displacement
+ boundary conditions for the axial displacement test.
+ """
+
+ def __init__(self):
+ """
+ Constructor.
+ """
+ return
+
+
+ def run(self):
+ """
+ Generate the database.
+ """
+ # Domain
+ x = numpy.arange(-4000.0, 4000.1, 1000.0)
+ y = numpy.arange(-4000.0, 4000.1, 1000.0)
+ z = numpy.arange(-6000.0, 0.1, 1000.0)
+ nxpts = x.shape[0]
+ nypts = y.shape[0]
+ nzpts = z.shape[0]
+
+ xx = numpy.tile(x, (1,nypts*nzpts))
+ yy = numpy.tile(y, (nxpts,nzpts)).transpose()
+ zz = numpy.tile(z, (nxpts*nypts,1)).transpose()
+
+ xyz = numpy.zeros( (nxpts*nypts*nzpts, 3), dtype=numpy.float64)
+ xyz[:,0] = numpy.ravel(xx)
+ xyz[:,1] = numpy.ravel(yy)
+ xyz[:,2] = numpy.ravel(zz)
+
+ from axialdisp_soln import AnalyticalSoln
+ soln = AnalyticalSoln()
+ disp = soln.displacement(xyz)
+
+ from spatialdata.geocoords.CSCart import CSCart
+ cs = CSCart()
+ cs.inventory.spaceDim = 3
+ cs._configure()
+ data = {'points': xyz,
+ 'coordsys': cs,
+ 'data_dim': 3,
+ 'values': [{'name': "displacement-x",
+ 'units': "m",
+ 'data': numpy.ravel(disp[:,0])},
+ {'name': "displacement-y",
+ 'units': "m",
+ 'data': numpy.ravel(disp[:,1])},
+ {'name': "displacement-z",
+ 'units': "m",
+ 'data': numpy.ravel(disp[:,2])}]}
+
+ from spatialdata.spatialdb.SimpleIOAscii import SimpleIOAscii
+ io = SimpleIOAscii()
+ io.inventory.filename = "axial_disp.spatialdb"
+ io._configure()
+ io.write(data)
+ return
+
+
+# ======================================================================
+if __name__ == "__main__":
+ app = GenerateDB()
+ app.run()
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_soln.py
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_soln.py (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/axialdisp_soln.py 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,107 @@
+#!/usr/bin/env python
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+## @file tests/2d/quad4/axialdisp_soln.py
+##
+## @brief Analytical solution to axial displacement problem.
+
+## 3-D axial extension test with linear hexahedral cells.
+
+import numpy
+
+# Physical properties
+p_density = 2500.0
+p_vs = 3000.0
+p_vp = 5291.502622129181
+
+p_mu = p_density*p_vs**2
+p_lambda = p_density*p_vp**2 - 2*p_mu
+
+# Uniform stress field (plane strain)
+sxx = 1.0e+7
+syy = 0.0
+szz = 0.0
+sxy = 0.0
+syz = 0.0
+sxz = 0.0
+
+# Uniform strain field
+exx = 1.0/(2*p_mu) * (sxx - p_lambda/(3*p_lambda+2*p_mu) * (sxx+syy+szz))
+eyy = 1.0/(2*p_mu) * (syy - p_lambda/(3*p_lambda+2*p_mu) * (sxx+syy+szz))
+ezz = 1.0/(2*p_mu) * (szz - p_lambda/(3*p_lambda+2*p_mu) * (sxx+syy+szz))
+
+exy = 1.0/(2*p_mu) * (sxy)
+eyz = 1.0/(2*p_mu) * (syz)
+exz = 1.0/(2*p_mu) * (sxz)
+
+#print exx,eyy,exy,ezz
+#print -exx*p_lambda/(p_lambda+2*p_mu)
+
+# ----------------------------------------------------------------------
+class AnalyticalSoln(object):
+ """
+ Analytical solution to axial/shear displacement problem.
+ """
+
+ def __init__(self):
+ return
+
+
+ def displacement(self, locs):
+ """
+ Compute displacement field at locations.
+ """
+ (npts, dim) = locs.shape
+ disp = numpy.zeros( (npts, 3), dtype=numpy.float64)
+ disp[:,0] = exx*locs[:,0] + exy*locs[:,1] + exz*locs[:,2]
+ disp[:,1] = eyy*locs[:,1] + exy*locs[:,0] + eyz*locs[:,2]
+ disp[:,2] = ezz*locs[:,2] + exz*locs[:,0] + eyz*locs[:,1]
+ return disp
+
+
+ def strain(self, locs):
+ """
+ Compute strain field at locations.
+ """
+ (npts, dim) = locs.shape
+ strain = numpy.zeros( (npts, 6), dtype=numpy.float64)
+ strain[:,0] = exx
+ strain[:,1] = eyy
+ strain[:,2] = ezz
+ strain[:,3] = exy
+ strain[:,4] = eyz
+ strain[:,5] = exz
+ return strain
+
+
+ def stress(self, locs):
+ """
+ Compute stress field at locations.
+ """
+ (npts, dim) = locs.shape
+ stress = numpy.zeros( (npts, 6), dtype=numpy.float64)
+ stress[:,0] = sxx
+ stress[:,1] = syy
+ stress[:,2] = szz
+ stress[:,3] = sxy
+ stress[:,4] = syz
+ stress[:,5] = sxz
+ return stress
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/geometry.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/geometry.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/geometry.jou 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,45 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 8 km x 8 km x 6 km
+# -4 km <= x <= 4 km
+# -4 km <= y <= 4 km
+# -6 km <= z <= 0 km
+reset
+brick x 8000 y 8000 z 6000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -3000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane xplane offset +2000
+surface 8 name "fault_surface_b"
+create planar surface with plane zplane offset -2000
+surface 9 name "material_interface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+webcut volume 1 with plane surface fault_surface_b
+webcut volume 1 with plane surface material_interface
+webcut volume 5 with plane surface material_interface
+webcut volume 6 with plane surface material_interface
+volume 1 name "elastic_xpos"
+volume 6 name "elastic_mid"
+volume 5 name "elastic_xneg"
+volume 7 name "visco_xpos"
+volume 9 name "visco_mid"
+volume 8 name "visco_xneg"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/matprops.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/matprops.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/matprops.spatialdb 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = density vs vp
+ value-units = kg/m**3 m/s m/s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2500.0 3000.0 5291.502622129181
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/mesh.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/mesh.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/mesh.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/mesh.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/mesh.jou 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,110 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 1000
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+volume all scheme tetmesh
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 5 6
+block 1 name "elastic"
+block 2 volume 7 9 8
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create nodesets for faults
+# ----------------------------------------------------------------------
+group "fault" add node in surface fault_surface
+group "fault" add node in fault_surface at A
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+group "fault_b" add node in surface fault_surface_b
+group "fault_b" add node in surface fault_surface_b at A
+nodeset 11 group fault_b
+nodeset 11 name "fault_b"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 31
+group "face_xpos" add node in surface 39
+nodeset 20 group face_xpos
+nodeset 20 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 41
+group "face_xneg" add node in surface 49
+nodeset 21 group face_xneg
+nodeset 21 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 32
+group "face_ypos" add node in surface 44
+group "face_ypos" add node in surface 54
+group "face_ypos" add node in surface 38
+group "face_ypos" add node in surface 46
+group "face_ypos" add node in surface 56
+nodeset 22 group face_ypos
+nodeset 22 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 34
+group "face_yneg" add node in surface 42
+group "face_yneg" add node in surface 52
+group "face_yneg" add node in surface 36
+group "face_yneg" add node in surface 48
+group "face_yneg" add node in surface 58
+nodeset 23 group face_yneg
+nodeset 23 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 18
+group "face_zpos" add node in surface 27
+group "face_zpos" add node in surface 23
+nodeset 24 group face_zpos
+nodeset 24 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 17
+group "face_zneg" add node in surface 28
+group "face_zneg" add node in surface 21
+nodeset 25 group face_zneg
+nodeset 25 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 26 group face_zneg_nofault
+nodeset 26 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "mesh.exo" dimension 3 overwrite
+
+
+
Added: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/testpylith.py
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/testpylith.py (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/testpylith.py 2010-09-03 00:33:01 UTC (rev 17169)
@@ -0,0 +1,51 @@
+#!/usr/bin/env python
+#
+# ======================================================================
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ======================================================================
+#
+
+__requires__ = "PyLith"
+
+import unittest
+
+def suite():
+ """
+ Create test suite.
+ """
+ suite = unittest.TestSuite()
+
+ from TestLgDeformRigidBody import TestRigidBody
+ suite.addTest(unittest.makeSuite(TestRigidBody))
+
+ from TestLgDeformTraction import TestTraction
+ suite.addTest(unittest.makeSuite(TestTraction))
+
+ return suite
+
+
+def main():
+ """
+ Run test suite.
+ """
+ unittest.TextTestRunner(verbosity=2).run(suite())
+ return
+
+
+# ----------------------------------------------------------------------
+if __name__ == '__main__':
+ main()
+
+
+# End of file
Property changes on: short/3D/PyLith/branches/v1.5-stable/tests_auto/3dnew/tet4/testpylith.py
___________________________________________________________________
Name: svn:executable
+ *
More information about the CIG-COMMITS
mailing list