From 593e26200f8d8d1c3eee3501584c686faf44880c Mon Sep 17 00:00:00 2001 From: anne-glerum Date: Tue, 21 Sep 2021 14:56:17 +0200 Subject: [PATCH 1/5] Add stress invariant pp --- .../visualization/stress_second_invariant.h | 65 +++++++++ .../visualization/stress_second_invariant.cc | 127 ++++++++++++++++++ 2 files changed, 192 insertions(+) create mode 100644 include/aspect/postprocess/visualization/stress_second_invariant.h create mode 100644 source/postprocess/visualization/stress_second_invariant.cc diff --git a/include/aspect/postprocess/visualization/stress_second_invariant.h b/include/aspect/postprocess/visualization/stress_second_invariant.h new file mode 100644 index 00000000000..acd5557bbe5 --- /dev/null +++ b/include/aspect/postprocess/visualization/stress_second_invariant.h @@ -0,0 +1,65 @@ +/* + Copyright (C) 2011 - 2020 by the authors of the ASPECT code. + + This file is part of ASPECT. + + ASPECT is free software; you can redistribute it and/or modify + it under the terms of the GNU General Public License as published by + the Free Software Foundation; either version 2, or (at your option) + any later version. + + ASPECT is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + GNU General Public License for more details. + + You should have received a copy of the GNU General Public License + along with ASPECT; see the file LICENSE. If not see + . +*/ + + +#ifndef _aspect_postprocess_visualization_stress_second_invariant_h +#define _aspect_postprocess_visualization_stress_second_invariant_h + +#include +#include + +#include + + +namespace aspect +{ + namespace Postprocess + { + namespace VisualizationPostprocessors + { + /** + * A class derived from DataPostprocessor that takes an output vector + * and computes a variable that represents the difference between the + * second moment invariant of the deviatoric stress and the yield stress. + * + * The member functions are all implementations of those declared in the + * base class. See there for their meaning. + */ + template + class StressSecondInvariant + : public DataPostprocessorScalar, + public SimulatorAccess, + public Interface + { + public: + /** + * Constructor. + */ + StressSecondInvariant (); + + void + evaluate_vector_field(const DataPostprocessorInputs::Vector &input_data, + std::vector> &computed_quantities) const override; + }; + } + } +} + +#endif diff --git a/source/postprocess/visualization/stress_second_invariant.cc b/source/postprocess/visualization/stress_second_invariant.cc new file mode 100644 index 00000000000..652c82e2a09 --- /dev/null +++ b/source/postprocess/visualization/stress_second_invariant.cc @@ -0,0 +1,127 @@ +/* + Copyright (C) 2011 - 2021 by the authors of the ASPECT code. + + This file is part of ASPECT. + + ASPECT is free software; you can redistribute it and/or modify + it under the terms of the GNU General Public License as published by + the Free Software Foundation; either version 2, or (at your option) + any later version. + + ASPECT is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + GNU General Public License for more details. + + You should have received a copy of the GNU General Public License + along with ASPECT; see the file LICENSE. If not see + . +*/ + + +#include + +namespace aspect +{ + namespace Postprocess + { + namespace VisualizationPostprocessors + { + template + StressSecondInvariant:: + StressSecondInvariant () + : + DataPostprocessorScalar ("stress_second_invariant", + update_values | update_gradients | update_quadrature_points) + {} + + + + template + void + StressSecondInvariant:: + evaluate_vector_field(const DataPostprocessorInputs::Vector &input_data, + std::vector> &computed_quantities) const + { + const unsigned int n_quadrature_points = input_data.solution_values.size(); + Assert(computed_quantities.size() == n_quadrature_points, ExcInternalError()); + Assert(computed_quantities[0].size() == 1, ExcInternalError()); + Assert(input_data.solution_values[0].size() == this->introspection().n_components, ExcInternalError()); + Assert(input_data.solution_gradients[0].size() == this->introspection().n_components, ExcInternalError()); + + MaterialModel::MaterialModelInputs in(input_data, + this->introspection()); + MaterialModel::MaterialModelOutputs out(n_quadrature_points, + this->n_compositional_fields()); + + this->get_material_model().create_additional_named_outputs(out); + + // Get the viscosity + this->get_material_model().evaluate(in, out); + + for (unsigned int q = 0; q < n_quadrature_points; ++q) + { + const SymmetricTensor<2, dim> strain_rate = in.strain_rate[q]; + const SymmetricTensor<2, dim> deviatoric_strain_rate = (this->get_material_model().is_compressible() + ? strain_rate - 1. / 3 * trace(strain_rate) * unit_symmetric_tensor() + : strain_rate); + + const double eta = out.viscosities[q]; + + // Compressive stress is positive in geoscience applications + // TODO what sign to use here? In the material model we don't add + // the minus sign for the comparison with the yield stress. + SymmetricTensor<2, dim> stress = -2. * eta * deviatoric_strain_rate + + in.pressure[q] * unit_symmetric_tensor(); + + // Add elastic stresses if existent + if (this->get_parameters().enable_elasticity == true) + { + stress[0][0] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_xx")]; + stress[1][1] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_yy")]; + stress[0][1] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_xy")]; + + if (dim == 3) + { + stress[2][2] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_zz")]; + stress[0][2] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_xz")]; + stress[1][2] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_yz")]; + } + } + + // Compute the deviatoric stress after adding the elastic stresses that might have added + // diagonal components + const SymmetricTensor<2, dim> deviatoric_stress = deviator(stress); + + // Compute the second moment invariant of the deviatoric stress + const double stress_invariant = std::sqrt(std::fabs(second_invariant(deviatoric_stress))); + + // Output the second stress invariant + computed_quantities[q](0) = stress_invariant; + } + + // average the values if requested + const auto &viz = this->get_postprocess_manager().template get_matching_postprocessor>(); + if (!viz.output_pointwise_stress_and_strain()) + average_quantities(computed_quantities); + } + } + } +} + + +// explicit instantiations +namespace aspect +{ + namespace Postprocess + { + namespace VisualizationPostprocessors + { + ASPECT_REGISTER_VISUALIZATION_POSTPROCESSOR(StressSecondInvariant, + "stress second invariant", + "A visualization output object that generates output " + "the second moment invariant " + "of the deviatoric stress tensor.") + } + } +} From f4bce3b0dd5ccea7e8c5da1bdd5bfead9a367565 Mon Sep 17 00:00:00 2001 From: anne-glerum Date: Thu, 30 Sep 2021 11:31:08 +0200 Subject: [PATCH 2/5] Clean up --- .../visualization/stress_second_invariant.h | 8 +++-- .../visualization/stress_second_invariant.cc | 35 +++++++++++-------- 2 files changed, 25 insertions(+), 18 deletions(-) diff --git a/include/aspect/postprocess/visualization/stress_second_invariant.h b/include/aspect/postprocess/visualization/stress_second_invariant.h index acd5557bbe5..458e7f4b1e2 100644 --- a/include/aspect/postprocess/visualization/stress_second_invariant.h +++ b/include/aspect/postprocess/visualization/stress_second_invariant.h @@ -1,5 +1,5 @@ /* - Copyright (C) 2011 - 2020 by the authors of the ASPECT code. + Copyright (C) 2011 - 2021 by the authors of the ASPECT code. This file is part of ASPECT. @@ -19,6 +19,7 @@ */ + #ifndef _aspect_postprocess_visualization_stress_second_invariant_h #define _aspect_postprocess_visualization_stress_second_invariant_h @@ -28,6 +29,7 @@ #include + namespace aspect { namespace Postprocess @@ -36,8 +38,8 @@ namespace aspect { /** * A class derived from DataPostprocessor that takes an output vector - * and computes a variable that represents the difference between the - * second moment invariant of the deviatoric stress and the yield stress. + * and computes a variable that represents the second moment invariant + * of the deviatoric stress. * * The member functions are all implementations of those declared in the * base class. See there for their meaning. diff --git a/source/postprocess/visualization/stress_second_invariant.cc b/source/postprocess/visualization/stress_second_invariant.cc index 652c82e2a09..f0c1f6ec265 100644 --- a/source/postprocess/visualization/stress_second_invariant.cc +++ b/source/postprocess/visualization/stress_second_invariant.cc @@ -19,8 +19,11 @@ */ + #include + + namespace aspect { namespace Postprocess @@ -49,14 +52,17 @@ namespace aspect Assert(input_data.solution_values[0].size() == this->introspection().n_components, ExcInternalError()); Assert(input_data.solution_gradients[0].size() == this->introspection().n_components, ExcInternalError()); + // Create the material model inputs and outputs to + // retrieve the current viscosity. MaterialModel::MaterialModelInputs in(input_data, - this->introspection()); + this->introspection(), + /*compute_strain_rate = */ true); + + in.requested_properties = MaterialModel::MaterialProperties::viscosity; + MaterialModel::MaterialModelOutputs out(n_quadrature_points, this->n_compositional_fields()); - this->get_material_model().create_additional_named_outputs(out); - - // Get the viscosity this->get_material_model().evaluate(in, out); for (unsigned int q = 0; q < n_quadrature_points; ++q) @@ -68,13 +74,11 @@ namespace aspect const double eta = out.viscosities[q]; - // Compressive stress is positive in geoscience applications - // TODO what sign to use here? In the material model we don't add - // the minus sign for the comparison with the yield stress. + // Compressive stress is positive in geoscience applications. SymmetricTensor<2, dim> stress = -2. * eta * deviatoric_strain_rate + in.pressure[q] * unit_symmetric_tensor(); - // Add elastic stresses if existent + // Add elastic stresses if existent. if (this->get_parameters().enable_elasticity == true) { stress[0][0] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_xx")]; @@ -89,15 +93,16 @@ namespace aspect } } - // Compute the deviatoric stress after adding the elastic stresses that might have added - // diagonal components + // Compute the deviatoric stress tensor after elastic stresses were added. const SymmetricTensor<2, dim> deviatoric_stress = deviator(stress); // Compute the second moment invariant of the deviatoric stress + // in the same way as the second moment invariant of the deviatoric + // strain rate is computed in the viscoplastic material model. + // TODO check that this is valid for the compressible case. const double stress_invariant = std::sqrt(std::fabs(second_invariant(deviatoric_stress))); - // Output the second stress invariant - computed_quantities[q](0) = stress_invariant; + computed_quantities[q](0) = stress_invariant; } // average the values if requested @@ -110,6 +115,7 @@ namespace aspect } + // explicit instantiations namespace aspect { @@ -119,9 +125,8 @@ namespace aspect { ASPECT_REGISTER_VISUALIZATION_POSTPROCESSOR(StressSecondInvariant, "stress second invariant", - "A visualization output object that generates output " - "the second moment invariant " - "of the deviatoric stress tensor.") + "A visualization output object that outputs " + "the second moment invariant of the deviatoric stress tensor.") } } } From 2c5636164e523d56edf29f238cd1aeead7418fca Mon Sep 17 00:00:00 2001 From: anne-glerum Date: Thu, 30 Sep 2021 11:51:51 +0200 Subject: [PATCH 3/5] Add pp to existing test --- tests/visualization_stress.prm | 4 +- .../solution/solution-00000.0000.gnuplot | 518 +++++++++--------- .../solution/solution-00000.0000.gnuplot | 518 +++++++++--------- 3 files changed, 521 insertions(+), 519 deletions(-) diff --git a/tests/visualization_stress.prm b/tests/visualization_stress.prm index e8c11bdd794..2fb1e669acf 100644 --- a/tests/visualization_stress.prm +++ b/tests/visualization_stress.prm @@ -6,6 +6,8 @@ # the stress is then simply this tensor times -2*eta plus the pressure and we set # eta=42 below. compared to the visualization_shear_stress test, the test # also includes a nonzero gravity so that we get a nonconstant pressure +# +# the second invariant of the deviatoric stress will be 42. set Dimension = 2 set CFL number = 1.0 @@ -104,6 +106,6 @@ subsection Postprocess set Interpolate output = false set Point-wise stress and strain = true set Output format = gnuplot - set List of output variables = stress + set List of output variables = stress, stress second invariant end end diff --git a/tests/visualization_stress/solution/solution-00000.0000.gnuplot b/tests/visualization_stress/solution/solution-00000.0000.gnuplot index 2ae6fd0a7ab..8c1989888d4 100644 --- a/tests/visualization_stress/solution/solution-00000.0000.gnuplot +++ b/tests/visualization_stress/solution/solution-00000.0000.gnuplot @@ -1,455 +1,455 @@ # This file was generated by the deal.II library. - - +# Date = 2021/9/30 +# Time = 11:46:01 # # For a description of the GNUPLOT format see the GNUPLOT manual. # -#

