From 624483ff390cc89299f6b0654c730dfd83dbd369 Mon Sep 17 00:00:00 2001 From: bangerth Date: Sun, 13 Aug 2006 04:01:40 +0000 Subject: [PATCH] Add more tests that check whether projection on weird meshes with hanging node constraints and face_orientation==false does any harm to the approximation qualities of an FE. git-svn-id: https://svn.dealii.org/trunk@13692 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/deal.II/project_abf_03.cc | 30 ++ tests/deal.II/project_common.cc | 25 +- tests/deal.II/project_dgp_03.cc | 29 ++ tests/deal.II/project_dgp_03/cmp/generic | 178 ++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 178 ++++++++++++ tests/deal.II/project_dgp_monomial_03.cc | 29 ++ .../project_dgp_monomial_03/cmp/generic | 178 ++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 178 ++++++++++++ tests/deal.II/project_dgp_nonparametric_03.cc | 29 ++ tests/deal.II/project_dgq_03.cc | 29 ++ tests/deal.II/project_dgq_03/cmp/generic | 274 ++++++++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 178 ++++++++++++ tests/deal.II/project_nedelec_03.cc | 31 ++ tests/deal.II/project_nedelec_03/cmp/generic | 71 +++++ tests/deal.II/project_q_03.cc | 29 ++ tests/deal.II/project_q_03/cmp/generic | 204 +++++++++++++ tests/deal.II/project_q_hierarchical_03.cc | 29 ++ .../project_q_hierarchical_03/cmp/generic | 204 +++++++++++++ .../cmp/x86_64-unknown-linux-gnu+gcc4.1 | 148 ++++++++++ tests/deal.II/project_q_system_03.cc | 31 ++ tests/deal.II/project_q_system_03/cmp/generic | 92 ++++++ tests/deal.II/project_rt_03.cc | 33 +++ tests/deal.II/project_rt_03/cmp/generic | 50 ++++ 23 files changed, 2249 insertions(+), 8 deletions(-) create mode 100644 tests/deal.II/project_abf_03.cc create mode 100644 tests/deal.II/project_dgp_03.cc create mode 100644 tests/deal.II/project_dgp_03/cmp/generic create mode 100644 tests/deal.II/project_dgp_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_dgp_monomial_03.cc create mode 100644 tests/deal.II/project_dgp_monomial_03/cmp/generic create mode 100644 tests/deal.II/project_dgp_monomial_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_dgp_nonparametric_03.cc create mode 100644 tests/deal.II/project_dgq_03.cc create mode 100644 tests/deal.II/project_dgq_03/cmp/generic create mode 100644 tests/deal.II/project_dgq_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_nedelec_03.cc create mode 100644 tests/deal.II/project_nedelec_03/cmp/generic create mode 100644 tests/deal.II/project_q_03.cc create mode 100644 tests/deal.II/project_q_03/cmp/generic create mode 100644 tests/deal.II/project_q_hierarchical_03.cc create mode 100644 tests/deal.II/project_q_hierarchical_03/cmp/generic create mode 100644 tests/deal.II/project_q_hierarchical_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 create mode 100644 tests/deal.II/project_q_system_03.cc create mode 100644 tests/deal.II/project_q_system_03/cmp/generic create mode 100644 tests/deal.II/project_rt_03.cc create mode 100644 tests/deal.II/project_rt_03/cmp/generic diff --git a/tests/deal.II/project_abf_03.cc b/tests/deal.II/project_abf_03.cc new file mode 100644 index 0000000000..46af5b7a4c --- /dev/null +++ b/tests/deal.II/project_abf_03.cc @@ -0,0 +1,30 @@ +//---------------------------- project_abf_03.cc --------------------------- +// $Id: project_abf_03.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_03.cc --------------------------- + + +// check that VectorTools::project works for ABF elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_abf_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + if (dim > 1) + for (unsigned int p=0; p<6-dim; ++p) + test_with_wrong_face_orientation (FE_ABF(p), p); +} diff --git a/tests/deal.II/project_common.cc b/tests/deal.II/project_common.cc index 77be2216e4..f5ee57aff3 100644 --- a/tests/deal.II/project_common.cc +++ b/tests/deal.II/project_common.cc @@ -184,22 +184,31 @@ void test_with_hanging_nodes (const FiniteElement &fe, // 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 +// +// the mesh we use is the 7 cells of the hyperball mesh in 3d, with each of +// the cells refined in turn. that then makes 7 meshes with 14 active cells +// each. this also cycles through all possibilities of coarser or finer cell +// having face_orientation==false template void test_with_wrong_face_orientation (const FiniteElement &fe, const unsigned int p, const unsigned int order_difference = 0) { - if (dim < 3) + 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); + for (unsigned int i=0; i<7; ++i) + { + Triangulation triangulation; + GridGenerator::hyper_ball (triangulation); + typename Triangulation::active_cell_iterator + cell = triangulation.begin_active(); + std::advance (cell, i); + cell->set_refine_flag (); + triangulation.execute_coarsening_and_refinement (); - do_project (triangulation, fe, p, order_difference); + do_project (triangulation, fe, p, order_difference); + } } diff --git a/tests/deal.II/project_dgp_03.cc b/tests/deal.II/project_dgp_03.cc new file mode 100644 index 0000000000..a0551f9d3a --- /dev/null +++ b/tests/deal.II/project_dgp_03.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_03.cc --------------------------- +// $Id: project_dgp_03.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_03.cc --------------------------- + + +// check that VectorTools::project works for DGP elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_dgp_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_wrong_face_orientation (FE_DGP(p), p); +} diff --git a/tests/deal.II/project_dgp_03/cmp/generic b/tests/deal.II/project_dgp_03/cmp/generic new file mode 100644 index 0000000000..adfd75df87 --- /dev/null +++ b/tests/deal.II/project_dgp_03/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_DGP<1>(0), P_0, rel. error=0 +DEAL:cg::Starting value 0.736 +DEAL:cg::Convergence step 1 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 2 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 1 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 2 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 1 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 1 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 1 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 1 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_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgp_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..6e360be741 --- /dev/null +++ b/tests/deal.II/project_dgp_03/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_03.cc b/tests/deal.II/project_dgp_monomial_03.cc new file mode 100644 index 0000000000..7a897eca93 --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_03.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_monomial_03.cc --------------------------- +// $Id: project_dgp_monomial_03.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_03.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_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_wrong_face_orientation (FE_DGPMonomial(p), p); +} diff --git a/tests/deal.II/project_dgp_monomial_03/cmp/generic b/tests/deal.II/project_dgp_monomial_03/cmp/generic new file mode 100644 index 0000000000..8a3ee90e5b --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_03/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_monomial_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgp_monomial_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..f0ec21fe78 --- /dev/null +++ b/tests/deal.II/project_dgp_monomial_03/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_03.cc b/tests/deal.II/project_dgp_nonparametric_03.cc new file mode 100644 index 0000000000..bddae02e62 --- /dev/null +++ b/tests/deal.II/project_dgp_nonparametric_03.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgp_nonparametric_03.cc --------------------------- +// $Id: project_dgp_nonparametric_03.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_03.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_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_wrong_face_orientation (FE_DGPNonparametric(p), p); +} diff --git a/tests/deal.II/project_dgq_03.cc b/tests/deal.II/project_dgq_03.cc new file mode 100644 index 0000000000..b5f2f3d008 --- /dev/null +++ b/tests/deal.II/project_dgq_03.cc @@ -0,0 +1,29 @@ +//---------------------------- project_dgq_03.cc --------------------------- +// $Id: project_dgq_03.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_03.cc --------------------------- + + +// check that VectorTools::project works for DGQ elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_dgq_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=0; p<6-dim; ++p) + test_with_wrong_face_orientation (FE_DGQ(p), p); +} diff --git a/tests/deal.II/project_dgq_03/cmp/generic b/tests/deal.II/project_dgq_03/cmp/generic new file mode 100644 index 0000000000..2fe4eccff4 --- /dev/null +++ b/tests/deal.II/project_dgq_03/cmp/generic @@ -0,0 +1,274 @@ + +DEAL::n_dofs=14 +DEAL:cg::Starting value 3.59 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 3.77 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0872 +DEAL:cg::Starting value 4.97 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0871 +DEAL::n_dofs=14 +DEAL:cg::Starting value 3.35 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 3.72 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.101 +DEAL:cg::Starting value 4.82 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0883 +DEAL::n_dofs=14 +DEAL:cg::Starting value 3.35 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 3.47 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0745 +DEAL:cg::Starting value 4.61 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0716 +DEAL::n_dofs=14 +DEAL:cg::Starting value 3.35 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 3.26 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0697 +DEAL:cg::Starting value 4.32 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0632 +DEAL::n_dofs=14 +DEAL:cg::Starting value 3.35 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 3.62 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0844 +DEAL:cg::Starting value 4.76 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0797 +DEAL::n_dofs=14 +DEAL:cg::Starting value 3.35 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 3.68 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0917 +DEAL:cg::Starting value 4.79 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0834 +DEAL::n_dofs=14 +DEAL:cg::Starting value 3.35 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 3.37 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_1, rel. error=0.0716 +DEAL:cg::Starting value 4.48 +DEAL:cg::Convergence step 2 value 0 +DEAL::FE_DGQ<3>(0), P_2, rel. error=0.0670 +DEAL::n_dofs=112 +DEAL:cg::Starting value 1.32 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.41 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.87 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00962 +DEAL:cg::Starting value 1.94 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00993 +DEAL::n_dofs=112 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.39 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.81 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00950 +DEAL:cg::Starting value 1.91 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=0.0102 +DEAL::n_dofs=112 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.29 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.73 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00751 +DEAL:cg::Starting value 1.79 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00767 +DEAL::n_dofs=112 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.62 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00712 +DEAL:cg::Starting value 1.63 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00699 +DEAL::n_dofs=112 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.79 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00825 +DEAL:cg::Starting value 1.87 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00868 +DEAL::n_dofs=112 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.37 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.80 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00877 +DEAL:cg::Starting value 1.89 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00934 +DEAL::n_dofs=112 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 1.25 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 1.68 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_2, rel. error=0.00725 +DEAL:cg::Starting value 1.72 +DEAL:cg::Convergence step 26 value 0 +DEAL::FE_DGQ<3>(1), P_3, rel. error=0.00727 +DEAL::n_dofs=378 +DEAL:cg::Starting value 1.30 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.38 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.77 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.83 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00166 +DEAL:cg::Starting value 1.96 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00164 +DEAL::n_dofs=378 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.36 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.72 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.80 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00169 +DEAL:cg::Starting value 1.92 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00161 +DEAL::n_dofs=378 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.27 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.64 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.70 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00126 +DEAL:cg::Starting value 1.82 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00124 +DEAL::n_dofs=378 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.19 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.53 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.55 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00119 +DEAL:cg::Starting value 1.65 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00116 +DEAL::n_dofs=378 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.32 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.70 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.77 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00142 +DEAL:cg::Starting value 1.89 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00138 +DEAL::n_dofs=378 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.34 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.71 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.79 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00153 +DEAL:cg::Starting value 1.90 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00147 +DEAL::n_dofs=378 +DEAL:cg::Starting value 1.21 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.23 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.59 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.63 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_3, rel. error=0.00121 +DEAL:cg::Starting value 1.75 +DEAL:cg::Convergence step 25 value 0 +DEAL::FE_DGQ<3>(2), P_4, rel. error=0.00119 diff --git a/tests/deal.II/project_dgq_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgq_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..802c43f920 --- /dev/null +++ b/tests/deal.II/project_dgq_03/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_03.cc b/tests/deal.II/project_nedelec_03.cc new file mode 100644 index 0000000000..bea7bdda22 --- /dev/null +++ b/tests/deal.II/project_nedelec_03.cc @@ -0,0 +1,31 @@ +//---------------------------- project_nedelec_03.cc --------------------------- +// $Id: project_nedelec_03.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_03.cc --------------------------- + + +// check that VectorTools::project works for Nedelec elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_nedelec_03/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_wrong_face_orientation (FE_Nedelec(p), p, 1); +} diff --git a/tests/deal.II/project_nedelec_03/cmp/generic b/tests/deal.II/project_nedelec_03/cmp/generic new file mode 100644 index 0000000000..d5ceb9b0fa --- /dev/null +++ b/tests/deal.II/project_nedelec_03/cmp/generic @@ -0,0 +1,71 @@ + +DEAL::n_dofs=86 +DEAL:cg::Starting value 10.6 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 15.1 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.0356 +DEAL:cg::Starting value 21.3 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0328 +DEAL::n_dofs=86 +DEAL:cg::Starting value 10.8 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 15.4 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.0345 +DEAL:cg::Starting value 21.8 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0312 +DEAL::n_dofs=86 +DEAL:cg::Starting value 10.6 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 14.9 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.0279 +DEAL:cg::Starting value 21.0 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0256 +DEAL::n_dofs=86 +DEAL:cg::Starting value 10.8 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 15.3 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.0316 +DEAL:cg::Starting value 21.7 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0293 +DEAL::n_dofs=86 +DEAL:cg::Starting value 10.6 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 15.0 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.0294 +DEAL:cg::Starting value 21.1 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0272 +DEAL::n_dofs=86 +DEAL:cg::Starting value 10.7 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 15.2 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.0323 +DEAL:cg::Starting value 21.4 +DEAL:cg::Convergence step 24 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0291 +DEAL::n_dofs=86 +DEAL:cg::Starting value 10.7 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 15.0 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.0299 +DEAL:cg::Starting value 21.2 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0277 diff --git a/tests/deal.II/project_q_03.cc b/tests/deal.II/project_q_03.cc new file mode 100644 index 0000000000..c7d7ac9724 --- /dev/null +++ b/tests/deal.II/project_q_03.cc @@ -0,0 +1,29 @@ +//---------------------------- project_q_03.cc --------------------------- +// $Id: project_q_03.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_03.cc --------------------------- + + +// check that VectorTools::project works for Q elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_q_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<6-dim; ++p) + test_with_wrong_face_orientation (FE_Q(p), p); +} diff --git a/tests/deal.II/project_q_03/cmp/generic b/tests/deal.II/project_q_03/cmp/generic new file mode 100644 index 0000000000..6e3ce60ae9 --- /dev/null +++ b/tests/deal.II/project_q_03/cmp/generic @@ -0,0 +1,204 @@ + +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.36 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.47 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.28 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.0200 +DEAL:cg::Starting value 3.36 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.0213 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.45 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.23 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.0210 +DEAL:cg::Starting value 3.33 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.0229 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.39 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.18 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.0147 +DEAL:cg::Starting value 3.25 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.0160 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.35 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.13 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.0154 +DEAL:cg::Starting value 3.18 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.0156 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.42 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.21 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.0162 +DEAL:cg::Starting value 3.30 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.0175 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.44 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.22 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.0193 +DEAL:cg::Starting value 3.31 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.0206 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.37 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.16 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_2, rel. error=0.0158 +DEAL:cg::Starting value 3.22 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q<3>(1), P_3, rel. error=0.0165 +DEAL::n_dofs=195 +DEAL:cg::Starting value 1.44 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.54 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.99 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 2.06 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=0.00312 +DEAL:cg::Starting value 2.20 +DEAL:cg::Convergence step 23 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=0.00305 +DEAL::n_dofs=194 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 21 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.52 +DEAL:cg::Convergence step 21 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.93 +DEAL:cg::Convergence step 21 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 2.02 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=0.00316 +DEAL:cg::Starting value 2.16 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=0.00295 +DEAL::n_dofs=194 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 21 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.42 +DEAL:cg::Convergence step 21 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.85 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.91 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=0.00234 +DEAL:cg::Starting value 2.05 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=0.00227 +DEAL::n_dofs=194 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.34 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.74 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.76 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=0.00211 +DEAL:cg::Starting value 1.89 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=0.00204 +DEAL::n_dofs=194 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.48 +DEAL:cg::Convergence step 21 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.90 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.99 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=0.00267 +DEAL:cg::Starting value 2.13 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=0.00253 +DEAL::n_dofs=194 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.50 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.92 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 2.01 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=0.00288 +DEAL:cg::Starting value 2.14 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=0.00270 +DEAL::n_dofs=194 +DEAL:cg::Starting value 1.35 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 1.38 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 1.80 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 1.84 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_3, rel. error=0.00221 +DEAL:cg::Starting value 1.98 +DEAL:cg::Convergence step 22 value 0 +DEAL::FE_Q<3>(2), P_4, rel. error=0.00214 diff --git a/tests/deal.II/project_q_hierarchical_03.cc b/tests/deal.II/project_q_hierarchical_03.cc new file mode 100644 index 0000000000..3dde6e3774 --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_03.cc @@ -0,0 +1,29 @@ +//---------------------------- project_q_hierarchical_03.cc --------------------------- +// $Id: project_q_hierarchical_03.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_03.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_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<6-dim; ++p) + test_with_wrong_face_orientation (FE_Q_Hierarchical(p), p); +} diff --git a/tests/deal.II/project_q_hierarchical_03/cmp/generic b/tests/deal.II/project_q_hierarchical_03/cmp/generic new file mode 100644 index 0000000000..5a3c966224 --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_03/cmp/generic @@ -0,0 +1,204 @@ + +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.36 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.47 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.28 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=0.0200 +DEAL:cg::Starting value 3.36 +DEAL:cg::Convergence step 16 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=0.0213 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.45 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.23 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=0.0210 +DEAL:cg::Starting value 3.33 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=0.0229 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.39 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.18 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=0.0147 +DEAL:cg::Starting value 3.25 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=0.0160 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.35 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.13 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=0.0154 +DEAL:cg::Starting value 3.18 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=0.0156 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.42 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.21 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=0.0162 +DEAL:cg::Starting value 3.30 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=0.0175 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.44 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.22 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=0.0193 +DEAL:cg::Starting value 3.31 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=0.0206 +DEAL::n_dofs=35 +DEAL:cg::Starting value 2.29 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0 +DEAL:cg::Starting value 2.37 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0 +DEAL:cg::Starting value 3.16 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=0.0158 +DEAL:cg::Starting value 3.22 +DEAL:cg::Convergence step 17 value 0 +DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=0.0165 +DEAL::n_dofs=195 +DEAL:cg::Starting value 4.87 +DEAL:cg::Convergence step 140 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 5.09 +DEAL:cg::Convergence step 140 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 6.59 +DEAL:cg::Convergence step 140 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 6.75 +DEAL:cg::Convergence step 139 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=0.00640 +DEAL:cg::Starting value 7.23 +DEAL:cg::Convergence step 140 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=0.00592 +DEAL::n_dofs=194 +DEAL:cg::Starting value 4.66 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 5.04 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 6.45 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 6.66 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=0.00630 +DEAL:cg::Starting value 7.11 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=0.00563 +DEAL::n_dofs=194 +DEAL:cg::Starting value 4.66 +DEAL:cg::Convergence step 149 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 4.83 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 6.28 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 6.41 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=0.00499 +DEAL:cg::Starting value 6.88 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=0.00466 +DEAL::n_dofs=194 +DEAL:cg::Starting value 4.66 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 4.68 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 6.08 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 6.14 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=0.00467 +DEAL:cg::Starting value 6.59 +DEAL:cg::Convergence step 153 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=0.00437 +DEAL::n_dofs=194 +DEAL:cg::Starting value 4.66 +DEAL:cg::Convergence step 153 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 4.95 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 6.39 +DEAL:cg::Convergence step 153 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 6.57 +DEAL:cg::Convergence step 153 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=0.00550 +DEAL:cg::Starting value 7.04 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=0.00499 +DEAL::n_dofs=194 +DEAL:cg::Starting value 4.66 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 5.00 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 6.42 +DEAL:cg::Convergence step 153 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 6.62 +DEAL:cg::Convergence step 153 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=0.00584 +DEAL:cg::Starting value 7.08 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=0.00524 +DEAL::n_dofs=194 +DEAL:cg::Starting value 4.66 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0 +DEAL:cg::Starting value 4.76 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0 +DEAL:cg::Starting value 6.18 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0 +DEAL:cg::Starting value 6.29 +DEAL:cg::Convergence step 151 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=0.00480 +DEAL:cg::Starting value 6.75 +DEAL:cg::Convergence step 152 value 0 +DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=0.00450 diff --git a/tests/deal.II/project_q_hierarchical_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_q_hierarchical_03/cmp/x86_64-unknown-linux-gnu+gcc4.1 new file mode 100644 index 0000000000..0002f1b2b5 --- /dev/null +++ b/tests/deal.II/project_q_hierarchical_03/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_03.cc b/tests/deal.II/project_q_system_03.cc new file mode 100644 index 0000000000..23067aa4ef --- /dev/null +++ b/tests/deal.II/project_q_system_03.cc @@ -0,0 +1,31 @@ +//---------------------------- project_q_system_03.cc --------------------------- +// $Id: project_q_system_03.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_03.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_03/output"; + + +#include "project_common.cc" + + +template +void test () +{ + for (unsigned int p=1; p<5-dim; ++p) + test_with_wrong_face_orientation (FESystem(FE_Q(p), 1, + FE_DGQ(p+1), 1), + p); +} diff --git a/tests/deal.II/project_q_system_03/cmp/generic b/tests/deal.II/project_q_system_03/cmp/generic new file mode 100644 index 0000000000..c763391215 --- /dev/null +++ b/tests/deal.II/project_q_system_03/cmp/generic @@ -0,0 +1,92 @@ + +DEAL::n_dofs=413 +DEAL:cg::Starting value 3.06 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 3.62 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 4.93 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=0.00269 +DEAL:cg::Starting value 5.50 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=0.00261 +DEAL::n_dofs=413 +DEAL:cg::Starting value 2.93 +DEAL:cg::Convergence step 26 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 3.53 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 4.77 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=0.00273 +DEAL:cg::Starting value 5.33 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=0.00270 +DEAL::n_dofs=413 +DEAL:cg::Starting value 2.93 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 3.42 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 4.67 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=0.00220 +DEAL:cg::Starting value 5.19 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=0.00223 +DEAL::n_dofs=413 +DEAL:cg::Starting value 2.93 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 3.34 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 4.57 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=0.00238 +DEAL:cg::Starting value 5.04 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=0.00226 +DEAL::n_dofs=413 +DEAL:cg::Starting value 2.93 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 3.48 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 4.74 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=0.00230 +DEAL:cg::Starting value 5.28 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=0.00228 +DEAL::n_dofs=413 +DEAL:cg::Starting value 2.93 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 3.51 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 4.75 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=0.00263 +DEAL:cg::Starting value 5.31 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=0.00255 +DEAL::n_dofs=413 +DEAL:cg::Starting value 2.93 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0 +DEAL:cg::Starting value 3.38 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0 +DEAL:cg::Starting value 4.62 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=0.00239 +DEAL:cg::Starting value 5.12 +DEAL:cg::Convergence step 27 value 0 +DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=0.00234 diff --git a/tests/deal.II/project_rt_03.cc b/tests/deal.II/project_rt_03.cc new file mode 100644 index 0000000000..cce90f97f3 --- /dev/null +++ b/tests/deal.II/project_rt_03.cc @@ -0,0 +1,33 @@ +//---------------------------- project_rt_03.cc --------------------------- +// $Id: project_rt_03.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_03.cc --------------------------- + + +// check that VectorTools::project works for RT elements correctly on +// a uniformly refined mesh for functions of degree q + +char logname[] = "project_rt_03/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_wrong_face_orientation (FE_RaviartThomas(p), p+1, 1); +} diff --git a/tests/deal.II/project_rt_03/cmp/generic b/tests/deal.II/project_rt_03/cmp/generic new file mode 100644 index 0000000000..af7f5e1406 --- /dev/null +++ b/tests/deal.II/project_rt_03/cmp/generic @@ -0,0 +1,50 @@ + +DEAL::n_dofs=144 +DEAL:cg::Starting value 5.65 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_RaviartThomas<2>(0), P_0, rel. error=0 +DEAL:cg::Starting value 11.4 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_RaviartThomas<2>(0), P_1, rel. error=0.0129 +DEAL:cg::Starting value 17.4 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_RaviartThomas<2>(0), P_2, rel. error=0.0220 +DEAL::n_dofs=544 +DEAL:cg::Starting value 5.83 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_RaviartThomas<2>(1), P_0, rel. error=0 +DEAL:cg::Starting value 11.8 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_RaviartThomas<2>(1), P_1, rel. error=0 +DEAL:cg::Starting value 18.1 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_RaviartThomas<2>(1), P_2, rel. error=0.000294 +DEAL:cg::Starting value 24.7 +DEAL:cg::Convergence step 13 value 0 +DEAL::FE_RaviartThomas<2>(1), P_3, rel. error=0.000688 +DEAL::n_dofs=1200 +DEAL:cg::Starting value 5.83 +DEAL:cg::Convergence step 14 value 0 +DEAL::FE_RaviartThomas<2>(2), P_0, rel. error=0 +DEAL:cg::Starting value 11.8 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_RaviartThomas<2>(2), P_1, rel. error=0 +DEAL:cg::Starting value 18.1 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_RaviartThomas<2>(2), P_2, rel. error=0 +DEAL:cg::Starting value 24.7 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_RaviartThomas<2>(2), P_3, rel. error=9.05e-06 +DEAL:cg::Starting value 31.6 +DEAL:cg::Convergence step 15 value 0 +DEAL::FE_RaviartThomas<2>(2), P_4, rel. error=2.65e-05 +DEAL::n_dofs=1728 +DEAL:cg::Starting value 44.2 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_RaviartThomas<3>(0), P_0, rel. error=0 +DEAL:cg::Starting value 89.0 +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_RaviartThomas<3>(0), P_1, rel. error=0.0313 +DEAL:cg::Starting value 135. +DEAL:cg::Convergence step 8 value 0 +DEAL::FE_RaviartThomas<3>(0), P_2, rel. error=0.0542 -- 2.39.5