From: bangerth Date: Sun, 3 Sep 2006 20:22:52 +0000 (+0000) Subject: Make sure we understand the output and fix it X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=35b93027445c6894e56f3389ce8948a3a51882a2;p=dealii-svn.git Make sure we understand the output and fix it git-svn-id: https://svn.dealii.org/trunk@13820 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/bits/mesh_3d_18.cc b/tests/bits/mesh_3d_18.cc index 1aa8c3cf38..df40e020de 100644 --- a/tests/bits/mesh_3d_18.cc +++ b/tests/bits/mesh_3d_18.cc @@ -47,36 +47,49 @@ void test_with_wrong_face_orientation () ++cell; ++cell; - deallog << "Cell=" << cell << std::endl; - deallog << "Neighbor=" << cell->neighbor(3) << std::endl; + deallog << "cell=" << cell << std::endl; + deallog << "cell->neighbor(3)=" << cell->neighbor(3) << std::endl; + deallog << "cell->face_orientation(3)=" + << (cell->face_orientation(3) ? "true" : "false") + << std::endl; const Triangulation<3>::active_cell_iterator neighbor_child = cell->neighbor_child_on_subface (3, 1); - deallog << "Neighbor_child(3,1)=" << neighbor_child << std::endl; - for (unsigned int i=0; i<6; ++i) - { - deallog << "Neighbor_child neighbors="; - if (neighbor_child->at_boundary(i)) - deallog << "(at boundary)"; - else - deallog << neighbor_child->neighbor(i); - deallog << std::endl; - } - - deallog << "Neighbor_of_neighbor=" << cell->neighbor_of_neighbor(3) << std::endl; - - deallog << "Face=" << cell->face(3) << std::endl; - for (unsigned int i=0; i<4; ++i) - deallog << "Face_child=" << cell->face(3)->child(i) << std::endl; - - for (unsigned int i=0; i<4; ++i) - deallog << "Neighbor_face=" << cell->neighbor(3)->face(i) << std::endl; + deallog << "cell->neighbor_child_on_subface(3,1)=" << neighbor_child << std::endl; + deallog << "cell->neighbor_child_on_subface(3,1)->neighbor(5)=" + << neighbor_child->neighbor(5) << std::endl; + deallog << "cell->neighbor_child_on_subface(3,1)->face_orientation(5)=" + << (neighbor_child->face_orientation(5) ? "true" : "false") + << std::endl; + + deallog << "cell->face(3)=" << cell->face(3) << std::endl; for (unsigned int i=0; i<4; ++i) - deallog << "Neighbor_child_face=" << neighbor_child->face(i) << std::endl; + deallog << "cell->face(3)->child(" << i << ")=" + << cell->face(3)->child(i) << std::endl; - Assert (cell->face(3)->child(1) == + for (unsigned int i=0; i<6; ++i) + deallog << "cell->neighbor(3)->face(" << i << ")=" + << cell->neighbor(3)->face(i) << std::endl; + + for (unsigned int i=0; i<6; ++i) + deallog << "cell->neighbor_child_on_subface(3,1)->face(" << i << ")=" + << cell->neighbor_child_on_subface(3,1)->face(i) << std::endl; + + // The following assertion was originally + // in make_hanging_node_constraints with an + // equality sign and triggered for the mesh + // and cell here at hand. however, if one + // carefully reads the comment for + // CellAccessor::neighbor_child_on_subface, + // one realizes that it isn't true if + // cell->face_orientation()==false for the + // face we are presently considering. the + // original assertion therefore + // failed. let's make sure we get the + // status we expect. + Assert (cell->face(3)->child(1) != neighbor_child->face(cell->neighbor_of_neighbor(3)), ExcInternalError()); } diff --git a/tests/bits/mesh_3d_18/cmp/generic b/tests/bits/mesh_3d_18/cmp/generic index 57873fae27..7a8a5a89aa 100644 --- a/tests/bits/mesh_3d_18/cmp/generic +++ b/tests/bits/mesh_3d_18/cmp/generic @@ -1,274 +1,25 @@ -DEAL::n_dofs=14 -DEAL:cg::Starting value 3.59 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_0, rel. error=0 -DEAL:cg::Starting value 3.77 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0872 -DEAL:cg::Starting value 4.97 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0871 -DEAL::n_dofs=14 -DEAL:cg::Starting value 3.35 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_0, rel. error=0 -DEAL:cg::Starting value 3.72 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_1, rel. error=0.101 -DEAL:cg::Starting value 4.82 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0883 -DEAL::n_dofs=14 -DEAL:cg::Starting value 3.35 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_0, rel. error=0 -DEAL:cg::Starting value 3.47 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0745 -DEAL:cg::Starting value 4.61 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0716 -DEAL::n_dofs=14 -DEAL:cg::Starting value 3.35 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_0, rel. error=0 -DEAL:cg::Starting value 3.26 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0697 -DEAL:cg::Starting value 4.32 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0632 -DEAL::n_dofs=14 -DEAL:cg::Starting value 3.35 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_0, rel. error=0 -DEAL:cg::Starting value 3.62 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0844 -DEAL:cg::Starting value 4.76 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0797 -DEAL::n_dofs=14 -DEAL:cg::Starting value 3.35 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_0, rel. error=0 -DEAL:cg::Starting value 3.68 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0917 -DEAL:cg::Starting value 4.79 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0834 -DEAL::n_dofs=14 -DEAL:cg::Starting value 3.35 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_0, rel. error=0 -DEAL:cg::Starting value 3.37 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0716 -DEAL:cg::Starting value 4.48 -DEAL:cg::Convergence step 1 value 0 -DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0670 -DEAL::n_dofs=112 -DEAL:cg::Starting value 1.32 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_0, rel. error=0 -DEAL:cg::Starting value 1.41 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_1, rel. error=0 -DEAL:cg::Starting value 1.87 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00962 -DEAL:cg::Starting value 1.94 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00993 -DEAL::n_dofs=112 -DEAL:cg::Starting value 1.23 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_0, rel. error=0 -DEAL:cg::Starting value 1.39 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_1, rel. error=0 -DEAL:cg::Starting value 1.81 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00950 -DEAL:cg::Starting value 1.91 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_3, rel. error=0.0102 -DEAL::n_dofs=112 -DEAL:cg::Starting value 1.23 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_0, rel. error=0 -DEAL:cg::Starting value 1.29 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_1, rel. error=0 -DEAL:cg::Starting value 1.73 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00751 -DEAL:cg::Starting value 1.79 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00767 -DEAL::n_dofs=112 -DEAL:cg::Starting value 1.23 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_0, rel. error=0 -DEAL:cg::Starting value 1.21 -DEAL:cg::Convergence step 20 value 0 -DEAL::FE_DGQ<3>(1), P_1, rel. error=0 -DEAL:cg::Starting value 1.62 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00712 -DEAL:cg::Starting value 1.63 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00699 -DEAL::n_dofs=112 -DEAL:cg::Starting value 1.23 -DEAL:cg::Convergence step 20 value 0 -DEAL::FE_DGQ<3>(1), P_0, rel. error=0 -DEAL:cg::Starting value 1.35 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_1, rel. error=0 -DEAL:cg::Starting value 1.79 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00825 -DEAL:cg::Starting value 1.87 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00868 -DEAL::n_dofs=112 -DEAL:cg::Starting value 1.23 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_0, rel. error=0 -DEAL:cg::Starting value 1.37 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_1, rel. error=0 -DEAL:cg::Starting value 1.80 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00877 -DEAL:cg::Starting value 1.89 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00934 -DEAL::n_dofs=112 -DEAL:cg::Starting value 1.23 -DEAL:cg::Convergence step 20 value 0 -DEAL::FE_DGQ<3>(1), P_0, rel. error=0 -DEAL:cg::Starting value 1.25 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_1, rel. error=0 -DEAL:cg::Starting value 1.68 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00725 -DEAL:cg::Starting value 1.72 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00727 -DEAL::n_dofs=378 -DEAL:cg::Starting value 1.30 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_0, rel. error=0 -DEAL:cg::Starting value 1.38 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_1, rel. error=0 -DEAL:cg::Starting value 1.77 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_2, rel. error=0 -DEAL:cg::Starting value 1.83 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00166 -DEAL:cg::Starting value 1.96 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00164 -DEAL::n_dofs=378 -DEAL:cg::Starting value 1.21 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_0, rel. error=0 -DEAL:cg::Starting value 1.36 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_1, rel. error=0 -DEAL:cg::Starting value 1.72 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_2, rel. error=0 -DEAL:cg::Starting value 1.80 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00169 -DEAL:cg::Starting value 1.92 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00161 -DEAL::n_dofs=378 -DEAL:cg::Starting value 1.21 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_0, rel. error=0 -DEAL:cg::Starting value 1.27 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_1, rel. error=0 -DEAL:cg::Starting value 1.64 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_2, rel. error=0 -DEAL:cg::Starting value 1.70 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00126 -DEAL:cg::Starting value 1.82 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00124 -DEAL::n_dofs=378 -DEAL:cg::Starting value 1.21 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_0, rel. error=0 -DEAL:cg::Starting value 1.19 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_1, rel. error=0 -DEAL:cg::Starting value 1.53 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_2, rel. error=0 -DEAL:cg::Starting value 1.55 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00119 -DEAL:cg::Starting value 1.65 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00116 -DEAL::n_dofs=378 -DEAL:cg::Starting value 1.21 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_0, rel. error=0 -DEAL:cg::Starting value 1.32 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_1, rel. error=0 -DEAL:cg::Starting value 1.70 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_2, rel. error=0 -DEAL:cg::Starting value 1.77 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00142 -DEAL:cg::Starting value 1.89 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00138 -DEAL::n_dofs=378 -DEAL:cg::Starting value 1.21 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_0, rel. error=0 -DEAL:cg::Starting value 1.34 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_1, rel. error=0 -DEAL:cg::Starting value 1.71 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_2, rel. error=0 -DEAL:cg::Starting value 1.79 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00153 -DEAL:cg::Starting value 1.90 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00147 -DEAL::n_dofs=378 -DEAL:cg::Starting value 1.21 -DEAL:cg::Convergence step 18 value 0 -DEAL::FE_DGQ<3>(2), P_0, rel. error=0 -DEAL:cg::Starting value 1.23 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_1, rel. error=0 -DEAL:cg::Starting value 1.59 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_2, rel. error=0 -DEAL:cg::Starting value 1.63 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00121 -DEAL:cg::Starting value 1.75 -DEAL:cg::Convergence step 19 value 0 -DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00119 +DEAL::cell=0.3 +DEAL::cell->neighbor(3)=0.0 +DEAL::cell->face_orientation(3)=false +DEAL::cell->neighbor_child_on_subface(3,1)=1.6 +DEAL::cell->neighbor_child_on_subface(3,1)->neighbor(5)=0.3 +DEAL::cell->neighbor_child_on_subface(3,1)->face_orientation(5)=true +DEAL::cell->face(3)=4 +DEAL::cell->face(3)->child(0)=40 +DEAL::cell->face(3)->child(1)=41 +DEAL::cell->face(3)->child(2)=42 +DEAL::cell->face(3)->child(3)=43 +DEAL::cell->neighbor(3)->face(0)=1 +DEAL::cell->neighbor(3)->face(1)=3 +DEAL::cell->neighbor(3)->face(2)=0 +DEAL::cell->neighbor(3)->face(3)=5 +DEAL::cell->neighbor(3)->face(4)=2 +DEAL::cell->neighbor(3)->face(5)=4 +DEAL::cell->neighbor_child_on_subface(3,1)->face(0)=31 +DEAL::cell->neighbor_child_on_subface(3,1)->face(1)=51 +DEAL::cell->neighbor_child_on_subface(3,1)->face(2)=53 +DEAL::cell->neighbor_child_on_subface(3,1)->face(3)=45 +DEAL::cell->neighbor_child_on_subface(3,1)->face(4)=58 +DEAL::cell->neighbor_child_on_subface(3,1)->face(5)=42 +DEAL::OK