]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Add new tests that check projection onto meshes with hanging node constraints and...
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 13 Aug 2006 03:47:11 +0000 (03:47 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 13 Aug 2006 03:47:11 +0000 (03:47 +0000)
git-svn-id: https://svn.dealii.org/trunk@13691 0785d39b-7218-0410-832d-ea1e28bc413d

23 files changed:
tests/deal.II/project_abf_02.cc [new file with mode: 0644]
tests/deal.II/project_common.cc
tests/deal.II/project_dgp_02.cc [new file with mode: 0644]
tests/deal.II/project_dgp_02/cmp/generic [new file with mode: 0644]
tests/deal.II/project_dgp_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 [new file with mode: 0644]
tests/deal.II/project_dgp_monomial_02.cc [new file with mode: 0644]
tests/deal.II/project_dgp_monomial_02/cmp/generic [new file with mode: 0644]
tests/deal.II/project_dgp_monomial_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 [new file with mode: 0644]
tests/deal.II/project_dgp_nonparametric_02.cc [new file with mode: 0644]
tests/deal.II/project_dgq_02.cc [new file with mode: 0644]
tests/deal.II/project_dgq_02/cmp/generic [new file with mode: 0644]
tests/deal.II/project_dgq_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 [new file with mode: 0644]
tests/deal.II/project_nedelec_02.cc [new file with mode: 0644]
tests/deal.II/project_nedelec_02/cmp/generic [new file with mode: 0644]
tests/deal.II/project_q_02.cc [new file with mode: 0644]
tests/deal.II/project_q_02/cmp/generic [new file with mode: 0644]
tests/deal.II/project_q_hierarchical_01/cmp/generic
tests/deal.II/project_q_hierarchical_02.cc [new file with mode: 0644]
tests/deal.II/project_q_hierarchical_02/cmp/generic [new file with mode: 0644]
tests/deal.II/project_q_hierarchical_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 [new file with mode: 0644]
tests/deal.II/project_q_system_02.cc [new file with mode: 0644]
tests/deal.II/project_q_system_02/cmp/generic [new file with mode: 0644]
tests/deal.II/project_rt_02.cc [new file with mode: 0644]

diff --git a/tests/deal.II/project_abf_02.cc b/tests/deal.II/project_abf_02.cc
new file mode 100644 (file)
index 0000000..426af62
--- /dev/null
@@ -0,0 +1,30 @@
+//----------------------------  project_abf_02.cc  ---------------------------
+//    $Id: project_abf_02.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_02.cc  ---------------------------
+
+
+// check that VectorTools::project works for ABF elements correctly on
+// a uniformly refined mesh for functions of degree q
+
+char logname[] = "project_abf_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  if (dim > 1)
+    for (unsigned int p=0; p<6-dim; ++p)
+      test_with_hanging_nodes (FE_ABF<dim>(p), p);
+}
index 2b01bfac8017084db0db3abf79bddc5833ef8447..77be2216e4c38756746db4046f947a98edc51a41 100644 (file)
@@ -93,26 +93,20 @@ class F :  public Function<dim>
 
 
 
-// check the given element of polynomial order p. the last parameter, if
-// given, denotes a gap in convergence order; for example, the Nedelec element
-// of polynomial degree p has normal components of degree p-1 and therefore
-// can only represent polynomials of degree p-1 exactly. the gap is then 1.
 template <int dim>
