]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Add tests for FE_Q_Hierarchical and FE_DGPMonomial.
authorhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 12 Mar 2004 10:01:51 +0000 (10:01 +0000)
committerhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 12 Mar 2004 10:01:51 +0000 (10:01 +0000)
git-svn-id: https://svn.dealii.org/trunk@8719 0785d39b-7218-0410-832d-ea1e28bc413d

tests/fe/shapes.cc

index a2f123f9369c9a27b8dc4363a28a5bfeff4d72e6..c6c65d39ed56eafb671974d78832ed71da418f1b 100644 (file)
@@ -12,7 +12,9 @@
 #include <dofs/dof_accessor.h>
 #include <grid/grid_generator.h>
 #include <fe/fe_q.h>
+#include <fe/fe_q_hierarchical.h>
 #include <fe/fe_dgp.h>
+#include <fe/fe_dgp_monomial.h>
 #include <fe/fe_dgp_nonparametric.h>
 #include <fe/fe_dgq.h>
 #include <fe/fe_nedelec.h>
@@ -416,6 +418,39 @@ void plot_FE_Q_shape_functions()
 }
 
 
+template<int dim>
+void plot_FE_Q_Hierarchical_shape_functions()
+{
+  MappingQ1<dim> m;
+
+  FE_Q_Hierarchical<dim> q1(1);
+  plot_shape_functions(m, q1, "QHierarchical1");
+  plot_face_shape_functions(m, q1, "QHierarchical1");
+  test_compute_functions(m, q1, "QHierarchical1");
+
+  FE_Q_Hierarchical<dim> q2(2);
+  plot_shape_functions(m, q2, "QHierarchical2");
+  plot_face_shape_functions(m, q2, "QHierarchical2");
+  test_compute_functions(m, q2, "QHierarchical2");
+
+                                  // skip the following tests to
+                                  // reduce run-time
+  if (dim < 3)
+    {
+      FE_Q_Hierarchical<dim> q3(3);
+      plot_shape_functions(m, q3, "QHierarchical3");
+      plot_face_shape_functions(m, q3, "QHierarchical3");
+      test_compute_functions(m, q3, "QHierarchical3");
+
+      FE_Q_Hierarchical<dim> q4(4);
+      plot_shape_functions(m, q4, "QHierarchical4");
+      plot_face_shape_functions(m, q4, "QHierarchical4");
+      test_compute_functions(m, q4, "QHierarchical4");
+    }
+}
+
+
+
 template<int dim>
 void plot_FE_DGQ_shape_functions()
 {
@@ -496,6 +531,31 @@ void plot_FE_DGP_shape_functions()
 }
 
 
+template<int dim>
+void plot_FE_DGPMonomial_shape_functions()
+{
+  MappingQ1<dim> m;
+
+  FE_DGPMonomial<dim> p1(1);
+  plot_shape_functions(m, p1, "DGPMonomial1");
+  plot_face_shape_functions(m, p1, "DGPMonomial1");
+  test_compute_functions(m, p1, "DGPMonomial1");
+
+  FE_DGPMonomial<dim> p2(2);
+  plot_shape_functions(m, p2, "DGPMonomial2");
+  plot_face_shape_functions(m, p2, "DGPMonomial2");
+  test_compute_functions(m, p2, "DGPMonomial2");
+
+  if (dim<3)
+    {
+      FE_DGPMonomial<dim> p3(3);
+      plot_shape_functions(m, p3, "DGPMonomial3");
+      plot_face_shape_functions(m, p3, "DGPMonomial3");
+      test_compute_functions(m, p3, "DGPMonomial3");
+    }
+}
+
+
 template<int dim>
 void plot_FE_DGPNonparametric_shape_functions()
 {
@@ -601,10 +661,18 @@ main()
   plot_FE_Q_shape_functions<2>();
   plot_FE_Q_shape_functions<3>();
 
+  plot_FE_Q_Hierarchical_shape_functions<1>();
+  plot_FE_Q_Hierarchical_shape_functions<2>();
+  plot_FE_Q_Hierarchical_shape_functions<3>();
+
   plot_FE_DGP_shape_functions<1>();
   plot_FE_DGP_shape_functions<2>();
   plot_FE_DGP_shape_functions<3>();
   
+  plot_FE_DGPMonomial_shape_functions<1>();
+  plot_FE_DGPMonomial_shape_functions<2>();
+  plot_FE_DGPMonomial_shape_functions<3>();
+
   plot_FE_DGPNonparametric_shape_functions<1>();
   plot_FE_DGPNonparametric_shape_functions<2>();
   plot_FE_DGPNonparametric_shape_functions<3>();

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.