--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2005 - 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// test hp_constraints when neither_element_dominates.
+// 2D and 3D hp-refinement.
+// Check continuity across the face by evaluating the
+// field from both sides.
+
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/utilities.h>
+#include <deal.II/base/index_set.h>
+#include <deal.II/base/conditional_ostream.h>
+#include <deal.II/base/function.h>
+#include <deal.II/base/function_parser.h>
+
+#include <deal.II/numerics/data_postprocessor.h>
+#include <deal.II/numerics/data_out.h>
+
+#include <deal.II/hp/dof_handler.h>
+#include <deal.II/hp/fe_values.h>
+#include <deal.II/hp/q_collection.h>
+#include <deal.II/hp/fe_collection.h>
+
+
+#include <deal.II/dofs/dof_renumbering.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_nothing.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/fe/fe_tools.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/lac/sparsity_tools.h>
+#include <deal.II/lac/compressed_simple_sparsity_pattern.h>
+#include <deal.II/hp/fe_values.h>
+
+#include <deal.II/lac/petsc_parallel_vector.h>
+#include <deal.II/lac/petsc_parallel_sparse_matrix.h>
+#include <deal.II/lac/petsc_solver.h>
+#include <deal.II/lac/petsc_precondition.h>
+
+#include <deal.II/lac/slepc_solver.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <deal.II/grid/grid_refinement.h>
+#include <deal.II/numerics/error_estimator.h>
+
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_in.h>
+#include <deal.II/grid/grid_tools.h>
+#include <deal.II/lac/vector.h>
+
+#include "../tests.h"
+
+#include <fstream>
+#include <iostream>
+
+//#define DEBUG_OUTPUT_VTK
+unsigned int counter=0;
+
+const double eps = 1e-10;
+
+using namespace dealii;
+
+template <int dim>
+struct less_than_key
+{
+ inline bool operator() (const std::pair<Point<dim>,Vector<double>> &pair1, const std::pair<Point<dim>,Vector<double>> &pair2)
+ {
+ const double precision = 1e-3;
+ const Point<dim> &p1 = pair1.first;
+ const Point<dim> &p2 = pair2.first;
+
+ for (unsigned int d = 0; d < dim; d++)
+ {
+ const bool is_equal = ( std::abs(p1[d] - p2[d]) < precision );
+ if (!is_equal)
+ return p1[d] < p2[d];
+ }
+
+ return false; //< They are equal
+ }
+};
+
+// ------------
+// |0 |0 | |
+// |-----| 2 |
+// |1 |1 | |
+// ------------
+template<int dim>
+void test2cells(const FiniteElement<dim> &fe_0,
+ const FiniteElement<dim> &fe_1,
+ const FiniteElement<dim> &fe_2,
+ const FiniteElement<dim> &fe_common)
+{
+ const unsigned int n_comp = fe_0.n_components ();
+ Triangulation<dim> triangulation;
+ {
+ Point<dim> p1,p2;
+ for (unsigned int d=0; d<dim; d++)
+ p2[d] = 1.0;
+ p1[0] = -1.0;
+ std::vector<unsigned int> repetitoins(dim,1);
+ repetitoins[0]=2;
+ GridGenerator::subdivided_hyper_rectangle(triangulation,
+ repetitoins,
+ p1,
+ p2);
+ }
+
+ hp::DoFHandler<dim> dof_handler(triangulation);
+
+ hp::FECollection<dim> fe_collection;
+ fe_collection.push_back(fe_0);
+ fe_collection.push_back(fe_1);
+ fe_collection.push_back(fe_2);
+ fe_collection.push_back(fe_common);
+
+ const QIterated<dim-1> quad_formula(QTrapez<1>(),2);
+
+ hp::QCollection< dim-1 > q_face_collection;
+ q_face_collection.push_back(quad_formula);
+ q_face_collection.push_back(quad_formula);
+ q_face_collection.push_back(quad_formula);
+ q_face_collection.push_back(quad_formula);
+
+ deallog << "hp-ref: "<<fe_0.get_name()<< " vs "<<fe_1.get_name()<<" vs " << fe_2.get_name() <<" with "<< fe_common.get_name()<<std::endl;
+
+ // refine once:
+ triangulation.begin_active()->set_refine_flag();
+ triangulation.execute_coarsening_and_refinement();
+
+ for (typename hp::DoFHandler<dim>::active_cell_iterator
+ cell = dof_handler.begin_active();
+ cell != dof_handler.end(); ++cell)
+ {
+ const Point<dim> ¢er = cell->center();
+ if (center[0] > 0.0 ) // right cell
+ cell->set_active_fe_index(2);
+ else if (center[1] > 0.5)
+ cell->set_active_fe_index(0);
+ else
+ cell->set_active_fe_index(1);
+ }
+ dof_handler.distribute_dofs(fe_collection);
+
+ ConstraintMatrix constraints;
+ constraints.clear();
+ dealii::DoFTools::make_hanging_node_constraints (dof_handler, constraints);
+ constraints.close ();
+
+ constraints.print(deallog.get_file_stream());
+
+#ifdef DEBUG_OUTPUT_VTK
+ // output to check if all is good:
+ counter++;
+ std::vector<Vector<double>> shape_functions;
+ std::vector<std::string> names;
+ for (unsigned int s=0; s < dof_handler.n_dofs(); s++)
+ {
+ Vector<double> shape_function;
+ shape_function.reinit(dof_handler.n_dofs());
+ shape_function[s] = 1.0;
+
+ // if the dof is constrained, first output unconstrained vector
+ if (constraints.is_constrained(s))
+ {
+ names.push_back(std::string("UN_") +
+ dealii::Utilities::int_to_string(s,2));
+ shape_functions.push_back(shape_function);
+ }
+
+ names.push_back(std::string("N_") +
+ dealii::Utilities::int_to_string(s,2));
+
+ // make continuous/constrain:
+ constraints.distribute(shape_function);
+ shape_functions.push_back(shape_function);
+ }
+
+ DataOut<dim,hp::DoFHandler<dim>> data_out;
+ data_out.attach_dof_handler (dof_handler);
+
+ // get material ids:
+ Vector<float> fe_index(triangulation.n_active_cells());
+ typename hp::DoFHandler<dim>::active_cell_iterator
+ cell = dof_handler.begin_active (),
+ endc = dof_handler.end ();
+ for (unsigned int index=0; cell!=endc; ++cell,++index)
+ {
+ fe_index[index] = cell->active_fe_index();
+ }
+ data_out.add_data_vector(fe_index, "fe_index");
+
+ for (unsigned int i = 0; i < shape_functions.size(); i++)
+ data_out.add_data_vector (shape_functions[i], names[i]);
+
+ data_out.build_patches(0);
+
+ std::string filename = "shape_functions_"+dealii::Utilities::int_to_string(counter,1)
+ +"_"+dealii::Utilities::int_to_string(dim)+"D.vtu";
+ std::ofstream output (filename.c_str ());
+ data_out.write_vtu (output);
+#endif
+
+ // fill some vector
+ Vector<double> solution(dof_handler.n_dofs());
+ for (unsigned int dof=0; dof < dof_handler.n_dofs(); dof++)
+ solution[dof] = 1.5*(dof%7) + 0.5*dim + 2.0*(dof%3);
+
+ constraints.distribute(solution);
+
+
+ // evaluate field at the interface:
+ hp::FEFaceValues< dim > fe_face_values_hp(fe_collection,
+ q_face_collection,
+ update_values | update_quadrature_points);
+
+ std::vector<unsigned int> local_face_dof_indices;
+ std::vector<Vector<double> > values;
+
+ std::vector<std::pair<Point<dim>,Vector<double>> > pairs_point_value;
+
+ for (typename hp::DoFHandler<dim>::active_cell_iterator
+ cell = dof_handler.begin_active();
+ cell != dof_handler.end(); ++cell)
+ {
+ const unsigned int fe_index = cell->active_fe_index();
+ local_face_dof_indices.resize(fe_collection[fe_index].dofs_per_face);
+ for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
+ if (std::abs(cell->face(f)->center()[0]) < 0.1)
+ {
+ // deallog << "cell="<<cell<<" face="<<f<<std::endl;
+ fe_face_values_hp.reinit(cell,f);
+ const FEFaceValues<dim> &fe_face_values = fe_face_values_hp.get_present_fe_values();
+
+ const unsigned int n_q_points = fe_face_values.n_quadrature_points;
+ values.resize (n_q_points,
+ Vector<double>(n_comp));
+
+ fe_face_values.get_function_values (solution,values);
+
+ const std::vector<dealii::Point<dim> > &q_points = fe_face_values.get_quadrature_points();
+
+ for (unsigned int q = 0; q < n_q_points; q++)
+ {
+ // since our face is [0,1]^{dim-1}, the quadrature rule
+ // will coincide with quadrature points on mother face.
+ // Use that to limit output of sub-faces at the same quadrature
+ // points only.
+ Point<dim-1> qpt;
+ for (unsigned int d=0; d < dim-1; d++)
+ qpt[d] = q_points[q][d+1];
+
+ unsigned int q_found = 0;
+ for (; q_found < quad_formula.size(); q_found++)
+ if (quad_formula.point(q_found).distance(qpt) < 1e-5)
+ break;
+
+ if (q_found < quad_formula.size())
+ {
+ pairs_point_value.push_back(
+ std::make_pair(q_points[q],
+ values[q]));
+// deallog << "@"<<q_points[q]<<" u = {"<<values[q][0];
+// for (unsigned int c = 1; c < n_comp; c++)
+// deallog << ","<<values[q][c];
+// deallog << "}"<<std::endl;
+ }
+ }
+ }
+ }
+
+ // sort
+ std::sort(pairs_point_value.begin(),
+ pairs_point_value.end(),
+ less_than_key<dim>());
+ for (unsigned int p=0; p < pairs_point_value.size(); p++)
+ {
+ const Point<dim> &pt = pairs_point_value[p].first;
+ const Vector<double> &val = pairs_point_value[p].second;
+
+ Assert(val.size() == n_comp,
+ ExcInternalError());
+ deallog << "@"<<pt<<" u = {"<<val[0];
+ for (unsigned int c = 1; c < n_comp; c++)
+ deallog << ","<<val[c];
+ deallog << "}"<<std::endl;
+ }
+
+ dof_handler.clear();
+}
+int main (int argc,char **argv)
+{
+ std::ofstream logfile ("output");
+ deallog << std::setprecision(4);
+ deallog << std::fixed;
+ deallog.attach(logfile);
+ deallog.depth_console(0);
+ deallog.threshold_double(1.e-10);
+
+
+ try
+ {
+ {
+ const unsigned int dim = 2;
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ }
+
+ {
+ const unsigned int dim = 3;
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(2),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(2),1),
+ FESystem<dim>(FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1,
+ FE_Q<dim>(1),1)
+ );
+ }
+ }
+ catch (std::exception &exc)
+ {
+ std::cerr << std::endl << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+ std::cerr << "Exception on processing: " << std::endl
+ << exc.what() << std::endl
+ << "Aborting!" << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+
+ return 1;
+ }
+ catch (...)
+ {
+ std::cerr << std::endl << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+ std::cerr << "Unknown exception!" << std::endl
+ << "Aborting!" << std::endl
+ << "----------------------------------------------------"
+ << std::endl;
+ return 1;
+ };
+}
--- /dev/null
+
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+ 6 24: 0.500000
+ 6 51: 0.500000
+ 20 30: 0.500000
+ 20 33: 0.500000
+ 25 23: 0.750000
+ 25 50: 0.250000
+ 27 23: 0.250000
+ 27 33: 0.500000
+ 27 50: 0.250000
+ 43 29: 0.500000
+ 43 32: 0.500000
+ 46 22: 0.500000
+ 46 49: 0.500000
+ 47 23: 0.500000
+ 47 50: 0.500000
+ 48 24: 0.500000
+ 48 51: 0.500000
+ 52 22: 0.250000
+ 52 49: 0.750000
+ 53 22: 0.250000
+ 53 32: 0.500000
+ 53 49: 0.250000
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.5000 u = {3.7500,7.2500,4.7500}
+DEAL::@0.0000 0.5000 u = {3.7500,7.2500,4.7500}
+DEAL::@0.0000 0.5000 u = {3.7500,7.2500,4.7500}
+DEAL::@0.0000 1.0000 u = {3.0000,6.5000,4.0000}
+DEAL::@0.0000 1.0000 u = {3.0000,6.5000,4.0000}
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(2)-FE_Q<2>(2)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+ 6 32: 0.500000
+ 6 59: 0.500000
+ 7 33: 0.500000
+ 7 60: 0.500000
+ 8 34: 0.500000
+ 8 61: 0.500000
+ 30 40: 0.500000
+ 30 43: 0.500000
+ 35 33: 0.750000
+ 35 60: 0.250000
+ 37 33: 0.250000
+ 37 43: 0.500000
+ 37 60: 0.250000
+ 53 39: 0.500000
+ 53 42: 0.500000
+ 56 32: 0.500000
+ 56 59: 0.500000
+ 57 33: 0.500000
+ 57 60: 0.500000
+ 58 34: 0.500000
+ 58 61: 0.500000
+ 62 32: 0.250000
+ 62 59: 0.750000
+ 63 32: 0.250000
+ 63 42: 0.500000
+ 63 59: 0.250000
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.5000 u = {10.2500,7.7500,11.2500}
+DEAL::@0.0000 0.5000 u = {10.2500,7.7500,11.2500}
+DEAL::@0.0000 0.5000 u = {10.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 u = {9.5000,7.0000,10.5000}
+DEAL::@0.0000 1.0000 u = {9.5000,7.0000,10.5000}
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(2)-FE_Q<2>(2)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+ 6 32: 0.500000
+ 6 54: 0.500000
+ 7 33: 0.500000
+ 7 55: 0.500000
+ 8 34: 0.500000
+ 8 56: 0.500000
+ 30 40: 0.500000
+ 30 43: 0.500000
+ 35 33: 0.750000
+ 35 55: 0.250000
+ 37 33: 0.250000
+ 37 43: 0.500000
+ 37 55: 0.250000
+ 51 32: 0.500000
+ 51 54: 0.500000
+ 52 33: 0.500000
+ 52 55: 0.500000
+ 53 34: 0.500000
+ 53 56: 0.500000
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.5000 u = {9.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 u = {9.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 u = {9.7500,10.2500,8.5000}
+DEAL::@0.0000 1.0000 u = {8.5000,12.0000,5.0000}
+DEAL::@0.0000 1.0000 u = {8.5000,12.0000,5.0000}
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+ 6 24: 0.500000
+ 6 50: 0.500000
+ 25 22: 0.750000
+ 25 48: 0.250000
+ 45 22: 0.500000
+ 45 48: 0.500000
+ 46 23: 0.500000
+ 46 49: 0.500000
+ 47 24: 0.500000
+ 47 50: 0.500000
+ 51 22: 0.250000
+ 51 48: 0.750000
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.5000 u = {7.2500,5.5000,6.0000}
+DEAL::@0.0000 0.5000 u = {7.2500,5.5000,6.0000}
+DEAL::@0.0000 0.5000 u = {7.2500,5.5000,6.0000}
+DEAL::@0.0000 1.0000 u = {10.0000,3.0000,6.5000}
+DEAL::@0.0000 1.0000 u = {10.0000,3.0000,6.5000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+ 12 64: 0.500000
+ 12 130: 0.500000
+ 16 173: 0.500000
+ 16 215: 0.500000
+ 20 64: 0.500000
+ 20 173: 0.500000
+ 22 130: 0.500000
+ 22 215: 0.500000
+ 24 64: 0.250000
+ 24 130: 0.250000
+ 24 173: 0.250000
+ 24 215: 0.250000
+ 46 83: 0.500000
+ 46 86: 0.500000
+ 50 95: 0.500000
+ 50 98: 0.500000
+ 53 83: 0.500000
+ 53 95: 0.500000
+ 54 86: 0.500000
+ 54 98: 0.500000
+ 58 83: 0.250000
+ 58 86: 0.250000
+ 58 95: 0.250000
+ 58 98: 0.250000
+ 65 62: 0.500000
+ 65 171: 0.500000
+ 66 63: 0.500000
+ 66 172: 0.500000
+ 67 64: 0.500000
+ 67 173: 0.500000
+ 68 63: 0.750000
+ 68 129: 0.250000
+ 70 63: 0.250000
+ 70 86: 0.500000
+ 70 129: 0.250000
+ 71 63: 0.375000
+ 71 129: 0.125000
+ 71 172: 0.375000
+ 71 214: 0.125000
+ 73 63: 0.125000
+ 73 98: 0.500000
+ 73 129: 0.125000
+ 73 172: 0.125000
+ 73 214: 0.125000
+ 74 63: 0.750000
+ 74 172: 0.250000
+ 75 63: 0.375000
+ 75 129: 0.375000
+ 75 172: 0.125000
+ 75 214: 0.125000
+ 76 63: 0.562500
+ 76 129: 0.187500
+ 76 172: 0.187500
+ 76 214: 0.0625000
+ 78 63: 0.187500
+ 78 86: 0.250000
+ 78 98: 0.250000
+ 78 129: 0.187500
+ 78 172: 0.0625000
+ 78 214: 0.0625000
+ 108 82: 0.500000
+ 108 85: 0.500000
+ 112 94: 0.500000
+ 112 97: 0.500000
+ 114 82: 0.500000
+ 114 94: 0.500000
+ 115 85: 0.500000
+ 115 97: 0.500000
+ 120 82: 0.250000
+ 120 85: 0.250000
+ 120 94: 0.250000
+ 120 97: 0.250000
+ 125 62: 0.500000
+ 125 128: 0.500000
+ 126 63: 0.500000
+ 126 129: 0.500000
+ 127 64: 0.500000
+ 127 130: 0.500000
+ 131 62: 0.250000
+ 131 128: 0.250000
+ 131 171: 0.250000
+ 131 213: 0.250000
+ 132 63: 0.250000
+ 132 129: 0.250000
+ 132 172: 0.250000
+ 132 214: 0.250000
+ 133 64: 0.250000
+ 133 130: 0.250000
+ 133 173: 0.250000
+ 133 215: 0.250000
+ 134 128: 0.500000
+ 134 213: 0.500000
+ 135 129: 0.500000
+ 135 214: 0.500000
+ 136 130: 0.500000
+ 136 215: 0.500000
+ 137 62: 0.250000
+ 137 128: 0.750000
+ 138 62: 0.250000
+ 138 85: 0.500000
+ 138 128: 0.250000
+ 140 62: 0.125000
+ 140 128: 0.375000
+ 140 171: 0.125000
+ 140 213: 0.375000
+ 141 62: 0.125000
+ 141 97: 0.500000
+ 141 128: 0.125000
+ 141 171: 0.125000
+ 141 213: 0.125000
+ 143 62: 0.375000
+ 143 128: 0.375000
+ 143 171: 0.125000
+ 143 213: 0.125000
+ 144 128: 0.750000
+ 144 213: 0.250000
+ 145 62: 0.187500
+ 145 128: 0.562500
+ 145 171: 0.0625000
+ 145 213: 0.187500
+ 146 62: 0.187500
+ 146 85: 0.250000
+ 146 97: 0.250000
+ 146 128: 0.187500
+ 146 171: 0.0625000
+ 146 213: 0.0625000
+ 160 185: 0.500000
+ 160 188: 0.500000
+ 163 95: 0.500000
+ 163 185: 0.500000
+ 164 98: 0.500000
+ 164 188: 0.500000
+ 168 95: 0.250000
+ 168 98: 0.250000
+ 168 185: 0.250000
+ 168 188: 0.250000
+ 174 172: 0.750000
+ 174 214: 0.250000
+ 176 172: 0.250000
+ 176 188: 0.500000
+ 176 214: 0.250000
+ 177 63: 0.250000
+ 177 172: 0.750000
+ 178 63: 0.125000
+ 178 129: 0.125000
+ 178 172: 0.375000
+ 178 214: 0.375000
+ 179 63: 0.187500
+ 179 129: 0.0625000
+ 179 172: 0.562500
+ 179 214: 0.187500
+ 181 63: 0.0625000
+ 181 98: 0.250000
+ 181 129: 0.0625000
+ 181 172: 0.187500
+ 181 188: 0.250000
+ 181 214: 0.187500
+ 198 184: 0.500000
+ 198 187: 0.500000
+ 200 94: 0.500000
+ 200 184: 0.500000
+ 201 97: 0.500000
+ 201 187: 0.500000
+ 206 94: 0.250000
+ 206 97: 0.250000
+ 206 184: 0.250000
+ 206 187: 0.250000
+ 210 171: 0.500000
+ 210 213: 0.500000
+ 211 172: 0.500000
+ 211 214: 0.500000
+ 212 173: 0.500000
+ 212 215: 0.500000
+ 216 171: 0.250000
+ 216 213: 0.750000
+ 217 171: 0.250000
+ 217 187: 0.500000
+ 217 213: 0.250000
+ 219 62: 0.125000
+ 219 128: 0.125000
+ 219 171: 0.375000
+ 219 213: 0.375000
+ 220 128: 0.250000
+ 220 213: 0.750000
+ 221 62: 0.0625000
+ 221 128: 0.187500
+ 221 171: 0.187500
+ 221 213: 0.562500
+ 222 62: 0.0625000
+ 222 97: 0.250000
+ 222 128: 0.0625000
+ 222 171: 0.187500
+ 222 187: 0.250000
+ 222 213: 0.187500
+DEAL::@0.0000 0.0000 0.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.0000 0.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.2500,5.5000,9.0000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.2500,5.5000,9.0000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.2500,5.5000,9.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {11.5000,3.7500,7.2500}
+DEAL::@0.0000 0.5000 0.0000 u = {11.5000,3.7500,7.2500}
+DEAL::@0.0000 0.5000 0.0000 u = {11.5000,3.7500,7.2500}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {8.5000,6.0000,9.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {8.5000,6.0000,9.5000}
+DEAL::@0.0000 1.0000 0.5000 u = {7.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {7.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {7.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 1.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(2)-FE_Q<3>(2)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+ 12 100: 0.500000
+ 12 166: 0.500000
+ 13 101: 0.500000
+ 13 167: 0.500000
+ 14 102: 0.500000
+ 14 168: 0.500000
+ 24 209: 0.500000
+ 24 251: 0.500000
+ 25 210: 0.500000
+ 25 252: 0.500000
+ 26 211: 0.500000
+ 26 253: 0.500000
+ 36 100: 0.500000
+ 36 209: 0.500000
+ 37 101: 0.500000
+ 37 210: 0.500000
+ 38 102: 0.500000
+ 38 211: 0.500000
+ 42 166: 0.500000
+ 42 251: 0.500000
+ 43 167: 0.500000
+ 43 252: 0.500000
+ 44 168: 0.500000
+ 44 253: 0.500000
+ 48 100: 0.250000
+ 48 166: 0.250000
+ 48 209: 0.250000
+ 48 251: 0.250000
+ 49 101: 0.250000
+ 49 167: 0.250000
+ 49 210: 0.250000
+ 49 252: 0.250000
+ 50 102: 0.250000
+ 50 168: 0.250000
+ 50 211: 0.250000
+ 50 253: 0.250000
+ 84 121: 0.500000
+ 84 124: 0.500000
+ 88 133: 0.500000
+ 88 136: 0.500000
+ 91 121: 0.500000
+ 91 133: 0.500000
+ 92 124: 0.500000
+ 92 136: 0.500000
+ 96 121: 0.250000
+ 96 124: 0.250000
+ 96 133: 0.250000
+ 96 136: 0.250000
+ 103 100: 0.500000
+ 103 209: 0.500000
+ 104 101: 0.500000
+ 104 210: 0.500000
+ 105 102: 0.500000
+ 105 211: 0.500000
+ 106 101: 0.750000
+ 106 167: 0.250000
+ 108 101: 0.250000
+ 108 124: 0.500000
+ 108 167: 0.250000
+ 109 101: 0.375000
+ 109 167: 0.125000
+ 109 210: 0.375000
+ 109 252: 0.125000
+ 111 101: 0.125000
+ 111 136: 0.500000
+ 111 167: 0.125000
+ 111 210: 0.125000
+ 111 252: 0.125000
+ 112 101: 0.750000
+ 112 210: 0.250000
+ 113 101: 0.375000
+ 113 167: 0.375000
+ 113 210: 0.125000
+ 113 252: 0.125000
+ 114 101: 0.562500
+ 114 167: 0.187500
+ 114 210: 0.187500
+ 114 252: 0.0625000
+ 116 101: 0.187500
+ 116 124: 0.250000
+ 116 136: 0.250000
+ 116 167: 0.187500
+ 116 210: 0.0625000
+ 116 252: 0.0625000
+ 146 120: 0.500000
+ 146 123: 0.500000
+ 150 132: 0.500000
+ 150 135: 0.500000
+ 152 120: 0.500000
+ 152 132: 0.500000
+ 153 123: 0.500000
+ 153 135: 0.500000
+ 158 120: 0.250000
+ 158 123: 0.250000
+ 158 132: 0.250000
+ 158 135: 0.250000
+ 163 100: 0.500000
+ 163 166: 0.500000
+ 164 101: 0.500000
+ 164 167: 0.500000
+ 165 102: 0.500000
+ 165 168: 0.500000
+ 169 100: 0.250000
+ 169 166: 0.250000
+ 169 209: 0.250000
+ 169 251: 0.250000
+ 170 101: 0.250000
+ 170 167: 0.250000
+ 170 210: 0.250000
+ 170 252: 0.250000
+ 171 102: 0.250000
+ 171 168: 0.250000
+ 171 211: 0.250000
+ 171 253: 0.250000
+ 172 166: 0.500000
+ 172 251: 0.500000
+ 173 167: 0.500000
+ 173 252: 0.500000
+ 174 168: 0.500000
+ 174 253: 0.500000
+ 175 100: 0.250000
+ 175 166: 0.750000
+ 176 100: 0.250000
+ 176 123: 0.500000
+ 176 166: 0.250000
+ 178 100: 0.125000
+ 178 166: 0.375000
+ 178 209: 0.125000
+ 178 251: 0.375000
+ 179 100: 0.125000
+ 179 135: 0.500000
+ 179 166: 0.125000
+ 179 209: 0.125000
+ 179 251: 0.125000
+ 181 100: 0.375000
+ 181 166: 0.375000
+ 181 209: 0.125000
+ 181 251: 0.125000
+ 182 166: 0.750000
+ 182 251: 0.250000
+ 183 100: 0.187500
+ 183 166: 0.562500
+ 183 209: 0.0625000
+ 183 251: 0.187500
+ 184 100: 0.187500
+ 184 123: 0.250000
+ 184 135: 0.250000
+ 184 166: 0.187500
+ 184 209: 0.0625000
+ 184 251: 0.0625000
+ 198 223: 0.500000
+ 198 226: 0.500000
+ 201 133: 0.500000
+ 201 223: 0.500000
+ 202 136: 0.500000
+ 202 226: 0.500000
+ 206 133: 0.250000
+ 206 136: 0.250000
+ 206 223: 0.250000
+ 206 226: 0.250000
+ 212 210: 0.750000
+ 212 252: 0.250000
+ 214 210: 0.250000
+ 214 226: 0.500000
+ 214 252: 0.250000
+ 215 101: 0.250000
+ 215 210: 0.750000
+ 216 101: 0.125000
+ 216 167: 0.125000
+ 216 210: 0.375000
+ 216 252: 0.375000
+ 217 101: 0.187500
+ 217 167: 0.0625000
+ 217 210: 0.562500
+ 217 252: 0.187500
+ 219 101: 0.0625000
+ 219 136: 0.250000
+ 219 167: 0.0625000
+ 219 210: 0.187500
+ 219 226: 0.250000
+ 219 252: 0.187500
+ 236 222: 0.500000
+ 236 225: 0.500000
+ 238 132: 0.500000
+ 238 222: 0.500000
+ 239 135: 0.500000
+ 239 225: 0.500000
+ 244 132: 0.250000
+ 244 135: 0.250000
+ 244 222: 0.250000
+ 244 225: 0.250000
+ 248 209: 0.500000
+ 248 251: 0.500000
+ 249 210: 0.500000
+ 249 252: 0.500000
+ 250 211: 0.500000
+ 250 253: 0.500000
+ 254 209: 0.250000
+ 254 251: 0.750000
+ 255 209: 0.250000
+ 255 225: 0.500000
+ 255 251: 0.250000
+ 257 100: 0.125000
+ 257 166: 0.125000
+ 257 209: 0.375000
+ 257 251: 0.375000
+ 258 166: 0.250000
+ 258 251: 0.750000
+ 259 100: 0.0625000
+ 259 166: 0.187500
+ 259 209: 0.187500
+ 259 251: 0.562500
+ 260 100: 0.0625000
+ 260 135: 0.250000
+ 260 166: 0.0625000
+ 260 209: 0.187500
+ 260 225: 0.250000
+ 260 251: 0.187500
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.5000,5.7500,6.2500}
+DEAL::@0.0000 0.0000 0.5000 u = {10.5000,5.7500,6.2500}
+DEAL::@0.0000 0.0000 0.5000 u = {10.5000,5.7500,6.2500}
+DEAL::@0.0000 0.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {8.7500,12.2500,4.5000}
+DEAL::@0.0000 0.5000 0.0000 u = {8.7500,12.2500,4.5000}
+DEAL::@0.0000 0.5000 0.0000 u = {8.7500,12.2500,4.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {11.0000,14.5000,1.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {11.0000,14.5000,1.5000}
+DEAL::@0.0000 1.0000 0.5000 u = {12.7500,8.0000,3.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {12.7500,8.0000,3.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {12.7500,8.0000,3.2500}
+DEAL::@0.0000 1.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 1.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(2)-FE_Q<3>(2)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+ 12 100: 0.500000
+ 12 147: 0.500000
+ 13 101: 0.500000
+ 13 148: 0.500000
+ 14 102: 0.500000
+ 14 149: 0.500000
+ 24 176: 0.500000
+ 24 204: 0.500000
+ 25 177: 0.500000
+ 25 205: 0.500000
+ 26 178: 0.500000
+ 26 206: 0.500000
+ 36 100: 0.500000
+ 36 176: 0.500000
+ 37 101: 0.500000
+ 37 177: 0.500000
+ 38 102: 0.500000
+ 38 178: 0.500000
+ 42 147: 0.500000
+ 42 204: 0.500000
+ 43 148: 0.500000
+ 43 205: 0.500000
+ 44 149: 0.500000
+ 44 206: 0.500000
+ 48 100: 0.250000
+ 48 147: 0.250000
+ 48 176: 0.250000
+ 48 204: 0.250000
+ 49 101: 0.250000
+ 49 148: 0.250000
+ 49 177: 0.250000
+ 49 205: 0.250000
+ 50 102: 0.250000
+ 50 149: 0.250000
+ 50 178: 0.250000
+ 50 206: 0.250000
+ 84 121: 0.500000
+ 84 124: 0.500000
+ 88 133: 0.500000
+ 88 136: 0.500000
+ 91 121: 0.500000
+ 91 133: 0.500000
+ 92 124: 0.500000
+ 92 136: 0.500000
+ 96 121: 0.250000
+ 96 124: 0.250000
+ 96 133: 0.250000
+ 96 136: 0.250000
+ 103 100: 0.500000
+ 103 176: 0.500000
+ 104 101: 0.500000
+ 104 177: 0.500000
+ 105 102: 0.500000
+ 105 178: 0.500000
+ 106 101: 0.750000
+ 106 148: 0.250000
+ 108 101: 0.250000
+ 108 124: 0.500000
+ 108 148: 0.250000
+ 109 101: 0.375000
+ 109 148: 0.125000
+ 109 177: 0.375000
+ 109 205: 0.125000
+ 111 101: 0.125000
+ 111 136: 0.500000
+ 111 148: 0.125000
+ 111 177: 0.125000
+ 111 205: 0.125000
+ 112 101: 0.750000
+ 112 177: 0.250000
+ 113 101: 0.375000
+ 113 148: 0.375000
+ 113 177: 0.125000
+ 113 205: 0.125000
+ 114 101: 0.562500
+ 114 148: 0.187500
+ 114 177: 0.187500
+ 114 205: 0.0625000
+ 116 101: 0.187500
+ 116 124: 0.250000
+ 116 136: 0.250000
+ 116 148: 0.187500
+ 116 177: 0.0625000
+ 116 205: 0.0625000
+ 144 100: 0.500000
+ 144 147: 0.500000
+ 145 101: 0.500000
+ 145 148: 0.500000
+ 146 102: 0.500000
+ 146 149: 0.500000
+ 150 100: 0.250000
+ 150 147: 0.250000
+ 150 176: 0.250000
+ 150 204: 0.250000
+ 151 101: 0.250000
+ 151 148: 0.250000
+ 151 177: 0.250000
+ 151 205: 0.250000
+ 152 102: 0.250000
+ 152 149: 0.250000
+ 152 178: 0.250000
+ 152 206: 0.250000
+ 153 147: 0.500000
+ 153 204: 0.500000
+ 154 148: 0.500000
+ 154 205: 0.500000
+ 155 149: 0.500000
+ 155 206: 0.500000
+ 165 190: 0.500000
+ 165 193: 0.500000
+ 168 133: 0.500000
+ 168 190: 0.500000
+ 169 136: 0.500000
+ 169 193: 0.500000
+ 173 133: 0.250000
+ 173 136: 0.250000
+ 173 190: 0.250000
+ 173 193: 0.250000
+ 179 177: 0.750000
+ 179 205: 0.250000
+ 181 177: 0.250000
+ 181 193: 0.500000
+ 181 205: 0.250000
+ 182 101: 0.250000
+ 182 177: 0.750000
+ 183 101: 0.125000
+ 183 148: 0.125000
+ 183 177: 0.375000
+ 183 205: 0.375000
+ 184 101: 0.187500
+ 184 148: 0.0625000
+ 184 177: 0.562500
+ 184 205: 0.187500
+ 186 101: 0.0625000
+ 186 136: 0.250000
+ 186 148: 0.0625000
+ 186 177: 0.187500
+ 186 193: 0.250000
+ 186 205: 0.187500
+ 201 176: 0.500000
+ 201 204: 0.500000
+ 202 177: 0.500000
+ 202 205: 0.500000
+ 203 178: 0.500000
+ 203 206: 0.500000
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {6.7500,7.2500,7.7500}
+DEAL::@0.0000 0.0000 0.5000 u = {6.7500,7.2500,7.7500}
+DEAL::@0.0000 0.0000 0.5000 u = {6.7500,7.2500,7.7500}
+DEAL::@0.0000 0.0000 1.0000 u = {7.0000,4.5000,8.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {7.0000,4.5000,8.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.0000,7.5000,8.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.0000,7.5000,8.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.0000,7.5000,8.0000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 1.0000 u = {5.0000,5.5000,9.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {5.0000,5.5000,9.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {5.0000,5.5000,9.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {1.5000,5.0000,8.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {1.5000,5.0000,8.5000}
+DEAL::@0.0000 1.0000 0.5000 u = {2.2500,5.7500,9.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {2.2500,5.7500,9.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {2.2500,5.7500,9.2500}
+DEAL::@0.0000 1.0000 1.0000 u = {3.0000,6.5000,10.0000}
+DEAL::@0.0000 1.0000 1.0000 u = {3.0000,6.5000,10.0000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+ 12 60: 0.500000
+ 12 122: 0.500000
+ 16 161: 0.500000
+ 16 200: 0.500000
+ 20 60: 0.500000
+ 20 161: 0.500000
+ 22 122: 0.500000
+ 22 200: 0.500000
+ 24 60: 0.250000
+ 24 122: 0.250000
+ 24 161: 0.250000
+ 24 200: 0.250000
+ 61 58: 0.500000
+ 61 159: 0.500000
+ 62 59: 0.500000
+ 62 160: 0.500000
+ 63 60: 0.500000
+ 63 161: 0.500000
+ 64 58: 0.750000
+ 64 120: 0.250000
+ 66 58: 0.375000
+ 66 120: 0.125000
+ 66 159: 0.375000
+ 66 198: 0.125000
+ 68 58: 0.750000
+ 68 159: 0.250000
+ 69 58: 0.562500
+ 69 120: 0.187500
+ 69 159: 0.187500
+ 69 198: 0.0625000
+ 117 58: 0.500000
+ 117 120: 0.500000
+ 118 59: 0.500000
+ 118 121: 0.500000
+ 119 60: 0.500000
+ 119 122: 0.500000
+ 123 58: 0.250000
+ 123 120: 0.250000
+ 123 159: 0.250000
+ 123 198: 0.250000
+ 124 59: 0.250000
+ 124 121: 0.250000
+ 124 160: 0.250000
+ 124 199: 0.250000
+ 125 60: 0.250000
+ 125 122: 0.250000
+ 125 161: 0.250000
+ 125 200: 0.250000
+ 126 120: 0.500000
+ 126 198: 0.500000
+ 127 121: 0.500000
+ 127 199: 0.500000
+ 128 122: 0.500000
+ 128 200: 0.500000
+ 129 58: 0.250000
+ 129 120: 0.750000
+ 132 58: 0.125000
+ 132 120: 0.375000
+ 132 159: 0.125000
+ 132 198: 0.375000
+ 135 58: 0.375000
+ 135 120: 0.375000
+ 135 159: 0.125000
+ 135 198: 0.125000
+ 136 120: 0.750000
+ 136 198: 0.250000
+ 137 58: 0.187500
+ 137 120: 0.562500
+ 137 159: 0.0625000
+ 137 198: 0.187500
+ 162 159: 0.750000
+ 162 198: 0.250000
+ 164 58: 0.250000
+ 164 159: 0.750000
+ 165 58: 0.187500
+ 165 120: 0.0625000
+ 165 159: 0.562500
+ 165 198: 0.187500
+ 195 159: 0.500000
+ 195 198: 0.500000
+ 196 160: 0.500000
+ 196 199: 0.500000
+ 197 161: 0.500000
+ 197 200: 0.500000
+ 201 159: 0.250000
+ 201 198: 0.750000
+ 204 58: 0.125000
+ 204 120: 0.125000
+ 204 159: 0.375000
+ 204 198: 0.375000
+ 205 120: 0.250000
+ 205 198: 0.750000
+ 206 58: 0.0625000
+ 206 120: 0.187500
+ 206 159: 0.187500
+ 206 198: 0.562500
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {7.7500,11.2500,6.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {7.7500,11.2500,6.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {7.7500,11.2500,6.5000}
+DEAL::@0.0000 0.0000 1.0000 u = {9.0000,12.5000,5.5000}
+DEAL::@0.0000 0.0000 1.0000 u = {9.0000,12.5000,5.5000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.7500,8.2500,8.7500}
+DEAL::@0.0000 0.5000 0.0000 u = {4.7500,8.2500,8.7500}
+DEAL::@0.0000 0.5000 0.0000 u = {4.7500,8.2500,8.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 1.0000 u = {6.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.7500,10.2500,8.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {3.0000,6.5000,10.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {3.0000,6.5000,10.0000}
+DEAL::@0.0000 1.0000 0.5000 u = {3.7500,7.2500,10.7500}
+DEAL::@0.0000 1.0000 0.5000 u = {3.7500,7.2500,10.7500}
+DEAL::@0.0000 1.0000 0.5000 u = {3.7500,7.2500,10.7500}
+DEAL::@0.0000 1.0000 1.0000 u = {4.5000,8.0000,11.5000}
+DEAL::@0.0000 1.0000 1.0000 u = {4.5000,8.0000,11.5000}
projects / dealii.git / commitdiff