]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Plenty of new tests
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 11 Aug 2006 17:40:23 +0000 (17:40 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 11 Aug 2006 17:40:23 +0000 (17:40 +0000)
git-svn-id: https://svn.dealii.org/trunk@13675 0785d39b-7218-0410-832d-ea1e28bc413d

tests/deal.II/project_common.cc
tests/deal.II/project_dgp_01.cc
tests/deal.II/project_dgp_monomial_01.cc [new file with mode: 0644]
tests/deal.II/project_dgp_monomial_01/cmp/generic [new file with mode: 0644]
tests/deal.II/project_dgp_nonparametric_01.cc [new file with mode: 0644]
tests/deal.II/project_nedelec_01.cc [new file with mode: 0644]
tests/deal.II/project_nedelec_01/cmp/generic [new file with mode: 0644]
tests/deal.II/project_q_hierarchical_01.cc [new file with mode: 0644]
tests/deal.II/project_q_hierarchical_01/cmp/generic [new file with mode: 0644]
tests/deal.II/project_q_system_01.cc [new file with mode: 0644]
tests/deal.II/project_q_system_01/cmp/generic [new file with mode: 0644]

index faebd59cf692873b977b43459d736e1cc45d5675..12615f8a552f16b048e4331084495a7b5dbeeeda 100644 (file)
@@ -110,7 +110,7 @@ void test_no_hanging_nodes (const FiniteElement<dim> &fe,
   ConstraintMatrix constraints;
   constraints.close ();
 
-  Vector<double> interpolant (dof_handler.n_dofs());
+  Vector<double> projection (dof_handler.n_dofs());
   Vector<float>  error (triangulation.n_active_cells());
   for (unsigned int q=0; q<=p+2; ++q)
     {
@@ -119,21 +119,21 @@ void test_no_hanging_nodes (const FiniteElement<dim> &fe,
                            constraints,
                            QGauss<dim>(p+2),
                            F<dim> (q, fe.n_components()),
-                           interpolant);
+                           projection);
       
                                       // then compute the interpolation error
       VectorTools::integrate_difference (dof_handler,
-                                        interpolant,
+                                        projection,
                                         F<dim> (q, fe.n_components()),
                                         error,
                                         QGauss<dim>(std::max(p,q)+1),
                                         VectorTools::L2_norm);
       deallog << fe.get_name() << ", P_" << q
-             << ", rel. error=" << error.l2_norm() / interpolant.l2_norm()
+             << ", rel. error=" << error.l2_norm() / projection.l2_norm()
              << std::endl;
          
       if (q<=p)
-       Assert (error.l2_norm() < 1e-12*interpolant.l2_norm(),
+       Assert (error.l2_norm() <= 1e-12*projection.l2_norm(),
                ExcInternalError());
     }
 }
index fa9b0881d4ffaf8d8b3d8a86b133d9a30f7b305c..8fa4a9b23f9bf216d28132737f6931f8f8ce7d44 100644 (file)
@@ -12,7 +12,7 @@
 //----------------------------  project_dgp_01.cc  ---------------------------
 
 
-// check that VectorTools::project works for Q elements correctly on
+// check that VectorTools::project works for DGP elements correctly on
 // a uniformly refined mesh for functions of degree q
 
 char logname[] = "project_dgp_01/output";
diff --git a/tests/deal.II/project_dgp_monomial_01.cc b/tests/deal.II/project_dgp_monomial_01.cc
new file mode 100644 (file)
index 0000000..1a01629
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_dgp_monomial_01.cc  ---------------------------
+//    $Id: project_dgp_monomial_01.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_01.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_01/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=0; p<6-dim; ++p)
+    test_no_hanging_nodes (FE_DGPMonomial<dim>(p), p);
+}
diff --git a/tests/deal.II/project_dgp_monomial_01/cmp/generic b/tests/deal.II/project_dgp_monomial_01/cmp/generic
new file mode 100644 (file)
index 0000000..8a3ee90
--- /dev/null
@@ -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_nonparametric_01.cc b/tests/deal.II/project_dgp_nonparametric_01.cc
new file mode 100644 (file)
index 0000000..f1d0d4b
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_dgp_nonparametric_01.cc  ---------------------------
+//    $Id: project_dgp_nonparametric_01.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_01.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_01/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=0; p<6-dim; ++p)
+    test_no_hanging_nodes (FE_DGPNonparametric<dim>(p), p);
+}
diff --git a/tests/deal.II/project_nedelec_01.cc b/tests/deal.II/project_nedelec_01.cc
new file mode 100644 (file)
index 0000000..12682fb
--- /dev/null
@@ -0,0 +1,31 @@
+//----------------------------  project_nedelec_01.cc  ---------------------------
+//    $Id: project_nedelec_01.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_01.cc  ---------------------------
+
+
+// check that VectorTools::project works for Nedelec elements correctly on
+// a uniformly refined mesh for functions of degree q
+
+char logname[] = "project_nedelec_01/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  if (dim > 1)
+                                    // only p=1 implemented at present
+    for (unsigned int p=1; p<2; ++p)
+      test_no_hanging_nodes (FE_Nedelec<dim>(p), p);
+}
diff --git a/tests/deal.II/project_nedelec_01/cmp/generic b/tests/deal.II/project_nedelec_01/cmp/generic
new file mode 100644 (file)
index 0000000..6c31f57
--- /dev/null
@@ -0,0 +1,148 @@
+
+DEAL::n_dofs=9
+DEAL:cg::Starting value 0.342
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q<1>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.709
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q<1>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.12
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q<1>(1), P_2, rel. error=0.000336
+DEAL:cg::Starting value 1.57
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q<1>(1), P_3, rel. error=0.000737
+DEAL::n_dofs=17
+DEAL:cg::Starting value 0.262
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Q<1>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.544
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Q<1>(2), P_1, rel. error=0
+DEAL:cg::Starting value 0.869
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(2), P_2, rel. error=0
+DEAL:cg::Starting value 1.23
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(2), P_3, rel. error=1.24e-05
+DEAL:cg::Starting value 1.63
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(2), P_4, rel. error=3.41e-05
+DEAL::n_dofs=25
+DEAL:cg::Starting value 0.206
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(3), P_0, rel. error=0
+DEAL:cg::Starting value 0.429
+DEAL:cg::Convergence step 14 value 0
+DEAL::FE_Q<1>(3), P_1, rel. error=0
+DEAL:cg::Starting value 0.685
+DEAL:cg::Convergence step 14 value 0
+DEAL::FE_Q<1>(3), P_2, rel. error=0
+DEAL:cg::Starting value 0.974
+DEAL:cg::Convergence step 14 value 0
+DEAL::FE_Q<1>(3), P_3, rel. error=0
+DEAL:cg::Starting value 1.29
+DEAL:cg::Convergence step 14 value 0
+DEAL::FE_Q<1>(3), P_4, rel. error=1.77e-07
+DEAL:cg::Starting value 1.64
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(3), P_5, rel. error=6.01e-07
+DEAL::n_dofs=33
+DEAL:cg::Starting value 0.191
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_0, rel. error=0
+DEAL:cg::Starting value 0.398
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_1, rel. error=0
+DEAL:cg::Starting value 0.638
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<1>(4), P_2, rel. error=0
+DEAL:cg::Starting value 0.908
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<1>(4), P_3, rel. error=0
+DEAL:cg::Starting value 1.21
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<1>(4), P_4, rel. error=0
+DEAL:cg::Starting value 1.53
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<1>(4), P_5, rel. error=6.36e-09
+DEAL:cg::Starting value 1.88
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<1>(4), P_6, rel. error=2.50e-08
+DEAL::n_dofs=81
+DEAL:cg::Starting value 0.352
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<2>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.717
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<2>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.10
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<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<2>(1), P_3, rel. error=0.000198
+DEAL::n_dofs=289
+DEAL:cg::Starting value 0.206
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<2>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.421
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_1, rel. error=0
+DEAL:cg::Starting value 0.654
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_2, rel. error=0
+DEAL:cg::Starting value 0.905
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_3, rel. error=2.41e-06
+DEAL:cg::Starting value 1.17
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_4, rel. error=6.79e-06
+DEAL::n_dofs=625
+DEAL:cg::Starting value 0.127
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Q<2>(3), P_0, rel. error=0
+DEAL:cg::Starting value 0.261
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Q<2>(3), P_1, rel. error=0
+DEAL:cg::Starting value 0.406
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_2, rel. error=0
+DEAL:cg::Starting value 0.562
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_3, rel. error=0
+DEAL:cg::Starting value 0.730
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_4, rel. error=2.90e-08
+DEAL:cg::Starting value 0.908
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_5, rel. error=1.00e-07
+DEAL::n_dofs=729
+DEAL:cg::Starting value 0.241
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q<3>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.488
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q<3>(1), P_1, rel. error=0
+DEAL:cg::Starting value 0.744
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q<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<3>(1), P_3, rel. error=5.72e-05
+DEAL::n_dofs=4913
+DEAL:cg::Starting value 0.108
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q<3>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.219
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(2), P_1, rel. error=0
+DEAL:cg::Starting value 0.337
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(2), P_2, rel. error=0
+DEAL:cg::Starting value 0.461
+DEAL:cg::Convergence step 24 value 0
+DEAL::FE_Q<3>(2), P_3, rel. error=5.06e-07
+DEAL:cg::Starting value 0.592
+DEAL:cg::Convergence step 24 value 0
+DEAL::FE_Q<3>(2), P_4, rel. error=1.44e-06
diff --git a/tests/deal.II/project_q_hierarchical_01.cc b/tests/deal.II/project_q_hierarchical_01.cc
new file mode 100644 (file)
index 0000000..a703b9f
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_q_hierarchical_01.cc  ---------------------------
+//    $Id: project_q_hierarchical_01.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_01.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_01/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=1; p<6-dim; ++p)
+    test_no_hanging_nodes (FE_Q_Hierarchical<dim>(p), p);
+}
diff --git a/tests/deal.II/project_q_hierarchical_01/cmp/generic b/tests/deal.II/project_q_hierarchical_01/cmp/generic
new file mode 100644 (file)
index 0000000..15a4b99
--- /dev/null
@@ -0,0 +1,148 @@
+
+DEAL::n_dofs=9
+DEAL:cg::Starting value 0.342
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q_Hierarchical<1>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.709
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q_Hierarchical<1>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.12
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q_Hierarchical<1>(1), P_2, rel. error=0.000336
+DEAL:cg::Starting value 1.57
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_Q_Hierarchical<1>(1), P_3, rel. error=0.000737
+DEAL::n_dofs=17
+DEAL:cg::Starting value 0.416
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q_Hierarchical<1>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.862
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q_Hierarchical<1>(2), P_1, rel. error=0
+DEAL:cg::Starting value 1.37
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q_Hierarchical<1>(2), P_2, rel. error=0
+DEAL:cg::Starting value 1.93
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q_Hierarchical<1>(2), P_3, rel. error=1.64e-05
+DEAL:cg::Starting value 2.53
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q_Hierarchical<1>(2), P_4, rel. error=4.45e-05
+DEAL::n_dofs=25
+DEAL:cg::Starting value 0.416
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q_Hierarchical<1>(3), P_0, rel. error=0
+DEAL:cg::Starting value 0.862
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q_Hierarchical<1>(3), P_1, rel. error=0
+DEAL:cg::Starting value 1.37
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q_Hierarchical<1>(3), P_2, rel. error=0
+DEAL:cg::Starting value 1.93
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Q_Hierarchical<1>(3), P_3, rel. error=0
+DEAL:cg::Starting value 2.53
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Q_Hierarchical<1>(3), P_4, rel. error=2.75e-07
+DEAL:cg::Starting value 3.18
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Q_Hierarchical<1>(3), P_5, rel. error=9.15e-07
+DEAL::n_dofs=33
+DEAL:cg::Starting value 0.503
+DEAL:cg::Convergence step 33 value 0
+DEAL::FE_Q_Hierarchical<1>(4), P_0, rel. error=0
+DEAL:cg::Starting value 1.04
+DEAL:cg::Convergence step 33 value 0
+DEAL::FE_Q_Hierarchical<1>(4), P_1, rel. error=0
+DEAL:cg::Starting value 1.66
+DEAL:cg::Convergence step 33 value 0
+DEAL::FE_Q_Hierarchical<1>(4), P_2, rel. error=0
+DEAL:cg::Starting value 2.35
+DEAL:cg::Convergence step 34 value 0
+DEAL::FE_Q_Hierarchical<1>(4), P_3, rel. error=0
+DEAL:cg::Starting value 3.09
+DEAL:cg::Convergence step 34 value 0
+DEAL::FE_Q_Hierarchical<1>(4), P_4, rel. error=0
+DEAL:cg::Starting value 3.89
+DEAL:cg::Convergence step 34 value 0
+DEAL::FE_Q_Hierarchical<1>(4), P_5, rel. error=1.10e-08
+DEAL:cg::Starting value 4.74
+DEAL:cg::Convergence step 34 value 0
+DEAL::FE_Q_Hierarchical<1>(4), P_6, rel. error=4.22e-08
+DEAL::n_dofs=81
+DEAL:cg::Starting value 0.352
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q_Hierarchical<2>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.717
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q_Hierarchical<2>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.10
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q_Hierarchical<2>(1), P_2, rel. error=8.80e-05
+DEAL:cg::Starting value 1.51
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q_Hierarchical<2>(1), P_3, rel. error=0.000198
+DEAL::n_dofs=289
+DEAL:cg::Starting value 0.518
+DEAL:cg::Convergence step 110 value 0
+DEAL::FE_Q_Hierarchical<2>(2), P_0, rel. error=0
+DEAL:cg::Starting value 1.06
+DEAL:cg::Convergence step 111 value 0
+DEAL::FE_Q_Hierarchical<2>(2), P_1, rel. error=0
+DEAL:cg::Starting value 1.64
+DEAL:cg::Convergence step 112 value 0
+DEAL::FE_Q_Hierarchical<2>(2), P_2, rel. error=0
+DEAL:cg::Starting value 2.25
+DEAL:cg::Convergence step 112 value 0
+DEAL::FE_Q_Hierarchical<2>(2), P_3, rel. error=4.42e-06
+DEAL:cg::Starting value 2.90
+DEAL:cg::Convergence step 114 value 0
+DEAL::FE_Q_Hierarchical<2>(2), P_4, rel. error=1.22e-05
+DEAL::n_dofs=625
+DEAL:cg::Starting value 0.518
+DEAL:cg::Convergence step 116 value 0
+DEAL::FE_Q_Hierarchical<2>(3), P_0, rel. error=0
+DEAL:cg::Starting value 1.06
+DEAL:cg::Convergence step 117 value 0
+DEAL::FE_Q_Hierarchical<2>(3), P_1, rel. error=0
+DEAL:cg::Starting value 1.64
+DEAL:cg::Convergence step 118 value 0
+DEAL::FE_Q_Hierarchical<2>(3), P_2, rel. error=0
+DEAL:cg::Starting value 2.25
+DEAL:cg::Convergence step 118 value 0
+DEAL::FE_Q_Hierarchical<2>(3), P_3, rel. error=0
+DEAL:cg::Starting value 2.90
+DEAL:cg::Convergence step 119 value 0
+DEAL::FE_Q_Hierarchical<2>(3), P_4, rel. error=7.58e-08
+DEAL:cg::Starting value 3.57
+DEAL:cg::Convergence step 123 value 0
+DEAL::FE_Q_Hierarchical<2>(3), P_5, rel. error=2.57e-07
+DEAL::n_dofs=729
+DEAL:cg::Starting value 0.241
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q_Hierarchical<3>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.488
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q_Hierarchical<3>(1), P_1, rel. error=0
+DEAL:cg::Starting value 0.744
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q_Hierarchical<3>(1), P_2, rel. error=2.51e-05
+DEAL:cg::Starting value 1.01
+DEAL:cg::Convergence step 22 value 0
+DEAL::FE_Q_Hierarchical<3>(1), P_3, rel. error=5.72e-05
+DEAL::n_dofs=4913
+DEAL:cg::Starting value 0.431
+DEAL:cg::Convergence step 393 value 0
+DEAL::FE_Q_Hierarchical<3>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.874
+DEAL:cg::Convergence step 401 value 0
+DEAL::FE_Q_Hierarchical<3>(2), P_1, rel. error=0
+DEAL:cg::Starting value 1.34
+DEAL:cg::Convergence step 407 value 0
+DEAL::FE_Q_Hierarchical<3>(2), P_2, rel. error=0
+DEAL:cg::Starting value 1.82
+DEAL:cg::Convergence step 411 value 0
+DEAL::FE_Q_Hierarchical<3>(2), P_3, rel. error=1.27e-06
+DEAL:cg::Starting value 2.32
+DEAL:cg::Convergence step 415 value 0
+DEAL::FE_Q_Hierarchical<3>(2), P_4, rel. error=3.57e-06
diff --git a/tests/deal.II/project_q_system_01.cc b/tests/deal.II/project_q_system_01.cc
new file mode 100644 (file)
index 0000000..f2b5661
--- /dev/null
@@ -0,0 +1,31 @@
+//----------------------------  project_q_system_01.cc  ---------------------------
+//    $Id: project_q_system_01.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_01.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_01/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=1; p<5-dim; ++p)
+    test_no_hanging_nodes (FESystem<dim>(FE_Q<dim>(p), 1,
+                                        FE_DGQ<dim>(p+1), 1),
+                          p);
+}
diff --git a/tests/deal.II/project_q_system_01/cmp/generic b/tests/deal.II/project_q_system_01/cmp/generic
new file mode 100644 (file)
index 0000000..dd07c4e
--- /dev/null
@@ -0,0 +1,91 @@
+
+DEAL::n_dofs=33
+DEAL:cg::Starting value 0.606
+DEAL:cg::Convergence step 11 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_0, rel. error=0
+DEAL:cg::Starting value 1.23
+DEAL:cg::Convergence step 11 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_1, rel. error=0
+DEAL:cg::Starting value 1.91
+DEAL:cg::Convergence step 11 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_2, rel. error=0.000109
+DEAL:cg::Starting value 2.62
+DEAL:cg::Convergence step 11 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_3, rel. error=0.000256
+DEAL::n_dofs=49
+DEAL:cg::Starting value 0.474
+DEAL:cg::Convergence step 13 value 0
+DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_0, rel. error=0
+DEAL:cg::Starting value 0.967
+DEAL:cg::Convergence step 13 value 0
+DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_1, rel. error=0
+DEAL:cg::Starting value 1.50
+DEAL:cg::Convergence step 13 value 0
+DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_2, rel. error=0
+DEAL:cg::Starting value 2.07
+DEAL:cg::Convergence step 14 value 0
+DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_3, rel. error=4.87e-06
+DEAL:cg::Starting value 2.67
+DEAL:cg::Convergence step 14 value 0
+DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_4, rel. error=1.41e-05
+DEAL::n_dofs=65
+DEAL:cg::Starting value 0.429
+DEAL:cg::Convergence step 14 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_0, rel. error=0
+DEAL:cg::Starting value 0.872
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_1, rel. error=0
+DEAL:cg::Starting value 1.35
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_2, rel. error=0
+DEAL:cg::Starting value 1.86
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_3, rel. error=0
+DEAL:cg::Starting value 2.39
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_4, rel. error=7.70e-08
+DEAL:cg::Starting value 2.96
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_5, rel. error=2.71e-07
+DEAL::n_dofs=657
+DEAL:cg::Starting value 0.470
+DEAL:cg::Convergence step 17 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_0, rel. error=0
+DEAL:cg::Starting value 0.955
+DEAL:cg::Convergence step 17 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_1, rel. error=0
+DEAL:cg::Starting value 1.46
+DEAL:cg::Convergence step 17 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_2, rel. error=2.06e-05
+DEAL:cg::Starting value 1.99
+DEAL:cg::Convergence step 17 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_3, rel. error=4.87e-05
+DEAL::n_dofs=1313
+DEAL:cg::Starting value 0.284
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_0, rel. error=0
+DEAL:cg::Starting value 0.576
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_1, rel. error=0
+DEAL:cg::Starting value 0.887
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_2, rel. error=0
+DEAL:cg::Starting value 1.22
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_3, rel. error=7.98e-07
+DEAL:cg::Starting value 1.56
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_4, rel. error=2.32e-06
+DEAL::n_dofs=14553
+DEAL:cg::Starting value 0.279
+DEAL:cg::Convergence step 22 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0
+DEAL:cg::Starting value 0.564
+DEAL:cg::Convergence step 23 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0
+DEAL:cg::Starting value 0.860
+DEAL:cg::Convergence step 23 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=3.98e-06
+DEAL:cg::Starting value 1.16
+DEAL:cg::Convergence step 23 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=9.45e-06

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.