]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Test on non local dofs fe_data.
authorZhaowei Liu <zhaowei.liu@glasgow.ac.uk>
Mon, 27 Jul 2020 06:51:41 +0000 (08:51 +0200)
committerLuca Heltai <luca.heltai@sissa.it>
Thu, 30 Jul 2020 05:26:57 +0000 (07:26 +0200)
tests/non_local/fe_data_non_local_test_01.cc [new file with mode: 0644]
tests/non_local/fe_data_non_local_test_01.output [new file with mode: 0644]

diff --git a/tests/non_local/fe_data_non_local_test_01.cc b/tests/non_local/fe_data_non_local_test_01.cc
new file mode 100644 (file)
index 0000000..1f30d50
--- /dev/null
@@ -0,0 +1,66 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2020 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+// This is the test function for FiniteElementData with non-local dofs, where
+// the size of dofs_per_object is equal to dim+2. The test function output the
+// total number of dofs per cell and the number of non-local dofs per cell.
+
+
+#include <deal.II/fe/fe_base.h>
+
+#include <iostream>
+
+#include "../tests.h"
+
+template <int dim>
+void
+test(const std::vector<std::vector<unsigned int>> &dof_per_object_vecs)
+{
+  for (auto dof_per_object : dof_per_object_vecs)
+    {
+      FiniteElementData<dim> fe_data(dof_per_object,
+                                     1,
+                                     1,
+                                     FiniteElementData<dim>::H2);
+      deallog << "total number of dofs per cell is "
+              << fe_data.n_dofs_per_cell()
+              << " and number of non-local dofs per cell is "
+              << fe_data.n_non_local_dofs_per_cell() << std::endl;
+    }
+}
+
+int
+main()
+{
+  initlog();
+
+  deallog << "2d elements: " << std::endl;
+
+  std::vector<unsigned int> dof_per_object1{0, 0, 0, 4};
+  std::vector<unsigned int> dof_per_object2{0, 0, 4, 10};
+  std::vector<unsigned int> dof_per_object3{0, 0, 0, 4};
+  std::vector<unsigned int> dof_per_object4{1, 1, 4, 4};
+
+  test<2>({dof_per_object1, dof_per_object2, dof_per_object3, dof_per_object4});
+
+  deallog << "3d elements: " << std::endl;
+
+  std::vector<unsigned int> dof_per_object5{0, 0, 0, 0, 8};
+  std::vector<unsigned int> dof_per_object6{0, 0, 4, 1, 10};
+  std::vector<unsigned int> dof_per_object7{0, 0, 0, 2, 4};
+  std::vector<unsigned int> dof_per_object8{1, 1, 4, 3, 4};
+
+  test<3>({dof_per_object5, dof_per_object6, dof_per_object7, dof_per_object8});
+}
\ No newline at end of file
diff --git a/tests/non_local/fe_data_non_local_test_01.output b/tests/non_local/fe_data_non_local_test_01.output
new file mode 100644 (file)
index 0000000..861c116
--- /dev/null
@@ -0,0 +1,11 @@
+
+DEAL::2d elements: 
+DEAL::total number of dofs per cell is 4 and number of non-local dofs per cell is 4
+DEAL::total number of dofs per cell is 14 and number of non-local dofs per cell is 10
+DEAL::total number of dofs per cell is 4 and number of non-local dofs per cell is 4
+DEAL::total number of dofs per cell is 16 and number of non-local dofs per cell is 4
+DEAL::3d elements: 
+DEAL::total number of dofs per cell is 8 and number of non-local dofs per cell is 8
+DEAL::total number of dofs per cell is 35 and number of non-local dofs per cell is 10
+DEAL::total number of dofs per cell is 6 and number of non-local dofs per cell is 4
+DEAL::total number of dofs per cell is 51 and number of non-local dofs per cell is 4

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.