]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
projects / dealii-svn.git / commitdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | patch | inline | side by side (parent: c279175)
New test that currently fails.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 15 Dec 2010 14:57:09 +0000 (14:57 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 15 Dec 2010 14:57:09 +0000 (14:57 +0000)
git-svn-id: https://svn.dealii.org/trunk@22965 0785d39b-7218-0410-832d-ea1e28bc413d

tests/codim_one/normal_vectors_01.cc [new file with mode: 0644] patch | blob
tests/codim_one/normal_vectors_01/cmp/generic [new file with mode: 0644] patch | blob

diff --git a/tests/codim_one/normal_vectors_01.cc b/tests/codim_one/normal_vectors_01.cc
new file mode 100644 (file)
index 0000000..8d9a96b
--- /dev/null
+++ b/tests/codim_one/normal_vectors_01.cc
@@ -0,0 +1,103 @@
+//----------------------------  normal_vectors_01.cc  ---------------------------
+//    $Id: bem.cc 22693 2010-11-11 20:11:47Z kanschat $
+//    Version: $Name$
+//
+//    Copyright (C) 2010 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  normal_vectors_01.cc  ---------------------------
+
+/*
+  Asking for normal vectors in the codim-1 case led to aborts.
+*/
+
+#include "../tests.h"
+
+#include <base/quadrature_lib.h>
+
+#include <grid/tria.h>
+#include <grid/grid_generator.h>
+#include <grid/tria_boundary_lib.h>
+#include <grid/grid_out.h>
+#include <grid/grid_tools.h>
+
+#include <fe/fe_q.h>
+#include <fe/mapping_q.h>
+#include <fe/fe_values.h>
+#include <dofs/dof_handler.h>
+#include <dofs/dof_accessor.h>
+
+template <int dim>
+void test ()
+{
+  Triangulation<dim-1,dim> tria;
+  std::map<typename Triangulation<dim-1,dim>::cell_iterator,
+    typename Triangulation<dim,dim>::face_iterator> surface_to_volume_mapping;
+
+  HyperBallBoundary<dim> boundary_description;
+  Triangulation<dim> volume_mesh;
+  GridGenerator::half_hyper_ball(volume_mesh);
+  
+  volume_mesh.set_boundary (1, boundary_description);
+  volume_mesh.set_boundary (0, boundary_description);
+  volume_mesh.refine_global (1);
+  
+  static HyperBallBoundary<dim-1,dim> surface_description;
+  tria.set_boundary (1, surface_description);
+  tria.set_boundary (0, surface_description);
+  
+  std::set<unsigned char> boundary_ids;
+  boundary_ids.insert(0);
+  
+  GridTools::extract_boundary_mesh (volume_mesh, tria,
+                                   surface_to_volume_mapping,
+                                   boundary_ids);
+  
+  
+  tria.refine_global(4);
+                                  // END OF SURFACE CONSTRUCTION      
+  
+  FE_Q<dim-1,dim> fe (1);
+  DoFHandler<dim-1,dim> dh(tria);
+  MappingQ<dim-1,dim> mapping(1,true);
+  
+  dh.distribute_dofs (fe);
+  
+  QGauss<dim-1>  quadrature_formula(2);
+  QGauss<dim-2>  face_quadrature_formula(2);
+  
+  FEValues<dim-1,dim> fe_values (mapping,fe, QGauss<dim-1>(2), 
+                                update_values   | update_cell_normal_vectors |
+                                update_gradients |
+                                update_quadrature_points | update_JxW_values);
+  
+  
+                                  // Neumann BC
+  FEFaceValues<dim-1,dim> fe_face_values (mapping, fe, QGauss<dim-2>(2), 
+                                         update_values         | update_quadrature_points  |
+                                         update_normal_vectors | update_JxW_values);
+  
+  typename DoFHandler<dim-1,dim>::active_cell_iterator
+    cell = dh.begin_active();
+  
+  fe_values.reinit (cell);  
+  fe_face_values.reinit (cell,0);
+}
+
+
+
+int main ()
+{
+  std::ofstream logfile("normal_vectors_01/output");
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+
+  test<2> ();
+  test<3> ();
+  
+  return 0;
+}
diff --git a/tests/codim_one/normal_vectors_01/cmp/generic b/tests/codim_one/normal_vectors_01/cmp/generic
new file mode 100644 (file)
index 0000000..d8b5693
--- /dev/null
+++ b/tests/codim_one/normal_vectors_01/cmp/generic
@@ -0,0 +1,293 @@
+
+DEAL::Testing hyper_cube in dim: 3...
+DEAL::The last column should be 0
+DEAL::Vol faces                Surf cell               Distance
+DEAL::
+DEAL::-0.577350 -0.577350 -0.577350            -0.577350 -0.577350 -0.577350                   0
+DEAL::0.00000 -0.577350 -0.577350              0.00000 -0.577350 -0.577350                     0
+DEAL::-0.577350 0.00000 -0.577350              -0.577350 0.00000 -0.577350                     0
+DEAL::0.00000 0.00000 -0.577350                0.00000 0.00000 -0.577350                       0
+DEAL::
+DEAL::0.00000 -0.577350 -0.577350              0.00000 -0.577350 -0.577350                     0
+DEAL::0.577350 -0.577350 -0.577350             0.577350 -0.577350 -0.577350                    0
+DEAL::0.00000 0.00000 -0.577350                0.00000 0.00000 -0.577350                       0
+DEAL::0.577350 0.00000 -0.577350               0.577350 0.00000 -0.577350                      0
+DEAL::
+DEAL::-0.577350 0.00000 -0.577350              -0.577350 0.00000 -0.577350                     0
+DEAL::0.00000 0.00000 -0.577350                0.00000 0.00000 -0.577350                       0
+DEAL::-0.577350 0.577350 -0.577350             -0.577350 0.577350 -0.577350                    0
+DEAL::0.00000 0.577350 -0.577350               0.00000 0.577350 -0.577350                      0
+DEAL::
+DEAL::0.00000 0.00000 -0.577350                0.00000 0.00000 -0.577350                       0
+DEAL::0.577350 0.00000 -0.577350               0.577350 0.00000 -0.577350                      0
+DEAL::0.00000 0.577350 -0.577350               0.00000 0.577350 -0.577350                      0
+DEAL::0.577350 0.577350 -0.577350              0.577350 0.577350 -0.577350                     0
+DEAL::
+DEAL::0.577350 -0.577350 -0.577350             0.577350 -0.577350 -0.577350                    0
+DEAL::0.577350 -0.577350 0.00000               0.577350 -0.577350 0.00000                      0
+DEAL::0.577350 0.00000 -0.577350               0.577350 0.00000 -0.577350                      0
+DEAL::0.577350 0.00000 0.00000         0.577350 0.00000 0.00000                        0
+DEAL::
+DEAL::0.577350 -0.577350 0.00000               0.577350 -0.577350 0.00000                      0
+DEAL::0.577350 -0.577350 0.577350              0.577350 -0.577350 0.577350                     0
+DEAL::0.577350 0.00000 0.00000         0.577350 0.00000 0.00000                        0
+DEAL::0.577350 0.00000 0.577350                0.577350 0.00000 0.577350                       0
+DEAL::
+DEAL::0.577350 0.00000 -0.577350               0.577350 0.00000 -0.577350                      0
+DEAL::0.577350 0.00000 0.00000         0.577350 0.00000 0.00000                        0
+DEAL::0.577350 0.577350 -0.577350              0.577350 0.577350 -0.577350                     0
+DEAL::0.577350 0.577350 0.00000                0.577350 0.577350 0.00000                       0
+DEAL::
+DEAL::0.577350 0.00000 0.00000         0.577350 0.00000 0.00000                        0
+DEAL::0.577350 0.00000 0.577350                0.577350 0.00000 0.577350                       0
+DEAL::0.577350 0.577350 0.00000                0.577350 0.577350 0.00000                       0
+DEAL::0.577350 0.577350 0.577350               0.577350 0.577350 0.577350                      0
+DEAL::
+DEAL::-0.577350 -0.577350 0.577350             -0.577350 -0.577350 0.577350                    0
+DEAL::-0.577350 0.00000 0.577350               -0.577350 0.00000 0.577350                      0
+DEAL::0.00000 -0.577350 0.577350               0.00000 -0.577350 0.577350                      0
+DEAL::0.00000 0.00000 0.577350         0.00000 0.00000 0.577350                        0
+DEAL::
+DEAL::-0.577350 0.00000 0.577350               -0.577350 0.00000 0.577350                      0
+DEAL::-0.577350 0.577350 0.577350              -0.577350 0.577350 0.577350                     0
+DEAL::0.00000 0.00000 0.577350         0.00000 0.00000 0.577350                        0
+DEAL::0.00000 0.577350 0.577350                0.00000 0.577350 0.577350                       0
+DEAL::
+DEAL::0.00000 -0.577350 0.577350               0.00000 -0.577350 0.577350                      0
+DEAL::0.00000 0.00000 0.577350         0.00000 0.00000 0.577350                        0
+DEAL::0.577350 -0.577350 0.577350              0.577350 -0.577350 0.577350                     0
+DEAL::0.577350 0.00000 0.577350                0.577350 0.00000 0.577350                       0
+DEAL::
+DEAL::0.00000 0.00000 0.577350         0.00000 0.00000 0.577350                        0
+DEAL::0.00000 0.577350 0.577350                0.00000 0.577350 0.577350                       0
+DEAL::0.577350 0.00000 0.577350                0.577350 0.00000 0.577350                       0
+DEAL::0.577350 0.577350 0.577350               0.577350 0.577350 0.577350                      0
+DEAL::
+DEAL::-0.577350 -0.577350 -0.577350            -0.577350 -0.577350 -0.577350                   0
+DEAL::-0.577350 0.00000 -0.577350              -0.577350 0.00000 -0.577350                     0
+DEAL::-0.577350 -0.577350 0.00000              -0.577350 -0.577350 0.00000                     0
+DEAL::-0.577350 0.00000 0.00000                -0.577350 0.00000 0.00000                       0
+DEAL::
+DEAL::-0.577350 0.00000 -0.577350              -0.577350 0.00000 -0.577350                     0
+DEAL::-0.577350 0.577350 -0.577350             -0.577350 0.577350 -0.577350                    0
+DEAL::-0.577350 0.00000 0.00000                -0.577350 0.00000 0.00000                       0
+DEAL::-0.577350 0.577350 0.00000               -0.577350 0.577350 0.00000                      0
+DEAL::
+DEAL::-0.577350 -0.577350 0.00000              -0.577350 -0.577350 0.00000                     0
+DEAL::-0.577350 0.00000 0.00000                -0.577350 0.00000 0.00000                       0
+DEAL::-0.577350 -0.577350 0.577350             -0.577350 -0.577350 0.577350                    0
+DEAL::-0.577350 0.00000 0.577350               -0.577350 0.00000 0.577350                      0
+DEAL::
+DEAL::-0.577350 0.00000 0.00000                -0.577350 0.00000 0.00000                       0
+DEAL::-0.577350 0.577350 0.00000               -0.577350 0.577350 0.00000                      0
+DEAL::-0.577350 0.00000 0.577350               -0.577350 0.00000 0.577350                      0
+DEAL::-0.577350 0.577350 0.577350              -0.577350 0.577350 0.577350                     0
+DEAL::
+DEAL::-0.577350 -0.577350 -0.577350            -0.577350 -0.577350 -0.577350                   0
+DEAL::-0.577350 -0.577350 0.00000              -0.577350 -0.577350 0.00000                     0
+DEAL::0.00000 -0.577350 -0.577350              0.00000 -0.577350 -0.577350                     0
+DEAL::0.00000 -0.577350 0.00000                0.00000 -0.577350 0.00000                       0
+DEAL::
+DEAL::-0.577350 -0.577350 0.00000              -0.577350 -0.577350 0.00000                     0
+DEAL::-0.577350 -0.577350 0.577350             -0.577350 -0.577350 0.577350                    0
+DEAL::0.00000 -0.577350 0.00000                0.00000 -0.577350 0.00000                       0
+DEAL::0.00000 -0.577350 0.577350               0.00000 -0.577350 0.577350                      0
+DEAL::
+DEAL::0.00000 -0.577350 -0.577350              0.00000 -0.577350 -0.577350                     0
+DEAL::0.00000 -0.577350 0.00000                0.00000 -0.577350 0.00000                       0
+DEAL::0.577350 -0.577350 -0.577350             0.577350 -0.577350 -0.577350                    0
+DEAL::0.577350 -0.577350 0.00000               0.577350 -0.577350 0.00000                      0
+DEAL::
+DEAL::0.00000 -0.577350 0.00000                0.00000 -0.577350 0.00000                       0
+DEAL::0.00000 -0.577350 0.577350               0.00000 -0.577350 0.577350                      0
+DEAL::0.577350 -0.577350 0.00000               0.577350 -0.577350 0.00000                      0
+DEAL::0.577350 -0.577350 0.577350              0.577350 -0.577350 0.577350                     0
+DEAL::
+DEAL::-0.577350 0.577350 -0.577350             -0.577350 0.577350 -0.577350                    0
+DEAL::0.00000 0.577350 -0.577350               0.00000 0.577350 -0.577350                      0
+DEAL::-0.577350 0.577350 0.00000               -0.577350 0.577350 0.00000                      0
+DEAL::0.00000 0.577350 0.00000         0.00000 0.577350 0.00000                        0
+DEAL::
+DEAL::0.00000 0.577350 -0.577350               0.00000 0.577350 -0.577350                      0
+DEAL::0.577350 0.577350 -0.577350              0.577350 0.577350 -0.577350                     0
+DEAL::0.00000 0.577350 0.00000         0.00000 0.577350 0.00000                        0
+DEAL::0.577350 0.577350 0.00000                0.577350 0.577350 0.00000                       0
+DEAL::
+DEAL::-0.577350 0.577350 0.00000               -0.577350 0.577350 0.00000                      0
+DEAL::0.00000 0.577350 0.00000         0.00000 0.577350 0.00000                        0
+DEAL::-0.577350 0.577350 0.577350              -0.577350 0.577350 0.577350                     0
+DEAL::0.00000 0.577350 0.577350                0.00000 0.577350 0.577350                       0
+DEAL::
+DEAL::0.00000 0.577350 0.00000         0.00000 0.577350 0.00000                        0
+DEAL::0.577350 0.577350 0.00000                0.577350 0.577350 0.00000                       0
+DEAL::0.00000 0.577350 0.577350                0.00000 0.577350 0.577350                       0
+DEAL::0.577350 0.577350 0.577350               0.577350 0.577350 0.577350                      0
+DEAL::
+-0.577350 -0.577350 -0.577350 1 0
+0.00000 -0.577350 -0.577350 1 0
+0.00000 0.00000 -0.577350 1 0
+-0.577350 0.00000 -0.577350 1 0
+-0.577350 -0.577350 -0.577350 1 0
+
+
+0.00000 -0.577350 -0.577350 1 0
+0.577350 -0.577350 -0.577350 1 0
+0.577350 0.00000 -0.577350 1 0
+0.00000 0.00000 -0.577350 1 0
+0.00000 -0.577350 -0.577350 1 0
+
+
+-0.577350 0.00000 -0.577350 1 0
+0.00000 0.00000 -0.577350 1 0
+0.00000 0.577350 -0.577350 1 0
+-0.577350 0.577350 -0.577350 1 0
+-0.577350 0.00000 -0.577350 1 0
+
+
+0.00000 0.00000 -0.577350 1 0
+0.577350 0.00000 -0.577350 1 0
+0.577350 0.577350 -0.577350 1 0
+0.00000 0.577350 -0.577350 1 0
+0.00000 0.00000 -0.577350 1 0
+
+
+0.577350 -0.577350 -0.577350 1 0
+0.577350 -0.577350 0.00000 1 0
+0.577350 0.00000 0.00000 1 0
+0.577350 0.00000 -0.577350 1 0
+0.577350 -0.577350 -0.577350 1 0
+
+
+0.577350 -0.577350 0.00000 1 0
+0.577350 -0.577350 0.577350 1 0
+0.577350 0.00000 0.577350 1 0
+0.577350 0.00000 0.00000 1 0
+0.577350 -0.577350 0.00000 1 0
+
+
+0.577350 0.00000 -0.577350 1 0
+0.577350 0.00000 0.00000 1 0
+0.577350 0.577350 0.00000 1 0
+0.577350 0.577350 -0.577350 1 0
+0.577350 0.00000 -0.577350 1 0
+
+
+0.577350 0.00000 0.00000 1 0
+0.577350 0.00000 0.577350 1 0
+0.577350 0.577350 0.577350 1 0
+0.577350 0.577350 0.00000 1 0
+0.577350 0.00000 0.00000 1 0
+
+
+-0.577350 -0.577350 0.577350 1 0
+-0.577350 0.00000 0.577350 1 0
+0.00000 0.00000 0.577350 1 0
+0.00000 -0.577350 0.577350 1 0
+-0.577350 -0.577350 0.577350 1 0
+
+
+-0.577350 0.00000 0.577350 1 0
+-0.577350 0.577350 0.577350 1 0
+0.00000 0.577350 0.577350 1 0
+0.00000 0.00000 0.577350 1 0
+-0.577350 0.00000 0.577350 1 0
+
+
+0.00000 -0.577350 0.577350 1 0
+0.00000 0.00000 0.577350 1 0
+0.577350 0.00000 0.577350 1 0
+0.577350 -0.577350 0.577350 1 0
+0.00000 -0.577350 0.577350 1 0
+
+
+0.00000 0.00000 0.577350 1 0
+0.00000 0.577350 0.577350 1 0
+0.577350 0.577350 0.577350 1 0
+0.577350 0.00000 0.577350 1 0
+0.00000 0.00000 0.577350 1 0
+
+
+-0.577350 -0.577350 -0.577350 1 0
+-0.577350 0.00000 -0.577350 1 0
+-0.577350 0.00000 0.00000 1 0
+-0.577350 -0.577350 0.00000 1 0
+-0.577350 -0.577350 -0.577350 1 0
+
+
+-0.577350 0.00000 -0.577350 1 0
+-0.577350 0.577350 -0.577350 1 0
+-0.577350 0.577350 0.00000 1 0
+-0.577350 0.00000 0.00000 1 0
+-0.577350 0.00000 -0.577350 1 0
+
+
+-0.577350 -0.577350 0.00000 1 0
+-0.577350 0.00000 0.00000 1 0
+-0.577350 0.00000 0.577350 1 0
+-0.577350 -0.577350 0.577350 1 0
+-0.577350 -0.577350 0.00000 1 0
+
+
+-0.577350 0.00000 0.00000 1 0
+-0.577350 0.577350 0.00000 1 0
+-0.577350 0.577350 0.577350 1 0
+-0.577350 0.00000 0.577350 1 0
+-0.577350 0.00000 0.00000 1 0
+
+
+-0.577350 -0.577350 -0.577350 1 0
+-0.577350 -0.577350 0.00000 1 0
+0.00000 -0.577350 0.00000 1 0
+0.00000 -0.577350 -0.577350 1 0
+-0.577350 -0.577350 -0.577350 1 0
+
+
+-0.577350 -0.577350 0.00000 1 0
+-0.577350 -0.577350 0.577350 1 0
+0.00000 -0.577350 0.577350 1 0
+0.00000 -0.577350 0.00000 1 0
+-0.577350 -0.577350 0.00000 1 0
+
+
+0.00000 -0.577350 -0.577350 1 0
+0.00000 -0.577350 0.00000 1 0
+0.577350 -0.577350 0.00000 1 0
+0.577350 -0.577350 -0.577350 1 0
+0.00000 -0.577350 -0.577350 1 0
+
+
+0.00000 -0.577350 0.00000 1 0
+0.00000 -0.577350 0.577350 1 0
+0.577350 -0.577350 0.577350 1 0
+0.577350 -0.577350 0.00000 1 0
+0.00000 -0.577350 0.00000 1 0
+
+
+-0.577350 0.577350 -0.577350 1 0
+0.00000 0.577350 -0.577350 1 0
+0.00000 0.577350 0.00000 1 0
+-0.577350 0.577350 0.00000 1 0
+-0.577350 0.577350 -0.577350 1 0
+
+
+0.00000 0.577350 -0.577350 1 0
+0.577350 0.577350 -0.577350 1 0
+0.577350 0.577350 0.00000 1 0
+0.00000 0.577350 0.00000 1 0
+0.00000 0.577350 -0.577350 1 0
+
+
+-0.577350 0.577350 0.00000 1 0
+0.00000 0.577350 0.00000 1 0
+0.00000 0.577350 0.577350 1 0
+-0.577350 0.577350 0.577350 1 0
+-0.577350 0.577350 0.00000 1 0
+
+
+0.00000 0.577350 0.00000 1 0
+0.577350 0.577350 0.00000 1 0
+0.577350 0.577350 0.577350 1 0
+0.00000 0.577350 0.577350 1 0
+0.00000 0.577350 0.00000 1 0
+
+
Git mirror of the historical svn repository

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.