]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Move test that now works.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 10 Aug 2011 22:57:51 +0000 (22:57 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 10 Aug 2011 22:57:51 +0000 (22:57 +0000)
git-svn-id: https://svn.dealii.org/trunk@24045 0785d39b-7218-0410-832d-ea1e28bc413d

tests/deal.II/no_flux_08.cc [moved from tests/fail/half_hyper_normals.cc with 79% similarity]
tests/deal.II/no_flux_08/cmp/generic [new file with mode: 0644]

similarity index 79%
rename from tests/fail/half_hyper_normals.cc
rename to tests/deal.II/no_flux_08.cc
index 5cd923b5a0bb82e9f5b62f7643991aba1e5ad2d5..21cdd93ac6fc1649ec7e11e6182531803c5aff80 100644 (file)
@@ -1,8 +1,8 @@
 //----------------------------  vectors_rhs_hp_02.cc  ---------------------------
 //    $Id: vectors_rhs_hp_02.cc 23710 2011-05-17 04:50:10Z bangerth $
-//    Version: $Name$ 
+//    Version: $Name$
 //
-//    Copyright (C) 2000, 2001, 2003, 2004, 2006, 2007 by the deal.II authors
+//    Copyright (C) 2000, 2001, 2003, 2004, 2006, 2007, 2011 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
 //----------------------------  vectors_rhs_hp_02.cc  ---------------------------
 
 
-// some normals seem to be wrong when using  GridGenerator::half_hyper_shell() in 2d:
+// we used to get this crash:
+//
 // An error occurred in line <4646> of file </w/heister/deal-trunk/deal.II/include/deal.II/numerics/vectors.templates.h> in function
 //     static void dealii::VectorTools::compute_no_normal_flux_constraints(const DH<dim, spacedim>&, unsigned int, const std::set<unsigned char>&, dealii::ConstraintMatrix&, const dealii::Mapping<dim, spacedim>&) [with int dim = 3, DH = dealii::DoFHandler, int spacedim = 3]
-// The violated condition was: 
+// The violated condition was:
 //     std::fabs (tangent.norm()-1) < 1e-12
 // The name and call sequence of the exception was:
 //     ExcInternalError()
+//
 // quarter_hyper_shell works, even though it is a very similar mesh.
-
+//
+// this was fixed with r24044 together with the no_flux_07 test that
+// reduces it to its essence
 
 #include "../tests.h"
 
@@ -37,39 +41,34 @@ template <int dim>
 void
 check ()
 {
-  Triangulation<dim> tr; 
+  Triangulation<dim> tr;
    GridGenerator::half_hyper_shell (tr,
                              Point<dim>(),
                              0.5,
-                             1,0); 
-/*   GridGenerator::quarter_hyper_shell (tr,
-                             Point<dim>(),
-                             0.5,
-                             1,0); */   
-//  tr.refine_global (1);
+                             1,0);
 
   ConstraintMatrix cm;
   MappingQ<dim> mapping(1);
-  
+
   FESystem<dim> fe(FE_Q<dim>(1),dim);
   DoFHandler<dim> dofh(tr);
-  
+
   dofh.distribute_dofs (fe);
-  
+
   std::set<unsigned char> no_normal_flux_boundaries;
   no_normal_flux_boundaries.insert (0);
   VectorTools::compute_no_normal_flux_constraints (dofh, 0, no_normal_flux_boundaries, cm, mapping);
-  
-  deallog << "ok" << std::endl;
+
+  cm.print (deallog.get_file_stream ());
 }
 
 
 
 int main ()
 {
-  std::ofstream logfile ("half_hyper_normals/output");
+  std::ofstream logfile ("no_flux_08/output");
   logfile.precision (4);
-  logfile.setf(std::ios::fixed);  
+  logfile.setf(std::ios::fixed);
   deallog.attach(logfile);
   deallog.depth_console (0);
 
diff --git a/tests/deal.II/no_flux_08/cmp/generic b/tests/deal.II/no_flux_08/cmp/generic
new file mode 100644 (file)
index 0000000..397f529
--- /dev/null
@@ -0,0 +1,49 @@
+
+    0 = 0
+    1 = 0
+    3 2:  0.7265
+    4 = 0
+    5 = 0
+    7 6:  0.7265
+    8 9:  0.3249
+    10 11:  0.3249
+    12 13:  -0.3249
+    14 15:  -0.3249
+    17 16:  -0.7265
+    19 18:  -0.7265
+    20 = 0
+    21 = 0
+    22 = 0
+    23 = 0
+    0 = 0
+    1 2:  -1.0000
+    3 = 0
+    4 5:  -1.0000
+    6 7:  -0.5273
+    6 8:  0.5273
+    9 10:  0.5273
+    9 11:  0.5273
+    12 = 0
+    13 14:  -1.0000
+    15 = 0
+    16 17:  -1.0000
+    18 19:  -0.5273
+    18 20:  0.5273
+    21 22:  0.5273
+    21 23:  0.5273
+    24 = 0
+    25 26:  -1.0000
+    27 = 0
+    28 29:  -1.0000
+    30 31:  -0.5273
+    30 32:  -0.5273
+    33 34:  -0.5273
+    33 35:  -0.5273
+    36 = 0
+    37 38:  -1.0000
+    39 40:  0.5273
+    39 41:  -0.5273
+    42 = 0
+    43 44:  -1.0000
+    45 46:  0.5273
+    45 47:  -0.5273

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.