From 2be7803984295e86eb3762bb71f6476e8b8b28fc Mon Sep 17 00:00:00 2001 From: hartmann Date: Fri, 12 Mar 2004 10:01:51 +0000 Subject: [PATCH] Add tests for FE_Q_Hierarchical and FE_DGPMonomial. git-svn-id: https://svn.dealii.org/trunk@8719 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/fe/shapes.cc | 68 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 68 insertions(+) diff --git a/tests/fe/shapes.cc b/tests/fe/shapes.cc index a2f123f936..c6c65d39ed 100644 --- a/tests/fe/shapes.cc +++ b/tests/fe/shapes.cc @@ -12,7 +12,9 @@ #include #include #include +#include #include +#include #include #include #include @@ -416,6 +418,39 @@ void plot_FE_Q_shape_functions() } +template +void plot_FE_Q_Hierarchical_shape_functions() +{ + MappingQ1 m; + + FE_Q_Hierarchical q1(1); + plot_shape_functions(m, q1, "QHierarchical1"); + plot_face_shape_functions(m, q1, "QHierarchical1"); + test_compute_functions(m, q1, "QHierarchical1"); + + FE_Q_Hierarchical 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 q3(3); + plot_shape_functions(m, q3, "QHierarchical3"); + plot_face_shape_functions(m, q3, "QHierarchical3"); + test_compute_functions(m, q3, "QHierarchical3"); + + FE_Q_Hierarchical q4(4); + plot_shape_functions(m, q4, "QHierarchical4"); + plot_face_shape_functions(m, q4, "QHierarchical4"); + test_compute_functions(m, q4, "QHierarchical4"); + } +} + + + template void plot_FE_DGQ_shape_functions() { @@ -496,6 +531,31 @@ void plot_FE_DGP_shape_functions() } +template +void plot_FE_DGPMonomial_shape_functions() +{ + MappingQ1 m; + + FE_DGPMonomial p1(1); + plot_shape_functions(m, p1, "DGPMonomial1"); + plot_face_shape_functions(m, p1, "DGPMonomial1"); + test_compute_functions(m, p1, "DGPMonomial1"); + + FE_DGPMonomial 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 p3(3); + plot_shape_functions(m, p3, "DGPMonomial3"); + plot_face_shape_functions(m, p3, "DGPMonomial3"); + test_compute_functions(m, p3, "DGPMonomial3"); + } +} + + template 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>(); -- 2.39.5