]> https://gitweb.dealii.org/ - dealii.git/commitdiff
New test.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sat, 10 Oct 2009 16:57:19 +0000 (16:57 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Sat, 10 Oct 2009 16:57:19 +0000 (16:57 +0000)
git-svn-id: https://svn.dealii.org/trunk@19806 0785d39b-7218-0410-832d-ea1e28bc413d

tests/bits/subdomain_ids_07.cc [new file with mode: 0644]
tests/bits/subdomain_ids_07/cmp/generic [new file with mode: 0644]

diff --git a/tests/bits/subdomain_ids_07.cc b/tests/bits/subdomain_ids_07.cc
new file mode 100644 (file)
index 0000000..9f6e9c5
--- /dev/null
@@ -0,0 +1,109 @@
+//----------------------------  subdomain_ids_07.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 2001, 2002, 2003, 2004, 2005, 2009 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.
+//
+//----------------------------  subdomain_ids_07.cc  ---------------------------
+
+// check DoFTools::dof_indices_with_subdomain_association
+
+
+#include "../tests.h"
+#include <base/logstream.h>
+#include <grid/tria.h>
+#include <grid/grid_generator.h>
+#include <grid/tria_iterator.h>
+#include <dofs/dof_handler.h>
+#include <dofs/dof_accessor.h>
+#include <dofs/dof_tools.h>
+#include <fe/fe_dgq.h>
+#include <fe/fe_q.h>
+#include <fe/fe_system.h>
+#include <dofs/dof_tools.h>
+
+#include <fstream>
+#include <algorithm>
+#include <cmath>
+
+
+std::ofstream logfile("subdomain_ids_07/output");
+
+
+template <int dim>
+void test ()
+{
+  deallog << dim << 'D' << std::endl;
+  Triangulation<dim> tria;
+  GridGenerator::hyper_cube(tria, -1, 1);
+  tria.refine_global (2);
+
+                                  // we now have a number of cells,
+                                  // flag them with some subdomain
+                                  // ids based on their position, in
+                                  // particular we take the quadrant
+                                  // (octant)
+  typename Triangulation<dim>::active_cell_iterator
+    cell = tria.begin_active (),
+    endc = tria.end ();
+  for (; cell!=endc; ++cell)
+    {
+      unsigned int subdomain = 0;
+      for (unsigned int d=0; d<dim; ++d)
+       if (cell->center()(d) > 0)
+         subdomain |= (1<<d);
+      Assert (subdomain < (1<<dim), ExcInternalError());
+
+      cell->set_subdomain_id (subdomain);
+    };
+
+                                  // distribute some degrees of freedom and
+                                  // output some information on them
+  FESystem<dim> fe(FE_Q<dim>(2),dim, FE_DGQ<dim>(1), 1);
+  DoFHandler<dim> dof_handler (tria);
+  dof_handler.distribute_dofs (fe);
+  deallog << dof_handler.n_dofs() << std::endl;
+
+  std::vector<unsigned int> subdomain_association (dof_handler.n_dofs());
+  DoFTools::get_subdomain_association (dof_handler,
+                                       subdomain_association);
+  for (unsigned int subdomain=0; subdomain<(1<<dim); ++subdomain)
+    {
+      const IndexSet index_set
+       = DoFTools::dof_indices_with_subdomain_association (dof_handler,
+                                                           subdomain);
+
+      deallog << "Index set is "
+             << (index_set.is_contiguous() ? "" : "not ")
+             << "contiguous."
+             << std::endl;
+
+      for (unsigned int i=0; i<dof_handler.n_dofs(); ++i)
+       Assert ((subdomain_association[i] == subdomain)
+               ==
+               (index_set.is_element(i) == true),
+               ExcInternalError());
+    }
+
+  deallog << "OK" << std::endl;
+}
+
+
+int main ()
+{
+  deallog << std::setprecision(4);
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  test<1> ();
+  test<2> ();
+  test<3> ();
+
+  return 0;
+}
diff --git a/tests/bits/subdomain_ids_07/cmp/generic b/tests/bits/subdomain_ids_07/cmp/generic
new file mode 100644 (file)
index 0000000..0810951
--- /dev/null
@@ -0,0 +1,24 @@
+
+DEAL::1D
+DEAL::17
+DEAL::Index set is contiguous.
+DEAL::Index set is contiguous.
+DEAL::OK
+DEAL::2D
+DEAL::226
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::OK
+DEAL::3D
+DEAL::2699
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+DEAL::Index set is not contiguous.
+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.