From 28cad03bb96d0bd95cab6022b6d7531f361626c5 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Tue, 18 May 2021 10:25:19 +0200 Subject: [PATCH] Remove second-order Nedelec and Raviart-Thomas elements --- tests/numerics/generalized_interpolation_02.cc | 15 +-------------- .../numerics/generalized_interpolation_02.output | 12 ------------ tests/numerics/generalized_interpolation_04.cc | 2 +- .../numerics/generalized_interpolation_04.output | 4 ++-- tests/numerics/generalized_interpolation_05.cc | 2 +- .../numerics/generalized_interpolation_05.output | 4 ++-- 6 files changed, 7 insertions(+), 32 deletions(-) diff --git a/tests/numerics/generalized_interpolation_02.cc b/tests/numerics/generalized_interpolation_02.cc index ab6073e8b1..ec5dd16470 100644 --- a/tests/numerics/generalized_interpolation_02.cc +++ b/tests/numerics/generalized_interpolation_02.cc @@ -48,11 +48,6 @@ main() 2, false, true); - test<2>(FE_RaviartThomas<2>(2), - Functions::ConstantFunction<2>(1.0, 2), - 3, - false, - true); test<3>(FE_RaviartThomas<3>(0), Functions::ConstantFunction<3>(1.0, 3), 1, @@ -74,11 +69,7 @@ main() 2, true, true); - test<2>(FE_RaviartThomas<2>(2), - Functions::ConstantFunction<2>(1.0, 2), - 3, - true, - true); + // lowest order RT in 3D does not contain constant 1 function on a // distorted mesh. test<3>(FE_RaviartThomas<3>(1), @@ -91,12 +82,8 @@ main() FE_Nedelec<2>(0), Functions::ConstantFunction<2>(1.0, 2), 1, false, true); test<2>( FE_Nedelec<2>(1), Functions::ConstantFunction<2>(1.0, 2), 2, false, true); - test<2>( - FE_Nedelec<2>(2), Functions::ConstantFunction<2>(1.0, 2), 3, false, true); test<2>( FE_Nedelec<2>(0), Functions::ConstantFunction<2>(1.0, 2), 1, true, true); test<2>( FE_Nedelec<2>(1), Functions::ConstantFunction<2>(1.0, 2), 2, true, true); - test<2>( - FE_Nedelec<2>(2), Functions::ConstantFunction<2>(1.0, 2), 3, true, true); } diff --git a/tests/numerics/generalized_interpolation_02.output b/tests/numerics/generalized_interpolation_02.output index 3ce63ec0b5..d36d72fe0d 100644 --- a/tests/numerics/generalized_interpolation_02.output +++ b/tests/numerics/generalized_interpolation_02.output @@ -23,9 +23,6 @@ DEAL::Check projection property: 5.55112e-17 DEAL::dim 2 FE_RaviartThomas<2>(1) DEAL::Function value at (0.0,0.0): 1.00000 1.00000 DEAL::Check projection property: 2.35922e-16 -DEAL::dim 2 FE_RaviartThomas<2>(2) -DEAL::Function value at (0.0,0.0): 1.00000 1.00000 -DEAL::Check projection property: 0.169271 DEAL::dim 3 FE_RaviartThomas<3>(0) DEAL::Function value at (0.0,0.0): 1.00000 1.00000 1.00000 DEAL::Check projection property: 1.11022e-16 @@ -38,9 +35,6 @@ DEAL::Check projection property: 1.11022e-16 DEAL::dim 2 FE_RaviartThomas<2>(1) DEAL::Function value at (0.0,0.0): 1.00000 1.00000 DEAL::Check projection property: 4.44089e-16 -DEAL::dim 2 FE_RaviartThomas<2>(2) -DEAL::Function value at (0.0,0.0): 1.00000 1.00000 -DEAL::Check projection property: 0.169442 DEAL::dim 3 FE_RaviartThomas<3>(1) DEAL::Function value at (0.0,0.0): 1.00000 1.00000 1.00000 DEAL::Check projection property: 2.42861e-16 @@ -50,15 +44,9 @@ DEAL::Check projection property: 2.77556e-17 DEAL::dim 2 FE_Nedelec<2>(1) DEAL::Function value at (0.0,0.0): 1.00000 1.00000 DEAL::Check projection property: 4.44089e-16 -DEAL::dim 2 FE_Nedelec<2>(2) -DEAL::Function value at (0.0,0.0): 1.00000 1.00000 -DEAL::Check projection property: 0.378501 DEAL::dim 2 FE_Nedelec<2>(0) DEAL::Function value at (0.0,0.0): 1.00000 1.00000 DEAL::Check projection property: 1.11022e-16 DEAL::dim 2 FE_Nedelec<2>(1) DEAL::Function value at (0.0,0.0): 1.00000 1.00000 DEAL::Check projection property: 4.44089e-16 -DEAL::dim 2 FE_Nedelec<2>(2) -DEAL::Function value at (0.0,0.0): 1.00000 1.00000 -DEAL::Check projection property: 0.379571 diff --git a/tests/numerics/generalized_interpolation_04.cc b/tests/numerics/generalized_interpolation_04.cc index 964f68ceb3..e4c0c27e78 100644 --- a/tests/numerics/generalized_interpolation_04.cc +++ b/tests/numerics/generalized_interpolation_04.cc @@ -36,7 +36,7 @@ main() 1, FE_Q(1), 1, - FE_Nedelec(2), + FE_Nedelec(1), 2, FESystem(FE_Q(1), dim), 1); diff --git a/tests/numerics/generalized_interpolation_04.output b/tests/numerics/generalized_interpolation_04.output index a1da03bd38..17a349df6f 100644 --- a/tests/numerics/generalized_interpolation_04.output +++ b/tests/numerics/generalized_interpolation_04.output @@ -1,3 +1,3 @@ -DEAL::dim 2 FESystem<2>[FE_RaviartThomas<2>(1)-FE_Q<2>(1)-FE_Nedelec<2>(2)^2-FESystem<2>[FE_Q<2>(1)^2]] -DEAL::Check projection property: 0.699901 +DEAL::dim 2 FESystem<2>[FE_RaviartThomas<2>(1)-FE_Q<2>(1)-FE_Nedelec<2>(1)^2-FESystem<2>[FE_Q<2>(1)^2]] +DEAL::Check projection property: 6.66134e-16 diff --git a/tests/numerics/generalized_interpolation_05.cc b/tests/numerics/generalized_interpolation_05.cc index 96e275261c..400da00946 100644 --- a/tests/numerics/generalized_interpolation_05.cc +++ b/tests/numerics/generalized_interpolation_05.cc @@ -36,7 +36,7 @@ main() 2, FE_Q(1), 1, - FESystem(FE_Nedelec(2), 2), + FESystem(FE_Nedelec(1), 2), 1, FESystem(FE_Q(1), dim), 1); diff --git a/tests/numerics/generalized_interpolation_05.output b/tests/numerics/generalized_interpolation_05.output index ffc98b78f6..8d316bec53 100644 --- a/tests/numerics/generalized_interpolation_05.output +++ b/tests/numerics/generalized_interpolation_05.output @@ -1,3 +1,3 @@ -DEAL::dim 2 FESystem<2>[FE_RaviartThomas<2>(1)^2-FE_Q<2>(1)-FESystem<2>[FE_Nedelec<2>(2)^2]-FESystem<2>[FE_Q<2>(1)^2]] -DEAL::Check projection property: 0.699901 +DEAL::dim 2 FESystem<2>[FE_RaviartThomas<2>(1)^2-FE_Q<2>(1)-FESystem<2>[FE_Nedelec<2>(1)^2]-FESystem<2>[FE_Q<2>(1)^2]] +DEAL::Check projection property: 6.66134e-16 -- 2.39.5