]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Make sure we understand the output and fix it
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 3 Sep 2006 20:22:52 +0000 (20:22 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 3 Sep 2006 20:22:52 +0000 (20:22 +0000)
git-svn-id: https://svn.dealii.org/trunk@13820 0785d39b-7218-0410-832d-ea1e28bc413d

tests/bits/mesh_3d_18.cc
tests/bits/mesh_3d_18/cmp/generic

index 1aa8c3cf382845df6cb2a5fa0b753e69947974a1..df40e020dea4d3c36a6be09a6194150b5c0f6757 100644 (file)
@@ -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());
 }
index 57873fae277e5f6b1067046a9cf1e0e18d6c707b..7a8a5a89aa511a953c7b546092af5d9884dd2075 100644 (file)
 
-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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.