From 73ef7d436cafff19a6394db298a5502cda87e79a Mon Sep 17 00:00:00 2001 From: bangerth Date: Sun, 13 Aug 2006 03:47:11 +0000 Subject: [PATCH] Add new tests that check projection onto meshes with hanging node constraints and make sure that each element satisfies its promised approximation order there. git-svn-id: https://svn.dealii.org/trunk@13691 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/deal.II/project_abf_02.cc | 30 +++ tests/deal.II/project_common.cc | 85 +++++++-- tests/deal.II/project_dgp_02.cc | 29 +++ tests/deal.II/project_dgp_02/cmp/generic | 178 ++++++++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 178 ++++++++++++++++++ tests/deal.II/project_dgp_monomial_02.cc | 29 +++ .../project_dgp_monomial_02/cmp/generic | 178 ++++++++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 178 ++++++++++++++++++ tests/deal.II/project_dgp_nonparametric_02.cc | 29 +++ tests/deal.II/project_dgq_02.cc | 29 +++ tests/deal.II/project_dgq_02/cmp/generic | 178 ++++++++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 178 ++++++++++++++++++ tests/deal.II/project_nedelec_02.cc | 31 +++ tests/deal.II/project_nedelec_02/cmp/generic | 21 +++ tests/deal.II/project_q_02.cc | 29 +++ tests/deal.II/project_q_02/cmp/generic | 148 +++++++++++++++ .../project_q_hierarchical_01/cmp/generic | 4 +- tests/deal.II/project_q_hierarchical_02.cc | 29 +++ .../project_q_hierarchical_02/cmp/generic | 148 +++++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 148 +++++++++++++++ tests/deal.II/project_q_system_02.cc | 31 +++ tests/deal.II/project_q_system_02/cmp/generic | 91 +++++++++ tests/deal.II/project_rt_02.cc | 33 ++++ 23 files changed, 1997 insertions(+), 15 deletions(-) create mode 100644 tests/deal.II/project_abf_02.cc create mode 100644 tests/deal.II/project_dgp_02.cc create mode 100644 tests/deal.II/project_dgp_02/cmp/generic create mode 100644 tests/deal.II/project_dgp_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_dgp_monomial_02.cc create mode 100644 tests/deal.II/project_dgp_monomial_02/cmp/generic create mode 100644 tests/deal.II/project_dgp_monomial_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_dgp_nonparametric_02.cc create mode 100644 tests/deal.II/project_dgq_02.cc create mode 100644 tests/deal.II/project_dgq_02/cmp/generic create mode 100644 tests/deal.II/project_dgq_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_nedelec_02.cc create mode 100644 tests/deal.II/project_nedelec_02/cmp/generic create mode 100644 tests/deal.II/project_q_02.cc create mode 100644 tests/deal.II/project_q_02/cmp/generic create mode 100644 tests/deal.II/project_q_hierarchical_02.cc create mode 100644 tests/deal.II/project_q_hierarchical_02/cmp/generic create mode 100644 tests/deal.II/project_q_hierarchical_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_q_system_02.cc create mode 100644 tests/deal.II/project_q_system_02/cmp/generic create mode 100644 tests/deal.II/project_rt_02.cc diff --git a/tests/deal.II/project_abf_02.cc b/tests/deal.II/project_abf_02.cc new file mode 100644 index 0000000000..426af62491 --- /dev/null +++ b/tests/deal.II/project_abf_02.cc @@ -0,0 +1,30 @@ +//---------------------------- project_abf_02.cc --------------------------- +// $Id: project_abf_02.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_abf_02.cc --------------------------- + + +// check that VectorTools::project works for ABF elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_abf_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + if (dim > 1) + for (unsigned int p=0; p<6-dim; ++p) + test_with_hanging_nodes (FE_ABF(p), p); +} diff --git a/tests/deal.II/project_common.cc b/tests/deal.II/project_common.cc index 2b01bfac80..77be2216e4 100644 --- a/tests/deal.II/project_common.cc +++ b/tests/deal.II/project_common.cc @@ -93,26 +93,20 @@ class F : public Function -// check the given element of polynomial order p. the last parameter, if -// given, denotes a gap in convergence order; for example, the Nedelec element -// of polynomial degree p has normal components of degree p-1 and therefore -// can only represent polynomials of degree p-1 exactly. the gap is then 1. template -void test_no_hanging_nodes (const FiniteElement &fe, - const unsigned int p, - const unsigned int order_difference = 0) -{ - Triangulation triangulation; - GridGenerator::hyper_cube (triangulation); - triangulation.refine_global (3); - +void do_project (const Triangulation &triangulation, + const FiniteElement &fe, + const unsigned int p, + const unsigned int order_difference) +{ DoFHandler dof_handler(triangulation); dof_handler.distribute_dofs (fe); deallog << "n_dofs=" << dof_handler.n_dofs() << std::endl; - // there are no constraints here ConstraintMatrix constraints; + DoFTools::make_hanging_node_constraints (dof_handler, + constraints); constraints.close (); Vector projection (dof_handler.n_dofs()); @@ -125,6 +119,10 @@ void test_no_hanging_nodes (const FiniteElement &fe, QGauss(p+2), F (q, fe.n_components()), projection); + // just to make sure it doesn't get + // forgotten: handle hanging node + // constraints + constraints.distribute (projection); // then compute the interpolation error VectorTools::integrate_difference (dof_handler, @@ -145,6 +143,67 @@ void test_no_hanging_nodes (const FiniteElement &fe, +// check the given element of polynomial order p. the last parameter, if +// given, denotes a gap in convergence order; for example, the Nedelec element +// of polynomial degree p has normal components of degree p-1 and therefore +// can only represent polynomials of degree p-1 exactly. the gap is then 1. +template +void test_no_hanging_nodes (const FiniteElement &fe, + const unsigned int p, + const unsigned int order_difference = 0) +{ + Triangulation triangulation; + GridGenerator::hyper_cube (triangulation); + triangulation.refine_global (3); + + do_project (triangulation, fe, p, order_difference); +} + + + +// same test as above, but this time with a mesh that has hanging nodes +template +void test_with_hanging_nodes (const FiniteElement &fe, + const unsigned int p, + const unsigned int order_difference = 0) +{ + Triangulation triangulation; + GridGenerator::hyper_cube (triangulation); + triangulation.refine_global (1); + triangulation.begin_active()->set_refine_flag (); + triangulation.execute_coarsening_and_refinement (); + triangulation.refine_global (1); + + do_project (triangulation, fe, p, order_difference); +} + + + +// test with a 3d grid that has cells with face_orientation==false and hanging +// nodes. this trips up all sorts of pieces of code, for example there was a +// crash when computing hanging node constraints on such faces (see +// bits/face_orientation_crash), and it triggers all sorts of other +// assumptions that may be hidden in places +template +void test_with_wrong_face_orientation (const FiniteElement &fe, + const unsigned int p, + const unsigned int order_difference = 0) +{ + if (dim < 3) + return; + + Triangulation triangulation; + GridGenerator::hyper_cube (triangulation); + triangulation.refine_global (1); + triangulation.begin_active()->set_refine_flag (); + triangulation.execute_coarsening_and_refinement (); + triangulation.refine_global (1); + + do_project (triangulation, fe, p, order_difference); +} + + + int main () { std::ofstream logfile(logname); diff --git a/tests/deal.II/project_dgp_02.cc b/tests/deal.II/project_dgp_02.cc new file mode 100644 index 0000000000..ff4dda5ef7 --- /dev/null +++ b/tests/deal.II/project_dgp_02.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_02.cc --------------------------- +// $Id: project_dgp_02.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_02.cc --------------------------- + + +// check that VectorTools::project works for DGP elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_dgp_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_hanging_nodes (FE_DGP(p), p); +} diff --git a/tests/deal.II/project_dgp_02/cmp/generic b/tests/deal.II/project_dgp_02/cmp/generic new file mode 100644 index 0000000000..9d7503b69e --- /dev/null +++ b/tests/deal.II/project_dgp_02/cmp/generic @@ -0,0 +1,178 @@ + +DEAL::n_dofs=6 +DEAL:cg::Starting value 0.433 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.967 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(0), P_1, rel. error=0.0244 +DEAL:cg::Starting value 1.60 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(0), P_2, rel. error=0.0487 +DEAL::n_dofs=12 +DEAL:cg::Starting value 0.433 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.968 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.60 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(1), P_2, rel. error=0.00142 +DEAL:cg::Starting value 2.33 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(1), P_3, rel. error=0.00405 +DEAL::n_dofs=18 +DEAL:cg::Starting value 0.433 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.968 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.60 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 2.33 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(2), P_3, rel. error=8.49e-05 +DEAL:cg::Starting value 3.12 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(2), P_4, rel. error=0.000309 +DEAL::n_dofs=24 +DEAL:cg::Starting value 0.433 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.968 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.60 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 2.33 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 3.12 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_4, rel. error=5.07e-06 +DEAL:cg::Starting value 3.98 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_5, rel. error=2.23e-05 +DEAL::n_dofs=30 +DEAL:cg::Starting value 0.433 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 0.968 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 1.60 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 2.33 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 3.12 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 3.98 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_5, rel. error=3.03e-07 +DEAL:cg::Starting value 4.91 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_6, rel. error=1.57e-06 +DEAL::n_dofs=28 +DEAL:cg::Starting value 0.676 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(0), P_0, rel. error=0 +DEAL:cg::Starting value 1.45 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<2>(0), P_1, rel. error=0.00998 +DEAL:cg::Starting value 2.31 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(0), P_2, rel. error=0.0192 +DEAL::n_dofs=84 +DEAL:cg::Starting value 0.676 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.45 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 2.32 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(1), P_2, rel. error=0.000642 +DEAL:cg::Starting value 3.27 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(1), P_3, rel. error=0.00167 +DEAL::n_dofs=168 +DEAL:cg::Starting value 0.676 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.45 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 2.32 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 3.27 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_3, rel. error=4.03e-05 +DEAL:cg::Starting value 4.28 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_4, rel. error=0.000131 +DEAL::n_dofs=280 +DEAL:cg::Starting value 0.676 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 1.45 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 2.32 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 3.27 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 4.28 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_4, rel. error=2.49e-06 +DEAL:cg::Starting value 5.36 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_5, rel. error=9.67e-06 +DEAL::n_dofs=120 +DEAL:cg::Starting value 0.708 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 1.48 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(0), P_1, rel. error=0.00433 +DEAL:cg::Starting value 2.30 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(0), P_2, rel. error=0.00820 +DEAL::n_dofs=480 +DEAL:cg::Starting value 0.708 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.48 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 2.31 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_2, rel. error=0.000288 +DEAL:cg::Starting value 3.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_3, rel. error=0.000729 +DEAL::n_dofs=1200 +DEAL:cg::Starting value 0.708 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.48 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 2.31 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 3.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_3, rel. error=1.85e-05 +DEAL:cg::Starting value 4.14 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_4, rel. error=5.78e-05 diff --git a/tests/deal.II/project_dgp_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgp_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..6e360be741 --- /dev/null +++ b/tests/deal.II/project_dgp_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 @@ -0,0 +1,178 @@ + +DEAL::n_dofs=8 +DEAL:cg::Starting value 0.354 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(0), P_1, rel. error=0.0123 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(0), P_2, rel. error=0.0203 +DEAL::n_dofs=16 +DEAL:cg::Starting value 0.354 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(1), P_2, rel. error=0.000370 +DEAL:cg::Starting value 1.68 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(1), P_3, rel. error=0.000834 +DEAL::n_dofs=24 +DEAL:cg::Starting value 0.354 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.68 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<1>(2), P_3, rel. error=1.10e-05 +DEAL:cg::Starting value 2.24 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(2), P_4, rel. error=3.12e-05 +DEAL::n_dofs=32 +DEAL:cg::Starting value 0.354 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 1.68 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 2.24 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_4, rel. error=3.25e-07 +DEAL:cg::Starting value 2.84 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(3), P_5, rel. error=1.11e-06 +DEAL::n_dofs=40 +DEAL:cg::Starting value 0.354 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 1.68 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 2.24 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 2.84 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_5, rel. error=9.64e-09 +DEAL:cg::Starting value 3.49 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<1>(4), P_6, rel. error=3.87e-08 +DEAL::n_dofs=64 +DEAL:cg::Starting value 0.375 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<2>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.767 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<2>(0), P_1, rel. error=0.00329 +DEAL:cg::Starting value 1.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(0), P_2, rel. error=0.00558 +DEAL::n_dofs=192 +DEAL:cg::Starting value 0.375 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGP<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.767 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(1), P_2, rel. error=0.000102 +DEAL:cg::Starting value 1.66 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(1), P_3, rel. error=0.000236 +DEAL::n_dofs=384 +DEAL:cg::Starting value 0.375 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.767 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.66 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_3, rel. error=3.10e-06 +DEAL:cg::Starting value 2.17 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(2), P_4, rel. error=9.01e-06 +DEAL::n_dofs=640 +DEAL:cg::Starting value 0.375 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.767 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 1.66 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 2.17 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_4, rel. error=9.38e-08 +DEAL:cg::Starting value 2.70 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<2>(3), P_5, rel. error=3.27e-07 +DEAL::n_dofs=512 +DEAL:cg::Starting value 0.265 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.539 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(0), P_1, rel. error=0.000979 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(0), P_2, rel. error=0.00168 +DEAL::n_dofs=2048 +DEAL:cg::Starting value 0.265 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.539 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_2, rel. error=3.07e-05 +DEAL:cg::Starting value 1.14 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(1), P_3, rel. error=7.18e-05 +DEAL::n_dofs=5120 +DEAL:cg::Starting value 0.265 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.539 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.14 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_3, rel. error=9.44e-07 +DEAL:cg::Starting value 1.47 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGP<3>(2), P_4, rel. error=2.77e-06 diff --git a/tests/deal.II/project_dgp_monomial_02.cc b/tests/deal.II/project_dgp_monomial_02.cc new file mode 100644 index 0000000000..3a55733740 --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_02.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_monomial_02.cc --------------------------- +// $Id: project_dgp_monomial_02.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_02.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_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_hanging_nodes (FE_DGPMonomial(p), p); +} diff --git a/tests/deal.II/project_dgp_monomial_02/cmp/generic b/tests/deal.II/project_dgp_monomial_02/cmp/generic new file mode 100644 index 0000000000..680c7eb486 --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_02/cmp/generic @@ -0,0 +1,178 @@ + +DEAL::n_dofs=6 +DEAL:cg::Starting value 0.433 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<1>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.967 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<1>(0), P_1, rel. error=0.0244 +DEAL:cg::Starting value 1.60 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<1>(0), P_2, rel. error=0.0487 +DEAL::n_dofs=12 +DEAL:cg::Starting value 0.484 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.09 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGPMonomial<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.81 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGPMonomial<1>(1), P_2, rel. error=0.00159 +DEAL:cg::Starting value 2.62 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGPMonomial<1>(1), P_3, rel. error=0.00451 +DEAL::n_dofs=18 +DEAL:cg::Starting value 0.505 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.14 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.90 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 2.76 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<1>(2), P_3, rel. error=9.83e-05 +DEAL:cg::Starting value 3.70 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<1>(2), P_4, rel. error=0.000362 +DEAL::n_dofs=24 +DEAL:cg::Starting value 0.517 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 1.17 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 1.95 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 2.84 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 3.82 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_4, rel. error=5.96e-06 +DEAL:cg::Starting value 4.88 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<1>(3), P_5, rel. error=2.64e-05 +DEAL::n_dofs=30 +DEAL:cg::Starting value 0.524 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 1.98 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 2.89 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 3.89 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 4.98 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGPMonomial<1>(4), P_5, rel. error=3.58e-07 +DEAL:cg::Starting value 6.14 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGPMonomial<1>(4), P_6, rel. error=1.85e-06 +DEAL::n_dofs=28 +DEAL:cg::Starting value 0.676 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<2>(0), P_0, rel. error=0 +DEAL:cg::Starting value 1.45 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGPMonomial<2>(0), P_1, rel. error=0.00998 +DEAL:cg::Starting value 2.31 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<2>(0), P_2, rel. error=0.0192 +DEAL::n_dofs=84 +DEAL:cg::Starting value 0.828 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.79 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 2.86 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_2, rel. error=0.000743 +DEAL:cg::Starting value 4.03 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGPMonomial<2>(1), P_3, rel. error=0.00198 +DEAL::n_dofs=168 +DEAL:cg::Starting value 0.903 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.96 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 3.15 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 4.44 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_3, rel. error=4.86e-05 +DEAL:cg::Starting value 5.83 +DEAL:cg::Convergence step 10 value 0 +DEAL::FE_DGPMonomial<2>(2), P_4, rel. error=0.000162 +DEAL::n_dofs=280 +DEAL:cg::Starting value 0.948 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_DGPMonomial<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 2.06 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGPMonomial<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 3.32 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 4.69 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 6.16 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_4, rel. error=3.09e-06 +DEAL:cg::Starting value 7.73 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_DGPMonomial<2>(3), P_5, rel. error=1.22e-05 +DEAL::n_dofs=120 +DEAL:cg::Starting value 0.708 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 1.48 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<3>(0), P_1, rel. error=0.00433 +DEAL:cg::Starting value 2.30 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGPMonomial<3>(0), P_2, rel. error=0.00820 +DEAL::n_dofs=480 +DEAL:cg::Starting value 0.936 +DEAL:cg::Convergence step 6 value 0 +DEAL::FE_DGPMonomial<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.96 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.08 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGPMonomial<3>(1), P_2, rel. error=0.000339 +DEAL:cg::Starting value 4.26 +DEAL:cg::Convergence step 7 value 0 +DEAL::FE_DGPMonomial<3>(1), P_3, rel. error=0.000889 +DEAL::n_dofs=1200 +DEAL:cg::Starting value 1.07 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_DGPMonomial<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 2.25 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_DGPMonomial<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 3.53 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_DGPMonomial<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 4.90 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGPMonomial<3>(2), P_3, rel. error=2.27e-05 +DEAL:cg::Starting value 6.35 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGPMonomial<3>(2), P_4, rel. error=7.36e-05 diff --git a/tests/deal.II/project_dgp_monomial_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgp_monomial_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..f0ec21fe78 --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 @@ -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 2 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 5 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 2 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_02.cc b/tests/deal.II/project_dgp_nonparametric_02.cc new file mode 100644 index 0000000000..f32cb616f6 --- /dev/null +++ b/tests/deal.II/project_dgp_nonparametric_02.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_nonparametric_02.cc --------------------------- +// $Id: project_dgp_nonparametric_02.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_02.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_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_hanging_nodes (FE_DGPNonparametric(p), p); +} diff --git a/tests/deal.II/project_dgq_02.cc b/tests/deal.II/project_dgq_02.cc new file mode 100644 index 0000000000..c685340163 --- /dev/null +++ b/tests/deal.II/project_dgq_02.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgq_02.cc --------------------------- +// $Id: project_dgq_02.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_dgq_02.cc --------------------------- + + +// check that VectorTools::project works for DGQ elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_dgq_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_hanging_nodes (FE_DGQ(p), p); +} diff --git a/tests/deal.II/project_dgq_02/cmp/generic b/tests/deal.II/project_dgq_02/cmp/generic new file mode 100644 index 0000000000..563a22f4d0 --- /dev/null +++ b/tests/deal.II/project_dgq_02/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_DGQ<1>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGQ<1>(0), P_1, rel. error=0.0123 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(0), P_2, rel. error=0.0203 +DEAL::n_dofs=16 +DEAL:cg::Starting value 0.250 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.520 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.833 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_2, rel. error=0.000261 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_3, rel. error=0.000586 +DEAL::n_dofs=24 +DEAL:cg::Starting value 0.250 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.520 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.833 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_3, rel. error=6.30e-06 +DEAL:cg::Starting value 1.57 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_4, rel. error=1.78e-05 +DEAL::n_dofs=32 +DEAL:cg::Starting value 0.198 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.411 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.659 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 0.940 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 1.25 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_4, rel. error=1.61e-07 +DEAL:cg::Starting value 1.59 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_5, rel. error=5.51e-07 +DEAL::n_dofs=40 +DEAL:cg::Starting value 0.188 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 0.391 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 0.627 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 0.894 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 1.51 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_5, rel. error=4.28e-09 +DEAL:cg::Starting value 1.85 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_6, rel. error=1.71e-08 +DEAL::n_dofs=64 +DEAL:cg::Starting value 0.375 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGQ<2>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.767 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGQ<2>(0), P_1, rel. error=0.00329 +DEAL:cg::Starting value 1.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<2>(0), P_2, rel. error=0.00558 +DEAL::n_dofs=256 +DEAL:cg::Starting value 0.188 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.383 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.598 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_2, rel. error=5.09e-05 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_3, rel. error=0.000117 +DEAL::n_dofs=576 +DEAL:cg::Starting value 0.187 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.383 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.598 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 0.830 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_3, rel. error=1.03e-06 +DEAL:cg::Starting value 1.08 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_4, rel. error=2.98e-06 +DEAL::n_dofs=1024 +DEAL:cg::Starting value 0.117 +DEAL:cg::Convergence step 14 value 0 +DEAL::FE_DGQ<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.240 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.374 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 0.520 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 0.676 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_4, rel. error=2.33e-08 +DEAL:cg::Starting value 0.843 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_5, rel. error=8.12e-08 +DEAL::n_dofs=512 +DEAL:cg::Starting value 0.265 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.539 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.000979 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.00168 +DEAL::n_dofs=4096 +DEAL:cg::Starting value 0.0938 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.190 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.294 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=1.08e-05 +DEAL:cg::Starting value 0.404 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=2.53e-05 +DEAL::n_dofs=13824 +DEAL:cg::Starting value 0.0938 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.190 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.294 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 0.403 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=1.81e-07 +DEAL:cg::Starting value 0.519 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=5.29e-07 diff --git a/tests/deal.II/project_dgq_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgq_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..802c43f920 --- /dev/null +++ b/tests/deal.II/project_dgq_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 @@ -0,0 +1,178 @@ + +DEAL::n_dofs=8 +DEAL:cg::Starting value 0.354 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGQ<1>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(0), P_1, rel. error=0.0123 +DEAL:cg::Starting value 1.18 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(0), P_2, rel. error=0.0203 +DEAL::n_dofs=16 +DEAL:cg::Starting value 0.250 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.520 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.833 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_2, rel. error=0.000261 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<1>(1), P_3, rel. error=0.000586 +DEAL::n_dofs=24 +DEAL:cg::Starting value 0.250 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.520 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.833 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_3, rel. error=6.30e-06 +DEAL:cg::Starting value 1.57 +DEAL:cg::Convergence step 3 value 0 +DEAL::FE_DGQ<1>(2), P_4, rel. error=1.78e-05 +DEAL::n_dofs=32 +DEAL:cg::Starting value 0.198 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.411 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.659 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 0.940 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 1.25 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_4, rel. error=1.61e-07 +DEAL:cg::Starting value 1.59 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<1>(3), P_5, rel. error=5.51e-07 +DEAL::n_dofs=40 +DEAL:cg::Starting value 0.188 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 0.391 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 0.627 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 0.894 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 1.51 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_5, rel. error=4.28e-09 +DEAL:cg::Starting value 1.85 +DEAL:cg::Convergence step 5 value 0 +DEAL::FE_DGQ<1>(4), P_6, rel. error=1.71e-08 +DEAL::n_dofs=64 +DEAL:cg::Starting value 0.375 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGQ<2>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.767 +DEAL:cg::Convergence step 1 value 0 +DEAL::FE_DGQ<2>(0), P_1, rel. error=0.00329 +DEAL:cg::Starting value 1.20 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<2>(0), P_2, rel. error=0.00558 +DEAL::n_dofs=256 +DEAL:cg::Starting value 0.188 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.383 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.598 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_2, rel. error=5.09e-05 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 4 value 0 +DEAL::FE_DGQ<2>(1), P_3, rel. error=0.000117 +DEAL::n_dofs=576 +DEAL:cg::Starting value 0.187 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.383 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.598 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 0.830 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_3, rel. error=1.03e-06 +DEAL:cg::Starting value 1.08 +DEAL:cg::Convergence step 9 value 0 +DEAL::FE_DGQ<2>(2), P_4, rel. error=2.98e-06 +DEAL::n_dofs=1024 +DEAL:cg::Starting value 0.117 +DEAL:cg::Convergence step 14 value 0 +DEAL::FE_DGQ<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.240 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.374 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 0.520 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 0.676 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_4, rel. error=2.33e-08 +DEAL:cg::Starting value 0.843 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_DGQ<2>(3), P_5, rel. error=8.12e-08 +DEAL::n_dofs=512 +DEAL:cg::Starting value 0.265 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.539 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.000979 +DEAL:cg::Starting value 0.831 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.00168 +DEAL::n_dofs=4096 +DEAL:cg::Starting value 0.0938 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.190 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 0.294 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=1.08e-05 +DEAL:cg::Starting value 0.404 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=2.53e-05 +DEAL::n_dofs=13824 +DEAL:cg::Starting value 0.0938 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.190 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.294 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 0.403 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=1.81e-07 +DEAL:cg::Starting value 0.519 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=5.29e-07 diff --git a/tests/deal.II/project_nedelec_02.cc b/tests/deal.II/project_nedelec_02.cc new file mode 100644 index 0000000000..c31b95ac11 --- /dev/null +++ b/tests/deal.II/project_nedelec_02.cc @@ -0,0 +1,31 @@ +//---------------------------- project_nedelec_02.cc --------------------------- +// $Id: project_nedelec_02.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_02.cc --------------------------- + + +// check that VectorTools::project works for Nedelec elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_nedelec_02/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_with_hanging_nodes (FE_Nedelec(p), p, 1); +} diff --git a/tests/deal.II/project_nedelec_02/cmp/generic b/tests/deal.II/project_nedelec_02/cmp/generic new file mode 100644 index 0000000000..3cdf4b0d00 --- /dev/null +++ b/tests/deal.II/project_nedelec_02/cmp/generic @@ -0,0 +1,21 @@ + +DEAL::n_dofs=72 +DEAL:cg::Starting value 5.43 +DEAL:cg::Convergence step 11 value 0 +DEAL::FE_Nedelec<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 11.0 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Nedelec<2>(1), P_1, rel. error=0.0217 +DEAL:cg::Starting value 16.7 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Nedelec<2>(1), P_2, rel. error=0.0391 +DEAL::n_dofs=576 +DEAL:cg::Starting value 6.85 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 13.9 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.00599 +DEAL:cg::Starting value 21.0 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0106 diff --git a/tests/deal.II/project_q_02.cc b/tests/deal.II/project_q_02.cc new file mode 100644 index 0000000000..87c55c44c1 --- /dev/null +++ b/tests/deal.II/project_q_02.cc @@ -0,0 +1,29 @@ +//---------------------------- project_q_02.cc --------------------------- +// $Id: project_q_02.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_02.cc --------------------------- + + +// check that VectorTools::project works for Q elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_q_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<6-dim; ++p) + test_with_hanging_nodes (FE_Q(p), p); +} diff --git a/tests/deal.II/project_q_02/cmp/generic b/tests/deal.II/project_q_02/cmp/generic new file mode 100644 index 0000000000..fedd7319ef --- /dev/null +++ b/tests/deal.II/project_q_02/cmp/generic @@ -0,0 +1,148 @@ + +DEAL::n_dofs=7 +DEAL:cg::Starting value 0.405 +DEAL:cg::Convergence step 7 value 0 +DEAL::FE_Q<1>(1), P_0, rel. error=0 +DEAL:cg::Starting value 0.884 +DEAL:cg::Convergence step 7 value 0 +DEAL::FE_Q<1>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.42 +DEAL:cg::Convergence step 7 value 0 +DEAL::FE_Q<1>(1), P_2, rel. error=0.00129 +DEAL:cg::Starting value 2.00 +DEAL:cg::Convergence step 7 value 0 +DEAL::FE_Q<1>(1), P_3, rel. error=0.00343 +DEAL::n_dofs=13 +DEAL:cg::Starting value 0.319 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Q<1>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.709 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Q<1>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.17 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Q<1>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.67 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Q<1>(2), P_3, rel. error=7.69e-05 +DEAL:cg::Starting value 2.22 +DEAL:cg::Convergence step 12 value 0 +DEAL::FE_Q<1>(2), P_4, rel. error=0.000263 +DEAL::n_dofs=19 +DEAL:cg::Starting value 0.251 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.559 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.923 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(3), P_2, rel. error=0 +DEAL:cg::Starting value 1.33 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(3), P_3, rel. error=0 +DEAL:cg::Starting value 1.78 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(3), P_4, rel. error=2.74e-06 +DEAL:cg::Starting value 2.26 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_Q<1>(3), P_5, rel. error=1.17e-05 +DEAL::n_dofs=25 +DEAL:cg::Starting value 0.234 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_0, rel. error=0 +DEAL:cg::Starting value 0.522 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_1, rel. error=0 +DEAL:cg::Starting value 0.863 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_2, rel. error=0 +DEAL:cg::Starting value 1.25 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_3, rel. error=0 +DEAL:cg::Starting value 1.67 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_4, rel. error=0 +DEAL:cg::Starting value 2.13 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_Q<1>(4), P_5, rel. error=1.71e-07 +DEAL:cg::Starting value 2.62 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<1>(4), P_6, rel. error=8.57e-07 +DEAL::n_dofs=41 +DEAL:cg::Starting value 0.597 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.26 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.97 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(1), P_2, rel. error=0.000512 +DEAL:cg::Starting value 2.69 +DEAL:cg::Convergence step 18 value 0 +DEAL::FE_Q<2>(1), P_3, rel. error=0.00127 +DEAL::n_dofs=141 +DEAL:cg::Starting value 0.367 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.785 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.25 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.74 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_3, rel. error=2.61e-05 +DEAL:cg::Starting value 2.25 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<2>(2), P_4, rel. error=7.92e-05 +DEAL::n_dofs=297 +DEAL:cg::Starting value 0.227 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_0, rel. error=0 +DEAL:cg::Starting value 0.488 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_1, rel. error=0 +DEAL:cg::Starting value 0.776 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_2, rel. error=0 +DEAL:cg::Starting value 1.09 +DEAL:cg::Convergence step 19 value 0 +DEAL::FE_Q<2>(3), P_3, rel. error=0 +DEAL:cg::Starting value 1.42 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_Q<2>(3), P_4, rel. error=7.31e-07 +DEAL:cg::Starting value 1.76 +DEAL:cg::Convergence step 20 value 0 +DEAL::FE_Q<2>(3), P_5, rel. error=2.79e-06 +DEAL::n_dofs=223 +DEAL:cg::Starting value 0.585 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.85 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.000203 +DEAL:cg::Starting value 2.49 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.000493 +DEAL::n_dofs=1375 +DEAL:cg::Starting value 0.282 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 0.588 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 0.913 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.25 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=8.13e-06 +DEAL:cg::Starting value 1.61 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=2.39e-05 diff --git a/tests/deal.II/project_q_hierarchical_01/cmp/generic b/tests/deal.II/project_q_hierarchical_01/cmp/generic index 15a4b99a24..d238b1ebc6 100644 --- a/tests/deal.II/project_q_hierarchical_01/cmp/generic +++ b/tests/deal.II/project_q_hierarchical_01/cmp/generic @@ -100,7 +100,7 @@ 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:cg::Convergence step 115 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 @@ -109,7 +109,7 @@ 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:cg::Convergence step 119 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 diff --git a/tests/deal.II/project_q_hierarchical_02.cc b/tests/deal.II/project_q_hierarchical_02.cc new file mode 100644 index 0000000000..0535b077c6 --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_02.cc @@ -0,0 +1,29 @@ +//---------------------------- project_q_hierarchical_02.cc --------------------------- +// $Id: project_q_hierarchical_02.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_02.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_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<6-dim; ++p) + test_with_hanging_nodes (FE_Q_Hierarchical(p), p); +} diff --git a/tests/deal.II/project_q_hierarchical_02/cmp/generic b/tests/deal.II/project_q_hierarchical_02/cmp/generic new file mode 100644 index 0000000000..15a4b99a24 --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_02/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_hierarchical_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_q_hierarchical_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..0002f1b2b5 --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 @@ -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 34 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 35 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 119 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_02.cc b/tests/deal.II/project_q_system_02.cc new file mode 100644 index 0000000000..66481a487f --- /dev/null +++ b/tests/deal.II/project_q_system_02.cc @@ -0,0 +1,31 @@ +//---------------------------- project_q_system_02.cc --------------------------- +// $Id: project_q_system_02.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_02.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_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<5-dim; ++p) + test_with_hanging_nodes (FESystem(FE_Q(p), 1, + FE_DGQ(p+1), 1), + p); +} diff --git a/tests/deal.II/project_q_system_02/cmp/generic b/tests/deal.II/project_q_system_02/cmp/generic new file mode 100644 index 0000000000..13b2c18ace --- /dev/null +++ b/tests/deal.II/project_q_system_02/cmp/generic @@ -0,0 +1,91 @@ + +DEAL::n_dofs=25 +DEAL:cg::Starting value 0.734 +DEAL:cg::Convergence step 10 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 1.56 +DEAL:cg::Convergence step 10 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 2.46 +DEAL:cg::Convergence step 10 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_2, rel. error=0.000418 +DEAL:cg::Starting value 3.42 +DEAL:cg::Convergence step 10 value 0 +DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_3, rel. error=0.00121 +DEAL::n_dofs=37 +DEAL:cg::Starting value 0.580 +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 1.24 +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.97 +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.76 +DEAL:cg::Convergence step 13 value 0 +DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_3, rel. error=2.98e-05 +DEAL:cg::Starting value 3.59 +DEAL:cg::Convergence step 13 value 0 +DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_4, rel. error=0.000108 +DEAL::n_dofs=49 +DEAL:cg::Starting value 0.524 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_0, rel. error=0 +DEAL:cg::Starting value 1.12 +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.77 +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 2.48 +DEAL:cg::Convergence step 14 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_3, rel. error=0 +DEAL:cg::Starting value 3.23 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_4, rel. error=1.17e-06 +DEAL:cg::Starting value 4.02 +DEAL:cg::Convergence step 15 value 0 +DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_5, rel. error=5.21e-06 +DEAL::n_dofs=293 +DEAL:cg::Starting value 0.821 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 1.72 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 2.68 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_2, rel. error=0.000125 +DEAL:cg::Starting value 3.66 +DEAL:cg::Convergence step 18 value 0 +DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_3, rel. error=0.000327 +DEAL::n_dofs=589 +DEAL:cg::Starting value 0.508 +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 1.07 +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 1.69 +DEAL:cg::Convergence step 19 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_2, rel. error=0 +DEAL:cg::Starting value 2.33 +DEAL:cg::Convergence step 19 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_3, rel. error=8.68e-06 +DEAL:cg::Starting value 3.01 +DEAL:cg::Convergence step 19 value 0 +DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_4, rel. error=2.75e-05 +DEAL::n_dofs=3463 +DEAL:cg::Starting value 0.695 +DEAL:cg::Convergence step 24 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 1.43 +DEAL:cg::Convergence step 24 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 2.19 +DEAL:cg::Convergence step 25 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=3.56e-05 +DEAL:cg::Starting value 2.96 +DEAL:cg::Convergence step 25 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=9.06e-05 diff --git a/tests/deal.II/project_rt_02.cc b/tests/deal.II/project_rt_02.cc new file mode 100644 index 0000000000..315de1af7e --- /dev/null +++ b/tests/deal.II/project_rt_02.cc @@ -0,0 +1,33 @@ +//---------------------------- project_rt_02.cc --------------------------- +// $Id: project_rt_02.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_rt_02.cc --------------------------- + + +// check that VectorTools::project works for RT elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_rt_02/output"; + + +#include "project_common.cc" + + +template +void test () +{ + if (dim != 1) + // this is interesting also in 3d, but is + // exceedingly slow there. limit to the + // case of RT(0) elements in 3d + for (unsigned int p=0; p<(dim == 2 ? 3 : 1); ++p) + test_with_hanging_nodes (FE_RaviartThomas(p), p+1, 1); +} -- 2.39.5