From: Martin Kronbichler Date: Tue, 17 Sep 2013 14:24:00 +0000 (+0000) Subject: Fix subface interpolation. X-Git-Tag: v8.1.0~789 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7b26d73402efc243994350dc0cb094db5698e9fb;p=dealii.git Fix subface interpolation. git-svn-id: https://svn.dealii.org/trunk@30765 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/fe/fe_face.cc b/deal.II/source/fe/fe_face.cc index a362e3cdec..641e20a684 100644 --- a/deal.II/source/fe/fe_face.cc +++ b/deal.II/source/fe/fe_face.cc @@ -436,7 +436,7 @@ get_subface_interpolation_matrix (const FiniteElement &x_source_fe // zero. const double eps = 2e-13*(this->degree+1)*(dim-1); - FullMatrix mass (face_quadrature.size(), this->dofs_per_face); + FullMatrix mass (face_quadrature.size(), source_fe->dofs_per_face); for (unsigned int k = 0; k < face_quadrature.size(); ++k) { @@ -446,18 +446,18 @@ get_subface_interpolation_matrix (const FiniteElement &x_source_fe GeometryInfo::child_to_cell_coordinates (face_quadrature.point(k), subface); - for (unsigned int j = 0; j < this->dofs_per_face; ++j) - mass (k , j) = this->poly_space.compute_value(j, p); + for (unsigned int j = 0; j < source_fe->dofs_per_face; ++j) + mass (k , j) = source_fe->poly_space.compute_value(j, p); } Householder H(mass); Vector v_in(face_quadrature.size()); - Vector v_out(this->dofs_per_face); + Vector v_out(source_fe->dofs_per_face); // compute the interpolation matrix by evaluating on the fine side and // then solving the least squares problem - for (unsigned int i=0; idofs_per_face; ++i) + for (unsigned int i=0; idofs_per_face; ++i) { for (unsigned int k = 0; k < face_quadrature.size(); ++k) { @@ -465,12 +465,12 @@ get_subface_interpolation_matrix (const FiniteElement &x_source_fe face_quadrature.point(k) : GeometryInfo::child_to_cell_coordinates (face_quadrature.point(k), subface); - v_in(k) = source_fe->poly_space.compute_value(i, p); + v_in(k) = this->poly_space.compute_value(i, p); } const double result = H.least_squares(v_out, v_in); Assert(result < 1e-12, FETools::ExcLeastSquaresError (result)); - for (unsigned int j = 0; j < this->dofs_per_face; ++j) + for (unsigned int j = 0; j < source_fe->dofs_per_face; ++j) { double matrix_entry = v_out(j); @@ -482,21 +482,9 @@ get_subface_interpolation_matrix (const FiniteElement &x_source_fe if (std::fabs (matrix_entry) < eps) matrix_entry = 0.0; - interpolation_matrix(i,j) = matrix_entry; + interpolation_matrix(j,i) = matrix_entry; } } - - // make sure that the row sum of each of the matrices is 1 at this - // point. this must be so since the shape functions sum up to 1 - for (unsigned int j=0; jdofs_per_face; ++j) - { - double sum = 0.; - - for (unsigned int i=0; idofs_per_face; ++i) - sum += interpolation_matrix(j,i); - - Assert (std::fabs(sum-1) < eps, ExcInternalError()); - } } else if (dynamic_cast *>(&x_source_fe) != 0) { diff --git a/tests/bits/faceq_face_interpolation.cc b/tests/bits/faceq_face_interpolation.cc index fe8f50e6dd..ca0d906738 100644 --- a/tests/bits/faceq_face_interpolation.cc +++ b/tests/bits/faceq_face_interpolation.cc @@ -19,7 +19,7 @@ #include // check -// FE_FaceQ::face_interpolation +// FE_FaceQ/P::face_interpolation std::string output_file_name = "faceq_face_interpolation/output"; @@ -125,12 +125,22 @@ void check (const unsigned int degree1, const unsigned int degree2) { - deallog << "Checking FE_FaceQ in " << dim << "d:" - << std::endl; + { + deallog << "Checking FE_FaceQ in " << dim << "d:" + << std::endl; + + FE_FaceQ fe1 (degree1), fe2 (degree2); + + check_this (fe1, fe2); + } + { + deallog << "Checking FE_FaceP in " << dim << "d:" + << std::endl; - FE_FaceQ fe1 (degree1), fe2 (degree2); + FE_FaceP fe1 (degree1), fe2 (degree2); - check_this (fe1, fe2); + check_this (fe1, fe2); + } } diff --git a/tests/bits/faceq_face_interpolation/cmp/generic b/tests/bits/faceq_face_interpolation/cmp/generic index 5b2aa67ac2..f477f4dd29 100644 --- a/tests/bits/faceq_face_interpolation/cmp/generic +++ b/tests/bits/faceq_face_interpolation/cmp/generic @@ -16,6 +16,23 @@ DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 DEAL::1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 @@ -28,6 +45,18 @@ DEAL::1.00000 0 1.00000 0 DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(1) DEAL::2.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::1.00000 1.00000 +DEAL::1.41421 +DEAL::1.00000 0 1.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 @@ -40,6 +69,18 @@ DEAL::1.00000 0 1.00000 1.00000 1.00000 0 DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(2) DEAL::3.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::1.00000 1.00000 +DEAL::1.73205 +DEAL::1.00000 0 1.00000 1.00000 1.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 @@ -52,6 +93,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(3) DEAL::4.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::2.00000 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(1) DEAL::1.00000 1.00000 @@ -64,6 +117,18 @@ DEAL::1.00000 0 1.00000 1.00000 1.00000 0 DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(2) DEAL::1.50000 1.00000 DEAL::1.00000 0 0.500000 0.500000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::1.00000 1.00000 +DEAL::1.41421 +DEAL::1.00000 0 1.00000 0 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::1.00000 1.00000 +DEAL::1.73205 +DEAL::1.00000 0 1.00000 1.00000 1.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 0 1.00000 0 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(1) DEAL::1.00000 1.00000 @@ -76,6 +141,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(3) DEAL::2.00000 1.00000 DEAL::1.00000 0 0.276393 0.723607 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::1.00000 1.00000 +DEAL::1.41421 +DEAL::1.00000 0 1.00000 0 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::2.00000 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 0 1.00000 0 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(2) vs. FE_FaceQ<2>(2) DEAL::1.00000 1.00000 @@ -88,6 +165,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 DEAL::FE_FaceQ<2>(2) vs. FE_FaceQ<2>(3) DEAL::1.60000 1.24721 DEAL::1.00000 0 0.800000 0.800000 0.323607 -0.123607 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::1.00000 1.00000 +DEAL::1.73205 +DEAL::1.00000 0 1.00000 1.00000 1.00000 0 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::2.00000 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 0 1.00000 1.00000 1.00000 0 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -105,6 +194,23 @@ DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 DEAL::1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -117,6 +223,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(1) DEAL::4.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.73205 +DEAL::1.00000 0 1.00000 1.00000 1.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -129,6 +247,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 1.00000 1.00000 0 1.00000 DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(2) DEAL::9.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::2.44949 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -141,6 +271,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00 DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(3) DEAL::16.0000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::3.16228 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(1) DEAL::1.00000 1.00000 @@ -153,6 +295,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 1.00000 1.00000 0 1.00000 DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(2) DEAL::2.25000 1.00000 DEAL::1.00000 0 0.500000 0 0 1.00000 0 0.500000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.73205 +DEAL::1.00000 0 1.00000 1.00000 1.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::2.44949 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 0 1.00000 1.00000 0 0 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(1) DEAL::1.00000 1.00000 @@ -165,6 +319,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00 DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(3) DEAL::4.00000 1.00000 DEAL::1.00000 0 0.276393 0 0 0.723607 0 0 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.73205 +DEAL::1.00000 0 1.00000 1.00000 1.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::3.16228 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 0 1.00000 1.00000 0 0 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(2) vs. FE_FaceQ<3>(2) DEAL::1.00000 1.00000 @@ -177,3 +343,15 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00 DEAL::FE_FaceQ<3>(2) vs. FE_FaceQ<3>(3) DEAL::2.56000 1.55554 DEAL::1.00000 0 0.800000 0 0.323607 0 0 0 0 0 -0.0988854 0.258885 0.0152786 0.104721 0 0 0 -0.123607 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::2.44949 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::3.16228 +DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 0 1.00000 0 1.00000 0 0 0 0 0 0 0 diff --git a/tests/bits/faceq_subface_interpolation.cc b/tests/bits/faceq_subface_interpolation.cc index c8bfa05d5e..ecf6cb3196 100644 --- a/tests/bits/faceq_subface_interpolation.cc +++ b/tests/bits/faceq_subface_interpolation.cc @@ -19,7 +19,7 @@ #include // check -// FE_FaceQ::subface_interpolation +// FE_FaceQ/P::subface_interpolation std::string output_file_name = "faceq_subface_interpolation/output"; @@ -128,12 +128,22 @@ void check (const unsigned int degree1, const unsigned int degree2) { - deallog << "Checking FE_FaceQ in " << dim << "d:" - << std::endl; + { + deallog << "Checking FE_FaceQ in " << dim << "d:" + << std::endl; + + FE_FaceQ fe1 (degree1), fe2 (degree2); + + check_this (fe1, fe2); + } + { + deallog << "Checking FE_FaceP in " << dim << "d:" + << std::endl; - FE_FaceQ fe1 (degree1), fe2 (degree2); + FE_FaceP fe1 (degree1), fe2 (degree2); - check_this (fe1, fe2); + check_this (fe1, fe2); + } } diff --git a/tests/bits/faceq_subface_interpolation/cmp/generic b/tests/bits/faceq_subface_interpolation/cmp/generic index df01dec611..5c0103af34 100644 --- a/tests/bits/faceq_subface_interpolation/cmp/generic +++ b/tests/bits/faceq_subface_interpolation/cmp/generic @@ -32,6 +32,39 @@ DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 DEAL::1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 @@ -55,6 +88,29 @@ DEAL::0.500000 0.500000 1.00000 0 DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(1) DEAL::2.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::3.73205 2.73205 +DEAL::2.82843 +DEAL::1.00000 1.73205 2.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::3.73205 2.73205 +DEAL::2.82843 +DEAL::1.00000 -1.73205 2.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 @@ -78,6 +134,29 @@ DEAL::0 0 0.750000 0.750000 1.00000 0 DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(2) DEAL::3.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::18.4542 9.74597 +DEAL::11.3578 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::18.4542 9.74597 +DEAL::11.3578 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(0) DEAL::1.00000 1.00000 @@ -101,6 +180,29 @@ DEAL::-0.125000 -0.125000 0.220492 0.934017 0.779508 -0.184017 1.00000 0 DEAL::FE_FaceQ<2>(0) vs. FE_FaceQ<2>(3) DEAL::4.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::97.4282 46.4066 +DEAL::54.2586 +DEAL::1.00000 29.1033 2.00000 7.74597 4.00000 0 8.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::97.4282 46.4066 +DEAL::54.2586 +DEAL::1.00000 -29.1033 2.00000 -7.74597 4.00000 0 8.00000 0 +DEAL::FE_FaceP<2>(0) vs. FE_FaceP<2>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(1) DEAL::1.50000 1.00000 @@ -124,6 +226,29 @@ DEAL::0 0 0.750000 0.750000 1.00000 0 DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(2) DEAL::2.25000 1.00000 DEAL::0.500000 0.500000 0.750000 0.250000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::3.73205 2.73205 +DEAL::2.82843 +DEAL::1.00000 1.73205 2.00000 0 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::18.4542 9.74597 +DEAL::11.3578 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(2) +DEAL::3.73205 2.73205 +DEAL::1.00000 1.73205 2.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::3.73205 2.73205 +DEAL::2.82843 +DEAL::1.00000 -1.73205 2.00000 0 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::18.4542 9.74597 +DEAL::11.3578 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(2) +DEAL::3.73205 2.73205 +DEAL::1.00000 -1.73205 2.00000 0 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(1) DEAL::1.50000 1.00000 @@ -147,6 +272,29 @@ DEAL::-0.125000 -0.125000 0.220492 0.934017 0.779508 -0.184017 1.00000 0 DEAL::FE_FaceQ<2>(1) vs. FE_FaceQ<2>(3) DEAL::3.00000 1.00000 DEAL::0.500000 0.500000 0.638197 0.361803 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::3.73205 2.73205 +DEAL::2.82843 +DEAL::1.00000 1.73205 2.00000 0 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::97.4282 46.4066 +DEAL::54.2586 +DEAL::1.00000 29.1033 2.00000 7.74597 4.00000 0 8.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(3) +DEAL::3.73205 2.73205 +DEAL::1.00000 1.73205 2.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(1) +DEAL::3.73205 2.73205 +DEAL::2.82843 +DEAL::1.00000 -1.73205 2.00000 0 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::97.4282 46.4066 +DEAL::54.2586 +DEAL::1.00000 -29.1033 2.00000 -7.74597 4.00000 0 8.00000 0 +DEAL::FE_FaceP<2>(1) vs. FE_FaceP<2>(3) +DEAL::3.73205 2.73205 +DEAL::1.00000 -1.73205 2.00000 0 DEAL::Checking FE_FaceQ in 2d: DEAL::FE_FaceQ<2>(2) vs. FE_FaceQ<2>(2) DEAL::1.75000 1.25000 @@ -170,6 +318,29 @@ DEAL::-0.125000 -0.125000 0.220492 0.934017 0.779508 -0.184017 1.00000 0 DEAL::FE_FaceQ<2>(2) vs. FE_FaceQ<2>(3) DEAL::2.40000 1.20000 DEAL::0 0 0.923607 0.923607 0.623607 -0.100000 +DEAL::Checking FE_FaceP in 2d: +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::18.4542 9.74597 +DEAL::11.3578 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::97.4282 46.4066 +DEAL::54.2586 +DEAL::1.00000 29.1033 2.00000 7.74597 4.00000 0 8.00000 0 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(3) +DEAL::18.4542 9.74597 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(2) +DEAL::18.4542 9.74597 +DEAL::11.3578 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 +DEAL::FE_FaceP<2>(3) vs. FE_FaceP<2>(3) +DEAL::97.4282 46.4066 +DEAL::54.2586 +DEAL::1.00000 -29.1033 2.00000 -7.74597 4.00000 0 8.00000 0 +DEAL::FE_FaceP<2>(2) vs. FE_FaceP<2>(3) +DEAL::18.4542 9.74597 +DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -235,6 +406,71 @@ DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 DEAL::1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -280,6 +516,51 @@ DEAL::0.250000 0.250000 0.500000 0 0.500000 0 1.00000 0 DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(1) DEAL::4.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(1) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -325,6 +606,51 @@ DEAL::0 0 0 0 0 0 0 0 0.562500 0.562500 0.750000 0 0 0 0.750000 0 1.00000 0 DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(2) DEAL::9.00000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(2) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(0) DEAL::1.00000 1.00000 @@ -370,6 +696,51 @@ DEAL::0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 DEAL::FE_FaceQ<3>(0) vs. FE_FaceQ<3>(3) DEAL::16.0000 1.00000 DEAL::1.00000 1.00000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(0) +DEAL::1.00000 1.00000 +DEAL::1.00000 +DEAL::1.00000 1.00000 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(0) vs. FE_FaceP<3>(3) +DEAL::1.00000 1.00000 +DEAL::1.00000 1.00000 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(1) DEAL::2.25000 1.00000 @@ -415,6 +786,51 @@ DEAL::0 0 0 0 0 0 0 0 0.562500 0.562500 0.750000 0 0 0 0.750000 0 1.00000 0 DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(2) DEAL::5.06250 1.00000 DEAL::0.250000 0.250000 0.375000 0.125000 0 0.500000 0.375000 0.125000 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(2) +DEAL::3.73205 4.46410 +DEAL::1.00000 1.73205 2.00000 2.00000 0 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(2) +DEAL::3.73205 4.46410 +DEAL::1.00000 1.73205 2.00000 2.00000 0 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(2) +DEAL::3.73205 4.46410 +DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(2) +DEAL::3.73205 4.46410 +DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(1) DEAL::2.25000 1.00000 @@ -460,6 +876,51 @@ DEAL::0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 DEAL::FE_FaceQ<3>(1) vs. FE_FaceQ<3>(3) DEAL::9.00000 1.00000 DEAL::0.250000 0.250000 0.319098 0.180902 0.0690983 0.430902 0.500000 0 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(3) +DEAL::3.73205 4.46410 +DEAL::1.00000 1.73205 2.00000 2.00000 0 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(3) +DEAL::3.73205 4.46410 +DEAL::1.00000 1.73205 2.00000 2.00000 0 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(3) +DEAL::3.73205 4.46410 +DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(1) +DEAL::3.73205 4.46410 +DEAL::3.87298 +DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(1) vs. FE_FaceP<3>(3) +DEAL::3.73205 4.46410 +DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 DEAL::Checking FE_FaceQ in 3d: DEAL::FE_FaceQ<3>(2) vs. FE_FaceQ<3>(2) DEAL::3.06250 1.56250 @@ -505,3 +966,48 @@ DEAL::0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 DEAL::FE_FaceQ<3>(2) vs. FE_FaceQ<3>(3) DEAL::5.76000 1.44000 DEAL::0 0 0 0 0 0 0 1.00000 0.923607 0.923607 0.162918 -0.0923607 -0.0176393 -0.0623607 0 0 0 0 +DEAL::Checking FE_FaceP in 3d: +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(3) +DEAL::18.4542 20.8805 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 0 0 3.46410 0 0 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(3) +DEAL::18.4542 20.8805 +DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 0 0 -3.46410 0 0 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(3) +DEAL::18.4542 20.8805 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 0 0 3.46410 0 0 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(2) +DEAL::18.4542 20.8805 +DEAL::17.4929 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 +DEAL::FE_FaceP<3>(3) vs. FE_FaceP<3>(3) +DEAL::97.4282 102.325 +DEAL::87.3499 +DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 +DEAL::FE_FaceP<3>(2) vs. FE_FaceP<3>(3) +DEAL::18.4542 20.8805 +DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 0 0 -3.46410 0 0 0