-void test_no_hanging_nodes (const FiniteElement<dim> &fe,
-                           const unsigned int        p,
-                           const unsigned int        order_difference = 0)
-{
-  Triangulation<dim>     triangulation;
-  GridGenerator::hyper_cube (triangulation);
-  triangulation.refine_global (3);
-
+void do_project (const Triangulation<dim> &triangulation,
+                const FiniteElement<dim> &fe,
+                const unsigned int        p,
+                const unsigned int        order_difference)
+{  
   DoFHandler<dim>        dof_handler(triangulation);
   dof_handler.distribute_dofs (fe);
 
   deallog << "n_dofs=" << dof_handler.n_dofs() << std::endl;
 
-                                  // there are no constraints here
   ConstraintMatrix constraints;
+  DoFTools::make_hanging_node_constraints (dof_handler,
+                                          constraints);
   constraints.close ();
 
   Vector<double> projection (dof_handler.n_dofs());
@@ -125,6 +119,10 @@ void test_no_hanging_nodes (const FiniteElement<dim> &fe,
                            QGauss<dim>(p+2),
                            F<dim> (q, fe.n_components()),
                            projection);
+                                      // just to make sure it doesn't get
+                                      // forgotten: handle hanging node
+                                      // constraints
+      constraints.distribute (projection);
       
                                       // then compute the interpolation error
       VectorTools::integrate_difference (dof_handler,
@@ -145,6 +143,67 @@ void test_no_hanging_nodes (const FiniteElement<dim> &fe,
 
 
 
+// check the given element of polynomial order p. the last parameter, if
+// given, denotes a gap in convergence order; for example, the Nedelec element
+// of polynomial degree p has normal components of degree p-1 and therefore
+// can only represent polynomials of degree p-1 exactly. the gap is then 1.
+template <int dim>
+void test_no_hanging_nodes (const FiniteElement<dim> &fe,
+                           const unsigned int        p,
+                           const unsigned int        order_difference = 0)
+{
+  Triangulation<dim>     triangulation;
+  GridGenerator::hyper_cube (triangulation);
+  triangulation.refine_global (3);
+
+  do_project (triangulation, fe, p, order_difference);
+}
+
+
+
+// same test as above, but this time with a mesh that has hanging nodes
+template <int dim>
+void test_with_hanging_nodes (const FiniteElement<dim> &fe,
+                             const unsigned int        p,
+                             const unsigned int        order_difference = 0)
+{
+  Triangulation<dim>     triangulation;
+  GridGenerator::hyper_cube (triangulation);
+  triangulation.refine_global (1);
+  triangulation.begin_active()->set_refine_flag ();
+  triangulation.execute_coarsening_and_refinement ();
+  triangulation.refine_global (1);
+  
+  do_project (triangulation, fe, p, order_difference);
+}
+
+
+
+// test with a 3d grid that has cells with face_orientation==false and hanging
+// nodes. this trips up all sorts of pieces of code, for example there was a
+// 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
+template <int dim>
+void test_with_wrong_face_orientation (const FiniteElement<dim> &fe,
+                                      const unsigned int        p,
+                                      const unsigned int        order_difference = 0)
+{
+  if (dim < 3)
+    return;
+  
+  Triangulation<dim>     triangulation;
+  GridGenerator::hyper_cube (triangulation);
+  triangulation.refine_global (1);
+  triangulation.begin_active()->set_refine_flag ();
+  triangulation.execute_coarsening_and_refinement ();
+  triangulation.refine_global (1);
+  
+  do_project (triangulation, fe, p, order_difference);
+}
+
+
+
 int main ()
 {
   std::ofstream logfile(logname);
diff --git a/tests/deal.II/project_dgp_02.cc b/tests/deal.II/project_dgp_02.cc
new file mode 100644 (file)
index 0000000..ff4dda5
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_dgp_02.cc  ---------------------------
+//    $Id: project_dgp_02.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_02.cc  ---------------------------
+
+
+// check that VectorTools::project works for DGP elements correctly on
+// a uniformly refined mesh for functions of degree q
+
+char logname[] = "project_dgp_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=0; p<6-dim; ++p)
+    test_with_hanging_nodes (FE_DGP<dim>(p), p);
+}
diff --git a/tests/deal.II/project_dgp_02/cmp/generic b/tests/deal.II/project_dgp_02/cmp/generic
new file mode 100644 (file)
index 0000000..9d7503b
--- /dev/null
@@ -0,0 +1,178 @@
+
+DEAL::n_dofs=6
+DEAL:cg::Starting value 0.433
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<1>(0), P_0, rel. error=0
+DEAL:cg::Starting value 0.967
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<1>(0), P_1, rel. error=0.0244
+DEAL:cg::Starting value 1.60
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<1>(0), P_2, rel. error=0.0487
+DEAL::n_dofs=12
+DEAL:cg::Starting value 0.433
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<1>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.968
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.60
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<1>(1), P_2, rel. error=0.00142
+DEAL:cg::Starting value 2.33
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(1), P_3, rel. error=0.00405
+DEAL::n_dofs=18
+DEAL:cg::Starting value 0.433
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<1>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.968
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(2), P_1, rel. error=0
+DEAL:cg::Starting value 1.60
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<1>(2), P_2, rel. error=0
+DEAL:cg::Starting value 2.33
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(2), P_3, rel. error=8.49e-05
+DEAL:cg::Starting value 3.12
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(2), P_4, rel. error=0.000309
+DEAL::n_dofs=24
+DEAL:cg::Starting value 0.433
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(3), P_0, rel. error=0
+DEAL:cg::Starting value 0.968
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(3), P_1, rel. error=0
+DEAL:cg::Starting value 1.60
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(3), P_2, rel. error=0
+DEAL:cg::Starting value 2.33
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(3), P_3, rel. error=0
+DEAL:cg::Starting value 3.12
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(3), P_4, rel. error=5.07e-06
+DEAL:cg::Starting value 3.98
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(3), P_5, rel. error=2.23e-05
+DEAL::n_dofs=30
+DEAL:cg::Starting value 0.433
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(4), P_0, rel. error=0
+DEAL:cg::Starting value 0.968
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(4), P_1, rel. error=0
+DEAL:cg::Starting value 1.60
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(4), P_2, rel. error=0
+DEAL:cg::Starting value 2.33
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(4), P_3, rel. error=0
+DEAL:cg::Starting value 3.12
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(4), P_4, rel. error=0
+DEAL:cg::Starting value 3.98
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(4), P_5, rel. error=3.03e-07
+DEAL:cg::Starting value 4.91
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<1>(4), P_6, rel. error=1.57e-06
+DEAL::n_dofs=28
+DEAL:cg::Starting value 0.676
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(0), P_0, rel. error=0
+DEAL:cg::Starting value 1.45
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<2>(0), P_1, rel. error=0.00998
+DEAL:cg::Starting value 2.31
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(0), P_2, rel. error=0.0192
+DEAL::n_dofs=84
+DEAL:cg::Starting value 0.676
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGP<2>(1), P_0, rel. error=0
+DEAL:cg::Starting value 1.45
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(1), P_1, rel. error=0
+DEAL:cg::Starting value 2.32
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(1), P_2, rel. error=0.000642
+DEAL:cg::Starting value 3.27
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(1), P_3, rel. error=0.00167
+DEAL::n_dofs=168
+DEAL:cg::Starting value 0.676
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(2), P_0, rel. error=0
+DEAL:cg::Starting value 1.45
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(2), P_1, rel. error=0
+DEAL:cg::Starting value 2.32
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(2), P_2, rel. error=0
+DEAL:cg::Starting value 3.27
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(2), P_3, rel. error=4.03e-05
+DEAL:cg::Starting value 4.28
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(2), P_4, rel. error=0.000131
+DEAL::n_dofs=280
+DEAL:cg::Starting value 0.676
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(3), P_0, rel. error=0
+DEAL:cg::Starting value 1.45
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(3), P_1, rel. error=0
+DEAL:cg::Starting value 2.32
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(3), P_2, rel. error=0
+DEAL:cg::Starting value 3.27
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(3), P_3, rel. error=0
+DEAL:cg::Starting value 4.28
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(3), P_4, rel. error=2.49e-06
+DEAL:cg::Starting value 5.36
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<2>(3), P_5, rel. error=9.67e-06
+DEAL::n_dofs=120
+DEAL:cg::Starting value 0.708
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(0), P_0, rel. error=0
+DEAL:cg::Starting value 1.48
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(0), P_1, rel. error=0.00433
+DEAL:cg::Starting value 2.30
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(0), P_2, rel. error=0.00820
+DEAL::n_dofs=480
+DEAL:cg::Starting value 0.708
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(1), P_0, rel. error=0
+DEAL:cg::Starting value 1.48
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(1), P_1, rel. error=0
+DEAL:cg::Starting value 2.31
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(1), P_2, rel. error=0.000288
+DEAL:cg::Starting value 3.20
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(1), P_3, rel. error=0.000729
+DEAL::n_dofs=1200
+DEAL:cg::Starting value 0.708
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(2), P_0, rel. error=0
+DEAL:cg::Starting value 1.48
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(2), P_1, rel. error=0
+DEAL:cg::Starting value 2.31
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(2), P_2, rel. error=0
+DEAL:cg::Starting value 3.20
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(2), P_3, rel. error=1.85e-05
+DEAL:cg::Starting value 4.14
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGP<3>(2), P_4, rel. error=5.78e-05
diff --git a/tests/deal.II/project_dgp_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgp_02/cmp/x86_64-unknown-linux-gnu+gcc4.1
new file mode 100644 (file)
index 0000000..6e360be
--- /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_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_02.cc b/tests/deal.II/project_dgp_monomial_02.cc
new file mode 100644 (file)
index 0000000..3a55733
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_dgp_monomial_02.cc  ---------------------------
+//    $Id: project_dgp_monomial_02.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_02.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_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=0; p<6-dim; ++p)
+    test_with_hanging_nodes (FE_DGPMonomial<dim>(p), p);
+}
diff --git a/tests/deal.II/project_dgp_monomial_02/cmp/generic b/tests/deal.II/project_dgp_monomial_02/cmp/generic
new file mode 100644 (file)
index 0000000..680c7eb
--- /dev/null
@@ -0,0 +1,178 @@
+
+DEAL::n_dofs=6
+DEAL:cg::Starting value 0.433
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGPMonomial<1>(0), P_0, rel. error=0
+DEAL:cg::Starting value 0.967
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGPMonomial<1>(0), P_1, rel. error=0.0244
+DEAL:cg::Starting value 1.60
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGPMonomial<1>(0), P_2, rel. error=0.0487
+DEAL::n_dofs=12
+DEAL:cg::Starting value 0.484
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGPMonomial<1>(1), P_0, rel. error=0
+DEAL:cg::Starting value 1.09
+DEAL:cg::Convergence step 3 value 0
+DEAL::FE_DGPMonomial<1>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.81
+DEAL:cg::Convergence step 3 value 0
+DEAL::FE_DGPMonomial<1>(1), P_2, rel. error=0.00159
+DEAL:cg::Starting value 2.62
+DEAL:cg::Convergence step 3 value 0
+DEAL::FE_DGPMonomial<1>(1), P_3, rel. error=0.00451
+DEAL::n_dofs=18
+DEAL:cg::Starting value 0.505
+DEAL:cg::Convergence step 4 value 0
+DEAL::FE_DGPMonomial<1>(2), P_0, rel. error=0
+DEAL:cg::Starting value 1.14
+DEAL:cg::Convergence step 4 value 0
+DEAL::FE_DGPMonomial<1>(2), P_1, rel. error=0
+DEAL:cg::Starting value 1.90
+DEAL:cg::Convergence step 5 value 0
+DEAL::FE_DGPMonomial<1>(2), P_2, rel. error=0
+DEAL:cg::Starting value 2.76
+DEAL:cg::Convergence step 5 value 0
+DEAL::FE_DGPMonomial<1>(2), P_3, rel. error=9.83e-05
+DEAL:cg::Starting value 3.70
+DEAL:cg::Convergence step 5 value 0
+DEAL::FE_DGPMonomial<1>(2), P_4, rel. error=0.000362
+DEAL::n_dofs=24
+DEAL:cg::Starting value 0.517
+DEAL:cg::Convergence step 6 value 0
+DEAL::FE_DGPMonomial<1>(3), P_0, rel. error=0
+DEAL:cg::Starting value 1.17
+DEAL:cg::Convergence step 6 value 0
+DEAL::FE_DGPMonomial<1>(3), P_1, rel. error=0
+DEAL:cg::Starting value 1.95
+DEAL:cg::Convergence step 6 value 0
+DEAL::FE_DGPMonomial<1>(3), P_2, rel. error=0
+DEAL:cg::Starting value 2.84
+DEAL:cg::Convergence step 6 value 0
+DEAL::FE_DGPMonomial<1>(3), P_3, rel. error=0
+DEAL:cg::Starting value 3.82
+DEAL:cg::Convergence step 6 value 0
+DEAL::FE_DGPMonomial<1>(3), P_4, rel. error=5.96e-06
+DEAL:cg::Starting value 4.88
+DEAL:cg::Convergence step 6 value 0
+DEAL::FE_DGPMonomial<1>(3), P_5, rel. error=2.64e-05
+DEAL::n_dofs=30
+DEAL:cg::Starting value 0.524
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_DGPMonomial<1>(4), P_0, rel. error=0
+DEAL:cg::Starting value 1.19
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_DGPMonomial<1>(4), P_1, rel. error=0
+DEAL:cg::Starting value 1.98
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_DGPMonomial<1>(4), P_2, rel. error=0
+DEAL:cg::Starting value 2.89
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_DGPMonomial<1>(4), P_3, rel. error=0
+DEAL:cg::Starting value 3.89
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_DGPMonomial<1>(4), P_4, rel. error=0
+DEAL:cg::Starting value 4.98
+DEAL:cg::Convergence step 8 value 0
+DEAL::FE_DGPMonomial<1>(4), P_5, rel. error=3.58e-07
+DEAL:cg::Starting value 6.14
+DEAL:cg::Convergence step 9 value 0
+DEAL::FE_DGPMonomial<1>(4), P_6, rel. error=1.85e-06
+DEAL::n_dofs=28
+DEAL:cg::Starting value 0.676
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGPMonomial<2>(0), P_0, rel. error=0
+DEAL:cg::Starting value 1.45
+DEAL:cg::Convergence step 1 value 0
+DEAL::FE_DGPMonomial<2>(0), P_1, rel. error=0.00998
+DEAL:cg::Starting value 2.31
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGPMonomial<2>(0), P_2, rel. error=0.0192
+DEAL::n_dofs=84
+DEAL:cg::Starting value 0.828
+DEAL:cg::Convergence step 4 value 0
+DEAL::FE_DGPMonomial<2>(1), P_0, rel. error=0
+DEAL:cg::Starting value 1.79
+DEAL:cg::Convergence step 4 value 0
+DEAL::FE_DGPMonomial<2>(1), P_1, rel. error=0
+DEAL:cg::Starting value 2.86
+DEAL:cg::Convergence step 4 value 0
+DEAL::FE_DGPMonomial<2>(1), P_2, rel. error=0.000743
+DEAL:cg::Starting value 4.03
+DEAL:cg::Convergence step 4 value 0
+DEAL::FE_DGPMonomial<2>(1), P_3, rel. error=0.00198
+DEAL::n_dofs=168
+DEAL:cg::Starting value 0.903
+DEAL:cg::Convergence step 10 value 0
+DEAL::FE_DGPMonomial<2>(2), P_0, rel. error=0
+DEAL:cg::Starting value 1.96
+DEAL:cg::Convergence step 10 value 0
+DEAL::FE_DGPMonomial<2>(2), P_1, rel. error=0
+DEAL:cg::Starting value 3.15
+DEAL:cg::Convergence step 10 value 0
+DEAL::FE_DGPMonomial<2>(2), P_2, rel. error=0
+DEAL:cg::Starting value 4.44
+DEAL:cg::Convergence step 10 value 0
+DEAL::FE_DGPMonomial<2>(2), P_3, rel. error=4.86e-05
+DEAL:cg::Starting value 5.83
+DEAL:cg::Convergence step 10 value 0
+DEAL::FE_DGPMonomial<2>(2), P_4, rel. error=0.000162
+DEAL::n_dofs=280
+DEAL:cg::Starting value 0.948
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_DGPMonomial<2>(3), P_0, rel. error=0
+DEAL:cg::Starting value 2.06
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_DGPMonomial<2>(3), P_1, rel. error=0
+DEAL:cg::Starting value 3.32
+DEAL:cg::Convergence step 20 value 0
+DEAL::FE_DGPMonomial<2>(3), P_2, rel. error=0
+DEAL:cg::Starting value 4.69
+DEAL:cg::Convergence step 20 value 0
+DEAL::FE_DGPMonomial<2>(3), P_3, rel. error=0
+DEAL:cg::Starting value 6.16
+DEAL:cg::Convergence step 20 value 0
+DEAL::FE_DGPMonomial<2>(3), P_4, rel. error=3.09e-06
+DEAL:cg::Starting value 7.73
+DEAL:cg::Convergence step 20 value 0
+DEAL::FE_DGPMonomial<2>(3), P_5, rel. error=1.22e-05
+DEAL::n_dofs=120
+DEAL:cg::Starting value 0.708
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGPMonomial<3>(0), P_0, rel. error=0
+DEAL:cg::Starting value 1.48
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGPMonomial<3>(0), P_1, rel. error=0.00433
+DEAL:cg::Starting value 2.30
+DEAL:cg::Convergence step 2 value 0
+DEAL::FE_DGPMonomial<3>(0), P_2, rel. error=0.00820
+DEAL::n_dofs=480
+DEAL:cg::Starting value 0.936
+DEAL:cg::Convergence step 6 value 0
+DEAL::FE_DGPMonomial<3>(1), P_0, rel. error=0
+DEAL:cg::Starting value 1.96
+DEAL:cg::Convergence step 5 value 0
+DEAL::FE_DGPMonomial<3>(1), P_1, rel. error=0
+DEAL:cg::Starting value 3.08
+DEAL:cg::Convergence step 5 value 0
+DEAL::FE_DGPMonomial<3>(1), P_2, rel. error=0.000339
+DEAL:cg::Starting value 4.26
+DEAL:cg::Convergence step 7 value 0
+DEAL::FE_DGPMonomial<3>(1), P_3, rel. error=0.000889
+DEAL::n_dofs=1200
+DEAL:cg::Starting value 1.07
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_DGPMonomial<3>(2), P_0, rel. error=0
+DEAL:cg::Starting value 2.25
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_DGPMonomial<3>(2), P_1, rel. error=0
+DEAL:cg::Starting value 3.53
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_DGPMonomial<3>(2), P_2, rel. error=0
+DEAL:cg::Starting value 4.90
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_DGPMonomial<3>(2), P_3, rel. error=2.27e-05
+DEAL:cg::Starting value 6.35
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_DGPMonomial<3>(2), P_4, rel. error=7.36e-05
diff --git a/tests/deal.II/project_dgp_monomial_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgp_monomial_02/cmp/x86_64-unknown-linux-gnu+gcc4.1
new file mode 100644 (file)
index 0000000..f0ec21f
--- /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 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_02.cc b/tests/deal.II/project_dgp_nonparametric_02.cc
new file mode 100644 (file)
index 0000000..f32cb61
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_dgp_nonparametric_02.cc  ---------------------------
+//    $Id: project_dgp_nonparametric_02.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_02.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_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=0; p<6-dim; ++p)
+    test_with_hanging_nodes (FE_DGPNonparametric<dim>(p), p);
+}
diff --git a/tests/deal.II/project_dgq_02.cc b/tests/deal.II/project_dgq_02.cc
new file mode 100644 (file)
index 0000000..c685340
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_dgq_02.cc  ---------------------------
+//    $Id: project_dgq_02.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_02.cc  ---------------------------
+
+
+// check that VectorTools::project works for DGQ elements correctly on
+// a uniformly refined mesh for functions of degree q
+
+char logname[] = "project_dgq_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=0; p<6-dim; ++p)
+    test_with_hanging_nodes (FE_DGQ<dim>(p), p);
+}
diff --git a/tests/deal.II/project_dgq_02/cmp/generic b/tests/deal.II/project_dgq_02/cmp/generic
new file mode 100644 (file)
index 0000000..563a22f
--- /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_DGQ<1>(0), P_0, rel. error=0
+DEAL:cg::Starting value 0.736
+DEAL:cg::Convergence step 1 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 1 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_dgq_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_dgq_02/cmp/x86_64-unknown-linux-gnu+gcc4.1
new file mode 100644 (file)
index 0000000..802c43f
--- /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_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_02.cc b/tests/deal.II/project_nedelec_02.cc
new file mode 100644 (file)
index 0000000..c31b95a
--- /dev/null
@@ -0,0 +1,31 @@
+//----------------------------  project_nedelec_02.cc  ---------------------------
+//    $Id: project_nedelec_02.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_02.cc  ---------------------------
+
+
+// check that VectorTools::project works for Nedelec elements correctly on
+// a uniformly refined mesh for functions of degree q
+
+char logname[] = "project_nedelec_02/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_with_hanging_nodes (FE_Nedelec<dim>(p), p, 1);
+}
diff --git a/tests/deal.II/project_nedelec_02/cmp/generic b/tests/deal.II/project_nedelec_02/cmp/generic
new file mode 100644 (file)
index 0000000..3cdf4b0
--- /dev/null
@@ -0,0 +1,21 @@
+
+DEAL::n_dofs=72
+DEAL:cg::Starting value 5.43
+DEAL:cg::Convergence step 11 value 0
+DEAL::FE_Nedelec<2>(1), P_0, rel. error=0
+DEAL:cg::Starting value 11.0
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Nedelec<2>(1), P_1, rel. error=0.0217
+DEAL:cg::Starting value 16.7
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Nedelec<2>(1), P_2, rel. error=0.0391
+DEAL::n_dofs=576
+DEAL:cg::Starting value 6.85
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Nedelec<3>(1), P_0, rel. error=0
+DEAL:cg::Starting value 13.9
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Nedelec<3>(1), P_1, rel. error=0.00599
+DEAL:cg::Starting value 21.0
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Nedelec<3>(1), P_2, rel. error=0.0106
diff --git a/tests/deal.II/project_q_02.cc b/tests/deal.II/project_q_02.cc
new file mode 100644 (file)
index 0000000..87c55c4
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_q_02.cc  ---------------------------
+//    $Id: project_q_02.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_02.cc  ---------------------------
+
+
+// check that VectorTools::project works for Q elements correctly on
+// a uniformly refined mesh for functions of degree q
+
+char logname[] = "project_q_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=1; p<6-dim; ++p)
+    test_with_hanging_nodes (FE_Q<dim>(p), p);
+}
diff --git a/tests/deal.II/project_q_02/cmp/generic b/tests/deal.II/project_q_02/cmp/generic
new file mode 100644 (file)
index 0000000..fedd731
--- /dev/null
@@ -0,0 +1,148 @@
+
+DEAL::n_dofs=7
+DEAL:cg::Starting value 0.405
+DEAL:cg::Convergence step 7 value 0
+DEAL::FE_Q<1>(1), P_0, rel. error=0
+DEAL:cg::Starting value 0.884
+DEAL:cg::Convergence step 7 value 0
+DEAL::FE_Q<1>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.42
+DEAL:cg::Convergence step 7 value 0
+DEAL::FE_Q<1>(1), P_2, rel. error=0.00129
+DEAL:cg::Starting value 2.00
+DEAL:cg::Convergence step 7 value 0
+DEAL::FE_Q<1>(1), P_3, rel. error=0.00343
+DEAL::n_dofs=13
+DEAL:cg::Starting value 0.319
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Q<1>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.709
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Q<1>(2), P_1, rel. error=0
+DEAL:cg::Starting value 1.17
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Q<1>(2), P_2, rel. error=0
+DEAL:cg::Starting value 1.67
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Q<1>(2), P_3, rel. error=7.69e-05
+DEAL:cg::Starting value 2.22
+DEAL:cg::Convergence step 12 value 0
+DEAL::FE_Q<1>(2), P_4, rel. error=0.000263
+DEAL::n_dofs=19
+DEAL:cg::Starting value 0.251
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(3), P_0, rel. error=0
+DEAL:cg::Starting value 0.559
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(3), P_1, rel. error=0
+DEAL:cg::Starting value 0.923
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(3), P_2, rel. error=0
+DEAL:cg::Starting value 1.33
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(3), P_3, rel. error=0
+DEAL:cg::Starting value 1.78
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(3), P_4, rel. error=2.74e-06
+DEAL:cg::Starting value 2.26
+DEAL:cg::Convergence step 13 value 0
+DEAL::FE_Q<1>(3), P_5, rel. error=1.17e-05
+DEAL::n_dofs=25
+DEAL:cg::Starting value 0.234
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_0, rel. error=0
+DEAL:cg::Starting value 0.522
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_1, rel. error=0
+DEAL:cg::Starting value 0.863
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_2, rel. error=0
+DEAL:cg::Starting value 1.25
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_3, rel. error=0
+DEAL:cg::Starting value 1.67
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_4, rel. error=0
+DEAL:cg::Starting value 2.13
+DEAL:cg::Convergence step 15 value 0
+DEAL::FE_Q<1>(4), P_5, rel. error=1.71e-07
+DEAL:cg::Starting value 2.62
+DEAL:cg::Convergence step 16 value 0
+DEAL::FE_Q<1>(4), P_6, rel. error=8.57e-07
+DEAL::n_dofs=41
+DEAL:cg::Starting value 0.597
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(1), P_0, rel. error=0
+DEAL:cg::Starting value 1.26
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.97
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(1), P_2, rel. error=0.000512
+DEAL:cg::Starting value 2.69
+DEAL:cg::Convergence step 18 value 0
+DEAL::FE_Q<2>(1), P_3, rel. error=0.00127
+DEAL::n_dofs=141
+DEAL:cg::Starting value 0.367
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.785
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_1, rel. error=0
+DEAL:cg::Starting value 1.25
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_2, rel. error=0
+DEAL:cg::Starting value 1.74
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_3, rel. error=2.61e-05
+DEAL:cg::Starting value 2.25
+DEAL:cg::Convergence step 17 value 0
+DEAL::FE_Q<2>(2), P_4, rel. error=7.92e-05
+DEAL::n_dofs=297
+DEAL:cg::Starting value 0.227
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_0, rel. error=0
+DEAL:cg::Starting value 0.488
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_1, rel. error=0
+DEAL:cg::Starting value 0.776
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_2, rel. error=0
+DEAL:cg::Starting value 1.09
+DEAL:cg::Convergence step 19 value 0
+DEAL::FE_Q<2>(3), P_3, rel. error=0
+DEAL:cg::Starting value 1.42
+DEAL:cg::Convergence step 20 value 0
+DEAL::FE_Q<2>(3), P_4, rel. error=7.31e-07
+DEAL:cg::Starting value 1.76
+DEAL:cg::Convergence step 20 value 0
+DEAL::FE_Q<2>(3), P_5, rel. error=2.79e-06
+DEAL::n_dofs=223
+DEAL:cg::Starting value 0.585
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(1), P_0, rel. error=0
+DEAL:cg::Starting value 1.21
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(1), P_1, rel. error=0
+DEAL:cg::Starting value 1.85
+DEAL:cg::Convergence step 24 value 0
+DEAL::FE_Q<3>(1), P_2, rel. error=0.000203
+DEAL:cg::Starting value 2.49
+DEAL:cg::Convergence step 24 value 0
+DEAL::FE_Q<3>(1), P_3, rel. error=0.000493
+DEAL::n_dofs=1375
+DEAL:cg::Starting value 0.282
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(2), P_0, rel. error=0
+DEAL:cg::Starting value 0.588
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(2), P_1, rel. error=0
+DEAL:cg::Starting value 0.913
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(2), P_2, rel. error=0
+DEAL:cg::Starting value 1.25
+DEAL:cg::Convergence step 23 value 0
+DEAL::FE_Q<3>(2), P_3, rel. error=8.13e-06
+DEAL:cg::Starting value 1.61
+DEAL:cg::Convergence step 24 value 0
+DEAL::FE_Q<3>(2), P_4, rel. error=2.39e-05
index 15a4b99a246010b6d716da7877441a53e1d3ff55..d238b1ebc6aa65bee14bd81ab691c053d857f65e 100644 (file)
@@ -100,7 +100,7 @@ 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:cg::Convergence step 115 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
@@ -109,7 +109,7 @@ 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: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
diff --git a/tests/deal.II/project_q_hierarchical_02.cc b/tests/deal.II/project_q_hierarchical_02.cc
new file mode 100644 (file)
index 0000000..0535b07
--- /dev/null
@@ -0,0 +1,29 @@
+//----------------------------  project_q_hierarchical_02.cc  ---------------------------
+//    $Id: project_q_hierarchical_02.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_02.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_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=1; p<6-dim; ++p)
+    test_with_hanging_nodes (FE_Q_Hierarchical<dim>(p), p);
+}
diff --git a/tests/deal.II/project_q_hierarchical_02/cmp/generic b/tests/deal.II/project_q_hierarchical_02/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_hierarchical_02/cmp/x86_64-unknown-linux-gnu+gcc4.1 b/tests/deal.II/project_q_hierarchical_02/cmp/x86_64-unknown-linux-gnu+gcc4.1
new file mode 100644 (file)
index 0000000..0002f1b
--- /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 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_02.cc b/tests/deal.II/project_q_system_02.cc
new file mode 100644 (file)
index 0000000..66481a4
--- /dev/null
@@ -0,0 +1,31 @@
+//----------------------------  project_q_system_02.cc  ---------------------------
+//    $Id: project_q_system_02.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_02.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_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+void test ()
+{
+  for (unsigned int p=1; p<5-dim; ++p)
+    test_with_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_02/cmp/generic b/tests/deal.II/project_q_system_02/cmp/generic
new file mode 100644 (file)
index 0000000..13b2c18
--- /dev/null
@@ -0,0 +1,91 @@
+
+DEAL::n_dofs=25
+DEAL:cg::Starting value 0.734
+DEAL:cg::Convergence step 10 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_0, rel. error=0
+DEAL:cg::Starting value 1.56
+DEAL:cg::Convergence step 10 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_1, rel. error=0
+DEAL:cg::Starting value 2.46
+DEAL:cg::Convergence step 10 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_2, rel. error=0.000418
+DEAL:cg::Starting value 3.42
+DEAL:cg::Convergence step 10 value 0
+DEAL::FESystem<1>[FE_Q<1>(1)-FE_DGQ<1>(2)], P_3, rel. error=0.00121
+DEAL::n_dofs=37
+DEAL:cg::Starting value 0.580
+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 1.24
+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.97
+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.76
+DEAL:cg::Convergence step 13 value 0
+DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_3, rel. error=2.98e-05
+DEAL:cg::Starting value 3.59
+DEAL:cg::Convergence step 13 value 0
+DEAL::FESystem<1>[FE_Q<1>(2)-FE_DGQ<1>(3)], P_4, rel. error=0.000108
+DEAL::n_dofs=49
+DEAL:cg::Starting value 0.524
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_0, rel. error=0
+DEAL:cg::Starting value 1.12
+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.77
+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 2.48
+DEAL:cg::Convergence step 14 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_3, rel. error=0
+DEAL:cg::Starting value 3.23
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_4, rel. error=1.17e-06
+DEAL:cg::Starting value 4.02
+DEAL:cg::Convergence step 15 value 0
+DEAL::FESystem<1>[FE_Q<1>(3)-FE_DGQ<1>(4)], P_5, rel. error=5.21e-06
+DEAL::n_dofs=293
+DEAL:cg::Starting value 0.821
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_0, rel. error=0
+DEAL:cg::Starting value 1.72
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_1, rel. error=0
+DEAL:cg::Starting value 2.68
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_2, rel. error=0.000125
+DEAL:cg::Starting value 3.66
+DEAL:cg::Convergence step 18 value 0
+DEAL::FESystem<2>[FE_Q<2>(1)-FE_DGQ<2>(2)], P_3, rel. error=0.000327
+DEAL::n_dofs=589
+DEAL:cg::Starting value 0.508
+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 1.07
+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 1.69
+DEAL:cg::Convergence step 19 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_2, rel. error=0
+DEAL:cg::Starting value 2.33
+DEAL:cg::Convergence step 19 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_3, rel. error=8.68e-06
+DEAL:cg::Starting value 3.01
+DEAL:cg::Convergence step 19 value 0
+DEAL::FESystem<2>[FE_Q<2>(2)-FE_DGQ<2>(3)], P_4, rel. error=2.75e-05
+DEAL::n_dofs=3463
+DEAL:cg::Starting value 0.695
+DEAL:cg::Convergence step 24 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_0, rel. error=0
+DEAL:cg::Starting value 1.43
+DEAL:cg::Convergence step 24 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_1, rel. error=0
+DEAL:cg::Starting value 2.19
+DEAL:cg::Convergence step 25 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_2, rel. error=3.56e-05
+DEAL:cg::Starting value 2.96
+DEAL:cg::Convergence step 25 value 0
+DEAL::FESystem<3>[FE_Q<3>(1)-FE_DGQ<3>(2)], P_3, rel. error=9.06e-05
diff --git a/tests/deal.II/project_rt_02.cc b/tests/deal.II/project_rt_02.cc
new file mode 100644 (file)
index 0000000..315de1a
--- /dev/null
@@ -0,0 +1,33 @@
+//----------------------------  project_rt_02.cc  ---------------------------
+//    $Id: project_rt_02.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_02.cc  ---------------------------
+
+
+// check that VectorTools::project works for RT elements correctly on
+// a uniformly refined mesh for functions of degree q
+
+char logname[] = "project_rt_02/output";
+
+
+#include "project_common.cc"
+
+
+template <int dim>
+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_hanging_nodes (FE_RaviartThomas<dim>(p), p+1, 1);
+}

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.