From: bangerth Date: Fri, 11 Aug 2006 17:40:23 +0000 (+0000) Subject: Plenty of new tests X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=05baf1ddaea5ffc19eb17e853074316549dc9814;p=dealii-svn.git Plenty of new tests git-svn-id: https://svn.dealii.org/trunk@13675 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/deal.II/project_common.cc b/tests/deal.II/project_common.cc index faebd59cf6..12615f8a55 100644 --- a/tests/deal.II/project_common.cc +++ b/tests/deal.II/project_common.cc @@ -110,7 +110,7 @@ void test_no_hanging_nodes (const FiniteElement &fe, ConstraintMatrix constraints; constraints.close (); - Vector interpolant (dof_handler.n_dofs()); + Vector projection (dof_handler.n_dofs()); Vector error (triangulation.n_active_cells()); for (unsigned int q=0; q<=p+2; ++q) { @@ -119,21 +119,21 @@ void test_no_hanging_nodes (const FiniteElement &fe, constraints, QGauss(p+2), F (q, fe.n_components()), - interpolant); + projection); // then compute the interpolation error VectorTools::integrate_difference (dof_handler, - interpolant, + projection, F (q, fe.n_components()), error, QGauss(std::max(p,q)+1), VectorTools::L2_norm); deallog << fe.get_name() << ", P_" << q - << ", rel. error=" << error.l2_norm() / interpolant.l2_norm() + << ", rel. error=" << error.l2_norm() / projection.l2_norm() << std::endl; if (q<=p) - Assert (error.l2_norm() < 1e-12*interpolant.l2_norm(), + Assert (error.l2_norm() <= 1e-12*projection.l2_norm(), ExcInternalError()); } } diff --git a/tests/deal.II/project_dgp_01.cc b/tests/deal.II/project_dgp_01.cc index fa9b0881d4..8fa4a9b23f 100644 --- a/tests/deal.II/project_dgp_01.cc +++ b/tests/deal.II/project_dgp_01.cc @@ -12,7 +12,7 @@ //---------------------------- project_dgp_01.cc --------------------------- -// check that VectorTools::project works for Q elements correctly on +// check that VectorTools::project works for DGP elements correctly on // a uniformly refined mesh for functions of degree q char logname[] = "project_dgp_01/output"; diff --git a/tests/deal.II/project_dgp_monomial_01.cc b/tests/deal.II/project_dgp_monomial_01.cc new file mode 100644 index 0000000000..1a0162987d --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_01.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_monomial_01.cc --------------------------- +// $Id: project_dgp_monomial_01.cc 12732 2006-03-28 23:15:45Z wolf $ +// Version: $Name$ +// +// Copyright (C) 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. +// +//---------------------------- project_dgp_monomial_01.cc --------------------------- + + +// check that VectorTools::project works for DGP_Monomial elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_dgp_monomial_01/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_no_hanging_nodes (FE_DGPMonomial(p), p); +} diff --git a/tests/deal.II/project_dgp_monomial_01/cmp/generic b/tests/deal.II/project_dgp_monomial_01/cmp/generic new file mode 100644 index 0000000000..8a3ee90e5b --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_01/cmp/generic @@ -0,0 +1,178 @@ + +DEAL::n_dofs=8 +DEAL:cg::Starting value 0.354 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<1>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<1>(0), P_1, rel. error=0.0123 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<1>(0), P_2, rel. error=0.0203 +DEAL::n_dofs=16 +DEAL:cg::Starting value 0.395 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.826 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGPMonomial<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.33 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGPMonomial<1>(1), P_2, rel. error=0.000402 +DEAL:cg::Starting value 1.89 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGPMonomial<1>(1), P_3, rel. error=0.000924 +DEAL::n_dofs=24 +DEAL:cg::Starting value 0.412 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.863 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.39 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.98 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<1>(2), P_3, rel. error=1.22e-05 +DEAL:cg::Starting value 2.64 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<1>(2), P_4, rel. error=3.54e-05 +DEAL::n_dofs=32 +DEAL:cg::Starting value 0.422 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.884 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.42 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 2.03 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 2.71 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_4, rel. error=3.68e-07 +DEAL:cg::Starting value 3.46 +DEAL:cg::Convergence step 7 value 0 +DEAL::FE_DGPMonomial<1>(3), P_5, rel. error=1.28e-06 +DEAL::n_dofs=40 +DEAL:cg::Starting value 0.428 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 0.897 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 1.44 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 2.07 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 2.76 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 3.52 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_5, rel. error=1.11e-08 +DEAL:cg::Starting value 4.33 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_6, rel. error=4.50e-08 +DEAL::n_dofs=64 +DEAL:cg::Starting value 0.375 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<2>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.767 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<2>(0), P_1, rel. error=0.00329 +DEAL:cg::Starting value 1.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<2>(0), P_2, rel. error=0.00558 +DEAL::n_dofs=192 +DEAL:cg::Starting value 0.459 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.942 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.47 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_2, rel. error=0.000112 +DEAL:cg::Starting value 2.05 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_3, rel. error=0.000266 +DEAL::n_dofs=384 +DEAL:cg::Starting value 0.501 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGPMonomial<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.03 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.61 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 2.25 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_3, rel. error=3.49e-06 +DEAL:cg::Starting value 2.93 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_4, rel. error=1.04e-05 +DEAL::n_dofs=640 +DEAL:cg::Starting value 0.526 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_DGPMonomial<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 1.08 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGPMonomial<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.70 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 2.36 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 3.09 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_4, rel. error=1.08e-07 +DEAL:cg::Starting value 3.86 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_5, rel. error=3.85e-07 +DEAL::n_dofs=512 +DEAL:cg::Starting value 0.265 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.539 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<3>(0), P_1, rel. error=0.000979 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<3>(0), P_2, rel. error=0.00168 +DEAL::n_dofs=2048 +DEAL:cg::Starting value 0.351 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.715 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.10 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<3>(1), P_2, rel. error=3.38e-05 +DEAL:cg::Starting value 1.52 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<3>(1), P_3, rel. error=8.15e-05 +DEAL::n_dofs=5120 +DEAL:cg::Starting value 0.400 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_DGPMonomial<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.816 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_DGPMonomial<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.26 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_DGPMonomial<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.74 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_DGPMonomial<3>(2), P_3, rel. error=1.07e-06 +DEAL:cg::Starting value 2.24 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_DGPMonomial<3>(2), P_4, rel. error=3.23e-06 diff --git a/tests/deal.II/project_dgp_nonparametric_01.cc b/tests/deal.II/project_dgp_nonparametric_01.cc new file mode 100644 index 0000000000..f1d0d4b3a8 --- /dev/null +++ b/tests/deal.II/project_dgp_nonparametric_01.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_nonparametric_01.cc --------------------------- +// $Id: project_dgp_nonparametric_01.cc 12732 2006-03-28 23:15:45Z wolf $ +// Version: $Name$ +// +// Copyright (C) 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. +// +//---------------------------- project_dgp_nonparametric_01.cc --------------------------- + + +// check that VectorTools::project works for DGPNonparametric elements +// correctly on a uniformly refined mesh for functions of degree q + +char logname[] = "project_dgp_nonparametric_01/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_no_hanging_nodes (FE_DGPNonparametric(p), p); +} diff --git a/tests/deal.II/project_nedelec_01.cc b/tests/deal.II/project_nedelec_01.cc new file mode 100644 index 0000000000..12682fbb0f --- /dev/null +++ b/tests/deal.II/project_nedelec_01.cc @@ -0,0 +1,31 @@ +//---------------------------- project_nedelec_01.cc --------------------------- +// $Id: project_nedelec_01.cc 12732 2006-03-28 23:15:45Z wolf $ +// Version: $Name$ +// +// Copyright (C) 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. +// +//---------------------------- project_nedelec_01.cc --------------------------- + + +// check that VectorTools::project works for Nedelec elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_nedelec_01/output"; + + +#include "project_common.cc" + + +template +void test () +{ + if (dim > 1) + // only p=1 implemented at present + for (unsigned int p=1; p<2; ++p) + test_no_hanging_nodes (FE_Nedelec(p), p); +} diff --git a/tests/deal.II/project_nedelec_01/cmp/generic b/tests/deal.II/project_nedelec_01/cmp/generic new file mode 100644 index 0000000000..6c31f57c16 --- /dev/null +++ b/tests/deal.II/project_nedelec_01/cmp/generic @@ -0,0 +1,148 @@ + +DEAL::n_dofs=9 +DEAL:cg::Starting value 0.342 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.709 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.12 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q<1>(1), P_2, rel. error=0.000336 +DEAL:cg::Starting value 1.57 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q<1>(1), P_3, rel. error=0.000737 +DEAL::n_dofs=17 +DEAL:cg::Starting value 0.262 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Q<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.544 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Q<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.869 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(2), P_3, rel. error=1.24e-05 +DEAL:cg::Starting value 1.63 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(2), P_4, rel. error=3.41e-05 +DEAL::n_dofs=25 +DEAL:cg::Starting value 0.206 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.429 +DEAL:cg::Convergence step 14 value 0 +DEAL::FE_Q<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.685 +DEAL:cg::Convergence step 14 value 0 +DEAL::FE_Q<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 0.974 +DEAL:cg::Convergence step 14 value 0 +DEAL::FE_Q<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 1.29 +DEAL:cg::Convergence step 14 value 0 +DEAL::FE_Q<1>(3), P_4, rel. error=1.77e-07 +DEAL:cg::Starting value 1.64 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(3), P_5, rel. error=6.01e-07 +DEAL::n_dofs=33 +DEAL:cg::Starting value 0.191 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 0.398 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 0.638 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 0.908 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 1.53 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<1>(4), P_5, rel. error=6.36e-09 +DEAL:cg::Starting value 1.88 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<1>(4), P_6, rel. error=2.50e-08 +DEAL::n_dofs=81 +DEAL:cg::Starting value 0.352 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.717 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.10 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(1), P_2, rel. error=8.80e-05 +DEAL:cg::Starting value 1.51 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(1), P_3, rel. error=0.000198 +DEAL::n_dofs=289 +DEAL:cg::Starting value 0.206 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.421 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.654 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 0.905 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_3, rel. error=2.41e-06 +DEAL:cg::Starting value 1.17 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_4, rel. error=6.79e-06 +DEAL::n_dofs=625 +DEAL:cg::Starting value 0.127 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Q<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.261 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Q<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.406 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 0.562 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 0.730 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_4, rel. error=2.90e-08 +DEAL:cg::Starting value 0.908 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_5, rel. error=1.00e-07 +DEAL::n_dofs=729 +DEAL:cg::Starting value 0.241 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.488 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.744 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=2.51e-05 +DEAL:cg::Starting value 1.01 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=5.72e-05 +DEAL::n_dofs=4913 +DEAL:cg::Starting value 0.108 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.219 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.337 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 0.461 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=5.06e-07 +DEAL:cg::Starting value 0.592 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=1.44e-06 diff --git a/tests/deal.II/project_q_hierarchical_01.cc b/tests/deal.II/project_q_hierarchical_01.cc new file mode 100644 index 0000000000..a703b9f7eb --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_01.cc @@ -0,0 +1,29 @@ +//---------------------------- project_q_hierarchical_01.cc --------------------------- +// $Id: project_q_hierarchical_01.cc 12732 2006-03-28 23:15:45Z wolf $ +// Version: $Name$ +// +// Copyright (C) 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. +// +//---------------------------- project_q_hierarchical_01.cc --------------------------- + + +// check that VectorTools::project works for QHierarchical elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_q_hierarchical_01/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<6-dim; ++p) + test_no_hanging_nodes (FE_Q_Hierarchical(p), p); +} diff --git a/tests/deal.II/project_q_hierarchical_01/cmp/generic b/tests/deal.II/project_q_hierarchical_01/cmp/generic new file mode 100644 index 0000000000..15a4b99a24 --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_01/cmp/generic @@ -0,0 +1,148 @@ + +DEAL::n_dofs=9 +DEAL:cg::Starting value 0.342 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q_Hierarchical<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.709 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q_Hierarchical<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.12 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q_Hierarchical<1>(1), P_2, rel. error=0.000336 +DEAL:cg::Starting value 1.57 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_Q_Hierarchical<1>(1), P_3, rel. error=0.000737 +DEAL::n_dofs=17 +DEAL:cg::Starting value 0.416 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q_Hierarchical<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.862 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q_Hierarchical<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.37 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q_Hierarchical<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.93 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q_Hierarchical<1>(2), P_3, rel. error=1.64e-05 +DEAL:cg::Starting value 2.53 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q_Hierarchical<1>(2), P_4, rel. error=4.45e-05 +DEAL::n_dofs=25 +DEAL:cg::Starting value 0.416 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.862 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.37 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 1.93 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Q_Hierarchical<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 2.53 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Q_Hierarchical<1>(3), P_4, rel. error=2.75e-07 +DEAL:cg::Starting value 3.18 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Q_Hierarchical<1>(3), P_5, rel. error=9.15e-07 +DEAL::n_dofs=33 +DEAL:cg::Starting value 0.503 +DEAL:cg::Convergence step 33 value 0 +DEAL::FE_Q_Hierarchical<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 1.04 +DEAL:cg::Convergence step 33 value 0 +DEAL::FE_Q_Hierarchical<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 1.66 +DEAL:cg::Convergence step 33 value 0 +DEAL::FE_Q_Hierarchical<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 2.35 +DEAL:cg::Convergence step 34 value 0 +DEAL::FE_Q_Hierarchical<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 3.09 +DEAL:cg::Convergence step 34 value 0 +DEAL::FE_Q_Hierarchical<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 3.89 +DEAL:cg::Convergence step 34 value 0 +DEAL::FE_Q_Hierarchical<1>(4), P_5, rel. error=1.10e-08 +DEAL:cg::Starting value 4.74 +DEAL:cg::Convergence step 34 value 0 +DEAL::FE_Q_Hierarchical<1>(4), P_6, rel. error=4.22e-08 +DEAL::n_dofs=81 +DEAL:cg::Starting value 0.352 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q_Hierarchical<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.717 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q_Hierarchical<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.10 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<2>(1), P_2, rel. error=8.80e-05 +DEAL:cg::Starting value 1.51 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<2>(1), P_3, rel. error=0.000198 +DEAL::n_dofs=289 +DEAL:cg::Starting value 0.518 +DEAL:cg::Convergence step 110 value 0 +DEAL::FE_Q_Hierarchical<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.06 +DEAL:cg::Convergence step 111 value 0 +DEAL::FE_Q_Hierarchical<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.64 +DEAL:cg::Convergence step 112 value 0 +DEAL::FE_Q_Hierarchical<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 2.25 +DEAL:cg::Convergence step 112 value 0 +DEAL::FE_Q_Hierarchical<2>(2), P_3, rel. error=4.42e-06 +DEAL:cg::Starting value 2.90 +DEAL:cg::Convergence step 114 value 0 +DEAL::FE_Q_Hierarchical<2>(2), P_4, rel. error=1.22e-05 +DEAL::n_dofs=625 +DEAL:cg::Starting value 0.518 +DEAL:cg::Convergence step 116 value 0 +DEAL::FE_Q_Hierarchical<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 1.06 +DEAL:cg::Convergence step 117 value 0 +DEAL::FE_Q_Hierarchical<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.64 +DEAL:cg::Convergence step 118 value 0 +DEAL::FE_Q_Hierarchical<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 2.25 +DEAL:cg::Convergence step 118 value 0 +DEAL::FE_Q_Hierarchical<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 2.90 +DEAL:cg::Convergence step 119 value 0 +DEAL::FE_Q_Hierarchical<2>(3), P_4, rel. error=7.58e-08 +DEAL:cg::Starting value 3.57 +DEAL:cg::Convergence step 123 value 0 +DEAL::FE_Q_Hierarchical<2>(3), P_5, rel. error=2.57e-07 +DEAL::n_dofs=729 +DEAL:cg::Starting value 0.241 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.488 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.744 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=2.51e-05 +DEAL:cg::Starting value 1.01 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=5.72e-05 +DEAL::n_dofs=4913 +DEAL:cg::Starting value 0.431 +DEAL:cg::Convergence step 393 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.874 +DEAL:cg::Convergence step 401 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.34 +DEAL:cg::Convergence step 407 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.82 +DEAL:cg::Convergence step 411 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=1.27e-06 +DEAL:cg::Starting value 2.32 +DEAL:cg::Convergence step 415 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=3.57e-06 diff --git a/tests/deal.II/project_q_system_01.cc b/tests/deal.II/project_q_system_01.cc new file mode 100644 index 0000000000..f2b5661726 --- /dev/null +++ b/tests/deal.II/project_q_system_01.cc @@ -0,0 +1,31 @@ +//---------------------------- project_q_system_01.cc --------------------------- +// $Id: project_q_system_01.cc 12732 2006-03-28 23:15:45Z wolf $ +// Version: $Name$ +// +// Copyright (C) 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. +// +//---------------------------- project_q_system_01.cc --------------------------- + + +// check that VectorTools::project works for FESystem(FE_Q) elements correctly +// on a uniformly refined mesh for functions of degree q + +char logname[] = "project_q_system_01/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<5-dim; ++p) + test_no_hanging_nodes (FESystem(FE_Q(p), 1, + FE_DGQ(p+1), 1), + p); +} diff --git a/tests/deal.II/project_q_system_01/cmp/generic b/tests/deal.II/project_q_system_01/cmp/generic new file mode 100644 index 0000000000..dd07c4e404 --- /dev/null +++ b/tests/deal.II/project_q_system_01/cmp/generic @@ -0,0 +1,91 @@ + +DEAL::n_dofs=33 +DEAL:cg::Starting value 0.606 +DEAL:cg::Convergence step 11 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 11 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 1.91 +DEAL:cg::Convergence step 11 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_2, rel. error=0.000109 +DEAL:cg::Starting value 2.62 +DEAL:cg::Convergence step 11 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_3, rel. error=0.000256 +DEAL::n_dofs=49 +DEAL:cg::Starting value 0.474 +DEAL:cg::Convergence step 13 value 0 +DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_0, rel. error=0 +DEAL:cg::Starting value 0.967 +DEAL:cg::Convergence step 13 value 0 +DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_1, rel. error=0 +DEAL:cg::Starting value 1.50 +DEAL:cg::Convergence step 13 value 0 +DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_2, rel. error=0 +DEAL:cg::Starting value 2.07 +DEAL:cg::Convergence step 14 value 0 +DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_3, rel. error=4.87e-06 +DEAL:cg::Starting value 2.67 +DEAL:cg::Convergence step 14 value 0 +DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_4, rel. error=1.41e-05 +DEAL::n_dofs=65 +DEAL:cg::Starting value 0.429 +DEAL:cg::Convergence step 14 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_0, rel. error=0 +DEAL:cg::Starting value 0.872 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_1, rel. error=0 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_2, rel. error=0 +DEAL:cg::Starting value 1.86 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_3, rel. error=0 +DEAL:cg::Starting value 2.39 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_4, rel. error=7.70e-08 +DEAL:cg::Starting value 2.96 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_5, rel. error=2.71e-07 +DEAL::n_dofs=657 +DEAL:cg::Starting value 0.470 +DEAL:cg::Convergence step 17 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 0.955 +DEAL:cg::Convergence step 17 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 1.46 +DEAL:cg::Convergence step 17 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_2, rel. error=2.06e-05 +DEAL:cg::Starting value 1.99 +DEAL:cg::Convergence step 17 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_3, rel. error=4.87e-05 +DEAL::n_dofs=1313 +DEAL:cg::Starting value 0.284 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_0, rel. error=0 +DEAL:cg::Starting value 0.576 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_1, rel. error=0 +DEAL:cg::Starting value 0.887 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_2, rel. error=0 +DEAL:cg::Starting value 1.22 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_3, rel. error=7.98e-07 +DEAL:cg::Starting value 1.56 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_4, rel. error=2.32e-06 +DEAL::n_dofs=14553 +DEAL:cg::Starting value 0.279 +DEAL:cg::Convergence step 22 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 0.564 +DEAL:cg::Convergence step 23 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 0.860 +DEAL:cg::Convergence step 23 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=3.98e-06 +DEAL:cg::Starting value 1.16 +DEAL:cg::Convergence step 23 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=9.45e-06