From dae6562052b3c5a2e32d1f32dba6cd9d8d419259 Mon Sep 17 00:00:00 2001 From: kronbichler Date: Wed, 30 Apr 2014 09:15:00 +0000 Subject: [PATCH] Fix FE_DGQArbitraryNodes::has_support_on_face git-svn-id: https://svn.dealii.org/trunk@32863 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/doc/news/changes.h | 7 + deal.II/source/fe/fe_dgq.cc | 13 +- tests/fe/fe_data_test.cc | 5 + tests/fe/fe_data_test.output | 348 +++++++++++++++++++++++++++-------- 4 files changed, 292 insertions(+), 81 deletions(-) diff --git a/deal.II/doc/news/changes.h b/deal.II/doc/news/changes.h index c442348d1b..d4be862392 100644 --- a/deal.II/doc/news/changes.h +++ b/deal.II/doc/news/changes.h @@ -148,6 +148,13 @@ inconvenience this causes.

Specific improvements

    +
  1. Fixed: FE_DGQArbitraryNodes::has_support_on_face was broken when + polynomials with support points not on the element boundary were used. + This is now fixed. +
    + (Martin Kronbichler, 2014/04/30) +
    2. +
  2. Fixed: parallel::distributed::Triangulation::load now has an additional parameter autopartition to control p4est's behavior of rebalancing triangulations between MPI nodes upon reading. It is diff --git a/deal.II/source/fe/fe_dgq.cc b/deal.II/source/fe/fe_dgq.cc index 8fe34b6c60..0ca4d1c021 100644 --- a/deal.II/source/fe/fe_dgq.cc +++ b/deal.II/source/fe/fe_dgq.cc @@ -696,10 +696,17 @@ FE_DGQ::has_support_on_face (const unsigned int shape_index, unsigned int n = this->degree+1; - // for DGQ(0) elements, the single - // shape functions is constant and + // for DGQ(0) elements or arbitrary node DGQ with support points not located + // at the element boundary, the single shape functions is constant and // therefore lives on the boundary - if (this->degree == 0) + bool support_points_on_boundary = true; + for (unsigned int d=0; dunit_support_points[0][d]) > 1e-13) + support_points_on_boundary = false; + for (unsigned int d=0; dunit_support_points.back()[d]-1.) > 1e-13) + support_points_on_boundary = false; + if (support_points_on_boundary == false) return true; unsigned int n2 = n*n; diff --git a/tests/fe/fe_data_test.cc b/tests/fe/fe_data_test.cc index a9e16506a1..d554c5cb35 100644 --- a/tests/fe/fe_data_test.cc +++ b/tests/fe/fe_data_test.cc @@ -21,6 +21,7 @@ #include #include +#include #include #include @@ -112,6 +113,10 @@ void test_fe_datas() deallog << (*fe_datas.rbegin())->get_name() << std::endl; fe_datas.push_back(new FE_DGQ (4)); deallog << (*fe_datas.rbegin())->get_name() << std::endl; + fe_datas.push_back(new FE_DGQArbitraryNodes (QGaussLobatto<1>(5))); + deallog << (*fe_datas.rbegin())->get_name() << std::endl; + fe_datas.push_back(new FE_DGQArbitraryNodes (QGauss<1>(3))); + deallog << (*fe_datas.rbegin())->get_name() << std::endl; fe_datas.push_back(new FE_DGP (1)); deallog << (*fe_datas.rbegin())->get_name() << std::endl; fe_datas.push_back(new FE_DGP (2)); diff --git a/tests/fe/fe_data_test.output b/tests/fe/fe_data_test.output index f99beab9f4..b3cdb29e7c 100644 --- a/tests/fe/fe_data_test.output +++ b/tests/fe/fe_data_test.output @@ -8,6 +8,8 @@ DEAL::FE_Q_Hierarchical<1>(4) DEAL::FE_DGQ<1>(1) DEAL::FE_DGQ<1>(2) DEAL::FE_DGQ<1>(4) +DEAL::FE_DGQArbitraryNodes<1>(QGaussLobatto(5)) +DEAL::FE_DGQArbitraryNodes<1>(QUnknownNodes(2)) DEAL::FE_DGP<1>(1) DEAL::FE_DGP<1>(2) DEAL::FESystem<1>[FE_Q<1>(2)^2] @@ -270,7 +272,61 @@ DEAL::face_to_cell_index: DEAL::support on face 0: 0 DEAL::support on face 1: 4 DEAL:: -DEAL::fe_data[9]:FE_DGP<1>(1) +DEAL::fe_data[9]:FE_DGQArbitraryNodes<1>(QGaussLobatto(5)) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=5 +DEAL::dofs_per_quad=0 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=5 +DEAL::first_hex_index=5 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=0 +DEAL::dofs_per_cell=5 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[5]->5 +DEAL::degree=4 +DEAL::conformity= L2 +DEAL::unit_support_points=5 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=5 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::support on face 0: 0 +DEAL::support on face 1: 4 +DEAL:: +DEAL::fe_data[10]:FE_DGQArbitraryNodes<1>(QUnknownNodes(2)) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=3 +DEAL::dofs_per_quad=0 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=3 +DEAL::first_hex_index=3 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=0 +DEAL::dofs_per_cell=3 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[3]->3 +DEAL::degree=2 +DEAL::conformity= L2 +DEAL::unit_support_points=3 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=3 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::support on face 0: 0 1 2 +DEAL::support on face 1: 0 1 2 +DEAL:: +DEAL::fe_data[11]:FE_DGP<1>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=0 @@ -297,7 +353,7 @@ DEAL::face_to_cell_index: DEAL::support on face 0: 0 1 DEAL::support on face 1: 0 1 DEAL:: -DEAL::fe_data[10]:FE_DGP<1>(2) +DEAL::fe_data[12]:FE_DGP<1>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=0 @@ -324,7 +380,7 @@ DEAL::face_to_cell_index: DEAL::support on face 0: 0 1 2 DEAL::support on face 1: 0 1 2 DEAL:: -DEAL::fe_data[11]:FESystem<1>[FE_Q<1>(2)^2] +DEAL::fe_data[13]:FESystem<1>[FE_Q<1>(2)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=0 @@ -351,7 +407,7 @@ DEAL::face_to_cell_index: 2 3 DEAL::support on face 0: 0 1 DEAL::support on face 1: 2 3 DEAL:: -DEAL::fe_data[12]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(2)] +DEAL::fe_data[14]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(2)] DEAL::dofs_per_vertex=3 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -378,7 +434,7 @@ DEAL::face_to_cell_index: 3 4 5 DEAL::support on face 0: 0 1 2 DEAL::support on face 1: 3 4 5 DEAL:: -DEAL::fe_data[13]:FE_FaceQ<1>(0) +DEAL::fe_data[15]:FE_FaceQ<1>(0) DEAL::dofs_per_vertex=1 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -405,7 +461,7 @@ DEAL::face_to_cell_index: 1 DEAL::support on face 0: 0 DEAL::support on face 1: 1 DEAL:: -DEAL::fe_data[14]:FE_FaceQ<1>(1) +DEAL::fe_data[16]:FE_FaceQ<1>(1) DEAL::dofs_per_vertex=1 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -432,7 +488,7 @@ DEAL::face_to_cell_index: 1 DEAL::support on face 0: 0 DEAL::support on face 1: 1 DEAL:: -DEAL::fe_data[15]:FE_FaceQ<1>(3) +DEAL::fe_data[17]:FE_FaceQ<1>(3) DEAL::dofs_per_vertex=1 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -459,7 +515,7 @@ DEAL::face_to_cell_index: 1 DEAL::support on face 0: 0 DEAL::support on face 1: 1 DEAL:: -DEAL::fe_data[16]:FE_FaceP<1>(0) +DEAL::fe_data[18]:FE_FaceP<1>(0) DEAL::dofs_per_vertex=1 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -486,7 +542,7 @@ DEAL::face_to_cell_index: 1 DEAL::support on face 0: 0 DEAL::support on face 1: 1 DEAL:: -DEAL::fe_data[17]:FE_FaceP<1>(1) +DEAL::fe_data[19]:FE_FaceP<1>(1) DEAL::dofs_per_vertex=1 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -513,7 +569,7 @@ DEAL::face_to_cell_index: 1 DEAL::support on face 0: 0 DEAL::support on face 1: 1 DEAL:: -DEAL::fe_data[18]:FE_FaceP<1>(3) +DEAL::fe_data[20]:FE_FaceP<1>(3) DEAL::dofs_per_vertex=1 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -540,7 +596,7 @@ DEAL::face_to_cell_index: 1 DEAL::support on face 0: 0 DEAL::support on face 1: 1 DEAL:: -DEAL::fe_data[19]:FESystem<1>[FE_Q<1>(3)^2] +DEAL::fe_data[21]:FESystem<1>[FE_Q<1>(3)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=4 DEAL::dofs_per_quad=0 @@ -567,7 +623,7 @@ DEAL::face_to_cell_index: 2 3 DEAL::support on face 0: 0 1 DEAL::support on face 1: 2 3 DEAL:: -DEAL::fe_data[20]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(3)] +DEAL::fe_data[22]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(3)] DEAL::dofs_per_vertex=3 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=0 @@ -594,7 +650,7 @@ DEAL::face_to_cell_index: 3 4 5 DEAL::support on face 0: 0 1 2 DEAL::support on face 1: 3 4 5 DEAL:: -DEAL::fe_data[21]:FESystem<1>[FE_Q<1>(4)^2] +DEAL::fe_data[23]:FESystem<1>[FE_Q<1>(4)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=6 DEAL::dofs_per_quad=0 @@ -621,7 +677,7 @@ DEAL::face_to_cell_index: 2 3 DEAL::support on face 0: 0 1 DEAL::support on face 1: 2 3 DEAL:: -DEAL::fe_data[22]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2] +DEAL::fe_data[24]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -648,7 +704,7 @@ DEAL::face_to_cell_index: 4 5 6 7 DEAL::support on face 0: 0 1 2 3 DEAL::support on face 1: 4 5 6 7 DEAL:: -DEAL::fe_data[23]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]-FESystem<1>[FE_DGQ<1>(1)^2]] +DEAL::fe_data[25]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]-FESystem<1>[FE_DGQ<1>(1)^2]] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=4 DEAL::dofs_per_quad=0 @@ -675,7 +731,7 @@ DEAL::face_to_cell_index: 2 3 DEAL::support on face 0: 0 1 4 6 DEAL::support on face 1: 2 3 5 7 DEAL:: -DEAL::fe_data[24]:FESystem<1>[FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]-FESystem<1>[FE_Q<1>(2)^2]-FESystem<1>[FE_DGQ<1>(2)^2]] +DEAL::fe_data[26]:FESystem<1>[FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]-FESystem<1>[FE_Q<1>(2)^2]-FESystem<1>[FE_DGQ<1>(2)^2]] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=9 DEAL::dofs_per_quad=0 @@ -702,7 +758,7 @@ DEAL::face_to_cell_index: 4 5 6 7 DEAL::support on face 0: 0 1 2 3 11 14 DEAL::support on face 1: 4 5 6 7 13 16 DEAL:: -DEAL::fe_data[25]:FESystem<1>[FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2]^2-FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]-FESystem<1>[FE_DGQ<1>(1)^2]]-FESystem<1>[FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]-FESystem<1>[FE_Q<1>(2)^2]-FESystem<1>[FE_DGQ<1>(2)^2]]^2] +DEAL::fe_data[27]:FESystem<1>[FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2]^2-FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]-FESystem<1>[FE_DGQ<1>(1)^2]]-FESystem<1>[FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]-FESystem<1>[FE_Q<1>(2)^2]-FESystem<1>[FE_DGQ<1>(2)^2]]^2] DEAL::dofs_per_vertex=18 DEAL::dofs_per_line=22 DEAL::dofs_per_quad=0 @@ -738,6 +794,8 @@ DEAL::FE_Q_Hierarchical<2>(4) DEAL::FE_DGQ<2>(1) DEAL::FE_DGQ<2>(2) DEAL::FE_DGQ<2>(4) +DEAL::FE_DGQArbitraryNodes<2>(QGaussLobatto(5)) +DEAL::FE_DGQArbitraryNodes<2>(QUnknownNodes(2)) DEAL::FE_DGP<2>(1) DEAL::FE_DGP<2>(2) DEAL::FESystem<2>[FE_Q<2>(2)^2] @@ -1053,7 +1111,69 @@ DEAL::support on face 1: 4 9 14 19 24 DEAL::support on face 2: 0 1 2 3 4 DEAL::support on face 3: 20 21 22 23 24 DEAL:: -DEAL::fe_data[9]:FE_DGP<2>(1) +DEAL::fe_data[9]:FE_DGQArbitraryNodes<2>(QGaussLobatto(5)) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=0 +DEAL::dofs_per_quad=25 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=0 +DEAL::first_hex_index=25 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=0 +DEAL::dofs_per_cell=25 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[25]->25 +DEAL::degree=4 +DEAL::conformity= L2 +DEAL::unit_support_points=25 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=25 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::support on face 0: 0 5 10 15 20 +DEAL::support on face 1: 4 9 14 19 24 +DEAL::support on face 2: 0 1 2 3 4 +DEAL::support on face 3: 20 21 22 23 24 +DEAL:: +DEAL::fe_data[10]:FE_DGQArbitraryNodes<2>(QUnknownNodes(2)) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=0 +DEAL::dofs_per_quad=9 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=0 +DEAL::first_hex_index=9 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=0 +DEAL::dofs_per_cell=9 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[9]->9 +DEAL::degree=2 +DEAL::conformity= L2 +DEAL::unit_support_points=9 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=9 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::support on face 0: 0 1 2 3 4 5 6 7 8 +DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 +DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 +DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 +DEAL:: +DEAL::fe_data[11]:FE_DGP<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=3 @@ -1084,7 +1204,7 @@ DEAL::support on face 1: 0 1 2 DEAL::support on face 2: 0 1 2 DEAL::support on face 3: 0 1 2 DEAL:: -DEAL::fe_data[10]:FE_DGP<2>(2) +DEAL::fe_data[12]:FE_DGP<2>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=6 @@ -1115,7 +1235,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 DEAL::support on face 2: 0 1 2 3 4 5 DEAL::support on face 3: 0 1 2 3 4 5 DEAL:: -DEAL::fe_data[11]:FESystem<2>[FE_Q<2>(2)^2] +DEAL::fe_data[13]:FESystem<2>[FE_Q<2>(2)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=2 @@ -1146,7 +1266,7 @@ DEAL::support on face 1: 2 3 6 7 10 11 DEAL::support on face 2: 0 1 2 3 12 13 DEAL::support on face 3: 4 5 6 7 14 15 DEAL:: -DEAL::fe_data[12]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(2)] +DEAL::fe_data[14]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(2)] DEAL::dofs_per_vertex=3 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=1 @@ -1177,7 +1297,7 @@ DEAL::support on face 1: 3 4 5 9 10 11 13 DEAL::support on face 2: 0 1 2 3 4 5 14 DEAL::support on face 3: 6 7 8 9 10 11 15 DEAL:: -DEAL::fe_data[13]:FE_FaceQ<2>(0) +DEAL::fe_data[15]:FE_FaceQ<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1208,7 +1328,7 @@ DEAL::support on face 1: 1 DEAL::support on face 2: 2 DEAL::support on face 3: 3 DEAL:: -DEAL::fe_data[14]:FE_FaceQ<2>(1) +DEAL::fe_data[16]:FE_FaceQ<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=0 @@ -1239,7 +1359,7 @@ DEAL::support on face 1: 2 3 DEAL::support on face 2: 4 5 DEAL::support on face 3: 6 7 DEAL:: -DEAL::fe_data[15]:FE_FaceQ<2>(3) +DEAL::fe_data[17]:FE_FaceQ<2>(3) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=4 DEAL::dofs_per_quad=0 @@ -1270,7 +1390,7 @@ DEAL::support on face 1: 4 5 6 7 DEAL::support on face 2: 8 9 10 11 DEAL::support on face 3: 12 13 14 15 DEAL:: -DEAL::fe_data[16]:FE_FaceP<2>(0) +DEAL::fe_data[18]:FE_FaceP<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1301,7 +1421,7 @@ DEAL::support on face 1: 1 DEAL::support on face 2: 2 DEAL::support on face 3: 3 DEAL:: -DEAL::fe_data[17]:FE_FaceP<2>(1) +DEAL::fe_data[19]:FE_FaceP<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=0 @@ -1332,7 +1452,7 @@ DEAL::support on face 1: 2 3 DEAL::support on face 2: 4 5 DEAL::support on face 3: 6 7 DEAL:: -DEAL::fe_data[18]:FE_FaceP<2>(3) +DEAL::fe_data[20]:FE_FaceP<2>(3) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=4 DEAL::dofs_per_quad=0 @@ -1363,7 +1483,7 @@ DEAL::support on face 1: 4 5 6 7 DEAL::support on face 2: 8 9 10 11 DEAL::support on face 3: 12 13 14 15 DEAL:: -DEAL::fe_data[19]:FE_DGRaviartThomas<2,2>(0) +DEAL::fe_data[21]:FE_DGRaviartThomas<2,2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -1394,7 +1514,7 @@ DEAL::support on face 1: 0 1 2 3 DEAL::support on face 2: 0 1 2 3 DEAL::support on face 3: 0 1 2 3 DEAL:: -DEAL::fe_data[20]:FE_DGRaviartThomas<2,2>(1) +DEAL::fe_data[22]:FE_DGRaviartThomas<2,2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=12 @@ -1425,7 +1545,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[21]:FE_DGNedelec<2,2>(0) +DEAL::fe_data[23]:FE_DGNedelec<2,2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -1456,7 +1576,7 @@ DEAL::support on face 1: 0 1 2 3 DEAL::support on face 2: 0 1 2 3 DEAL::support on face 3: 0 1 2 3 DEAL:: -DEAL::fe_data[22]:FE_DGNedelec<2,2>(1) +DEAL::fe_data[24]:FE_DGNedelec<2,2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=12 @@ -1487,7 +1607,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[23]:FE_RaviartThomas<2>(0) +DEAL::fe_data[25]:FE_RaviartThomas<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1518,7 +1638,7 @@ DEAL::support on face 1: 1 2 3 DEAL::support on face 2: 0 1 2 DEAL::support on face 3: 0 1 3 DEAL:: -DEAL::fe_data[24]:FE_RaviartThomas<2>(1) +DEAL::fe_data[26]:FE_RaviartThomas<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1549,7 +1669,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[25]:FE_RaviartThomas<2>(2) +DEAL::fe_data[27]:FE_RaviartThomas<2>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=12 @@ -1580,7 +1700,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 DEAL:: -DEAL::fe_data[26]:FESystem<2>[FE_RaviartThomas<2>(1)-FE_DGQ<2>(1)] +DEAL::fe_data[28]:FESystem<2>[FE_RaviartThomas<2>(1)-FE_DGQ<2>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=8 @@ -1611,7 +1731,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 13 15 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 14 15 DEAL:: -DEAL::fe_data[27]:FE_Nedelec<2>(0) +DEAL::fe_data[29]:FE_Nedelec<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1642,7 +1762,7 @@ DEAL::support on face 1: 1 2 3 DEAL::support on face 2: 0 1 2 DEAL::support on face 3: 0 1 3 DEAL:: -DEAL::fe_data[28]:FE_Nedelec<2>(1) +DEAL::fe_data[30]:FE_Nedelec<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1673,7 +1793,7 @@ DEAL::support on face 1: 2 3 4 5 6 7 DEAL::support on face 2: 0 1 2 3 4 5 DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[29]:FE_DGBDM<2,2>(1) +DEAL::fe_data[31]:FE_DGBDM<2,2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=8 @@ -1704,7 +1824,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 DEAL::support on face 2: 0 1 2 3 4 5 6 7 DEAL::support on face 3: 0 1 2 3 4 5 6 7 DEAL:: -DEAL::fe_data[30]:FE_DGBDM<2,2>(2) +DEAL::fe_data[32]:FE_DGBDM<2,2>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=14 @@ -1735,7 +1855,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL:: -DEAL::fe_data[31]:FE_BDM<2>(1) +DEAL::fe_data[33]:FE_BDM<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=0 @@ -1766,7 +1886,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 DEAL::support on face 2: 0 1 2 3 4 5 6 7 DEAL::support on face 3: 0 1 2 3 4 5 6 7 DEAL:: -DEAL::fe_data[32]:FE_BDM<2>(2) +DEAL::fe_data[34]:FE_BDM<2>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=2 @@ -1797,7 +1917,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL:: -DEAL::fe_data[33]:FE_RaviartThomasNodal<2>(0) +DEAL::fe_data[35]:FE_RaviartThomasNodal<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1828,7 +1948,7 @@ DEAL::support on face 1: 1 2 3 DEAL::support on face 2: 0 1 2 DEAL::support on face 3: 0 1 3 DEAL:: -DEAL::fe_data[34]:FE_RaviartThomasNodal<2>(1) +DEAL::fe_data[36]:FE_RaviartThomasNodal<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1859,7 +1979,7 @@ DEAL::support on face 1: 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 8 9 10 11 DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[35]:FESystem<2>[FE_RaviartThomasNodal<2>(1)-FE_DGQ<2>(1)] +DEAL::fe_data[37]:FESystem<2>[FE_RaviartThomasNodal<2>(1)-FE_DGQ<2>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=8 @@ -1890,7 +2010,7 @@ DEAL::support on face 1: 2 3 4 5 6 7 8 9 10 11 13 15 DEAL::support on face 2: 0 1 2 3 4 5 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 14 15 DEAL:: -DEAL::fe_data[36]:FESystem<2>[FE_Q<2>(3)^2] +DEAL::fe_data[38]:FESystem<2>[FE_Q<2>(3)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=4 DEAL::dofs_per_quad=8 @@ -1921,7 +2041,7 @@ DEAL::support on face 1: 2 3 6 7 12 13 14 15 DEAL::support on face 2: 0 1 2 3 16 17 18 19 DEAL::support on face 3: 4 5 6 7 20 21 22 23 DEAL:: -DEAL::fe_data[37]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(3)] +DEAL::fe_data[39]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(3)] DEAL::dofs_per_vertex=3 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1952,7 +2072,7 @@ DEAL::support on face 1: 3 4 5 9 10 11 14 15 DEAL::support on face 2: 0 1 2 3 4 5 16 17 DEAL::support on face 3: 6 7 8 9 10 11 18 19 DEAL:: -DEAL::fe_data[38]:FESystem<2>[FE_Q<2>(4)^2] +DEAL::fe_data[40]:FESystem<2>[FE_Q<2>(4)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=6 DEAL::dofs_per_quad=18 @@ -1983,7 +2103,7 @@ DEAL::support on face 1: 2 3 6 7 14 15 16 17 18 19 DEAL::support on face 2: 0 1 2 3 20 21 22 23 24 25 DEAL::support on face 3: 4 5 6 7 26 27 28 29 30 31 DEAL:: -DEAL::fe_data[39]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2] +DEAL::fe_data[41]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2014,7 +2134,7 @@ DEAL::support on face 1: 4 5 6 7 12 13 14 15 DEAL::support on face 2: 0 1 2 3 4 5 6 7 DEAL::support on face 3: 8 9 10 11 12 13 14 15 DEAL:: -DEAL::fe_data[40]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]] +DEAL::fe_data[42]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=8 @@ -2045,7 +2165,7 @@ DEAL::support on face 1: 2 3 6 7 9 11 13 15 DEAL::support on face 2: 0 1 2 3 8 9 12 13 DEAL::support on face 3: 4 5 6 7 10 11 14 15 DEAL:: -DEAL::fe_data[41]:FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]] +DEAL::fe_data[43]:FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=21 @@ -2076,7 +2196,7 @@ DEAL::support on face 1: 4 5 6 7 12 13 14 15 19 20 21 33 36 39 42 45 48 DEAL::support on face 2: 0 1 2 3 4 5 6 7 22 23 24 31 32 33 40 41 42 DEAL::support on face 3: 8 9 10 11 12 13 14 15 25 26 27 37 38 39 46 47 48 DEAL:: -DEAL::fe_data[42]:FESystem<2>[FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2]^2-FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]]-FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]]^2] +DEAL::fe_data[44]:FESystem<2>[FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2]^2-FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]]-FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]]^2] DEAL::dofs_per_vertex=18 DEAL::dofs_per_line=6 DEAL::dofs_per_quad=50 @@ -2116,6 +2236,8 @@ DEAL::FE_Q_Hierarchical<3>(4) DEAL::FE_DGQ<3>(1) DEAL::FE_DGQ<3>(2) DEAL::FE_DGQ<3>(4) +DEAL::FE_DGQArbitraryNodes<3>(QGaussLobatto(5)) +DEAL::FE_DGQArbitraryNodes<3>(QUnknownNodes(2)) DEAL::FE_DGP<3>(1) DEAL::FE_DGP<3>(2) DEAL::FESystem<3>[FE_Q<3>(2)^2] @@ -2460,7 +2582,77 @@ DEAL::support on face 3: 20 21 22 23 24 45 46 47 48 49 70 71 72 73 74 95 96 97 9 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 DEAL::support on face 5: 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 DEAL:: -DEAL::fe_data[9]:FE_DGP<3>(1) +DEAL::fe_data[9]:FE_DGQArbitraryNodes<3>(QGaussLobatto(5)) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=0 +DEAL::dofs_per_quad=0 +DEAL::dofs_per_hex=125 +DEAL::first_line_index=0 +DEAL::first_quad_index=0 +DEAL::first_hex_index=0 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=0 +DEAL::dofs_per_cell=125 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[125]->125 +DEAL::degree=4 +DEAL::conformity= L2 +DEAL::unit_support_points=125 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=125 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::support on face 0: 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 +DEAL::support on face 1: 4 9 14 19 24 29 34 39 44 49 54 59 64 69 74 79 84 89 94 99 104 109 114 119 124 +DEAL::support on face 2: 0 1 2 3 4 25 26 27 28 29 50 51 52 53 54 75 76 77 78 79 100 101 102 103 104 +DEAL::support on face 3: 20 21 22 23 24 45 46 47 48 49 70 71 72 73 74 95 96 97 98 99 120 121 122 123 124 +DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 +DEAL::support on face 5: 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 +DEAL:: +DEAL::fe_data[10]:FE_DGQArbitraryNodes<3>(QUnknownNodes(2)) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=0 +DEAL::dofs_per_quad=0 +DEAL::dofs_per_hex=27 +DEAL::first_line_index=0 +DEAL::first_quad_index=0 +DEAL::first_hex_index=0 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=0 +DEAL::dofs_per_cell=27 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[27]->27 +DEAL::degree=2 +DEAL::conformity= L2 +DEAL::unit_support_points=27 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=27 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::face_to_cell_index: +DEAL::support on face 0: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 +DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 +DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 +DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 +DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 +DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 +DEAL:: +DEAL::fe_data[11]:FE_DGP<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2495,7 +2687,7 @@ DEAL::support on face 3: 0 1 2 3 DEAL::support on face 4: 0 1 2 3 DEAL::support on face 5: 0 1 2 3 DEAL:: -DEAL::fe_data[10]:FE_DGP<3>(2) +DEAL::fe_data[12]:FE_DGP<3>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2530,7 +2722,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 DEAL:: -DEAL::fe_data[11]:FESystem<3>[FE_Q<3>(2)^2] +DEAL::fe_data[13]:FESystem<3>[FE_Q<3>(2)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=2 @@ -2565,7 +2757,7 @@ DEAL::support on face 3: 4 5 6 7 12 13 14 15 22 23 30 31 36 37 38 39 46 47 DEAL::support on face 4: 0 1 2 3 4 5 6 7 16 17 18 19 20 21 22 23 48 49 DEAL::support on face 5: 8 9 10 11 12 13 14 15 24 25 26 27 28 29 30 31 50 51 DEAL:: -DEAL::fe_data[12]:FESystem<3>[FE_Q<3>(1)^2-FE_Q<3>(2)] +DEAL::fe_data[14]:FESystem<3>[FE_Q<3>(1)^2-FE_Q<3>(2)] DEAL::dofs_per_vertex=3 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=1 @@ -2600,7 +2792,7 @@ DEAL::support on face 3: 6 7 8 9 10 11 18 19 20 21 22 23 27 31 34 35 39 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 24 25 26 27 40 DEAL::support on face 5: 12 13 14 15 16 17 18 19 20 21 22 23 28 29 30 31 41 DEAL:: -DEAL::fe_data[13]:FE_FaceQ<3>(0) +DEAL::fe_data[15]:FE_FaceQ<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=1 @@ -2635,7 +2827,7 @@ DEAL::support on face 3: 3 DEAL::support on face 4: 4 DEAL::support on face 5: 5 DEAL:: -DEAL::fe_data[14]:FE_FaceQ<3>(1) +DEAL::fe_data[16]:FE_FaceQ<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -2670,7 +2862,7 @@ DEAL::support on face 3: 12 13 14 15 DEAL::support on face 4: 16 17 18 19 DEAL::support on face 5: 20 21 22 23 DEAL:: -DEAL::fe_data[15]:FE_FaceQ<3>(3) +DEAL::fe_data[17]:FE_FaceQ<3>(3) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=16 @@ -2705,7 +2897,7 @@ DEAL::support on face 3: 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 DEAL::support on face 4: 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 DEAL::support on face 5: 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 DEAL:: -DEAL::fe_data[16]:FE_FaceP<3>(0) +DEAL::fe_data[18]:FE_FaceP<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=1 @@ -2740,7 +2932,7 @@ DEAL::support on face 3: 3 DEAL::support on face 4: 4 DEAL::support on face 5: 5 DEAL:: -DEAL::fe_data[17]:FE_FaceP<3>(1) +DEAL::fe_data[19]:FE_FaceP<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=3 @@ -2775,7 +2967,7 @@ DEAL::support on face 3: 9 10 11 DEAL::support on face 4: 12 13 14 DEAL::support on face 5: 15 16 17 DEAL:: -DEAL::fe_data[18]:FE_FaceP<3>(3) +DEAL::fe_data[20]:FE_FaceP<3>(3) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=10 @@ -2810,7 +3002,7 @@ DEAL::support on face 3: 30 31 32 33 34 35 36 37 38 39 DEAL::support on face 4: 40 41 42 43 44 45 46 47 48 49 DEAL::support on face 5: 50 51 52 53 54 55 56 57 58 59 DEAL:: -DEAL::fe_data[19]:FE_DGRaviartThomas<3,3>(0) +DEAL::fe_data[21]:FE_DGRaviartThomas<3,3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2845,7 +3037,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 DEAL::support on face 4: 0 1 2 3 4 5 DEAL::support on face 5: 0 1 2 3 4 5 DEAL:: -DEAL::fe_data[20]:FE_DGRaviartThomas<3,3>(1) +DEAL::fe_data[22]:FE_DGRaviartThomas<3,3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2880,7 +3072,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL:: -DEAL::fe_data[21]:FE_DGNedelec<3,3>(0) +DEAL::fe_data[23]:FE_DGNedelec<3,3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2915,7 +3107,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[22]:FE_DGNedelec<3,3>(1) +DEAL::fe_data[24]:FE_DGNedelec<3,3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2950,7 +3142,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 DEAL:: -DEAL::fe_data[23]:FE_RaviartThomas<3>(0) +DEAL::fe_data[25]:FE_RaviartThomas<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=1 @@ -2985,7 +3177,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 DEAL::support on face 4: 0 1 2 3 4 5 DEAL::support on face 5: 0 1 2 3 4 5 DEAL:: -DEAL::fe_data[24]:FE_RaviartThomas<3>(1) +DEAL::fe_data[26]:FE_RaviartThomas<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -3020,7 +3212,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL:: -DEAL::fe_data[25]:FE_RaviartThomas<3>(2) +DEAL::fe_data[27]:FE_RaviartThomas<3>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=9 @@ -3055,7 +3247,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 DEAL:: -DEAL::fe_data[26]:FESystem<3>[FE_RaviartThomas<3>(1)-FE_DGQ<3>(1)] +DEAL::fe_data[28]:FESystem<3>[FE_RaviartThomas<3>(1)-FE_DGQ<3>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -3090,7 +3282,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 40 41 42 43 DEAL:: -DEAL::fe_data[27]:FE_Nedelec<3>(0) +DEAL::fe_data[29]:FE_Nedelec<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -3125,7 +3317,7 @@ DEAL::support on face 3: 0 1 3 4 5 7 10 11 DEAL::support on face 4: 0 1 2 3 8 9 10 11 DEAL::support on face 5: 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[28]:FE_Nedelec<3>(1) +DEAL::fe_data[30]:FE_Nedelec<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -3160,7 +3352,7 @@ DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 14 15 20 21 22 23 26 27 30 31 36 DEAL::support on face 4: 0 1 2 3 4 5 6 7 16 17 18 19 20 21 22 23 24 25 28 29 34 35 38 39 40 41 42 43 DEAL::support on face 5: 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 28 29 34 35 38 39 44 45 46 47 DEAL:: -DEAL::fe_data[29]:FE_RaviartThomasNodal<3>(0) +DEAL::fe_data[31]:FE_RaviartThomasNodal<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=1 @@ -3195,7 +3387,7 @@ DEAL::support on face 3: 0 1 3 4 5 DEAL::support on face 4: 0 1 2 3 4 DEAL::support on face 5: 0 1 2 3 5 DEAL:: -DEAL::fe_data[30]:FE_RaviartThomasNodal<3>(1) +DEAL::fe_data[32]:FE_RaviartThomasNodal<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -3230,7 +3422,7 @@ DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL::support on face 5: 0 1 2 3 4 5 6 7 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL:: -DEAL::fe_data[31]:FESystem<3>[FE_RaviartThomasNodal<3>(1)-FE_DGQ<3>(1)] +DEAL::fe_data[33]:FESystem<3>[FE_RaviartThomasNodal<3>(1)-FE_DGQ<3>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -3265,7 +3457,7 @@ DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 DEAL::support on face 5: 0 1 2 3 4 5 6 7 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 40 41 42 43 DEAL:: -DEAL::fe_data[32]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2] +DEAL::fe_data[34]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -3300,7 +3492,7 @@ DEAL::support on face 3: 8 9 10 11 12 13 14 15 24 25 26 27 28 29 30 31 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 DEAL::support on face 5: 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 DEAL:: -DEAL::fe_data[33]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]] +DEAL::fe_data[35]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -3335,7 +3527,7 @@ DEAL::support on face 3: 4 5 6 7 12 13 14 15 18 19 22 23 26 27 30 31 DEAL::support on face 4: 0 1 2 3 4 5 6 7 16 17 18 19 24 25 26 27 DEAL::support on face 5: 8 9 10 11 12 13 14 15 20 21 22 23 28 29 30 31 DEAL:: -DEAL::fe_data[34]:FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]] +DEAL::fe_data[36]:FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=3 @@ -3370,7 +3562,7 @@ DEAL::support on face 3: 8 9 10 11 12 13 14 15 24 25 26 27 28 29 30 31 41 42 43 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 32 33 34 35 36 37 38 39 40 41 42 43 80 81 82 89 90 91 92 93 94 95 96 97 116 117 118 119 120 121 122 123 124 DEAL::support on face 5: 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 44 45 46 47 48 49 50 51 52 53 54 55 83 84 85 107 108 109 110 111 112 113 114 115 134 135 136 137 138 139 140 141 142 DEAL:: -DEAL::fe_data[35]:FESystem<3>[FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2]^2-FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]]-FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]]^2] +DEAL::fe_data[37]:FESystem<3>[FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2]^2-FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]]-FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]]^2] DEAL::dofs_per_vertex=18 DEAL::dofs_per_line=6 DEAL::dofs_per_quad=6 -- 2.39.5