From 4fc65610bb3846d805f6b20755a29ab1eedfec29 Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Sun, 2 Jul 2023 17:26:02 +0200 Subject: [PATCH] Add project test --- tests/simplex/simplex_project_cg.cc | 6 +- tests/simplex/simplex_project_cg.output | 82 +++++++++++++++---------- tests/simplex/simplex_project_dg.cc | 2 + tests/simplex/simplex_project_dg.output | 18 ++++++ 4 files changed, 74 insertions(+), 34 deletions(-) diff --git a/tests/simplex/simplex_project_cg.cc b/tests/simplex/simplex_project_cg.cc index e7d3ac6e93..314b5f0dbd 100644 --- a/tests/simplex/simplex_project_cg.cc +++ b/tests/simplex/simplex_project_cg.cc @@ -84,10 +84,10 @@ test(const unsigned int degree) { #ifdef SIMPLEX FE_SimplexP fe(degree); - QGaussSimplex quadrature(degree + 2); + QGaussSimplex quadrature(degree + 1); #else FE_Q fe(degree); - QGauss quadrature(degree + 2); + QGauss quadrature(degree + 1); #endif deallog << "FE = " << fe.get_name() << std::endl; @@ -199,7 +199,9 @@ main() test<2>(1); test<2>(2); + // test<3>(3); test<3>(1); test<3>(2); + // test<3>(3); } diff --git a/tests/simplex/simplex_project_cg.output b/tests/simplex/simplex_project_cg.output index 287c3090a0..b3608296a3 100644 --- a/tests/simplex/simplex_project_cg.output +++ b/tests/simplex/simplex_project_cg.output @@ -1,37 +1,55 @@ DEAL::FE = FE_SimplexP<2>(1) -DEAL::max error = 0.332362 -DEAL::ratio = 3.00877 -DEAL::max error = 0.100426 -DEAL::ratio = 3.30953 -DEAL::max error = 0.0271946 -DEAL::ratio = 3.69285 -DEAL::max error = 0.00747939 -DEAL::ratio = 3.63594 +DEAL::max error = 0.324312 +DEAL::ratio = 3.08345 +DEAL::max error = 0.0993124 +DEAL::ratio = 3.26557 +DEAL::max error = 0.0272879 +DEAL::ratio = 3.63943 +DEAL::max error = 0.00748538 +DEAL::ratio = 3.64549 DEAL::FE = FE_SimplexP<2>(2) -DEAL::max error = 0.0430400 -DEAL::ratio = 23.2342 -DEAL::max error = 0.0116055 -DEAL::ratio = 3.70859 -DEAL::max error = 0.00232162 -DEAL::ratio = 4.99887 -DEAL::max error = 0.000336599 -DEAL::ratio = 6.89729 +DEAL::max error = 0.0581409 +DEAL::ratio = 17.1996 +DEAL::max error = 0.0120766 +DEAL::ratio = 4.81435 +DEAL::max error = 0.00234218 +DEAL::ratio = 5.15614 +DEAL::max error = 0.000337327 +DEAL::ratio = 6.94335 +DEAL::FE = FE_SimplexP<3>(3) +DEAL::max error = 0.0812860 +DEAL::ratio = 12.3022 +DEAL::max error = 0.00669086 +DEAL::ratio = 12.1488 +DEAL::max error = 0.000438787 +DEAL::ratio = 15.2485 +DEAL::max error = 2.71681e-05 +DEAL::ratio = 16.1508 DEAL::FE = FE_SimplexP<3>(1) -DEAL::max error = 0.620494 -DEAL::ratio = 1.61162 -DEAL::max error = 0.259138 -DEAL::ratio = 2.39445 -DEAL::max error = 0.0619310 -DEAL::ratio = 4.18431 -DEAL::max error = 0.0151022 -DEAL::ratio = 4.10080 +DEAL::max error = 0.617024 +DEAL::ratio = 1.62068 +DEAL::max error = 0.261218 +DEAL::ratio = 2.36211 +DEAL::max error = 0.0620383 +DEAL::ratio = 4.21059 +DEAL::max error = 0.0151087 +DEAL::ratio = 4.10612 DEAL::FE = FE_SimplexP<3>(2) -DEAL::max error = 0.0550999 -DEAL::ratio = 18.1489 -DEAL::max error = 0.0466852 -DEAL::ratio = 1.18024 -DEAL::max error = 0.00502563 -DEAL::ratio = 9.28942 -DEAL::max error = 0.000810078 -DEAL::ratio = 6.20388 +DEAL::max error = 0.0563879 +DEAL::ratio = 17.7343 +DEAL::max error = 0.0452052 +DEAL::ratio = 1.24737 +DEAL::max error = 0.00506696 +DEAL::ratio = 8.92157 +DEAL::max error = 0.000811568 +DEAL::ratio = 6.24342 +DEAL::FE = FE_SimplexP<3>(3) +DEAL::max error = 0.0812860 +DEAL::ratio = 12.3022 +DEAL::max error = 0.00669086 +DEAL::ratio = 12.1488 +DEAL::max error = 0.000438787 +DEAL::ratio = 15.2485 +DEAL::max error = 2.71681e-05 +DEAL::ratio = 16.1508 diff --git a/tests/simplex/simplex_project_dg.cc b/tests/simplex/simplex_project_dg.cc index 1ec5d28a6f..ed5bac6a6d 100644 --- a/tests/simplex/simplex_project_dg.cc +++ b/tests/simplex/simplex_project_dg.cc @@ -105,7 +105,9 @@ main() test<2>(1); test<2>(2); + // test<2>(3); test<3>(1); test<3>(2); + // test<3>(3); } diff --git a/tests/simplex/simplex_project_dg.output b/tests/simplex/simplex_project_dg.output index 4bb911968d..9dae319d71 100644 --- a/tests/simplex/simplex_project_dg.output +++ b/tests/simplex/simplex_project_dg.output @@ -17,6 +17,15 @@ DEAL::max error = 0.000287343 DEAL::ratio = 7.91404 DEAL::max error = 3.60148e-05 DEAL::ratio = 7.97847 +DEAL::FE = FE_SimplexDGP<2>(3) +DEAL::max error = 0.00193694 +DEAL::ratio = 516.278 +DEAL::max error = 0.000124622 +DEAL::ratio = 15.5425 +DEAL::max error = 7.84536e-06 +DEAL::ratio = 15.8848 +DEAL::max error = 4.91221e-07 +DEAL::ratio = 15.9712 DEAL::FE = FE_SimplexDGP<3>(1) DEAL::max error = 0.208217 DEAL::ratio = 4.80269 @@ -35,3 +44,12 @@ DEAL::max error = 0.000842371 DEAL::ratio = 7.61471 DEAL::max error = 0.000106584 DEAL::ratio = 7.90336 +DEAL::FE = FE_SimplexDGP<3>(3) +DEAL::max error = 0.00566152 +DEAL::ratio = 176.631 +DEAL::max error = 0.000454735 +DEAL::ratio = 12.4501 +DEAL::max error = 3.01010e-05 +DEAL::ratio = 15.1069 +DEAL::max error = 1.90799e-06 +DEAL::ratio = 15.7763 -- 2.39.5