From: Denis Davydov Date: Wed, 2 Sep 2015 20:51:22 +0000 (+0200) Subject: test make_hanging_node_constraints() with hp-ref and neither_element_dominates X-Git-Tag: v8.4.0-rc2~484^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=599b12b1549383f2cbe995b6d3dd37ca0f3974c8;p=dealii.git test make_hanging_node_constraints() with hp-ref and neither_element_dominates --- diff --git a/tests/hp/hp_constraints_neither_dominate_02.cc b/tests/hp/hp_constraints_neither_dominate_02.cc new file mode 100644 index 0000000000..a9b6626cc4 --- /dev/null +++ b/tests/hp/hp_constraints_neither_dominate_02.cc @@ -0,0 +1,456 @@ +// --------------------------------------------------------------------- +// +// 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 +#include +#include +#include +#include +#include + +#include +#include + +#include +#include +#include +#include + + +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include +#include +#include + +#include + +#include + +#include +#include + +#include +#include +#include +#include + +#include "../tests.h" + +#include +#include + +//#define DEBUG_OUTPUT_VTK +unsigned int counter=0; + +const double eps = 1e-10; + +using namespace dealii; + +template +struct less_than_key +{ + inline bool operator() (const std::pair,Vector> &pair1, const std::pair,Vector> &pair2) + { + const double precision = 1e-3; + const Point &p1 = pair1.first; + const Point &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 +void test2cells(const FiniteElement &fe_0, + const FiniteElement &fe_1, + const FiniteElement &fe_2, + const FiniteElement &fe_common) +{ + const unsigned int n_comp = fe_0.n_components (); + Triangulation triangulation; + { + Point p1,p2; + for (unsigned int d=0; d repetitoins(dim,1); + repetitoins[0]=2; + GridGenerator::subdivided_hyper_rectangle(triangulation, + repetitoins, + p1, + p2); + } + + hp::DoFHandler dof_handler(triangulation); + + hp::FECollection 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 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: "<set_refine_flag(); + triangulation.execute_coarsening_and_refinement(); + + for (typename hp::DoFHandler::active_cell_iterator + cell = dof_handler.begin_active(); + cell != dof_handler.end(); ++cell) + { + const Point ¢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> shape_functions; + std::vector names; + for (unsigned int s=0; s < dof_handler.n_dofs(); s++) + { + Vector 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> data_out; + data_out.attach_dof_handler (dof_handler); + + // get material ids: + Vector fe_index(triangulation.n_active_cells()); + typename hp::DoFHandler::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 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 local_face_dof_indices; + std::vector > values; + + std::vector,Vector> > pairs_point_value; + + for (typename hp::DoFHandler::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::faces_per_cell; ++f) + if (std::abs(cell->face(f)->center()[0]) < 0.1) + { + // deallog << "cell="<()); + for (unsigned int p=0; p < pairs_point_value.size(); p++) + { + const Point &pt = pairs_point_value[p].first; + const Vector &val = pairs_point_value[p].second; + + Assert(val.size() == n_comp, + ExcInternalError()); + deallog << "@"<(FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(2),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1) + ); + test2cells(FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(2),1, + FE_Q(1),1), + FESystem(FE_Q(2),1, + FE_Q(2),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1) + ); + test2cells(FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(2),1, + FE_Q(1),1), + FESystem(FE_Q(2),1, + FE_Q(2),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1) + ); + test2cells(FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1) + ); + } + + { + const unsigned int dim = 3; + test2cells(FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(2),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1) + ); + test2cells(FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(2),1, + FE_Q(1),1), + FESystem(FE_Q(2),1, + FE_Q(2),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1) + ); + test2cells(FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(2),1, + FE_Q(1),1), + FESystem(FE_Q(2),1, + FE_Q(2),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(1),1) + ); + test2cells(FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(2),1, + FE_Q(1),1, + FE_Q(1),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(2),1), + FESystem(FE_Q(1),1, + FE_Q(1),1, + FE_Q(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; + }; +} diff --git a/tests/hp/hp_constraints_neither_dominate_02.output b/tests/hp/hp_constraints_neither_dominate_02.output new file mode 100644 index 0000000000..66d61ec0e2 --- /dev/null +++ b/tests/hp/hp_constraints_neither_dominate_02.output @@ -0,0 +1,869 @@ + +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}