From: Wolfgang Bangerth Date: Mon, 29 Oct 2007 03:44:59 +0000 (+0000) Subject: Add RT tests for hp. The corresponding functions in the RT element are currently... X-Git-Tag: v8.0.0~9709 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f383fa3ba30cee3c45e3c854d84c26f4e8788bc7;p=dealii.git Add RT tests for hp. The corresponding functions in the RT element are currently not implemented, so this isn't going to work. git-svn-id: https://svn.dealii.org/trunk@15384 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/fail/hp_constraints_rt_01.cc b/tests/fail/hp_constraints_rt_01.cc new file mode 100644 index 0000000000..b6cbcbd388 --- /dev/null +++ b/tests/fail/hp_constraints_rt_01.cc @@ -0,0 +1,34 @@ +//---------------------------- hp_constraints_rt_01.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2006, 2007 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. +// +//---------------------------- hp_constraints_rt_01.cc --------------------------- + + +// check that computation of hp constraints works for RT elements correctly + +char logname[] = "hp_constraints_rt_01/output"; + + +#include "../hp/hp_constraints_common.h" +#include + + +template +void test () +{ + if (dim == 1) + return; + + hp::FECollection fe; + for (unsigned int i=1; i<4; ++i) + fe.push_back (FE_RaviartThomas(i)); + test_no_hanging_nodes (fe); +} diff --git a/tests/fail/hp_constraints_rt_02.cc b/tests/fail/hp_constraints_rt_02.cc new file mode 100644 index 0000000000..d3e916e8e5 --- /dev/null +++ b/tests/fail/hp_constraints_rt_02.cc @@ -0,0 +1,34 @@ +//---------------------------- hp_constraints_rt_02.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2006, 2007 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. +// +//---------------------------- hp_constraints_rt_02.cc --------------------------- + + +// check that computation of hp constraints works for RT elements correctly + +char logname[] = "hp_constraints_rt_02/output"; + + +#include "../hp/hp_constraints_common.h" +#include + + +template +void test () +{ + if (dim == 1) + return; + + hp::FECollection fe; + for (unsigned int i=1; i<4; ++i) + fe.push_back (FE_RaviartThomas(i)); + test_with_hanging_nodes (fe); +} diff --git a/tests/fail/hp_constraints_rt_03.cc b/tests/fail/hp_constraints_rt_03.cc new file mode 100644 index 0000000000..d58263552f --- /dev/null +++ b/tests/fail/hp_constraints_rt_03.cc @@ -0,0 +1,34 @@ +//---------------------------- hp_constraints_rt_03.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2006, 2007 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. +// +//---------------------------- hp_constraints_rt_03.cc --------------------------- + + +// check that computation of hp constraints works for RT elements correctly + +char logname[] = "hp_constraints_rt_03/output"; + + +#include "../hp/hp_constraints_common.h" +#include + + +template +void test () +{ + if (dim == 1) + return; + + hp::FECollection fe; + for (unsigned int i=1; i<4; ++i) + fe.push_back (FE_RaviartThomas(i)); + test_with_wrong_face_orientation (fe); +} diff --git a/tests/fail/hp_constraints_rt_04.cc b/tests/fail/hp_constraints_rt_04.cc new file mode 100644 index 0000000000..1d6bfd7bc5 --- /dev/null +++ b/tests/fail/hp_constraints_rt_04.cc @@ -0,0 +1,34 @@ +//---------------------------- hp_constraints_rt_04.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2006, 2007 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. +// +//---------------------------- hp_constraints_rt_04.cc --------------------------- + + +// check that computation of hp constraints works for RT elements correctly + +char logname[] = "hp_constraints_rt_04/output"; + + +#include "../hp/hp_constraints_common.h" +#include + + +template +void test () +{ + if (dim == 1) + return; + + hp::FECollection fe; + for (unsigned int i=1; i<4; ++i) + fe.push_back (FE_RaviartThomas(i)); + test_with_2d_deformed_mesh (fe); +} diff --git a/tests/fail/hp_constraints_rt_05.cc b/tests/fail/hp_constraints_rt_05.cc new file mode 100644 index 0000000000..ed0b9f9083 --- /dev/null +++ b/tests/fail/hp_constraints_rt_05.cc @@ -0,0 +1,34 @@ +//---------------------------- hp_constraints_rt_05.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2006, 2007 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. +// +//---------------------------- hp_constraints_rt_05.cc --------------------------- + + +// check that computation of hp constraints works for RT elements correctly + +char logname[] = "hp_constraints_rt_05/output"; + + +#include "../hp/hp_constraints_common.h" +#include + + +template +void test () +{ + if (dim == 1) + return; + + hp::FECollection fe; + for (unsigned int i=1; i<4; ++i) + fe.push_back (FE_RaviartThomas(i)); + test_with_2d_deformed_refined_mesh (fe); +} diff --git a/tests/fail/hp_constraints_rt_06.cc b/tests/fail/hp_constraints_rt_06.cc new file mode 100644 index 0000000000..64980f8767 --- /dev/null +++ b/tests/fail/hp_constraints_rt_06.cc @@ -0,0 +1,39 @@ +//---------------------------- hp_constraints_rt_06.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2006, 2007 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. +// +//---------------------------- hp_constraints_rt_06.cc --------------------------- + + +// check that computation of hp constraints works for RT elements correctly + +char logname[] = "hp_constraints_rt_06/output"; + + +#include "../hp/hp_constraints_common.h" +#include + + +template +void test () +{ + if (dim == 1) + return; + + hp::FECollection fe; + std::vector degrees; + for (unsigned int i=1; i<7-dim; ++i) + { + fe.push_back (FE_RaviartThomas(i)); + degrees.push_back (i); + } + + test_interpolation (fe, degrees); +}