]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Add another test
authorkayser-herold <kayser-herold@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 27 Jul 2006 22:04:07 +0000 (22:04 +0000)
committerkayser-herold <kayser-herold@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 27 Jul 2006 22:04:07 +0000 (22:04 +0000)
git-svn-id: https://svn.dealii.org/trunk@13463 0785d39b-7218-0410-832d-ea1e28bc413d

tests/hp/hp_quad_dof_identities_q.cc [new file with mode: 0644]
tests/hp/hp_quad_dof_identities_q/cmp/generic [new file with mode: 0644]

diff --git a/tests/hp/hp_quad_dof_identities_q.cc b/tests/hp/hp_quad_dof_identities_q.cc
new file mode 100644 (file)
index 0000000..19f3dfd
--- /dev/null
@@ -0,0 +1,73 @@
+//----------------------------  hp_quad_dof_identities_q.cc  ---------------------------
+//    $Id: hp_quad_dof_identities_q.cc 12464 2006-02-23 01:13:17Z wolf $
+//    Version: $Name$ 
+//
+//    Copyright (C) 2005, 2006 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.
+//
+//----------------------------  hp_quad_dof_identities_q.cc  ---------------------------
+
+
+// check FE_Q::hp_quad_dof_identities
+
+
+#include <base/logstream.h>
+#include <fe/fe_collection.h>
+#include <fe/fe_q.h>
+
+#include <fstream>
+
+
+template <int dim>
+void test ()
+{
+  hp::FECollection<dim> fe_collection;
+  for (unsigned int i=1; i<8-dim; ++i)
+    fe_collection.push_back (FE_Q<dim>(i));
+
+  for (unsigned int i=0; i<fe_collection.size(); ++i)
+    for (unsigned int j=0; j<fe_collection.size(); ++j)
+      {
+       const std::vector<std::pair<unsigned int, unsigned int> >
+         identities = fe_collection[i].hp_quad_dof_identities (fe_collection[j]);
+
+       deallog << "Identities for "
+               << fe_collection[i].get_name() << " and "
+               << fe_collection[j].get_name() << ": "
+               << identities.size()
+               << std::endl;
+       
+       for (unsigned int k=0; k<identities.size(); ++k)
+         {
+           Assert (identities[k].first < fe_collection[i].dofs_per_quad,
+                   ExcInternalError());
+           Assert (identities[k].second < fe_collection[j].dofs_per_quad,
+                   ExcInternalError());
+           
+           deallog << identities[k].first << ' '
+                   << identities[k].second
+                   << std::endl;
+         }
+      }
+}
+
+
+
+int main ()
+{
+  std::ofstream logfile("hp_quad_dof_identities_q/output");
+  logfile.precision(2);
+  
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+  
+  test<2> ();
+  test<3> ();
+  
+  deallog << "OK" << std::endl;
+}
diff --git a/tests/hp/hp_quad_dof_identities_q/cmp/generic b/tests/hp/hp_quad_dof_identities_q/cmp/generic
new file mode 100644 (file)
index 0000000..2e68aec
--- /dev/null
@@ -0,0 +1,91 @@
+
+DEAL::Identities for FE_Q<2>(1) and FE_Q<2>(1): 0
+DEAL::Identities for FE_Q<2>(1) and FE_Q<2>(2): 0
+DEAL::Identities for FE_Q<2>(1) and FE_Q<2>(3): 0
+DEAL::Identities for FE_Q<2>(1) and FE_Q<2>(4): 0
+DEAL::Identities for FE_Q<2>(1) and FE_Q<2>(5): 0
+DEAL::Identities for FE_Q<2>(2) and FE_Q<2>(1): 0
+DEAL::Identities for FE_Q<2>(2) and FE_Q<2>(2): 1
+DEAL::0 0
+DEAL::Identities for FE_Q<2>(2) and FE_Q<2>(3): 0
+DEAL::Identities for FE_Q<2>(2) and FE_Q<2>(4): 1
+DEAL::0 4
+DEAL::Identities for FE_Q<2>(2) and FE_Q<2>(5): 0
+DEAL::Identities for FE_Q<2>(3) and FE_Q<2>(1): 0
+DEAL::Identities for FE_Q<2>(3) and FE_Q<2>(2): 0
+DEAL::Identities for FE_Q<2>(3) and FE_Q<2>(3): 4
+DEAL::0 0
+DEAL::1 1
+DEAL::2 2
+DEAL::3 3
+DEAL::Identities for FE_Q<2>(3) and FE_Q<2>(4): 0
+DEAL::Identities for FE_Q<2>(3) and FE_Q<2>(5): 0
+DEAL::Identities for FE_Q<2>(4) and FE_Q<2>(1): 0
+DEAL::Identities for FE_Q<2>(4) and FE_Q<2>(2): 1
+DEAL::4 0
+DEAL::Identities for FE_Q<2>(4) and FE_Q<2>(3): 0
+DEAL::Identities for FE_Q<2>(4) and FE_Q<2>(4): 9
+DEAL::0 0
+DEAL::1 1
+DEAL::2 2
+DEAL::3 3
+DEAL::4 4
+DEAL::5 5
+DEAL::6 6
+DEAL::7 7
+DEAL::8 8
+DEAL::Identities for FE_Q<2>(4) and FE_Q<2>(5): 0
+DEAL::Identities for FE_Q<2>(5) and FE_Q<2>(1): 0
+DEAL::Identities for FE_Q<2>(5) and FE_Q<2>(2): 0
+DEAL::Identities for FE_Q<2>(5) and FE_Q<2>(3): 0
+DEAL::Identities for FE_Q<2>(5) and FE_Q<2>(4): 0
+DEAL::Identities for FE_Q<2>(5) and FE_Q<2>(5): 16
+DEAL::0 0
+DEAL::1 1
+DEAL::2 2
+DEAL::3 3
+DEAL::4 4
+DEAL::5 5
+DEAL::6 6
+DEAL::7 7
+DEAL::8 8
+DEAL::9 9
+DEAL::10 10
+DEAL::11 11
+DEAL::12 12
+DEAL::13 13
+DEAL::14 14
+DEAL::15 15
+DEAL::Identities for FE_Q<3>(1) and FE_Q<3>(1): 0
+DEAL::Identities for FE_Q<3>(1) and FE_Q<3>(2): 0
+DEAL::Identities for FE_Q<3>(1) and FE_Q<3>(3): 0
+DEAL::Identities for FE_Q<3>(1) and FE_Q<3>(4): 0
+DEAL::Identities for FE_Q<3>(2) and FE_Q<3>(1): 0
+DEAL::Identities for FE_Q<3>(2) and FE_Q<3>(2): 1
+DEAL::0 0
+DEAL::Identities for FE_Q<3>(2) and FE_Q<3>(3): 0
+DEAL::Identities for FE_Q<3>(2) and FE_Q<3>(4): 1
+DEAL::0 4
+DEAL::Identities for FE_Q<3>(3) and FE_Q<3>(1): 0
+DEAL::Identities for FE_Q<3>(3) and FE_Q<3>(2): 0
+DEAL::Identities for FE_Q<3>(3) and FE_Q<3>(3): 4
+DEAL::0 0
+DEAL::1 1
+DEAL::2 2
+DEAL::3 3
+DEAL::Identities for FE_Q<3>(3) and FE_Q<3>(4): 0
+DEAL::Identities for FE_Q<3>(4) and FE_Q<3>(1): 0
+DEAL::Identities for FE_Q<3>(4) and FE_Q<3>(2): 1
+DEAL::4 0
+DEAL::Identities for FE_Q<3>(4) and FE_Q<3>(3): 0
+DEAL::Identities for FE_Q<3>(4) and FE_Q<3>(4): 9
+DEAL::0 0
+DEAL::1 1
+DEAL::2 2
+DEAL::3 3
+DEAL::4 4
+DEAL::5 5
+DEAL::6 6
+DEAL::7 7
+DEAL::8 8
+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.