-0 0 0 0 9.99962 1 9.99962 -42 -42 9.99962 -0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 +#

+0 0 0 0 9.99962 1 9.99962 -42 -42 9.99962 42 +0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 41.9999 -0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42.0001 -42.0001 8.75001 +0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 41.9999 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42.0001 -42.0001 8.75001 42.0001 -0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 -0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 +0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 41.9999 +0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42.0001 -42.0001 8.75001 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75011 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42.0001 -42.0001 8.75001 42.0001 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75011 42 -0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42 -42 8.75005 +0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 41.9999 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42 -42 8.75005 42 -0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -42 -42 7.50001 +0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 42 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -42 -42 7.50001 42 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42 -42 8.75005 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75009 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42 -42 8.75005 42 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75009 42 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50001 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50008 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50001 42 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50008 42 -0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 -0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 +0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 +0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75011 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.7501 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75011 42 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.7501 42 -0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 -0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 +0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 +0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.7501 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.7501 42 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75009 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.75009 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75009 42 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.75009 42 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50008 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50008 42 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 42 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.75009 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.75009 42 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 42 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 42 -0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -41.9999 -41.9999 7.50004 +0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 42 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -41.9999 -41.9999 7.50004 41.9999 -0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -42 -42 6.25003 +0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 42 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -42 -42 6.25003 42 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50004 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50007 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50004 42 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50007 42 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42.0001 -42.0001 6.25003 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25007 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42.0001 -42.0001 6.25003 42.0001 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25007 42 -0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -41.9999 -41.9999 6.25002 +0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 42 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -41.9999 -41.9999 6.25002 41.9999 -0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -42 -42 5.00008 +0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 42 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -42 -42 5.00008 42 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42 -42 6.25002 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25009 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42 -42 6.25002 42 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25009 42 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -42 -42 5.00008 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00004 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -42 -42 5.00008 42 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00004 42 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50007 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50007 42 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 42 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25007 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25004 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25007 42 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25004 42 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 42 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 42 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25004 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25004 42 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25009 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25003 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25009 42 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25003 42 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00004 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00003 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00004 42 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00003 42 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25003 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25003 42 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00003 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00004 -42 -42 5.00005 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00003 42 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00004 -42 -42 5.00005 42 -0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 -0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 +0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 +0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75007 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75007 42 -0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 -0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 +0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 +0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75007 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75007 42 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 42 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75008 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75008 42 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 42 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 42 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75008 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75008 42 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 42 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50007 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 42 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50007 42 -0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 -0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 +0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 +0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 41.9999 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42.0001 -42.0001 8.75016 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 42 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42.0001 -42.0001 8.75016 42.0001 -0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 -1 0 0 0 10.0006 0 10.0006 -42 -42 10.0006 +0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 41.9999 +1 0 0 0 10.0006 0 10.0006 -42 -42 10.0006 42 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42.0001 -42.0001 8.75016 -1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42.0001 -42.0001 8.75016 42.0001 +1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 41.9999 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42 -42 8.75013 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 42 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42 -42 8.75013 42 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50007 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50014 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50007 42 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50014 42 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42 -42 8.75013 -1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42 -42 8.75013 42 +1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 41.9999 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -42 -42 7.50014 -1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -42 -42 7.50014 42 +1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 42 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 42 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 42 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25008 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25008 42 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50008 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 42 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50008 42 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25008 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25005 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25008 42 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25005 42 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25009 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25009 42 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00006 -42 -42 5.00005 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00007 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00006 -42 -42 5.00005 42 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00007 42 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25009 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25004 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25009 42 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25004 42 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00007 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00006 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00007 42 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00006 42 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50008 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50011 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50008 42 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50011 42 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25005 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42.0001 -42.0001 6.25009 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25005 42 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42.0001 -42.0001 6.25009 42.0001 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -41.9999 -41.9999 7.50011 -1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -41.9999 -41.9999 7.50011 41.9999 +1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 42 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -42 -42 6.25009 -1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -42 -42 6.25009 42 +1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 42 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25004 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42 -42 6.25011 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25004 42 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42 -42 6.25011 42 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00006 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -42 -42 5.00002 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00006 42 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -42 -42 5.00002 42 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -41.9999 -41.9999 6.25011 -1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -41.9999 -41.9999 6.25011 41.9999 +1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 42 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -42 -42 5.00002 -1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -42 -42 5.00002 42 +1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 42 -0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -41.9999 -41.9999 5.00004 +0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 42 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -41.9999 -41.9999 5.00004 41.9999 -0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -42 -42 3.7501 +0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 42 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -42 -42 3.7501 42 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -41.9999 -41.9999 5.00004 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00008 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -41.9999 -41.9999 5.00004 41.9999 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00008 42 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -42.0001 -42.0001 3.7501 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -42.0001 -42.0001 3.7501 42.0001 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 42 -0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -41.9999 -41.9999 3.75004 +0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 42 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -41.9999 -41.9999 3.75004 41.9999 -0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42 -42 2.50005 +0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 42 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42 -42 2.50005 42 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -41.9999 -41.9999 3.75004 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -41.9999 -41.9999 3.75004 41.9999 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 42 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42 -42 2.50005 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -42 -42 2.50008 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42 -42 2.50005 42 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -42 -42 2.50008 42 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00008 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00005 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00008 42 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00005 42 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75006 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 42 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75006 42 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00005 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00004 -42 -42 5.00005 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00005 42 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00004 -42 -42 5.00005 42 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75006 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75006 -42 -42 3.75004 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75006 42 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75006 -42 -42 3.75004 42 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75011 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 42 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75011 42 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -42 -42 2.50008 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.50003 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -42 -42 2.50008 42 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.50003 42 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75011 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75006 -42 -42 3.75004 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75011 42 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75006 -42 -42 3.75004 42 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.50003 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.50003 42 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 41.9999 -0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42.0001 -42.0001 2.49999 +0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 42 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42.0001 -42.0001 2.49999 42.0001 -0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -41.9999 -41.9999 1.25 +0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 42 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -41.9999 -41.9999 1.25 41.9999 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42.0001 -42.0001 2.49999 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -41.9999 -41.9999 2.50004 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42.0001 -42.0001 2.49999 42.0001 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -41.9999 -41.9999 2.50004 41.9999 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -41.9999 -41.9999 1.25 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.25005 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -41.9999 -41.9999 1.25 41.9999 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.25005 42 -0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -42 -42 1.24998 +0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 42 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -42 -42 1.24998 42 -0 1 1 0 -0.000421673 1 -0.000421673 -42 -42 -0.000421673 -0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 +0 1 1 0 -0.000421673 1 -0.000421673 -42 -42 -0.000421673 42 +0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 41.9999 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -42 -42 1.24998 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.24998 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -42 -42 1.24998 42 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.24998 42 -0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 -0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 +0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 41.9999 +0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 42 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -41.9999 -41.9999 2.50004 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.5 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -41.9999 -41.9999 2.50004 41.9999 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.5 42 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42.0001 -42.0001 1.25005 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42.0001 -42.0001 1.25 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42.0001 -42.0001 1.25005 42.0001 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42.0001 -42.0001 1.25 42.0001 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.5 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.5 42 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 41.9999 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42.0001 -42.0001 1.25 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42.0001 -42.0001 1.25001 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42.0001 -42.0001 1.25 42.0001 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42.0001 -42.0001 1.25001 42.0001 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42 -42 1.24998 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42 -42 1.24999 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42 -42 1.24998 42 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42 -42 1.24999 42 -0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 -0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 +0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 42 +0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 42 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42 -42 1.24999 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42 -42 1.25001 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42 -42 1.24999 42 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42 -42 1.25001 42 -0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 -0.5 1 1 0 -4.75838e-09 0 -4.75838e-09 -42 -42 -1.72621e-07 +0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 42 +0.5 1 1 0 -4.7584e-09 0 -4.7584e-09 -42 -42 -1.72621e-07 42 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00006 -42 -42 5.00005 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00005 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00006 -42 -42 5.00005 42 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00005 42 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75002 -42 -42 3.75004 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.75001 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75002 -42 -42 3.75004 42 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.75001 42 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00005 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00002 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00005 42 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00002 42 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.75001 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75004 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.75001 42 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75004 42 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75002 -42 -42 3.75004 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.74997 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75002 -42 -42 3.75004 42 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.74997 42 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50002 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 41.9999 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50002 42 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.74997 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75005 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.74997 42 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75005 42 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50002 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -42 -42 2.49998 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50002 42 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -42 -42 2.49998 42 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00002 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -41.9999 -41.9999 5.00006 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00002 42 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -41.9999 -41.9999 5.00006 41.9999 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75004 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -42.0001 -42.0001 3.74998 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75004 42 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -42.0001 -42.0001 3.74998 42.0001 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -41.9999 -41.9999 5.00006 -1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -41.9999 -41.9999 5.00006 41.9999 +1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 42 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -42 -42 3.74998 -1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -42 -42 3.74998 42 +1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 42 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75005 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -41.9999 -41.9999 3.75004 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75005 42 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -41.9999 -41.9999 3.75004 41.9999 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -42 -42 2.49998 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42 -42 2.5 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -42 -42 2.49998 42 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42 -42 2.5 42 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -41.9999 -41.9999 3.75004 -1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -41.9999 -41.9999 3.75004 41.9999 +1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 42 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42 -42 2.5 -1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42 -42 2.5 42 +1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 42 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50005 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 41.9999 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50005 42 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42.0001 -42.0001 1.25001 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42.0001 -42.0001 1.25002 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42.0001 -42.0001 1.25001 42.0001 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42.0001 -42.0001 1.25002 42.0001 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50005 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -41.9999 -41.9999 2.50001 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50005 42 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -41.9999 -41.9999 2.50001 41.9999 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42.0001 -42.0001 1.25002 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42.0001 -42.0001 1.24998 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42.0001 -42.0001 1.25002 42.0001 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42.0001 -42.0001 1.24998 42.0001 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42 -42 1.25001 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42 -42 1.25003 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42 -42 1.25001 42 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42 -42 1.25003 42 -0.5 1 1 0 -4.75838e-09 0 -4.75838e-09 -42 -42 -1.72621e-07 -0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 +0.5 1 1 0 -4.7584e-09 0 -4.7584e-09 -42 -42 -1.72621e-07 42 +0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 42 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42 -42 1.25003 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42 -42 1.25005 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42 -42 1.25003 42 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42 -42 1.25005 42 -0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 -0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 +0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 42 +0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 42 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -41.9999 -41.9999 2.50001 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42.0001 -42.0001 2.50006 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -41.9999 -41.9999 2.50001 41.9999 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42.0001 -42.0001 2.50006 42.0001 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.24998 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -41.9999 -41.9999 1.25003 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.24998 42 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -41.9999 -41.9999 1.25003 41.9999 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42.0001 -42.0001 2.50006 -1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42.0001 -42.0001 2.50006 42.0001 +1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 42 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -41.9999 -41.9999 1.25003 -1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -41.9999 -41.9999 1.25003 41.9999 +1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 42 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.25005 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -42 -42 1.25005 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.25005 42 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -42 -42 1.25005 42 -0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 -0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 +0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 42 +0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 41.9999 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -42 -42 1.25005 -1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -42 -42 1.25005 42 +1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 42 -0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 -1 1 1 0 0.000421865 0 0.000421865 -42 -42 0.000421865 +0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 41.9999 +1 1 1 0 0.000421865 0 0.000421865 -42 -42 0.000421865 42 diff --git a/tests/visualization_stress_pointwise/solution/solution-00000.0000.gnuplot b/tests/visualization_stress_pointwise/solution/solution-00000.0000.gnuplot index 2ae6fd0a7ab..1639af82d4b 100644 --- a/tests/visualization_stress_pointwise/solution/solution-00000.0000.gnuplot +++ b/tests/visualization_stress_pointwise/solution/solution-00000.0000.gnuplot @@ -1,455 +1,455 @@ # This file was generated by the deal.II library. - - +# Date = 2021/9/30 +# Time = 11:49:29 # # For a description of the GNUPLOT format see the GNUPLOT manual. # -#

