From: Wolfgang Bangerth Date: Wed, 21 Dec 2016 17:40:02 +0000 (-0700) Subject: Provide the cell via DataOutRotation. Add tests. X-Git-Tag: v8.5.0-rc1~308^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9bcff3f5569c6ef739859e9259094a07382dffe4;p=dealii.git Provide the cell via DataOutRotation. Add tests. --- diff --git a/source/numerics/data_out_rotation.cc b/source/numerics/data_out_rotation.cc index 36c4aca01a..f05fb7b7ed 100644 --- a/source/numerics/data_out_rotation.cc +++ b/source/numerics/data_out_rotation.cc @@ -216,6 +216,12 @@ build_one_patch (const cell_iterator if (update_flags & update_quadrature_points) data.patch_values_scalar.evaluation_points = fe_patch_values.get_quadrature_points(); + const typename DoFHandlerType::active_cell_iterator dh_cell(&(*cell)->get_triangulation(), + (*cell)->level(), + (*cell)->index(), + this->dof_data[dataset]->dof_handler); + data.patch_values_scalar.template set_cell (dh_cell); + postprocessor-> evaluate_scalar_field(data.patch_values_scalar, data.postprocessed_values[dataset]); @@ -239,7 +245,12 @@ build_one_patch (const cell_iterator if (update_flags & update_quadrature_points) data.patch_values_system.evaluation_points = fe_patch_values.get_quadrature_points(); - std::vector > dummy_normals; + const typename DoFHandlerType::active_cell_iterator dh_cell(&(*cell)->get_triangulation(), + (*cell)->level(), + (*cell)->index(), + this->dof_data[dataset]->dof_handler); + data.patch_values_system.template set_cell (dh_cell); + postprocessor-> evaluate_vector_field(data.patch_values_system, data.postprocessed_values[dataset]); diff --git a/tests/numerics/data_out_rotation_postprocessor_scalar_02.cc b/tests/numerics/data_out_rotation_postprocessor_scalar_02.cc new file mode 100644 index 0000000000..fc7f2f13f0 --- /dev/null +++ b/tests/numerics/data_out_rotation_postprocessor_scalar_02.cc @@ -0,0 +1,125 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE at +// the top level of the deal.II distribution. +// +// --------------------------------------------------------------------- + + +// tests DataPostprocessor: access the cell we are currently working +// on for a scalar finite element field and DataOutRotation + + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include +#include + +#include + + +std::ofstream logfile("output"); + + +template +class MyPostprocessor : public DataPostprocessorScalar +{ +public: + MyPostprocessor () + : + DataPostprocessorScalar ("data", update_values) + {} + + virtual + void + evaluate_scalar_field (const DataPostprocessorInputs::Scalar &input_data, + std::vector > &computed_quantities) const + { + for (unsigned int q=0; q::cell_iterator + cell = input_data.template get_cell >(); + + Assert (input_data.solution_values[q] + == + double(cell->active_cell_index()), + ExcInternalError()); + + computed_quantities[q][0] = input_data.solution_values[q]; + } + } +}; + + + +template +void test () +{ + Triangulation triangulation; + FE_DGQ fe(1); + DoFHandler dof_handler(triangulation); + + GridGenerator::hyper_cube (triangulation, 0, 1); + triangulation.refine_global (1); + triangulation.begin_active()->set_refine_flag (); + triangulation.execute_coarsening_and_refinement (); + + dof_handler.distribute_dofs (fe); + + // create a vector with piecewise constants, and set the elements of + // the vector so that on each cell the field has values equal to the + // cell's active_cell_index + Vector solution (dof_handler.n_dofs()); + for (typename DoFHandler::active_cell_iterator cell = dof_handler.begin_active(); + cell != dof_handler.end(); ++cell) + { + std::vector local_dof_indices (cell->get_fe().dofs_per_cell); + cell->get_dof_indices (local_dof_indices); + for (unsigned int i=0; iactive_cell_index(); + } + + MyPostprocessor p; + DataOutRotation data_out; + data_out.attach_dof_handler (dof_handler); + data_out.add_data_vector (solution, p); + data_out.build_patches (4); + data_out.write_gnuplot (logfile); +} + + +int main () +{ + logfile << std::setprecision(2); + deallog << std::setprecision(2); + + test<1>(); + test<2>(); + + return 0; +} diff --git a/tests/numerics/data_out_rotation_postprocessor_scalar_02.output b/tests/numerics/data_out_rotation_postprocessor_scalar_02.output new file mode 100644 index 0000000000..dd7ae05a8a --- /dev/null +++ b/tests/numerics/data_out_rotation_postprocessor_scalar_02.output @@ -0,0 +1,1442 @@ +# This file was generated by the deal.II library. + + +# +# For a description of the GNUPLOT format see the GNUPLOT manual. +# +# +0.5 0 0 +1 0 0 + +3.1e-17 0.5 0 +6.1e-17 1 0 + + +3.1e-17 0.5 0 +6.1e-17 1 0 + +-0.5 6.1e-17 0 +-1 1.2e-16 0 + + +-0.5 6.1e-17 0 +-1 1.2e-16 0 + +-9.2e-17 -0.5 0 +-1.8e-16 -1 0 + + +-9.2e-17 -0.5 0 +-1.8e-16 -1 0 + +0.5 -1.2e-16 0 +1 -2.4e-16 0 + + +0 0 1 +0.25 0 1 + +0 0 1 +1.5e-17 0.25 1 + + +0 0 1 +1.5e-17 0.25 1 + +0 0 1 +-0.25 3.1e-17 1 + + +-0 0 1 +-0.25 3.1e-17 1 + +-0 0 1 +-4.6e-17 -0.25 1 + + +0 -0 1 +-4.6e-17 -0.25 1 + +0 -0 1 +0.25 -6.1e-17 1 + + +0.25 0 2 +0.5 0 2 + +1.5e-17 0.25 2 +3.1e-17 0.5 2 + + +1.5e-17 0.25 2 +3.1e-17 0.5 2 + +-0.25 3.1e-17 2 +-0.5 6.1e-17 2 + + +-0.25 3.1e-17 2 +-0.5 6.1e-17 2 + +-4.6e-17 -0.25 2 +-9.2e-17 -0.5 2 + + +-4.6e-17 -0.25 2 +-9.2e-17 -0.5 2 + +0.25 -6.1e-17 2 +0.5 -1.2e-16 2 + + +# This file was generated by the deal.II library. + + +# +# For a description of the GNUPLOT format see the GNUPLOT manual. +# +# +0 0.5 0 0 +0 1 0 0 + + +0 0.5 0 0 +0.5 0.5 0 0 + + +0 0.5 0 0 +0 3.1e-17 0.5 0 + + +0 1 0 0 +0.5 1 0 0 + + +0 1 0 0 +0 6.1e-17 1 0 + + +0.5 0.5 0 0 +0.5 1 0 0 + + +0.5 0.5 0 0 +0.5 3.1e-17 0.5 0 + + +0.5 1 0 0 +0.5 6.1e-17 1 0 + + +0 3.1e-17 0.5 0 +0 6.1e-17 1 0 + + +0 3.1e-17 0.5 0 +0.5 3.1e-17 0.5 0 + + +0 6.1e-17 1 0 +0.5 6.1e-17 1 0 + + +0.5 3.1e-17 0.5 0 +0.5 6.1e-17 1 0 + + +0 3.1e-17 0.5 0 +0 6.1e-17 1 0 + + +0 3.1e-17 0.5 0 +0.5 3.1e-17 0.5 0 + + +0 3.1e-17 0.5 0 +0 -0.5 6.1e-17 0 + + +0 6.1e-17 1 0 +0.5 6.1e-17 1 0 + + +0 6.1e-17 1 0 +0 -1 1.2e-16 0 + + +0.5 3.1e-17 0.5 0 +0.5 6.1e-17 1 0 + + +0.5 3.1e-17 0.5 0 +0.5 -0.5 6.1e-17 0 + + +0.5 6.1e-17 1 0 +0.5 -1 1.2e-16 0 + + +0 -0.5 6.1e-17 0 +0 -1 1.2e-16 0 + + +0 -0.5 6.1e-17 0 +0.5 -0.5 6.1e-17 0 + + +0 -1 1.2e-16 0 +0.5 -1 1.2e-16 0 + + +0.5 -0.5 6.1e-17 0 +0.5 -1 1.2e-16 0 + + +0 -0.5 6.1e-17 0 +0 -1 1.2e-16 0 + + +0 -0.5 6.1e-17 0 +0.5 -0.5 6.1e-17 0 + + +0 -0.5 6.1e-17 0 +0 -9.2e-17 -0.5 0 + + +0 -1 1.2e-16 0 +0.5 -1 1.2e-16 0 + + +0 -1 1.2e-16 0 +0 -1.8e-16 -1 0 + + +0.5 -0.5 6.1e-17 0 +0.5 -1 1.2e-16 0 + + +0.5 -0.5 6.1e-17 0 +0.5 -9.2e-17 -0.5 0 + + +0.5 -1 1.2e-16 0 +0.5 -1.8e-16 -1 0 + + +0 -9.2e-17 -0.5 0 +0 -1.8e-16 -1 0 + + +0 -9.2e-17 -0.5 0 +0.5 -9.2e-17 -0.5 0 + + +0 -1.8e-16 -1 0 +0.5 -1.8e-16 -1 0 + + +0.5 -9.2e-17 -0.5 0 +0.5 -1.8e-16 -1 0 + + +0 -9.2e-17 -0.5 0 +0 -1.8e-16 -1 0 + + +0 -9.2e-17 -0.5 0 +0.5 -9.2e-17 -0.5 0 + + +0 -9.2e-17 -0.5 0 +0 0.5 -1.2e-16 0 + + +0 -1.8e-16 -1 0 +0.5 -1.8e-16 -1 0 + + +0 -1.8e-16 -1 0 +0 1 -2.4e-16 0 + + +0.5 -9.2e-17 -0.5 0 +0.5 -1.8e-16 -1 0 + + +0.5 -9.2e-17 -0.5 0 +0.5 0.5 -1.2e-16 0 + + +0.5 -1.8e-16 -1 0 +0.5 1 -2.4e-16 0 + + +0 0.5 -1.2e-16 0 +0 1 -2.4e-16 0 + + +0 0.5 -1.2e-16 0 +0.5 0.5 -1.2e-16 0 + + +0 1 -2.4e-16 0 +0.5 1 -2.4e-16 0 + + +0.5 0.5 -1.2e-16 0 +0.5 1 -2.4e-16 0 + + +0.5 0 0 1 +0.5 0.5 0 1 + + +0.5 0 0 1 +1 0 0 1 + + +0.5 0 0 1 +0.5 0 0 1 + + +0.5 0.5 0 1 +1 0.5 0 1 + + +0.5 0.5 0 1 +0.5 3.1e-17 0.5 1 + + +1 0 0 1 +1 0.5 0 1 + + +1 0 0 1 +1 0 0 1 + + +1 0.5 0 1 +1 3.1e-17 0.5 1 + + +0.5 0 0 1 +0.5 3.1e-17 0.5 1 + + +0.5 0 0 1 +1 0 0 1 + + +0.5 3.1e-17 0.5 1 +1 3.1e-17 0.5 1 + + +1 0 0 1 +1 3.1e-17 0.5 1 + + +0.5 0 0 1 +0.5 3.1e-17 0.5 1 + + +0.5 0 0 1 +1 0 0 1 + + +0.5 0 0 1 +0.5 0 0 1 + + +0.5 3.1e-17 0.5 1 +1 3.1e-17 0.5 1 + + +0.5 3.1e-17 0.5 1 +0.5 -0.5 6.1e-17 1 + + +1 0 0 1 +1 3.1e-17 0.5 1 + + +1 0 0 1 +1 0 0 1 + + +1 3.1e-17 0.5 1 +1 -0.5 6.1e-17 1 + + +0.5 0 0 1 +0.5 -0.5 6.1e-17 1 + + +0.5 0 0 1 +1 0 0 1 + + +0.5 -0.5 6.1e-17 1 +1 -0.5 6.1e-17 1 + + +1 0 0 1 +1 -0.5 6.1e-17 1 + + +0.5 -0 0 1 +0.5 -0.5 6.1e-17 1 + + +0.5 -0 0 1 +1 -0 0 1 + + +0.5 -0 0 1 +0.5 -0 0 1 + + +0.5 -0.5 6.1e-17 1 +1 -0.5 6.1e-17 1 + + +0.5 -0.5 6.1e-17 1 +0.5 -9.2e-17 -0.5 1 + + +1 -0 0 1 +1 -0.5 6.1e-17 1 + + +1 -0 0 1 +1 -0 0 1 + + +1 -0.5 6.1e-17 1 +1 -9.2e-17 -0.5 1 + + +0.5 -0 0 1 +0.5 -9.2e-17 -0.5 1 + + +0.5 -0 0 1 +1 -0 0 1 + + +0.5 -9.2e-17 -0.5 1 +1 -9.2e-17 -0.5 1 + + +1 -0 0 1 +1 -9.2e-17 -0.5 1 + + +0.5 0 -0 1 +0.5 -9.2e-17 -0.5 1 + + +0.5 0 -0 1 +1 0 -0 1 + + +0.5 0 -0 1 +0.5 0 -0 1 + + +0.5 -9.2e-17 -0.5 1 +1 -9.2e-17 -0.5 1 + + +0.5 -9.2e-17 -0.5 1 +0.5 0.5 -1.2e-16 1 + + +1 0 -0 1 +1 -9.2e-17 -0.5 1 + + +1 0 -0 1 +1 0 -0 1 + + +1 -9.2e-17 -0.5 1 +1 0.5 -1.2e-16 1 + + +0.5 0 -0 1 +0.5 0.5 -1.2e-16 1 + + +0.5 0 -0 1 +1 0 -0 1 + + +0.5 0.5 -1.2e-16 1 +1 0.5 -1.2e-16 1 + + +1 0 -0 1 +1 0.5 -1.2e-16 1 + + +0.5 0.5 0 2 +0.5 1 0 2 + + +0.5 0.5 0 2 +1 0.5 0 2 + + +0.5 0.5 0 2 +0.5 3.1e-17 0.5 2 + + +0.5 1 0 2 +1 1 0 2 + + +0.5 1 0 2 +0.5 6.1e-17 1 2 + + +1 0.5 0 2 +1 1 0 2 + + +1 0.5 0 2 +1 3.1e-17 0.5 2 + + +1 1 0 2 +1 6.1e-17 1 2 + + +0.5 3.1e-17 0.5 2 +0.5 6.1e-17 1 2 + + +0.5 3.1e-17 0.5 2 +1 3.1e-17 0.5 2 + + +0.5 6.1e-17 1 2 +1 6.1e-17 1 2 + + +1 3.1e-17 0.5 2 +1 6.1e-17 1 2 + + +0.5 3.1e-17 0.5 2 +0.5 6.1e-17 1 2 + + +0.5 3.1e-17 0.5 2 +1 3.1e-17 0.5 2 + + +0.5 3.1e-17 0.5 2 +0.5 -0.5 6.1e-17 2 + + +0.5 6.1e-17 1 2 +1 6.1e-17 1 2 + + +0.5 6.1e-17 1 2 +0.5 -1 1.2e-16 2 + + +1 3.1e-17 0.5 2 +1 6.1e-17 1 2 + + +1 3.1e-17 0.5 2 +1 -0.5 6.1e-17 2 + + +1 6.1e-17 1 2 +1 -1 1.2e-16 2 + + +0.5 -0.5 6.1e-17 2 +0.5 -1 1.2e-16 2 + + +0.5 -0.5 6.1e-17 2 +1 -0.5 6.1e-17 2 + + +0.5 -1 1.2e-16 2 +1 -1 1.2e-16 2 + + +1 -0.5 6.1e-17 2 +1 -1 1.2e-16 2 + + +0.5 -0.5 6.1e-17 2 +0.5 -1 1.2e-16 2 + + +0.5 -0.5 6.1e-17 2 +1 -0.5 6.1e-17 2 + + +0.5 -0.5 6.1e-17 2 +0.5 -9.2e-17 -0.5 2 + + +0.5 -1 1.2e-16 2 +1 -1 1.2e-16 2 + + +0.5 -1 1.2e-16 2 +0.5 -1.8e-16 -1 2 + + +1 -0.5 6.1e-17 2 +1 -1 1.2e-16 2 + + +1 -0.5 6.1e-17 2 +1 -9.2e-17 -0.5 2 + + +1 -1 1.2e-16 2 +1 -1.8e-16 -1 2 + + +0.5 -9.2e-17 -0.5 2 +0.5 -1.8e-16 -1 2 + + +0.5 -9.2e-17 -0.5 2 +1 -9.2e-17 -0.5 2 + + +0.5 -1.8e-16 -1 2 +1 -1.8e-16 -1 2 + + +1 -9.2e-17 -0.5 2 +1 -1.8e-16 -1 2 + + +0.5 -9.2e-17 -0.5 2 +0.5 -1.8e-16 -1 2 + + +0.5 -9.2e-17 -0.5 2 +1 -9.2e-17 -0.5 2 + + +0.5 -9.2e-17 -0.5 2 +0.5 0.5 -1.2e-16 2 + + +0.5 -1.8e-16 -1 2 +1 -1.8e-16 -1 2 + + +0.5 -1.8e-16 -1 2 +0.5 1 -2.4e-16 2 + + +1 -9.2e-17 -0.5 2 +1 -1.8e-16 -1 2 + + +1 -9.2e-17 -0.5 2 +1 0.5 -1.2e-16 2 + + +1 -1.8e-16 -1 2 +1 1 -2.4e-16 2 + + +0.5 0.5 -1.2e-16 2 +0.5 1 -2.4e-16 2 + + +0.5 0.5 -1.2e-16 2 +1 0.5 -1.2e-16 2 + + +0.5 1 -2.4e-16 2 +1 1 -2.4e-16 2 + + +1 0.5 -1.2e-16 2 +1 1 -2.4e-16 2 + + +0 0 0 3 +0 0.25 0 3 + + +0 0 0 3 +0.25 0 0 3 + + +0 0 0 3 +0 0 0 3 + + +0 0.25 0 3 +0.25 0.25 0 3 + + +0 0.25 0 3 +0 1.5e-17 0.25 3 + + +0.25 0 0 3 +0.25 0.25 0 3 + + +0.25 0 0 3 +0.25 0 0 3 + + +0.25 0.25 0 3 +0.25 1.5e-17 0.25 3 + + +0 0 0 3 +0 1.5e-17 0.25 3 + + +0 0 0 3 +0.25 0 0 3 + + +0 1.5e-17 0.25 3 +0.25 1.5e-17 0.25 3 + + +0.25 0 0 3 +0.25 1.5e-17 0.25 3 + + +0 0 0 3 +0 1.5e-17 0.25 3 + + +0 0 0 3 +0.25 0 0 3 + + +0 0 0 3 +0 0 0 3 + + +0 1.5e-17 0.25 3 +0.25 1.5e-17 0.25 3 + + +0 1.5e-17 0.25 3 +0 -0.25 3.1e-17 3 + + +0.25 0 0 3 +0.25 1.5e-17 0.25 3 + + +0.25 0 0 3 +0.25 0 0 3 + + +0.25 1.5e-17 0.25 3 +0.25 -0.25 3.1e-17 3 + + +0 0 0 3 +0 -0.25 3.1e-17 3 + + +0 0 0 3 +0.25 0 0 3 + + +0 -0.25 3.1e-17 3 +0.25 -0.25 3.1e-17 3 + + +0.25 0 0 3 +0.25 -0.25 3.1e-17 3 + + +0 -0 0 3 +0 -0.25 3.1e-17 3 + + +0 -0 0 3 +0.25 -0 0 3 + + +0 -0 0 3 +0 -0 0 3 + + +0 -0.25 3.1e-17 3 +0.25 -0.25 3.1e-17 3 + + +0 -0.25 3.1e-17 3 +0 -4.6e-17 -0.25 3 + + +0.25 -0 0 3 +0.25 -0.25 3.1e-17 3 + + +0.25 -0 0 3 +0.25 -0 0 3 + + +0.25 -0.25 3.1e-17 3 +0.25 -4.6e-17 -0.25 3 + + +0 -0 0 3 +0 -4.6e-17 -0.25 3 + + +0 -0 0 3 +0.25 -0 0 3 + + +0 -4.6e-17 -0.25 3 +0.25 -4.6e-17 -0.25 3 + + +0.25 -0 0 3 +0.25 -4.6e-17 -0.25 3 + + +0 0 -0 3 +0 -4.6e-17 -0.25 3 + + +0 0 -0 3 +0.25 0 -0 3 + + +0 0 -0 3 +0 0 -0 3 + + +0 -4.6e-17 -0.25 3 +0.25 -4.6e-17 -0.25 3 + + +0 -4.6e-17 -0.25 3 +0 0.25 -6.1e-17 3 + + +0.25 0 -0 3 +0.25 -4.6e-17 -0.25 3 + + +0.25 0 -0 3 +0.25 0 -0 3 + + +0.25 -4.6e-17 -0.25 3 +0.25 0.25 -6.1e-17 3 + + +0 0 -0 3 +0 0.25 -6.1e-17 3 + + +0 0 -0 3 +0.25 0 -0 3 + + +0 0.25 -6.1e-17 3 +0.25 0.25 -6.1e-17 3 + + +0.25 0 -0 3 +0.25 0.25 -6.1e-17 3 + + +0 0.25 0 4 +0 0.5 0 4 + + +0 0.25 0 4 +0.25 0.25 0 4 + + +0 0.25 0 4 +0 1.5e-17 0.25 4 + + +0 0.5 0 4 +0.25 0.5 0 4 + + +0 0.5 0 4 +0 3.1e-17 0.5 4 + + +0.25 0.25 0 4 +0.25 0.5 0 4 + + +0.25 0.25 0 4 +0.25 1.5e-17 0.25 4 + + +0.25 0.5 0 4 +0.25 3.1e-17 0.5 4 + + +0 1.5e-17 0.25 4 +0 3.1e-17 0.5 4 + + +0 1.5e-17 0.25 4 +0.25 1.5e-17 0.25 4 + + +0 3.1e-17 0.5 4 +0.25 3.1e-17 0.5 4 + + +0.25 1.5e-17 0.25 4 +0.25 3.1e-17 0.5 4 + + +0 1.5e-17 0.25 4 +0 3.1e-17 0.5 4 + + +0 1.5e-17 0.25 4 +0.25 1.5e-17 0.25 4 + + +0 1.5e-17 0.25 4 +0 -0.25 3.1e-17 4 + + +0 3.1e-17 0.5 4 +0.25 3.1e-17 0.5 4 + + +0 3.1e-17 0.5 4 +0 -0.5 6.1e-17 4 + + +0.25 1.5e-17 0.25 4 +0.25 3.1e-17 0.5 4 + + +0.25 1.5e-17 0.25 4 +0.25 -0.25 3.1e-17 4 + + +0.25 3.1e-17 0.5 4 +0.25 -0.5 6.1e-17 4 + + +0 -0.25 3.1e-17 4 +0 -0.5 6.1e-17 4 + + +0 -0.25 3.1e-17 4 +0.25 -0.25 3.1e-17 4 + + +0 -0.5 6.1e-17 4 +0.25 -0.5 6.1e-17 4 + + +0.25 -0.25 3.1e-17 4 +0.25 -0.5 6.1e-17 4 + + +0 -0.25 3.1e-17 4 +0 -0.5 6.1e-17 4 + + +0 -0.25 3.1e-17 4 +0.25 -0.25 3.1e-17 4 + + +0 -0.25 3.1e-17 4 +0 -4.6e-17 -0.25 4 + + +0 -0.5 6.1e-17 4 +0.25 -0.5 6.1e-17 4 + + +0 -0.5 6.1e-17 4 +0 -9.2e-17 -0.5 4 + + +0.25 -0.25 3.1e-17 4 +0.25 -0.5 6.1e-17 4 + + +0.25 -0.25 3.1e-17 4 +0.25 -4.6e-17 -0.25 4 + + +0.25 -0.5 6.1e-17 4 +0.25 -9.2e-17 -0.5 4 + + +0 -4.6e-17 -0.25 4 +0 -9.2e-17 -0.5 4 + + +0 -4.6e-17 -0.25 4 +0.25 -4.6e-17 -0.25 4 + + +0 -9.2e-17 -0.5 4 +0.25 -9.2e-17 -0.5 4 + + +0.25 -4.6e-17 -0.25 4 +0.25 -9.2e-17 -0.5 4 + + +0 -4.6e-17 -0.25 4 +0 -9.2e-17 -0.5 4 + + +0 -4.6e-17 -0.25 4 +0.25 -4.6e-17 -0.25 4 + + +0 -4.6e-17 -0.25 4 +0 0.25 -6.1e-17 4 + + +0 -9.2e-17 -0.5 4 +0.25 -9.2e-17 -0.5 4 + + +0 -9.2e-17 -0.5 4 +0 0.5 -1.2e-16 4 + + +0.25 -4.6e-17 -0.25 4 +0.25 -9.2e-17 -0.5 4 + + +0.25 -4.6e-17 -0.25 4 +0.25 0.25 -6.1e-17 4 + + +0.25 -9.2e-17 -0.5 4 +0.25 0.5 -1.2e-16 4 + + +0 0.25 -6.1e-17 4 +0 0.5 -1.2e-16 4 + + +0 0.25 -6.1e-17 4 +0.25 0.25 -6.1e-17 4 + + +0 0.5 -1.2e-16 4 +0.25 0.5 -1.2e-16 4 + + +0.25 0.25 -6.1e-17 4 +0.25 0.5 -1.2e-16 4 + + +0.25 0 0 5 +0.25 0.25 0 5 + + +0.25 0 0 5 +0.5 0 0 5 + + +0.25 0 0 5 +0.25 0 0 5 + + +0.25 0.25 0 5 +0.5 0.25 0 5 + + +0.25 0.25 0 5 +0.25 1.5e-17 0.25 5 + + +0.5 0 0 5 +0.5 0.25 0 5 + + +0.5 0 0 5 +0.5 0 0 5 + + +0.5 0.25 0 5 +0.5 1.5e-17 0.25 5 + + +0.25 0 0 5 +0.25 1.5e-17 0.25 5 + + +0.25 0 0 5 +0.5 0 0 5 + + +0.25 1.5e-17 0.25 5 +0.5 1.5e-17 0.25 5 + + +0.5 0 0 5 +0.5 1.5e-17 0.25 5 + + +0.25 0 0 5 +0.25 1.5e-17 0.25 5 + + +0.25 0 0 5 +0.5 0 0 5 + + +0.25 0 0 5 +0.25 0 0 5 + + +0.25 1.5e-17 0.25 5 +0.5 1.5e-17 0.25 5 + + +0.25 1.5e-17 0.25 5 +0.25 -0.25 3.1e-17 5 + + +0.5 0 0 5 +0.5 1.5e-17 0.25 5 + + +0.5 0 0 5 +0.5 0 0 5 + + +0.5 1.5e-17 0.25 5 +0.5 -0.25 3.1e-17 5 + + +0.25 0 0 5 +0.25 -0.25 3.1e-17 5 + + +0.25 0 0 5 +0.5 0 0 5 + + +0.25 -0.25 3.1e-17 5 +0.5 -0.25 3.1e-17 5 + + +0.5 0 0 5 +0.5 -0.25 3.1e-17 5 + + +0.25 -0 0 5 +0.25 -0.25 3.1e-17 5 + + +0.25 -0 0 5 +0.5 -0 0 5 + + +0.25 -0 0 5 +0.25 -0 0 5 + + +0.25 -0.25 3.1e-17 5 +0.5 -0.25 3.1e-17 5 + + +0.25 -0.25 3.1e-17 5 +0.25 -4.6e-17 -0.25 5 + + +0.5 -0 0 5 +0.5 -0.25 3.1e-17 5 + + +0.5 -0 0 5 +0.5 -0 0 5 + + +0.5 -0.25 3.1e-17 5 +0.5 -4.6e-17 -0.25 5 + + +0.25 -0 0 5 +0.25 -4.6e-17 -0.25 5 + + +0.25 -0 0 5 +0.5 -0 0 5 + + +0.25 -4.6e-17 -0.25 5 +0.5 -4.6e-17 -0.25 5 + + +0.5 -0 0 5 +0.5 -4.6e-17 -0.25 5 + + +0.25 0 -0 5 +0.25 -4.6e-17 -0.25 5 + + +0.25 0 -0 5 +0.5 0 -0 5 + + +0.25 0 -0 5 +0.25 0 -0 5 + + +0.25 -4.6e-17 -0.25 5 +0.5 -4.6e-17 -0.25 5 + + +0.25 -4.6e-17 -0.25 5 +0.25 0.25 -6.1e-17 5 + + +0.5 0 -0 5 +0.5 -4.6e-17 -0.25 5 + + +0.5 0 -0 5 +0.5 0 -0 5 + + +0.5 -4.6e-17 -0.25 5 +0.5 0.25 -6.1e-17 5 + + +0.25 0 -0 5 +0.25 0.25 -6.1e-17 5 + + +0.25 0 -0 5 +0.5 0 -0 5 + + +0.25 0.25 -6.1e-17 5 +0.5 0.25 -6.1e-17 5 + + +0.5 0 -0 5 +0.5 0.25 -6.1e-17 5 + + +0.25 0.25 0 6 +0.25 0.5 0 6 + + +0.25 0.25 0 6 +0.5 0.25 0 6 + + +0.25 0.25 0 6 +0.25 1.5e-17 0.25 6 + + +0.25 0.5 0 6 +0.5 0.5 0 6 + + +0.25 0.5 0 6 +0.25 3.1e-17 0.5 6 + + +0.5 0.25 0 6 +0.5 0.5 0 6 + + +0.5 0.25 0 6 +0.5 1.5e-17 0.25 6 + + +0.5 0.5 0 6 +0.5 3.1e-17 0.5 6 + + +0.25 1.5e-17 0.25 6 +0.25 3.1e-17 0.5 6 + + +0.25 1.5e-17 0.25 6 +0.5 1.5e-17 0.25 6 + + +0.25 3.1e-17 0.5 6 +0.5 3.1e-17 0.5 6 + + +0.5 1.5e-17 0.25 6 +0.5 3.1e-17 0.5 6 + + +0.25 1.5e-17 0.25 6 +0.25 3.1e-17 0.5 6 + + +0.25 1.5e-17 0.25 6 +0.5 1.5e-17 0.25 6 + + +0.25 1.5e-17 0.25 6 +0.25 -0.25 3.1e-17 6 + + +0.25 3.1e-17 0.5 6 +0.5 3.1e-17 0.5 6 + + +0.25 3.1e-17 0.5 6 +0.25 -0.5 6.1e-17 6 + + +0.5 1.5e-17 0.25 6 +0.5 3.1e-17 0.5 6 + + +0.5 1.5e-17 0.25 6 +0.5 -0.25 3.1e-17 6 + + +0.5 3.1e-17 0.5 6 +0.5 -0.5 6.1e-17 6 + + +0.25 -0.25 3.1e-17 6 +0.25 -0.5 6.1e-17 6 + + +0.25 -0.25 3.1e-17 6 +0.5 -0.25 3.1e-17 6 + + +0.25 -0.5 6.1e-17 6 +0.5 -0.5 6.1e-17 6 + + +0.5 -0.25 3.1e-17 6 +0.5 -0.5 6.1e-17 6 + + +0.25 -0.25 3.1e-17 6 +0.25 -0.5 6.1e-17 6 + + +0.25 -0.25 3.1e-17 6 +0.5 -0.25 3.1e-17 6 + + +0.25 -0.25 3.1e-17 6 +0.25 -4.6e-17 -0.25 6 + + +0.25 -0.5 6.1e-17 6 +0.5 -0.5 6.1e-17 6 + + +0.25 -0.5 6.1e-17 6 +0.25 -9.2e-17 -0.5 6 + + +0.5 -0.25 3.1e-17 6 +0.5 -0.5 6.1e-17 6 + + +0.5 -0.25 3.1e-17 6 +0.5 -4.6e-17 -0.25 6 + + +0.5 -0.5 6.1e-17 6 +0.5 -9.2e-17 -0.5 6 + + +0.25 -4.6e-17 -0.25 6 +0.25 -9.2e-17 -0.5 6 + + +0.25 -4.6e-17 -0.25 6 +0.5 -4.6e-17 -0.25 6 + + +0.25 -9.2e-17 -0.5 6 +0.5 -9.2e-17 -0.5 6 + + +0.5 -4.6e-17 -0.25 6 +0.5 -9.2e-17 -0.5 6 + + +0.25 -4.6e-17 -0.25 6 +0.25 -9.2e-17 -0.5 6 + + +0.25 -4.6e-17 -0.25 6 +0.5 -4.6e-17 -0.25 6 + + +0.25 -4.6e-17 -0.25 6 +0.25 0.25 -6.1e-17 6 + + +0.25 -9.2e-17 -0.5 6 +0.5 -9.2e-17 -0.5 6 + + +0.25 -9.2e-17 -0.5 6 +0.25 0.5 -1.2e-16 6 + + +0.5 -4.6e-17 -0.25 6 +0.5 -9.2e-17 -0.5 6 + + +0.5 -4.6e-17 -0.25 6 +0.5 0.25 -6.1e-17 6 + + +0.5 -9.2e-17 -0.5 6 +0.5 0.5 -1.2e-16 6 + + +0.25 0.25 -6.1e-17 6 +0.25 0.5 -1.2e-16 6 + + +0.25 0.25 -6.1e-17 6 +0.5 0.25 -6.1e-17 6 + + +0.25 0.5 -1.2e-16 6 +0.5 0.5 -1.2e-16 6 + + +0.5 0.25 -6.1e-17 6 +0.5 0.5 -1.2e-16 6 + +