--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// solves a 2D Poisson equation with FE_DGQ elements and periodic boundary conditions
+
+#include "../tests.h"
+#include <deal.II/base/function.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/fe/fe_dgq.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/grid_tools.h>
+#include <deal.II/integrators/laplace.h>
+#include <deal.II/lac/vector.h>
+#include <deal.II/lac/compressed_sparsity_pattern.h>
+#include <deal.II/lac/solver_cg.h>
+#include <deal.II/lac/precondition.h>
+#include <deal.II/meshworker/dof_info.h>
+#include <deal.II/meshworker/integration_info.h>
+#include <deal.II/meshworker/assembler.h>
+#include <deal.II/meshworker/loop.h>
+#include <deal.II/numerics/data_out.h>
+#include <deal.II/numerics/vector_tools.h>
+
+
+#include <fstream>
+#include <iomanip>
+
+
+template <int dim>
+class MatrixIntegrator : public MeshWorker::LocalIntegrator<dim>
+{
+public:
+ void cell(MeshWorker::DoFInfo<dim> &dinfo,
+ typename MeshWorker::IntegrationInfo<dim> &info) const;
+ void boundary(MeshWorker::DoFInfo<dim> &dinfo,
+ typename MeshWorker::IntegrationInfo<dim> &info) const;
+ void face(MeshWorker::DoFInfo<dim> &dinfo1,
+ MeshWorker::DoFInfo<dim> &dinfo2,
+ typename MeshWorker::IntegrationInfo<dim> &info1,
+ typename MeshWorker::IntegrationInfo<dim> &info2) const;
+};
+
+template <int dim>
+void MatrixIntegrator<dim>
+::cell(MeshWorker::DoFInfo<dim> &dinfo,
+ typename MeshWorker::IntegrationInfo<dim> &info) const
+{
+ LocalIntegrators::Laplace::cell_matrix(dinfo.matrix(0,false).matrix,
+ info.fe_values());
+}
+
+template <int dim>
+void MatrixIntegrator<dim>
+::boundary(MeshWorker::DoFInfo<dim> &dinfo,
+ typename MeshWorker::IntegrationInfo<dim> &info) const
+{
+ const unsigned int deg = info.fe_values(0).get_fe().degree;
+ LocalIntegrators::Laplace
+ ::nitsche_matrix(dinfo.matrix(0,false).matrix, info.fe_values(0),
+ LocalIntegrators::Laplace::
+ compute_penalty(dinfo, dinfo, deg, deg));
+}
+
+template <int dim>
+void MatrixIntegrator<dim>
+::face(MeshWorker::DoFInfo<dim> &dinfo1,
+ MeshWorker::DoFInfo<dim> &dinfo2,
+ typename MeshWorker::IntegrationInfo<dim> &info1,
+ typename MeshWorker::IntegrationInfo<dim> &info2) const
+{
+ const unsigned int deg = info1.fe_values(0).get_fe().degree;
+ LocalIntegrators::Laplace
+ ::ip_matrix(dinfo1.matrix(0,false).matrix, dinfo1.matrix(0,true).matrix,
+ dinfo2.matrix(0,true).matrix, dinfo2.matrix(0,false).matrix,
+ info1.fe_values(0), info2.fe_values(0),
+ LocalIntegrators::Laplace::compute_penalty(dinfo1, dinfo2, deg, deg));
+}
+
+
+template <int dim>
+class Step4
+{
+public:
+ Step4 ();
+ void run ();
+
+private:
+ void make_grid ();
+ void setup_system();
+ void solve ();
+ void output_results (const unsigned int cycle) const;
+ void check_periodicity (const unsigned int cycle) const;
+
+
+ Triangulation<dim> triangulation;
+ FE_DGQ<dim> fe;
+ DoFHandler<dim> dof_handler;
+
+ SparsityPattern sparsity_pattern;
+ SparseMatrix<double> system_matrix;
+
+ Vector<double> solution;
+ Vector<double> system_rhs;
+};
+
+
+
+
+template <int dim>
+Step4<dim>::Step4 ()
+ :
+ fe (1),
+ dof_handler (triangulation)
+{}
+
+
+template <int dim>
+void Step4<dim>::make_grid ()
+{
+ GridGenerator::hyper_cube (triangulation, -1, 1, true);
+ typedef typename dealii::Triangulation<dim>::cell_iterator CellIteratorTria;
+ std::vector<dealii::GridTools::PeriodicFacePair<CellIteratorTria> > periodic_faces;
+ const unsigned int b_id1 = 2;
+ const unsigned int b_id2 = 3;
+ const unsigned int direction = 1;
+
+ dealii::GridTools::collect_periodic_faces (triangulation, b_id1, b_id2,
+ direction, periodic_faces, dealii::Tensor<1,dim>());
+ triangulation.add_periodicity(periodic_faces);
+ triangulation.refine_global (1);
+}
+
+
+
+template <int dim>
+void Step4<dim>::setup_system ()
+{
+ dof_handler.distribute_dofs (fe);
+
+ DynamicSparsityPattern dsp(dof_handler.n_dofs());
+ DoFTools::make_flux_sparsity_pattern (dof_handler, dsp);
+ sparsity_pattern.copy_from(dsp);
+ system_matrix.reinit (sparsity_pattern);
+
+ solution.reinit (dof_handler.n_dofs());
+ system_rhs.reinit (dof_handler.n_dofs());
+
+ MappingQGeneric<dim> mapping(1);
+ MeshWorker::IntegrationInfoBox<dim> info_box;
+ UpdateFlags update_flags = update_values | update_gradients;
+ info_box.add_update_flags_all(update_flags);
+ info_box.initialize(fe, mapping);
+
+ MeshWorker::DoFInfo<dim> dof_info(dof_handler);
+ MeshWorker::Assembler::MatrixSimple<SparseMatrix<double> > assembler;
+ assembler.initialize(system_matrix);
+ MatrixIntegrator<dim> integrator;
+ MeshWorker::integration_loop<dim, dim>(dof_handler.begin_active(),
+ dof_handler.end(),
+ dof_info, info_box,
+ integrator, assembler);
+
+ for (unsigned int i=0; i<system_rhs.size(); ++i)
+ system_rhs(i) = 0.01*i-0.000001*i*i;
+}
+
+
+
+template <int dim>
+void Step4<dim>::solve ()
+{
+ solution = 0;
+ SolverControl solver_control (1000, 1e-10, false, false);
+ SolverCG<> solver (solver_control);
+
+ PreconditionJacobi<SparseMatrix<double> > preconditioner;
+ preconditioner.initialize(system_matrix);
+
+ solver.solve (system_matrix, solution, system_rhs,
+ preconditioner);
+}
+
+
+
+template <int dim>
+void Step4<dim>::output_results (const unsigned int cycle) const
+{
+ std::string filename = "solution-"+dealii::Utilities::int_to_string(cycle,2);
+
+ dealii::DataOut<dim> data_out;
+ data_out.attach_dof_handler (dof_handler);
+ data_out.add_data_vector (solution, "u");
+
+ data_out.build_patches (fe.degree);
+ std::ofstream output ((filename + ".vtk").c_str());
+ data_out.write_vtk (output);
+}
+
+
+
+template <int dim>
+void Step4<dim>::check_periodicity (const unsigned int cycle) const
+{}
+
+template <>
+void Step4<2>::check_periodicity (const unsigned int cycle) const
+{
+ unsigned int n_points = 4;
+ for (unsigned int i = 0; i<cycle; i++)
+ n_points*=2;
+
+ //don't test exactly at the support points, since point_value is not stable there
+ const double eps = 1./(16.*n_points);
+
+ bool all_passed = true;
+
+ for (unsigned int i=1; i< n_points; i++)
+ {
+ Vector<double> value1(1);
+ Vector<double> value2(1);
+
+ Point<2> point1;
+ point1(0)=2*(1.*i/n_points+eps)-1;
+ point1(1)=-1.;
+ Point<2> point2;
+ point2(0)=2*(1.*i/n_points+eps)-1;
+ point2(1)=1.;
+
+ VectorTools::point_value (dof_handler, solution, point1, value1);
+ VectorTools::point_value (dof_handler, solution, point2, value2);
+
+ const double rel_error = std::abs((value2[0]-value1[0])/value1[0]);
+ const double rel_tol = 1./std::pow(2., cycle);
+
+ if (rel_error < rel_tol)
+ deallog << point1 << "\t pass" << std::endl;
+ else
+ {
+ deallog << point1 << "\t fail" << std::endl;
+ deallog << point1 << "\t" << value1[0] << "\t"
+ << value2[0] << "\t"
+ << rel_error << std::endl;
+ all_passed = false;
+ }
+ }
+ AssertThrow(all_passed, ExcInternalError());
+}
+
+
+
+template <int dim>
+void Step4<dim>::run()
+{
+ for (unsigned int cycle = 0; cycle < 4; ++cycle)
+ {
+ if (cycle == 0)
+ make_grid();
+ else
+ triangulation.refine_global(1);
+
+ setup_system();
+ solve();
+ //output_results(cycle);
+ deallog.push(Utilities::int_to_string(dof_handler.n_dofs(),5));
+ check_periodicity(cycle);
+ deallog.pop();
+ }
+}
+
+
+int main (int argc, char **argv)
+{
+ std::ofstream logfile("output");
+ deallog.attach(logfile);
+ deallog.threshold_double(1.e-10);
+
+ try
+ {
+ Step4<2> test;
+ test.run();
+ }
+ catch (std::exception &exc)
+ {
+ deallog << std::endl << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+ deallog << "Exception on processing: " << std::endl
+ << exc.what() << std::endl
+ << "Aborting!" << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+
+ return 1;
+ }
+ catch (...)
+ {
+ deallog << std::endl << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+ deallog << "Unknown exception!" << std::endl
+ << "Aborting!" << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+ return 1;
+ };
+}
--- /dev/null
+
+DEAL:00016::-0.468750 -1.00000 pass
+DEAL:00016::0.0312500 -1.00000 pass
+DEAL:00016::0.531250 -1.00000 pass
+DEAL:00064::-0.734375 -1.00000 pass
+DEAL:00064::-0.484375 -1.00000 pass
+DEAL:00064::-0.234375 -1.00000 pass
+DEAL:00064::0.0156250 -1.00000 pass
+DEAL:00064::0.265625 -1.00000 pass
+DEAL:00064::0.515625 -1.00000 pass
+DEAL:00064::0.765625 -1.00000 pass
+DEAL:00256::-0.867188 -1.00000 pass
+DEAL:00256::-0.742188 -1.00000 pass
+DEAL:00256::-0.617188 -1.00000 pass
+DEAL:00256::-0.492188 -1.00000 pass
+DEAL:00256::-0.367188 -1.00000 pass
+DEAL:00256::-0.242188 -1.00000 pass
+DEAL:00256::-0.117188 -1.00000 pass
+DEAL:00256::0.00781250 -1.00000 pass
+DEAL:00256::0.132812 -1.00000 pass
+DEAL:00256::0.257812 -1.00000 pass
+DEAL:00256::0.382812 -1.00000 pass
+DEAL:00256::0.507812 -1.00000 pass
+DEAL:00256::0.632812 -1.00000 pass
+DEAL:00256::0.757812 -1.00000 pass
+DEAL:00256::0.882812 -1.00000 pass
+DEAL:01024::-0.933594 -1.00000 pass
+DEAL:01024::-0.871094 -1.00000 pass
+DEAL:01024::-0.808594 -1.00000 pass
+DEAL:01024::-0.746094 -1.00000 pass
+DEAL:01024::-0.683594 -1.00000 pass
+DEAL:01024::-0.621094 -1.00000 pass
+DEAL:01024::-0.558594 -1.00000 pass
+DEAL:01024::-0.496094 -1.00000 pass
+DEAL:01024::-0.433594 -1.00000 pass
+DEAL:01024::-0.371094 -1.00000 pass
+DEAL:01024::-0.308594 -1.00000 pass
+DEAL:01024::-0.246094 -1.00000 pass
+DEAL:01024::-0.183594 -1.00000 pass
+DEAL:01024::-0.121094 -1.00000 pass
+DEAL:01024::-0.0585938 -1.00000 pass
+DEAL:01024::0.00390625 -1.00000 pass
+DEAL:01024::0.0664062 -1.00000 pass
+DEAL:01024::0.128906 -1.00000 pass
+DEAL:01024::0.191406 -1.00000 pass
+DEAL:01024::0.253906 -1.00000 pass
+DEAL:01024::0.316406 -1.00000 pass
+DEAL:01024::0.378906 -1.00000 pass
+DEAL:01024::0.441406 -1.00000 pass
+DEAL:01024::0.503906 -1.00000 pass
+DEAL:01024::0.566406 -1.00000 pass
+DEAL:01024::0.628906 -1.00000 pass
+DEAL:01024::0.691406 -1.00000 pass
+DEAL:01024::0.753906 -1.00000 pass
+DEAL:01024::0.816406 -1.00000 pass
+DEAL:01024::0.878906 -1.00000 pass
+DEAL:01024::0.941406 -1.00000 pass
dof_handler.n_dofs(),
dof_handler.n_locally_owned_dofs());
- SolverControl solver_control (dof_handler.n_dofs(), 1e-12);
+ SolverControl solver_control (dof_handler.n_dofs(), 1e-12, false, false);
PETScWrappers::SolverCG solver(solver_control, mpi_communicator);
PETScWrappers::PreconditionJacobi(system_matrix));
#endif
- pcout << " Solved in " << solver_control.last_step()
- << " iterations." << std::endl;
-
constraints.distribute (completely_distributed_solution);
locally_relevant_solution = completely_distributed_solution;
if (Utilities::MPI::this_mpi_process(mpi_communicator)==0)
{
- pcout << point1 << "\t" << get_real_assert_zero_imag(value1[0]) << std::endl;
if (std::abs(value2[0]-value1[0])>1e-8)
{
+ pcout << point1 << "\t" << "fail" << std::endl;
std::cout<<point1<< "\t" << value1[0] << std::endl;
std::cout<<point2<< "\t" << value2[0] << std::endl;
Assert(false, ExcInternalError());
}
+ else
+ {
+ pcout << point1 << "\t" << "pass" << std::endl;
+ }
}
}
}
if (Utilities::MPI::this_mpi_process(mpi_communicator)==0)
{
- pcout << point1 << "\t" << get_real_assert_zero_imag(value1[0]) << std::endl;
if (std::abs(value2[0]-value1[0])>1e-8)
{
+ pcout << point1 << "\t fail check 0" << std::endl;
std::cout<<point1<< "\t" << value1[0] << std::endl;
std::cout<<point2<< "\t" << value2[0] << std::endl;
Assert(false, ExcInternalError());
}
- pcout << point3 << "\t" << get_real_assert_zero_imag(value3[0]) << std::endl;
+ else
+ pcout << point1 << "\t pass check 0" << std::endl;
+
if (std::abs(value4[0]-value3[0])>1e-8)
{
+ pcout << point3 << "\t fail check 1" << std::endl;
std::cout<<point3<< "\t" << value3[0] << std::endl;
std::cout<<point4<< "\t" << value4[0] << std::endl;
Assert(false, ExcInternalError());
}
+ else
+ pcout << point3 << "\t pass check 1" << std::endl;
}
}
}
assemble_system ();
solve ();
//output_results (cycle);
+ deallog.push(Utilities::int_to_string(dof_handler.n_dofs(),5));
check_periodicity(cycle);
+ deallog.pop();
}
}
}
Testing for dim=2
Cycle 0:
-DEAL::Starting value 0.0139546
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.531250 0.00000 -0.00531794
+0.531250 0.00000 pass
Cycle 1:
-DEAL::Starting value 0.0223730
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.265625 0.00000 -0.00177028
-0.515625 0.00000 -0.00542131
-0.765625 0.00000 -0.00148613
+0.265625 0.00000 pass
+0.515625 0.00000 pass
+0.765625 0.00000 pass
Cycle 2:
-DEAL::Starting value 0.0270972
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.132812 0.00000 0.00130870
-0.257812 0.00000 -0.00146441
-0.382812 0.00000 -0.00443592
-0.507812 0.00000 -0.00544040
-0.632812 0.00000 -0.00415112
-0.757812 0.00000 -0.00167411
-0.882812 0.00000 0.000448975
+0.132812 0.00000 pass
+0.257812 0.00000 pass
+0.382812 0.00000 pass
+0.507812 0.00000 pass
+0.632812 0.00000 pass
+0.757812 0.00000 pass
+0.882812 0.00000 pass
Testing for dim=3
Cycle 0:
-DEAL::Starting value 0.00635044
-DEAL::Convergence step 7 value 0
- Solved in 7 iterations.
-0.531250 0.531250 0.00000 0.000730918
-0.531250 0.00000 0.531250 0.000460360
+0.531250 0.531250 0.00000 pass check 0
+0.531250 0.00000 0.531250 pass check 1
Cycle 1:
-DEAL::Starting value 0.0136523
-DEAL::Convergence step 7 value 0
- Solved in 7 iterations.
-0.265625 0.265625 0.00000 -7.01353e-05
-0.265625 0.00000 0.265625 -9.66620e-05
-0.265625 0.515625 0.00000 0.000111798
-0.265625 0.00000 0.515625 4.88063e-05
-0.265625 0.765625 0.00000 -6.75355e-05
-0.265625 0.00000 0.765625 -7.84425e-05
-0.515625 0.265625 0.00000 -0.000111934
-0.515625 0.00000 0.265625 -0.000242694
-0.515625 0.515625 0.00000 0.000800841
-0.515625 0.00000 0.515625 0.000505391
-0.515625 0.765625 0.00000 -0.000133271
-0.515625 0.00000 0.765625 -0.000148360
-0.765625 0.265625 0.00000 -8.91167e-05
-0.765625 0.00000 0.265625 -9.24145e-05
-0.765625 0.515625 0.00000 -7.45770e-05
-0.765625 0.00000 0.515625 -8.11949e-05
-0.765625 0.765625 0.00000 -8.35244e-05
-0.765625 0.00000 0.765625 -8.87256e-05
+0.265625 0.265625 0.00000 pass check 0
+0.265625 0.00000 0.265625 pass check 1
+0.265625 0.515625 0.00000 pass check 0
+0.265625 0.00000 0.515625 pass check 1
+0.265625 0.765625 0.00000 pass check 0
+0.265625 0.00000 0.765625 pass check 1
+0.515625 0.265625 0.00000 pass check 0
+0.515625 0.00000 0.265625 pass check 1
+0.515625 0.515625 0.00000 pass check 0
+0.515625 0.00000 0.515625 pass check 1
+0.515625 0.765625 0.00000 pass check 0
+0.515625 0.00000 0.765625 pass check 1
+0.765625 0.265625 0.00000 pass check 0
+0.765625 0.00000 0.265625 pass check 1
+0.765625 0.515625 0.00000 pass check 0
+0.765625 0.00000 0.515625 pass check 1
+0.765625 0.765625 0.00000 pass check 0
+0.765625 0.00000 0.765625 pass check 1
Cycle 2:
-DEAL::Starting value 0.0250422
-DEAL::Convergence step 7 value 0
- Solved in 7 iterations.
-0.132812 0.132812 0.00000 0.000646598
-0.132812 0.00000 0.132812 0.000661952
-0.132812 0.257812 0.00000 -1.72500e-06
-0.132812 0.00000 0.257812 2.56769e-05
-0.132812 0.382812 0.00000 -0.000451434
-0.132812 0.00000 0.382812 -0.000348091
-0.132812 0.507812 0.00000 -0.000513939
-0.132812 0.00000 0.507812 -0.000370142
-0.132812 0.632812 0.00000 -0.000305068
-0.132812 0.00000 0.632812 -0.000209183
-0.132812 0.757812 0.00000 -2.83548e-05
-0.132812 0.00000 0.757812 -3.09650e-05
-0.132812 0.882812 0.00000 0.000219467
-0.132812 0.00000 0.882812 0.000152889
-0.257812 0.132812 0.00000 -0.000165876
-0.257812 0.00000 0.132812 -0.000173604
-0.257812 0.257812 0.00000 -7.42803e-05
-0.257812 0.00000 0.257812 -9.44600e-05
-0.257812 0.382812 0.00000 1.16046e-05
-0.257812 0.00000 0.382812 -2.27690e-05
-0.257812 0.507812 0.00000 4.38914e-05
-0.257812 0.00000 0.507812 4.34583e-06
-0.257812 0.632812 0.00000 1.39282e-05
-0.257812 0.00000 0.632812 -1.65665e-05
-0.257812 0.757812 0.00000 -5.57917e-05
-0.257812 0.00000 0.757812 -6.88009e-05
-0.257812 0.882812 0.00000 -0.000130653
-0.257812 0.00000 0.882812 -0.000130155
-0.382812 0.132812 0.00000 -0.000995337
-0.382812 0.00000 0.132812 -0.00104505
-0.382812 0.257812 0.00000 -0.000164912
-0.382812 0.00000 0.257812 -0.000244631
-0.382812 0.382812 0.00000 0.000490431
-0.382812 0.00000 0.382812 0.000308958
-0.382812 0.507812 0.00000 0.000638078
-0.382812 0.00000 0.507812 0.000406714
-0.382812 0.632812 0.00000 0.000363883
-0.382812 0.00000 0.632812 0.000200367
-0.382812 0.757812 0.00000 -8.04884e-05
-0.382812 0.00000 0.757812 -0.000108485
-0.382812 0.882812 0.00000 -0.000527674
-0.382812 0.00000 0.882812 -0.000466275
-0.507812 0.132812 0.00000 -0.00125036
-0.507812 0.00000 0.132812 -0.00129986
-0.507812 0.257812 0.00000 -0.000203666
-0.507812 0.00000 0.257812 -0.000301637
-0.507812 0.382812 0.00000 0.000616785
-0.507812 0.00000 0.382812 0.000392221
-0.507812 0.507812 0.00000 0.000802115
-0.507812 0.00000 0.507812 0.000515750
-0.507812 0.632812 0.00000 0.000460227
-0.507812 0.00000 0.632812 0.000258189
-0.507812 0.757812 0.00000 -9.46937e-05
-0.507812 0.00000 0.757812 -0.000128448
-0.507812 0.882812 0.00000 -0.000656598
-0.507812 0.00000 0.882812 -0.000580311
-0.632812 0.132812 0.00000 -0.000846228
-0.632812 0.00000 0.132812 -0.000878904
-0.632812 0.257812 0.00000 -0.000169720
-0.632812 0.00000 0.257812 -0.000237495
-0.632812 0.382812 0.00000 0.000367537
-0.632812 0.00000 0.382812 0.000217482
-0.632812 0.507812 0.00000 0.000494364
-0.632812 0.00000 0.507812 0.000304766
-0.632812 0.632812 0.00000 0.000273516
-0.632812 0.00000 0.632812 0.000138834
-0.632812 0.757812 0.00000 -9.38392e-05
-0.632812 0.00000 0.757812 -0.000119091
-0.632812 0.882812 0.00000 -0.000470155
-0.632812 0.00000 0.882812 -0.000423658
-0.757812 0.132812 0.00000 -0.000128857
-0.757812 0.00000 0.132812 -0.000133514
-0.757812 0.257812 0.00000 -9.39660e-05
-0.757812 0.00000 0.257812 -0.000105571
-0.757812 0.382812 0.00000 -5.43154e-05
-0.757812 0.00000 0.382812 -7.19574e-05
-0.757812 0.507812 0.00000 -3.63902e-05
-0.757812 0.00000 0.507812 -5.47864e-05
-0.757812 0.632812 0.00000 -4.81446e-05
-0.757812 0.00000 0.632812 -6.35886e-05
-0.757812 0.757812 0.00000 -8.16497e-05
-0.757812 0.00000 0.757812 -9.07392e-05
-0.757812 0.882812 0.00000 -0.000121929
-0.757812 0.00000 0.882812 -0.000124171
-0.882812 0.132812 0.00000 0.000457188
-0.882812 0.00000 0.132812 0.000477415
-0.882812 0.257812 0.00000 -2.51848e-05
-0.882812 0.00000 0.257812 1.14568e-05
-0.882812 0.382812 0.00000 -0.000393924
-0.882812 0.00000 0.382812 -0.000301463
-0.882812 0.507812 0.00000 -0.000469132
-0.882812 0.00000 0.507812 -0.000344266
-0.882812 0.632812 0.00000 -0.000311403
-0.882812 0.00000 0.632812 -0.000228583
-0.882812 0.757812 0.00000 -6.69534e-05
-0.882812 0.00000 0.757812 -5.82098e-05
-0.882812 0.882812 0.00000 0.000172488
-0.882812 0.00000 0.882812 0.000134647
+0.132812 0.132812 0.00000 pass check 0
+0.132812 0.00000 0.132812 pass check 1
+0.132812 0.257812 0.00000 pass check 0
+0.132812 0.00000 0.257812 pass check 1
+0.132812 0.382812 0.00000 pass check 0
+0.132812 0.00000 0.382812 pass check 1
+0.132812 0.507812 0.00000 pass check 0
+0.132812 0.00000 0.507812 pass check 1
+0.132812 0.632812 0.00000 pass check 0
+0.132812 0.00000 0.632812 pass check 1
+0.132812 0.757812 0.00000 pass check 0
+0.132812 0.00000 0.757812 pass check 1
+0.132812 0.882812 0.00000 pass check 0
+0.132812 0.00000 0.882812 pass check 1
+0.257812 0.132812 0.00000 pass check 0
+0.257812 0.00000 0.132812 pass check 1
+0.257812 0.257812 0.00000 pass check 0
+0.257812 0.00000 0.257812 pass check 1
+0.257812 0.382812 0.00000 pass check 0
+0.257812 0.00000 0.382812 pass check 1
+0.257812 0.507812 0.00000 pass check 0
+0.257812 0.00000 0.507812 pass check 1
+0.257812 0.632812 0.00000 pass check 0
+0.257812 0.00000 0.632812 pass check 1
+0.257812 0.757812 0.00000 pass check 0
+0.257812 0.00000 0.757812 pass check 1
+0.257812 0.882812 0.00000 pass check 0
+0.257812 0.00000 0.882812 pass check 1
+0.382812 0.132812 0.00000 pass check 0
+0.382812 0.00000 0.132812 pass check 1
+0.382812 0.257812 0.00000 pass check 0
+0.382812 0.00000 0.257812 pass check 1
+0.382812 0.382812 0.00000 pass check 0
+0.382812 0.00000 0.382812 pass check 1
+0.382812 0.507812 0.00000 pass check 0
+0.382812 0.00000 0.507812 pass check 1
+0.382812 0.632812 0.00000 pass check 0
+0.382812 0.00000 0.632812 pass check 1
+0.382812 0.757812 0.00000 pass check 0
+0.382812 0.00000 0.757812 pass check 1
+0.382812 0.882812 0.00000 pass check 0
+0.382812 0.00000 0.882812 pass check 1
+0.507812 0.132812 0.00000 pass check 0
+0.507812 0.00000 0.132812 pass check 1
+0.507812 0.257812 0.00000 pass check 0
+0.507812 0.00000 0.257812 pass check 1
+0.507812 0.382812 0.00000 pass check 0
+0.507812 0.00000 0.382812 pass check 1
+0.507812 0.507812 0.00000 pass check 0
+0.507812 0.00000 0.507812 pass check 1
+0.507812 0.632812 0.00000 pass check 0
+0.507812 0.00000 0.632812 pass check 1
+0.507812 0.757812 0.00000 pass check 0
+0.507812 0.00000 0.757812 pass check 1
+0.507812 0.882812 0.00000 pass check 0
+0.507812 0.00000 0.882812 pass check 1
+0.632812 0.132812 0.00000 pass check 0
+0.632812 0.00000 0.132812 pass check 1
+0.632812 0.257812 0.00000 pass check 0
+0.632812 0.00000 0.257812 pass check 1
+0.632812 0.382812 0.00000 pass check 0
+0.632812 0.00000 0.382812 pass check 1
+0.632812 0.507812 0.00000 pass check 0
+0.632812 0.00000 0.507812 pass check 1
+0.632812 0.632812 0.00000 pass check 0
+0.632812 0.00000 0.632812 pass check 1
+0.632812 0.757812 0.00000 pass check 0
+0.632812 0.00000 0.757812 pass check 1
+0.632812 0.882812 0.00000 pass check 0
+0.632812 0.00000 0.882812 pass check 1
+0.757812 0.132812 0.00000 pass check 0
+0.757812 0.00000 0.132812 pass check 1
+0.757812 0.257812 0.00000 pass check 0
+0.757812 0.00000 0.257812 pass check 1
+0.757812 0.382812 0.00000 pass check 0
+0.757812 0.00000 0.382812 pass check 1
+0.757812 0.507812 0.00000 pass check 0
+0.757812 0.00000 0.507812 pass check 1
+0.757812 0.632812 0.00000 pass check 0
+0.757812 0.00000 0.632812 pass check 1
+0.757812 0.757812 0.00000 pass check 0
+0.757812 0.00000 0.757812 pass check 1
+0.757812 0.882812 0.00000 pass check 0
+0.757812 0.00000 0.882812 pass check 1
+0.882812 0.132812 0.00000 pass check 0
+0.882812 0.00000 0.132812 pass check 1
+0.882812 0.257812 0.00000 pass check 0
+0.882812 0.00000 0.257812 pass check 1
+0.882812 0.382812 0.00000 pass check 0
+0.882812 0.00000 0.382812 pass check 1
+0.882812 0.507812 0.00000 pass check 0
+0.882812 0.00000 0.507812 pass check 1
+0.882812 0.632812 0.00000 pass check 0
+0.882812 0.00000 0.632812 pass check 1
+0.882812 0.757812 0.00000 pass check 0
+0.882812 0.00000 0.757812 pass check 1
+0.882812 0.882812 0.00000 pass check 0
+0.882812 0.00000 0.882812 pass check 1
Testing for dim=2
Cycle 0:
-DEAL::Starting value 0.0140170
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.531250 0.00000 -0.00531794
+0.531250 0.00000 pass
Cycle 1:
-DEAL::Starting value 0.0219353
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.265625 0.00000 -0.00177028
-0.515625 0.00000 -0.00542131
-0.765625 0.00000 -0.00148613
+0.265625 0.00000 pass
+0.515625 0.00000 pass
+0.765625 0.00000 pass
Cycle 2:
-DEAL::Starting value 0.0268844
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.132812 0.00000 0.00130870
-0.257812 0.00000 -0.00146441
-0.382812 0.00000 -0.00443592
-0.507812 0.00000 -0.00544040
-0.632812 0.00000 -0.00415112
-0.757812 0.00000 -0.00167411
-0.882812 0.00000 0.000448975
+0.132812 0.00000 pass
+0.257812 0.00000 pass
+0.382812 0.00000 pass
+0.507812 0.00000 pass
+0.632812 0.00000 pass
+0.757812 0.00000 pass
+0.882812 0.00000 pass
Testing for dim=3
Cycle 0:
-DEAL::Starting value 0.00632368
-DEAL::Convergence step 7 value 0
- Solved in 7 iterations.
-0.531250 0.531250 0.00000 0.000730918
-0.531250 0.00000 0.531250 0.000460360
+0.531250 0.531250 0.00000 pass check 0
+0.531250 0.00000 0.531250 pass check 1
Cycle 1:
-DEAL::Starting value 0.0137083
-DEAL::Convergence step 7 value 0
- Solved in 7 iterations.
-0.265625 0.265625 0.00000 -7.01353e-05
-0.265625 0.00000 0.265625 -9.66620e-05
-0.265625 0.515625 0.00000 0.000111798
-0.265625 0.00000 0.515625 4.88063e-05
-0.265625 0.765625 0.00000 -6.75355e-05
-0.265625 0.00000 0.765625 -7.84425e-05
-0.515625 0.265625 0.00000 -0.000111934
-0.515625 0.00000 0.265625 -0.000242694
-0.515625 0.515625 0.00000 0.000800841
-0.515625 0.00000 0.515625 0.000505391
-0.515625 0.765625 0.00000 -0.000133271
-0.515625 0.00000 0.765625 -0.000148360
-0.765625 0.265625 0.00000 -8.91167e-05
-0.765625 0.00000 0.265625 -9.24145e-05
-0.765625 0.515625 0.00000 -7.45770e-05
-0.765625 0.00000 0.515625 -8.11949e-05
-0.765625 0.765625 0.00000 -8.35244e-05
-0.765625 0.00000 0.765625 -8.87256e-05
+0.265625 0.265625 0.00000 pass check 0
+0.265625 0.00000 0.265625 pass check 1
+0.265625 0.515625 0.00000 pass check 0
+0.265625 0.00000 0.515625 pass check 1
+0.265625 0.765625 0.00000 pass check 0
+0.265625 0.00000 0.765625 pass check 1
+0.515625 0.265625 0.00000 pass check 0
+0.515625 0.00000 0.265625 pass check 1
+0.515625 0.515625 0.00000 pass check 0
+0.515625 0.00000 0.515625 pass check 1
+0.515625 0.765625 0.00000 pass check 0
+0.515625 0.00000 0.765625 pass check 1
+0.765625 0.265625 0.00000 pass check 0
+0.765625 0.00000 0.265625 pass check 1
+0.765625 0.515625 0.00000 pass check 0
+0.765625 0.00000 0.515625 pass check 1
+0.765625 0.765625 0.00000 pass check 0
+0.765625 0.00000 0.765625 pass check 1
Cycle 2:
-DEAL::Starting value 0.0248658
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.132812 0.132812 0.00000 0.000646598
-0.132812 0.00000 0.132812 0.000661952
-0.132812 0.257812 0.00000 -1.72500e-06
-0.132812 0.00000 0.257812 2.56769e-05
-0.132812 0.382812 0.00000 -0.000451434
-0.132812 0.00000 0.382812 -0.000348091
-0.132812 0.507812 0.00000 -0.000513939
-0.132812 0.00000 0.507812 -0.000370142
-0.132812 0.632812 0.00000 -0.000305068
-0.132812 0.00000 0.632812 -0.000209183
-0.132812 0.757812 0.00000 -2.83548e-05
-0.132812 0.00000 0.757812 -3.09650e-05
-0.132812 0.882812 0.00000 0.000219467
-0.132812 0.00000 0.882812 0.000152889
-0.257812 0.132812 0.00000 -0.000165876
-0.257812 0.00000 0.132812 -0.000173604
-0.257812 0.257812 0.00000 -7.42803e-05
-0.257812 0.00000 0.257812 -9.44600e-05
-0.257812 0.382812 0.00000 1.16046e-05
-0.257812 0.00000 0.382812 -2.27690e-05
-0.257812 0.507812 0.00000 4.38914e-05
-0.257812 0.00000 0.507812 4.34583e-06
-0.257812 0.632812 0.00000 1.39282e-05
-0.257812 0.00000 0.632812 -1.65665e-05
-0.257812 0.757812 0.00000 -5.57917e-05
-0.257812 0.00000 0.757812 -6.88009e-05
-0.257812 0.882812 0.00000 -0.000130653
-0.257812 0.00000 0.882812 -0.000130155
-0.382812 0.132812 0.00000 -0.000995337
-0.382812 0.00000 0.132812 -0.00104505
-0.382812 0.257812 0.00000 -0.000164912
-0.382812 0.00000 0.257812 -0.000244631
-0.382812 0.382812 0.00000 0.000490431
-0.382812 0.00000 0.382812 0.000308958
-0.382812 0.507812 0.00000 0.000638078
-0.382812 0.00000 0.507812 0.000406714
-0.382812 0.632812 0.00000 0.000363883
-0.382812 0.00000 0.632812 0.000200367
-0.382812 0.757812 0.00000 -8.04884e-05
-0.382812 0.00000 0.757812 -0.000108485
-0.382812 0.882812 0.00000 -0.000527674
-0.382812 0.00000 0.882812 -0.000466275
-0.507812 0.132812 0.00000 -0.00125036
-0.507812 0.00000 0.132812 -0.00129986
-0.507812 0.257812 0.00000 -0.000203666
-0.507812 0.00000 0.257812 -0.000301637
-0.507812 0.382812 0.00000 0.000616785
-0.507812 0.00000 0.382812 0.000392221
-0.507812 0.507812 0.00000 0.000802115
-0.507812 0.00000 0.507812 0.000515750
-0.507812 0.632812 0.00000 0.000460227
-0.507812 0.00000 0.632812 0.000258189
-0.507812 0.757812 0.00000 -9.46937e-05
-0.507812 0.00000 0.757812 -0.000128448
-0.507812 0.882812 0.00000 -0.000656598
-0.507812 0.00000 0.882812 -0.000580311
-0.632812 0.132812 0.00000 -0.000846228
-0.632812 0.00000 0.132812 -0.000878904
-0.632812 0.257812 0.00000 -0.000169720
-0.632812 0.00000 0.257812 -0.000237495
-0.632812 0.382812 0.00000 0.000367537
-0.632812 0.00000 0.382812 0.000217482
-0.632812 0.507812 0.00000 0.000494364
-0.632812 0.00000 0.507812 0.000304766
-0.632812 0.632812 0.00000 0.000273516
-0.632812 0.00000 0.632812 0.000138834
-0.632812 0.757812 0.00000 -9.38392e-05
-0.632812 0.00000 0.757812 -0.000119091
-0.632812 0.882812 0.00000 -0.000470155
-0.632812 0.00000 0.882812 -0.000423658
-0.757812 0.132812 0.00000 -0.000128857
-0.757812 0.00000 0.132812 -0.000133514
-0.757812 0.257812 0.00000 -9.39660e-05
-0.757812 0.00000 0.257812 -0.000105571
-0.757812 0.382812 0.00000 -5.43154e-05
-0.757812 0.00000 0.382812 -7.19574e-05
-0.757812 0.507812 0.00000 -3.63902e-05
-0.757812 0.00000 0.507812 -5.47864e-05
-0.757812 0.632812 0.00000 -4.81446e-05
-0.757812 0.00000 0.632812 -6.35886e-05
-0.757812 0.757812 0.00000 -8.16497e-05
-0.757812 0.00000 0.757812 -9.07392e-05
-0.757812 0.882812 0.00000 -0.000121929
-0.757812 0.00000 0.882812 -0.000124171
-0.882812 0.132812 0.00000 0.000457188
-0.882812 0.00000 0.132812 0.000477415
-0.882812 0.257812 0.00000 -2.51848e-05
-0.882812 0.00000 0.257812 1.14568e-05
-0.882812 0.382812 0.00000 -0.000393924
-0.882812 0.00000 0.382812 -0.000301463
-0.882812 0.507812 0.00000 -0.000469132
-0.882812 0.00000 0.507812 -0.000344266
-0.882812 0.632812 0.00000 -0.000311403
-0.882812 0.00000 0.632812 -0.000228583
-0.882812 0.757812 0.00000 -6.69534e-05
-0.882812 0.00000 0.757812 -5.82098e-05
-0.882812 0.882812 0.00000 0.000172488
-0.882812 0.00000 0.882812 0.000134647
+0.132812 0.132812 0.00000 pass check 0
+0.132812 0.00000 0.132812 pass check 1
+0.132812 0.257812 0.00000 pass check 0
+0.132812 0.00000 0.257812 pass check 1
+0.132812 0.382812 0.00000 pass check 0
+0.132812 0.00000 0.382812 pass check 1
+0.132812 0.507812 0.00000 pass check 0
+0.132812 0.00000 0.507812 pass check 1
+0.132812 0.632812 0.00000 pass check 0
+0.132812 0.00000 0.632812 pass check 1
+0.132812 0.757812 0.00000 pass check 0
+0.132812 0.00000 0.757812 pass check 1
+0.132812 0.882812 0.00000 pass check 0
+0.132812 0.00000 0.882812 pass check 1
+0.257812 0.132812 0.00000 pass check 0
+0.257812 0.00000 0.132812 pass check 1
+0.257812 0.257812 0.00000 pass check 0
+0.257812 0.00000 0.257812 pass check 1
+0.257812 0.382812 0.00000 pass check 0
+0.257812 0.00000 0.382812 pass check 1
+0.257812 0.507812 0.00000 pass check 0
+0.257812 0.00000 0.507812 pass check 1
+0.257812 0.632812 0.00000 pass check 0
+0.257812 0.00000 0.632812 pass check 1
+0.257812 0.757812 0.00000 pass check 0
+0.257812 0.00000 0.757812 pass check 1
+0.257812 0.882812 0.00000 pass check 0
+0.257812 0.00000 0.882812 pass check 1
+0.382812 0.132812 0.00000 pass check 0
+0.382812 0.00000 0.132812 pass check 1
+0.382812 0.257812 0.00000 pass check 0
+0.382812 0.00000 0.257812 pass check 1
+0.382812 0.382812 0.00000 pass check 0
+0.382812 0.00000 0.382812 pass check 1
+0.382812 0.507812 0.00000 pass check 0
+0.382812 0.00000 0.507812 pass check 1
+0.382812 0.632812 0.00000 pass check 0
+0.382812 0.00000 0.632812 pass check 1
+0.382812 0.757812 0.00000 pass check 0
+0.382812 0.00000 0.757812 pass check 1
+0.382812 0.882812 0.00000 pass check 0
+0.382812 0.00000 0.882812 pass check 1
+0.507812 0.132812 0.00000 pass check 0
+0.507812 0.00000 0.132812 pass check 1
+0.507812 0.257812 0.00000 pass check 0
+0.507812 0.00000 0.257812 pass check 1
+0.507812 0.382812 0.00000 pass check 0
+0.507812 0.00000 0.382812 pass check 1
+0.507812 0.507812 0.00000 pass check 0
+0.507812 0.00000 0.507812 pass check 1
+0.507812 0.632812 0.00000 pass check 0
+0.507812 0.00000 0.632812 pass check 1
+0.507812 0.757812 0.00000 pass check 0
+0.507812 0.00000 0.757812 pass check 1
+0.507812 0.882812 0.00000 pass check 0
+0.507812 0.00000 0.882812 pass check 1
+0.632812 0.132812 0.00000 pass check 0
+0.632812 0.00000 0.132812 pass check 1
+0.632812 0.257812 0.00000 pass check 0
+0.632812 0.00000 0.257812 pass check 1
+0.632812 0.382812 0.00000 pass check 0
+0.632812 0.00000 0.382812 pass check 1
+0.632812 0.507812 0.00000 pass check 0
+0.632812 0.00000 0.507812 pass check 1
+0.632812 0.632812 0.00000 pass check 0
+0.632812 0.00000 0.632812 pass check 1
+0.632812 0.757812 0.00000 pass check 0
+0.632812 0.00000 0.757812 pass check 1
+0.632812 0.882812 0.00000 pass check 0
+0.632812 0.00000 0.882812 pass check 1
+0.757812 0.132812 0.00000 pass check 0
+0.757812 0.00000 0.132812 pass check 1
+0.757812 0.257812 0.00000 pass check 0
+0.757812 0.00000 0.257812 pass check 1
+0.757812 0.382812 0.00000 pass check 0
+0.757812 0.00000 0.382812 pass check 1
+0.757812 0.507812 0.00000 pass check 0
+0.757812 0.00000 0.507812 pass check 1
+0.757812 0.632812 0.00000 pass check 0
+0.757812 0.00000 0.632812 pass check 1
+0.757812 0.757812 0.00000 pass check 0
+0.757812 0.00000 0.757812 pass check 1
+0.757812 0.882812 0.00000 pass check 0
+0.757812 0.00000 0.882812 pass check 1
+0.882812 0.132812 0.00000 pass check 0
+0.882812 0.00000 0.132812 pass check 1
+0.882812 0.257812 0.00000 pass check 0
+0.882812 0.00000 0.257812 pass check 1
+0.882812 0.382812 0.00000 pass check 0
+0.882812 0.00000 0.382812 pass check 1
+0.882812 0.507812 0.00000 pass check 0
+0.882812 0.00000 0.507812 pass check 1
+0.882812 0.632812 0.00000 pass check 0
+0.882812 0.00000 0.632812 pass check 1
+0.882812 0.757812 0.00000 pass check 0
+0.882812 0.00000 0.757812 pass check 1
+0.882812 0.882812 0.00000 pass check 0
+0.882812 0.00000 0.882812 pass check 1
Testing for dim=2
Cycle 0:
-DEAL::Starting value 0.0139631
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.531250 0.00000 -0.00531794
+0.531250 0.00000 pass
Cycle 1:
-DEAL::Starting value 0.0222274
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.265625 0.00000 -0.00177028
-0.515625 0.00000 -0.00542131
-0.765625 0.00000 -0.00148613
+0.265625 0.00000 pass
+0.515625 0.00000 pass
+0.765625 0.00000 pass
Cycle 2:
-DEAL::Starting value 0.0268966
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.132812 0.00000 0.00130870
-0.257812 0.00000 -0.00146441
-0.382812 0.00000 -0.00443592
-0.507812 0.00000 -0.00544040
-0.632812 0.00000 -0.00415112
-0.757812 0.00000 -0.00167411
-0.882812 0.00000 0.000448975
+0.132812 0.00000 pass
+0.257812 0.00000 pass
+0.382812 0.00000 pass
+0.507812 0.00000 pass
+0.632812 0.00000 pass
+0.757812 0.00000 pass
+0.882812 0.00000 pass
Testing for dim=3
Cycle 0:
-DEAL::Starting value 0.00630053
-DEAL::Convergence step 7 value 0
- Solved in 7 iterations.
-0.531250 0.531250 0.00000 0.000730918
-0.531250 0.00000 0.531250 0.000460360
+0.531250 0.531250 0.00000 pass check 0
+0.531250 0.00000 0.531250 pass check 1
Cycle 1:
-DEAL::Starting value 0.0135106
-DEAL::Convergence step 7 value 0
- Solved in 7 iterations.
-0.265625 0.265625 0.00000 -7.01353e-05
-0.265625 0.00000 0.265625 -9.66620e-05
-0.265625 0.515625 0.00000 0.000111798
-0.265625 0.00000 0.515625 4.88063e-05
-0.265625 0.765625 0.00000 -6.75355e-05
-0.265625 0.00000 0.765625 -7.84425e-05
-0.515625 0.265625 0.00000 -0.000111934
-0.515625 0.00000 0.265625 -0.000242694
-0.515625 0.515625 0.00000 0.000800841
-0.515625 0.00000 0.515625 0.000505391
-0.515625 0.765625 0.00000 -0.000133271
-0.515625 0.00000 0.765625 -0.000148360
-0.765625 0.265625 0.00000 -8.91167e-05
-0.765625 0.00000 0.265625 -9.24145e-05
-0.765625 0.515625 0.00000 -7.45770e-05
-0.765625 0.00000 0.515625 -8.11949e-05
-0.765625 0.765625 0.00000 -8.35244e-05
-0.765625 0.00000 0.765625 -8.87256e-05
+0.265625 0.265625 0.00000 pass check 0
+0.265625 0.00000 0.265625 pass check 1
+0.265625 0.515625 0.00000 pass check 0
+0.265625 0.00000 0.515625 pass check 1
+0.265625 0.765625 0.00000 pass check 0
+0.265625 0.00000 0.765625 pass check 1
+0.515625 0.265625 0.00000 pass check 0
+0.515625 0.00000 0.265625 pass check 1
+0.515625 0.515625 0.00000 pass check 0
+0.515625 0.00000 0.515625 pass check 1
+0.515625 0.765625 0.00000 pass check 0
+0.515625 0.00000 0.765625 pass check 1
+0.765625 0.265625 0.00000 pass check 0
+0.765625 0.00000 0.265625 pass check 1
+0.765625 0.515625 0.00000 pass check 0
+0.765625 0.00000 0.515625 pass check 1
+0.765625 0.765625 0.00000 pass check 0
+0.765625 0.00000 0.765625 pass check 1
Cycle 2:
-DEAL::Starting value 0.0248824
-DEAL::Convergence step 8 value 0
- Solved in 8 iterations.
-0.132812 0.132812 0.00000 0.000646598
-0.132812 0.00000 0.132812 0.000661952
-0.132812 0.257812 0.00000 -1.72500e-06
-0.132812 0.00000 0.257812 2.56769e-05
-0.132812 0.382812 0.00000 -0.000451434
-0.132812 0.00000 0.382812 -0.000348091
-0.132812 0.507812 0.00000 -0.000513939
-0.132812 0.00000 0.507812 -0.000370142
-0.132812 0.632812 0.00000 -0.000305068
-0.132812 0.00000 0.632812 -0.000209183
-0.132812 0.757812 0.00000 -2.83548e-05
-0.132812 0.00000 0.757812 -3.09650e-05
-0.132812 0.882812 0.00000 0.000219467
-0.132812 0.00000 0.882812 0.000152889
-0.257812 0.132812 0.00000 -0.000165876
-0.257812 0.00000 0.132812 -0.000173604
-0.257812 0.257812 0.00000 -7.42803e-05
-0.257812 0.00000 0.257812 -9.44600e-05
-0.257812 0.382812 0.00000 1.16046e-05
-0.257812 0.00000 0.382812 -2.27690e-05
-0.257812 0.507812 0.00000 4.38914e-05
-0.257812 0.00000 0.507812 4.34583e-06
-0.257812 0.632812 0.00000 1.39282e-05
-0.257812 0.00000 0.632812 -1.65665e-05
-0.257812 0.757812 0.00000 -5.57917e-05
-0.257812 0.00000 0.757812 -6.88009e-05
-0.257812 0.882812 0.00000 -0.000130653
-0.257812 0.00000 0.882812 -0.000130155
-0.382812 0.132812 0.00000 -0.000995337
-0.382812 0.00000 0.132812 -0.00104505
-0.382812 0.257812 0.00000 -0.000164912
-0.382812 0.00000 0.257812 -0.000244631
-0.382812 0.382812 0.00000 0.000490431
-0.382812 0.00000 0.382812 0.000308958
-0.382812 0.507812 0.00000 0.000638078
-0.382812 0.00000 0.507812 0.000406714
-0.382812 0.632812 0.00000 0.000363883
-0.382812 0.00000 0.632812 0.000200367
-0.382812 0.757812 0.00000 -8.04884e-05
-0.382812 0.00000 0.757812 -0.000108485
-0.382812 0.882812 0.00000 -0.000527674
-0.382812 0.00000 0.882812 -0.000466275
-0.507812 0.132812 0.00000 -0.00125036
-0.507812 0.00000 0.132812 -0.00129986
-0.507812 0.257812 0.00000 -0.000203666
-0.507812 0.00000 0.257812 -0.000301637
-0.507812 0.382812 0.00000 0.000616785
-0.507812 0.00000 0.382812 0.000392221
-0.507812 0.507812 0.00000 0.000802115
-0.507812 0.00000 0.507812 0.000515750
-0.507812 0.632812 0.00000 0.000460227
-0.507812 0.00000 0.632812 0.000258189
-0.507812 0.757812 0.00000 -9.46937e-05
-0.507812 0.00000 0.757812 -0.000128448
-0.507812 0.882812 0.00000 -0.000656598
-0.507812 0.00000 0.882812 -0.000580311
-0.632812 0.132812 0.00000 -0.000846228
-0.632812 0.00000 0.132812 -0.000878904
-0.632812 0.257812 0.00000 -0.000169720
-0.632812 0.00000 0.257812 -0.000237495
-0.632812 0.382812 0.00000 0.000367537
-0.632812 0.00000 0.382812 0.000217482
-0.632812 0.507812 0.00000 0.000494364
-0.632812 0.00000 0.507812 0.000304766
-0.632812 0.632812 0.00000 0.000273516
-0.632812 0.00000 0.632812 0.000138834
-0.632812 0.757812 0.00000 -9.38392e-05
-0.632812 0.00000 0.757812 -0.000119091
-0.632812 0.882812 0.00000 -0.000470155
-0.632812 0.00000 0.882812 -0.000423658
-0.757812 0.132812 0.00000 -0.000128857
-0.757812 0.00000 0.132812 -0.000133514
-0.757812 0.257812 0.00000 -9.39660e-05
-0.757812 0.00000 0.257812 -0.000105571
-0.757812 0.382812 0.00000 -5.43154e-05
-0.757812 0.00000 0.382812 -7.19574e-05
-0.757812 0.507812 0.00000 -3.63902e-05
-0.757812 0.00000 0.507812 -5.47864e-05
-0.757812 0.632812 0.00000 -4.81446e-05
-0.757812 0.00000 0.632812 -6.35886e-05
-0.757812 0.757812 0.00000 -8.16497e-05
-0.757812 0.00000 0.757812 -9.07392e-05
-0.757812 0.882812 0.00000 -0.000121929
-0.757812 0.00000 0.882812 -0.000124171
-0.882812 0.132812 0.00000 0.000457188
-0.882812 0.00000 0.132812 0.000477415
-0.882812 0.257812 0.00000 -2.51848e-05
-0.882812 0.00000 0.257812 1.14568e-05
-0.882812 0.382812 0.00000 -0.000393924
-0.882812 0.00000 0.382812 -0.000301463
-0.882812 0.507812 0.00000 -0.000469132
-0.882812 0.00000 0.507812 -0.000344266
-0.882812 0.632812 0.00000 -0.000311403
-0.882812 0.00000 0.632812 -0.000228583
-0.882812 0.757812 0.00000 -6.69534e-05
-0.882812 0.00000 0.757812 -5.82098e-05
-0.882812 0.882812 0.00000 0.000172488
-0.882812 0.00000 0.882812 0.000134647
+0.132812 0.132812 0.00000 pass check 0
+0.132812 0.00000 0.132812 pass check 1
+0.132812 0.257812 0.00000 pass check 0
+0.132812 0.00000 0.257812 pass check 1
+0.132812 0.382812 0.00000 pass check 0
+0.132812 0.00000 0.382812 pass check 1
+0.132812 0.507812 0.00000 pass check 0
+0.132812 0.00000 0.507812 pass check 1
+0.132812 0.632812 0.00000 pass check 0
+0.132812 0.00000 0.632812 pass check 1
+0.132812 0.757812 0.00000 pass check 0
+0.132812 0.00000 0.757812 pass check 1
+0.132812 0.882812 0.00000 pass check 0
+0.132812 0.00000 0.882812 pass check 1
+0.257812 0.132812 0.00000 pass check 0
+0.257812 0.00000 0.132812 pass check 1
+0.257812 0.257812 0.00000 pass check 0
+0.257812 0.00000 0.257812 pass check 1
+0.257812 0.382812 0.00000 pass check 0
+0.257812 0.00000 0.382812 pass check 1
+0.257812 0.507812 0.00000 pass check 0
+0.257812 0.00000 0.507812 pass check 1
+0.257812 0.632812 0.00000 pass check 0
+0.257812 0.00000 0.632812 pass check 1
+0.257812 0.757812 0.00000 pass check 0
+0.257812 0.00000 0.757812 pass check 1
+0.257812 0.882812 0.00000 pass check 0
+0.257812 0.00000 0.882812 pass check 1
+0.382812 0.132812 0.00000 pass check 0
+0.382812 0.00000 0.132812 pass check 1
+0.382812 0.257812 0.00000 pass check 0
+0.382812 0.00000 0.257812 pass check 1
+0.382812 0.382812 0.00000 pass check 0
+0.382812 0.00000 0.382812 pass check 1
+0.382812 0.507812 0.00000 pass check 0
+0.382812 0.00000 0.507812 pass check 1
+0.382812 0.632812 0.00000 pass check 0
+0.382812 0.00000 0.632812 pass check 1
+0.382812 0.757812 0.00000 pass check 0
+0.382812 0.00000 0.757812 pass check 1
+0.382812 0.882812 0.00000 pass check 0
+0.382812 0.00000 0.882812 pass check 1
+0.507812 0.132812 0.00000 pass check 0
+0.507812 0.00000 0.132812 pass check 1
+0.507812 0.257812 0.00000 pass check 0
+0.507812 0.00000 0.257812 pass check 1
+0.507812 0.382812 0.00000 pass check 0
+0.507812 0.00000 0.382812 pass check 1
+0.507812 0.507812 0.00000 pass check 0
+0.507812 0.00000 0.507812 pass check 1
+0.507812 0.632812 0.00000 pass check 0
+0.507812 0.00000 0.632812 pass check 1
+0.507812 0.757812 0.00000 pass check 0
+0.507812 0.00000 0.757812 pass check 1
+0.507812 0.882812 0.00000 pass check 0
+0.507812 0.00000 0.882812 pass check 1
+0.632812 0.132812 0.00000 pass check 0
+0.632812 0.00000 0.132812 pass check 1
+0.632812 0.257812 0.00000 pass check 0
+0.632812 0.00000 0.257812 pass check 1
+0.632812 0.382812 0.00000 pass check 0
+0.632812 0.00000 0.382812 pass check 1
+0.632812 0.507812 0.00000 pass check 0
+0.632812 0.00000 0.507812 pass check 1
+0.632812 0.632812 0.00000 pass check 0
+0.632812 0.00000 0.632812 pass check 1
+0.632812 0.757812 0.00000 pass check 0
+0.632812 0.00000 0.757812 pass check 1
+0.632812 0.882812 0.00000 pass check 0
+0.632812 0.00000 0.882812 pass check 1
+0.757812 0.132812 0.00000 pass check 0
+0.757812 0.00000 0.132812 pass check 1
+0.757812 0.257812 0.00000 pass check 0
+0.757812 0.00000 0.257812 pass check 1
+0.757812 0.382812 0.00000 pass check 0
+0.757812 0.00000 0.382812 pass check 1
+0.757812 0.507812 0.00000 pass check 0
+0.757812 0.00000 0.507812 pass check 1
+0.757812 0.632812 0.00000 pass check 0
+0.757812 0.00000 0.632812 pass check 1
+0.757812 0.757812 0.00000 pass check 0
+0.757812 0.00000 0.757812 pass check 1
+0.757812 0.882812 0.00000 pass check 0
+0.757812 0.00000 0.882812 pass check 1
+0.882812 0.132812 0.00000 pass check 0
+0.882812 0.00000 0.132812 pass check 1
+0.882812 0.257812 0.00000 pass check 0
+0.882812 0.00000 0.257812 pass check 1
+0.882812 0.382812 0.00000 pass check 0
+0.882812 0.00000 0.382812 pass check 1
+0.882812 0.507812 0.00000 pass check 0
+0.882812 0.00000 0.507812 pass check 1
+0.882812 0.632812 0.00000 pass check 0
+0.882812 0.00000 0.632812 pass check 1
+0.882812 0.757812 0.00000 pass check 0
+0.882812 0.00000 0.757812 pass check 1
+0.882812 0.882812 0.00000 pass check 0
+0.882812 0.00000 0.882812 pass check 1