-0 0 0 0 9.99962 1 9.99962 -42 -42 9.99962 -0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 +#

+0 0 0 0 9.99962 1 9.99962 -42 -42 9.99962 42 +0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 41.9999 -0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42.0001 -42.0001 8.75001 +0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 41.9999 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42.0001 -42.0001 8.75001 42.0001 -0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 -0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 +0.125 0 0 0 10 0 10 -41.9999 -41.9999 10.0001 41.9999 +0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42.0001 -42.0001 8.75001 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75011 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42.0001 -42.0001 8.75001 42.0001 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75011 42 -0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42 -42 8.75005 +0 0.125 0.125 0 8.75005 1 8.7501 -41.9999 -41.9999 8.75005 41.9999 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.74999 -42 -42 8.75005 42 -0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -42 -42 7.50001 +0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 42 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -42 -42 7.50001 42 -0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42 -42 8.75005 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75009 +0.125 0.125 0.125 -5.82922e-08 8.75 0.125 8.75006 -42 -42 8.75005 42 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75 -42 -42 8.75009 42 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50001 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50008 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50001 42 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50008 42 -0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 -0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 +0.25 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 +0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75011 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.7501 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75011 42 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.7501 42 -0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 -0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 +0.375 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 +0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.7501 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.7501 42 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 -0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75009 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.75009 +0.25 0.125 0.125 -4.65763e-08 8.75006 0.25 8.75001 -42 -42 8.75009 42 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75009 -42 -42 8.75009 42 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50008 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50008 42 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 42 -0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.75009 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 +0.375 0.125 0.125 -2.93476e-08 8.75008 0.375 8.75013 -42 -42 8.75009 42 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 42 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 42 -0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -41.9999 -41.9999 7.50004 +0 0.25 0.25 0 7.50007 1 7.50003 -42 -42 7.50007 42 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50006 -41.9999 -41.9999 7.50004 41.9999 -0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -42 -42 6.25003 +0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 42 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -42 -42 6.25003 42 -0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50004 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50007 +0.125 0.25 0.25 -6.9154e-08 7.50005 0.125 7.50016 -42 -42 7.50004 42 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50006 -42 -42 7.50007 42 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42.0001 -42.0001 6.25003 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25007 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42.0001 -42.0001 6.25003 42.0001 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25007 42 -0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -41.9999 -41.9999 6.25002 +0 0.375 0.375 0 6.25007 1 6.25008 -42 -42 6.25007 42 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25 -41.9999 -41.9999 6.25002 41.9999 -0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -42 -42 5.00008 +0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 42 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -42 -42 5.00008 42 -0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42 -42 6.25002 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25009 +0.125 0.375 0.375 -4.42576e-08 6.25006 0.125 6.25015 -42 -42 6.25002 42 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25007 -42 -42 6.25009 42 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -42 -42 5.00008 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00004 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -42 -42 5.00008 42 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00004 42 -0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50007 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 +0.25 0.25 0.25 -5.00963e-08 7.5001 0.25 7.50007 -42 -42 7.50007 42 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50011 -42 -42 7.50004 42 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25007 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25004 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25007 42 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25004 42 -0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 +0.375 0.25 0.25 -1.37711e-08 7.50007 0.375 7.50012 -42 -42 7.50004 42 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50009 -42 -42 7.50008 42 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25004 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25004 42 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 -0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25009 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25003 +0.25 0.375 0.375 -1.2643e-08 6.25009 0.25 6.25006 -42 -42 6.25009 42 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25011 -42 -42 6.25003 42 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00004 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00003 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00004 42 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00003 42 -0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25003 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 +0.375 0.375 0.375 1.40865e-08 6.25005 0.375 6.25008 -42 -42 6.25003 42 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00003 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00004 -42 -42 5.00005 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00003 42 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00004 -42 -42 5.00005 42 -0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 -0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 +0.5 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 +0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75007 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75007 42 -0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 -0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 +0.625 0 0 0 10.0001 0 10.0001 -42 -42 10.0001 42 +0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75007 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75007 42 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 42 -0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75008 +0.5 0.125 0.125 1.24646e-10 8.75009 0.5 8.75009 -42 -42 8.75009 42 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75004 -42 -42 8.75008 42 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 42 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 42 -0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75008 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 +0.625 0.125 0.125 2.95971e-08 8.7501 0.625 8.75009 -42 -42 8.75008 42 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 42 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50007 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 42 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50007 42 -0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 -0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 +0.75 0 0 0 10.0001 0 10.0001 -41.9999 -41.9999 10.0001 41.9999 +0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 41.9999 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42.0001 -42.0001 8.75016 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75007 42 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42.0001 -42.0001 8.75016 42.0001 -0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 -1 0 0 0 10.0006 0 10.0006 -42 -42 10.0006 +0.875 0 0 0 10.0002 0 10.0002 -41.9999 -41.9999 10.0001 41.9999 +1 0 0 0 10.0006 0 10.0006 -42 -42 10.0006 42 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42.0001 -42.0001 8.75016 -1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42.0001 -42.0001 8.75016 42.0001 +1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 41.9999 -0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42 -42 8.75013 +0.75 0.125 0.125 4.68274e-08 8.75011 0.75 8.75017 -42 -42 8.75009 42 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75011 -42 -42 8.75013 42 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50007 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50014 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50007 42 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50014 42 -0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42 -42 8.75013 -1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 +0.875 0.125 0.125 5.85464e-08 8.75018 0.875 8.75019 -42 -42 8.75013 42 +1 0.125 0.125 0 8.75012 0 8.75007 -41.9999 -41.9999 8.75012 41.9999 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -42 -42 7.50014 -1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -42 -42 7.50014 42 +1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 42 -0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 +0.5 0.25 0.25 -3.36368e-11 7.50008 0.5 7.50006 -42 -42 7.50008 42 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50003 -42 -42 7.50011 42 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25008 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25008 42 -0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50008 +0.625 0.25 0.25 1.37039e-08 7.50008 0.625 7.50004 -42 -42 7.50011 42 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50008 -42 -42 7.50008 42 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25008 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25005 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25008 42 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25005 42 -0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25009 +0.5 0.375 0.375 9.4033e-12 6.25006 0.5 6.25006 -42 -42 6.25006 42 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25005 -42 -42 6.25009 42 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00006 -42 -42 5.00005 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00007 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00006 -42 -42 5.00005 42 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00007 42 -0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25009 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25004 +0.625 0.375 0.375 -1.40676e-08 6.25008 0.625 6.25001 -42 -42 6.25009 42 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25007 -42 -42 6.25004 42 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00007 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00006 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00007 42 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00006 42 -0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50008 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50011 +0.75 0.25 0.25 5.00295e-08 7.50005 0.75 7.50009 -42 -42 7.50008 42 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.49999 -42 -42 7.50011 42 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25005 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42.0001 -42.0001 6.25009 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25005 42 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42.0001 -42.0001 6.25009 42.0001 -0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -41.9999 -41.9999 7.50011 -1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 +0.875 0.25 0.25 6.9086e-08 7.5001 0.875 7.50009 -41.9999 -41.9999 7.50011 41.9999 +1 0.25 0.25 0 7.50008 0 7.50012 -42 -42 7.50008 42 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -42 -42 6.25009 -1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -42 -42 6.25009 42 +1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 42 -0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25004 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42 -42 6.25011 +0.75 0.375 0.375 1.26619e-08 6.25004 0.75 6.25005 -42 -42 6.25004 42 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.24998 -42 -42 6.25011 42 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00006 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -42 -42 5.00002 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00006 42 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -42 -42 5.00002 42 -0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -41.9999 -41.9999 6.25011 -1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 +0.875 0.375 0.375 4.42765e-08 6.25007 0.875 6.25013 -41.9999 -41.9999 6.25011 41.9999 +1 0.375 0.375 0 6.25006 0 6.25005 -42 -42 6.25006 42 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -42 -42 5.00002 -1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -42 -42 5.00002 42 +1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 42 -0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -41.9999 -41.9999 5.00004 +0 0.5 0.5 0 5.00001 1 5.00002 -42 -42 5.00001 42 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00001 -41.9999 -41.9999 5.00004 41.9999 -0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -42 -42 3.7501 +0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 42 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -42 -42 3.7501 42 -0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -41.9999 -41.9999 5.00004 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00008 +0.125 0.5 0.5 -3.24821e-08 5.00006 0.125 5.00016 -41.9999 -41.9999 5.00004 41.9999 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00003 -42 -42 5.00008 42 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -42.0001 -42.0001 3.7501 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -42.0001 -42.0001 3.7501 42.0001 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 42 -0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -41.9999 -41.9999 3.75004 +0 0.625 0.625 0 3.75002 1 3.75004 -42 -42 3.75002 42 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75003 -41.9999 -41.9999 3.75004 41.9999 -0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42 -42 2.50005 +0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 42 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42 -42 2.50005 42 -0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -41.9999 -41.9999 3.75004 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 +0.125 0.625 0.625 -4.01156e-08 3.75007 0.125 3.75018 -41.9999 -41.9999 3.75004 41.9999 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.74995 -42 -42 3.75003 42 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42 -42 2.50005 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -42 -42 2.50008 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42 -42 2.50005 42 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -42 -42 2.50008 42 -0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00008 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00005 +0.25 0.5 0.5 6.47671e-09 5.00007 0.25 5.00004 -42 -42 5.00008 42 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00008 -42 -42 5.00005 42 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75006 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 42 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75006 42 -0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00005 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00004 -42 -42 5.00005 +0.375 0.5 0.5 3.64313e-08 5.00004 0.375 5.00001 -42 -42 5.00005 42 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00004 -42 -42 5.00005 42 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75006 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75006 -42 -42 3.75004 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75006 42 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75006 -42 -42 3.75004 42 -0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75011 +0.25 0.625 0.625 9.45208e-10 3.75003 0.25 3.75006 -42 -42 3.75003 42 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74998 -42 -42 3.75011 42 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -42 -42 2.50008 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.50003 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -42 -42 2.50008 42 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.50003 42 -0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75011 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75006 -42 -42 3.75004 +0.375 0.625 0.625 -2.70802e-09 3.75002 0.375 3.74994 -42 -42 3.75011 42 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75006 -42 -42 3.75004 42 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.50003 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.50003 42 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 41.9999 -0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42.0001 -42.0001 2.49999 +0 0.75 0.75 0 2.49998 1 2.49998 -42 -42 2.49998 42 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50015 -42.0001 -42.0001 2.49999 42.0001 -0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -41.9999 -41.9999 1.25 +0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 42 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -41.9999 -41.9999 1.25 41.9999 -0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42.0001 -42.0001 2.49999 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -41.9999 -41.9999 2.50004 +0.125 0.75 0.75 5.65515e-09 2.50008 0.125 2.50007 -42.0001 -42.0001 2.49999 42.0001 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.49992 -41.9999 -41.9999 2.50004 41.9999 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -41.9999 -41.9999 1.25 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.25005 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -41.9999 -41.9999 1.25 41.9999 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.25005 42 -0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -42 -42 1.24998 +0 0.875 0.875 0 1.25022 1 1.25008 -42 -42 1.25022 42 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.2501 -42 -42 1.24998 42 -0 1 1 0 -0.000421673 1 -0.000421673 -42 -42 -0.000421673 -0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 +0 1 1 0 -0.000421673 1 -0.000421673 -42 -42 -0.000421673 42 +0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 41.9999 -0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -42 -42 1.24998 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.24998 +0.125 0.875 0.875 1.95172e-08 1.24993 0.125 1.24987 -42 -42 1.24998 42 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25005 -42 -42 1.24998 42 -0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 -0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 +0.125 1 1 0 5.77065e-05 0 5.77065e-05 -41.9999 -41.9999 3.4755e-05 41.9999 +0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 42 -0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -41.9999 -41.9999 2.50004 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.5 +0.25 0.75 0.75 -4.43965e-08 2.50001 0.25 2.50008 -41.9999 -41.9999 2.50004 41.9999 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.49997 -42 -42 2.5 42 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42.0001 -42.0001 1.25005 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42.0001 -42.0001 1.25 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42.0001 -42.0001 1.25005 42.0001 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42.0001 -42.0001 1.25 42.0001 -0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.5 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 +0.375 0.75 0.75 -5.74409e-08 2.50005 0.375 2.50005 -42 -42 2.5 42 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.50006 -41.9999 -41.9999 2.50003 41.9999 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42.0001 -42.0001 1.25 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42.0001 -42.0001 1.25001 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42.0001 -42.0001 1.25 42.0001 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42.0001 -42.0001 1.25001 42.0001 -0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42 -42 1.24998 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42 -42 1.24999 +0.25 0.875 0.875 -5.799e-08 1.25006 0.25 1.25009 -42 -42 1.24998 42 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.24999 -42 -42 1.24999 42 -0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 -0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 +0.25 1 1 0 -2.79584e-05 0 -2.79584e-05 -42 -42 -2.03641e-05 42 +0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 42 -0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42 -42 1.24999 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42 -42 1.25001 +0.375 0.875 0.875 -2.38002e-08 1.25 0.375 1.25017 -42 -42 1.24999 42 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.24988 -42 -42 1.25001 42 -0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 -0.5 1 1 0 -4.75838e-09 0 -4.75838e-09 -42 -42 -1.72621e-07 +0.375 1 1 0 1.91583e-05 0 1.91583e-05 -42 -42 -2.74365e-06 42 +0.5 1 1 0 -4.7584e-09 0 -4.7584e-09 -42 -42 -1.72621e-07 42 -0.5 0.5 0.5 -4.99415e-12 5.00005 0.5 5.00006 -42 -42 5.00005 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00005 +0.5 0.5 0.5 -4.99418e-12 5.00005 0.5 5.00006 -42 -42 5.00005 42 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00009 -42 -42 5.00005 42 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75002 -42 -42 3.75004 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.75001 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75002 -42 -42 3.75004 42 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.75001 42 -0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00005 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00002 +0.625 0.5 0.5 -3.64412e-08 5.00006 0.625 5.00002 -42 -42 5.00005 42 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00006 -42 -42 5.00002 42 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.75001 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75004 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.75001 42 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75004 42 -0.5 0.625 0.625 9.4032e-12 3.75004 0.5 3.75002 -42 -42 3.75004 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.74997 +0.5 0.625 0.625 9.40317e-12 3.75004 0.5 3.75002 -42 -42 3.75004 42 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.75014 -42 -42 3.74997 42 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50002 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 41.9999 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50002 42 -0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.74997 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75005 +0.625 0.625 0.625 2.72692e-09 3.75006 0.625 3.7501 -42 -42 3.74997 42 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75002 -42 -42 3.75005 42 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50002 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -42 -42 2.49998 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50002 42 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -42 -42 2.49998 42 -0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00002 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -41.9999 -41.9999 5.00006 +0.75 0.5 0.5 -6.48617e-09 5.00003 0.75 5.00007 -42 -42 5.00002 42 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 4.99994 -41.9999 -41.9999 5.00006 41.9999 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75004 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -42.0001 -42.0001 3.74998 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75004 42 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -42.0001 -42.0001 3.74998 42.0001 -0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -41.9999 -41.9999 5.00006 -1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 +0.875 0.5 0.5 3.24734e-08 5.00004 0.875 5.00009 -41.9999 -41.9999 5.00006 41.9999 +1 0.5 0.5 0 5.00009 0 5.00008 -42 -42 5.00009 42 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -42 -42 3.74998 -1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -42 -42 3.74998 42 +1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 42 -0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75005 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -41.9999 -41.9999 3.75004 +0.75 0.625 0.625 -9.2623e-10 3.75004 0.75 3.75012 -42 -42 3.75005 42 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.7499 -41.9999 -41.9999 3.75004 41.9999 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -42 -42 2.49998 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42 -42 2.5 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -42 -42 2.49998 42 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42 -42 2.5 42 -0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -41.9999 -41.9999 3.75004 -1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 +0.875 0.625 0.625 4.01345e-08 3.75001 0.875 3.75004 -41.9999 -41.9999 3.75004 41.9999 +1 0.625 0.625 0 3.75005 0 3.75004 -42 -42 3.75005 42 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42 -42 2.5 -1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42 -42 2.5 42 +1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 42 -0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50005 +0.5 0.75 0.75 -3.36369e-11 2.50003 0.5 2.49999 -41.9999 -41.9999 2.50003 41.9999 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.5 -42 -42 2.50005 42 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42.0001 -42.0001 1.25001 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42.0001 -42.0001 1.25002 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42.0001 -42.0001 1.25001 42.0001 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42.0001 -42.0001 1.25002 42.0001 -0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50005 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -41.9999 -41.9999 2.50001 +0.625 0.75 0.75 5.73737e-08 2.5 0.625 2.50008 -42 -42 2.50005 42 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.49997 -41.9999 -41.9999 2.50001 41.9999 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42.0001 -42.0001 1.25002 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42.0001 -42.0001 1.24998 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42.0001 -42.0001 1.25002 42.0001 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42.0001 -42.0001 1.24998 42.0001 -0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42 -42 1.25001 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42 -42 1.25003 +0.5 0.875 0.875 1.24646e-10 1.25001 0.5 1.25015 -42 -42 1.25001 42 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.24986 -42 -42 1.25003 42 -0.5 1 1 0 -4.75838e-09 0 -4.75838e-09 -42 -42 -1.72621e-07 -0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 +0.5 1 1 0 -4.7584e-09 0 -4.7584e-09 -42 -42 -1.72621e-07 42 +0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 42 -0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42 -42 1.25003 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42 -42 1.25005 +0.625 0.875 0.875 2.40497e-08 1.25002 0.625 1.25004 -42 -42 1.25003 42 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24993 -42 -42 1.25005 42 -0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 -0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 +0.625 1 1 0 -1.91746e-05 0 -1.91746e-05 -42 -42 2.39356e-06 42 +0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 42 -0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -41.9999 -41.9999 2.50001 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42.0001 -42.0001 2.50006 +0.75 0.75 0.75 4.43296e-08 2.50005 0.75 2.50013 -41.9999 -41.9999 2.50001 41.9999 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.49998 -42.0001 -42.0001 2.50006 42.0001 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.24998 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -41.9999 -41.9999 1.25003 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.24998 42 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -41.9999 -41.9999 1.25003 41.9999 -0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42.0001 -42.0001 2.50006 -1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 +0.875 0.75 0.75 -5.72309e-09 2.49997 0.875 2.4999 -42.0001 -42.0001 2.50006 42.0001 +1 0.75 0.75 0 2.50007 0 2.50007 -42 -42 2.50007 42 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -41.9999 -41.9999 1.25003 -1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -41.9999 -41.9999 1.25003 41.9999 +1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 42 -0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.25005 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -42 -42 1.25005 +0.75 0.875 0.875 5.82411e-08 1.24997 0.75 1.24997 -42 -42 1.25005 42 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.25015 -42 -42 1.25005 42 -0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 -0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 +0.75 1 1 0 2.79602e-05 0 2.79602e-05 -42 -42 2.00264e-05 42 +0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 41.9999 -0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -42 -42 1.25005 -1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 +0.875 0.875 0.875 -1.92631e-08 1.2501 0.875 1.24993 -42 -42 1.25005 42 +1 0.875 0.875 0 1.2498 0 1.24995 -42 -42 1.2498 42 -0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 -1 1 1 0 0.000421865 0 0.000421865 -42 -42 0.000421865 +0.875 1 1 0 -5.77834e-05 0 -5.77834e-05 -41.9999 -41.9999 -3.51986e-05 41.9999 +1 1 1 0 0.000421865 0 0.000421865 -42 -42 0.000421865 42 From f5dcc2bdee625bdec41d22c1175b315312cf597a Mon Sep 17 00:00:00 2001 From: anne-glerum Date: Tue, 12 Oct 2021 10:44:56 +0200 Subject: [PATCH 4/5] Fix elastic stresses --- .../visualization/stress_second_invariant.cc | 21 +++++++++++-------- 1 file changed, 12 insertions(+), 9 deletions(-) diff --git a/source/postprocess/visualization/stress_second_invariant.cc b/source/postprocess/visualization/stress_second_invariant.cc index f0c1f6ec265..92f46ed2468 100644 --- a/source/postprocess/visualization/stress_second_invariant.cc +++ b/source/postprocess/visualization/stress_second_invariant.cc @@ -67,16 +67,8 @@ namespace aspect for (unsigned int q = 0; q < n_quadrature_points; ++q) { - const SymmetricTensor<2, dim> strain_rate = in.strain_rate[q]; - const SymmetricTensor<2, dim> deviatoric_strain_rate = (this->get_material_model().is_compressible() - ? strain_rate - 1. / 3 * trace(strain_rate) * unit_symmetric_tensor() - : strain_rate); - - const double eta = out.viscosities[q]; - // Compressive stress is positive in geoscience applications. - SymmetricTensor<2, dim> stress = -2. * eta * deviatoric_strain_rate + - in.pressure[q] * unit_symmetric_tensor(); + SymmetricTensor<2, dim> stress = in.pressure[q] * unit_symmetric_tensor(); // Add elastic stresses if existent. if (this->get_parameters().enable_elasticity == true) @@ -92,6 +84,17 @@ namespace aspect stress[1][2] += in.composition[q][this->introspection().compositional_index_for_name("ve_stress_yz")]; } } + else + { + const SymmetricTensor<2, dim> strain_rate = in.strain_rate[q]; + const SymmetricTensor<2, dim> deviatoric_strain_rate = (this->get_material_model().is_compressible() + ? strain_rate - 1. / 3 * trace(strain_rate) * unit_symmetric_tensor() + : strain_rate); + + const double eta = out.viscosities[q]; + + stress += -2. * eta * deviatoric_strain_rate; + } // Compute the deviatoric stress tensor after elastic stresses were added. const SymmetricTensor<2, dim> deviatoric_stress = deviator(stress); From 8f0551954a9f0242a3da625bec57a9473b34dff4 Mon Sep 17 00:00:00 2001 From: anne-glerum Date: Fri, 22 Oct 2021 10:20:46 +0200 Subject: [PATCH 5/5] Add changelog --- doc/modules/changes/20211022_glerum | 4 ++++ 1 file changed, 4 insertions(+) create mode 100644 doc/modules/changes/20211022_glerum diff --git a/doc/modules/changes/20211022_glerum b/doc/modules/changes/20211022_glerum new file mode 100644 index 00000000000..aa295bc03f4 --- /dev/null +++ b/doc/modules/changes/20211022_glerum @@ -0,0 +1,4 @@ +New: There is now a postprocessor that computes the +second invariant of the deviatoric stress tensor. +
+(Anne Glerum, 2021/10/22)