From: kanschat Date: Mon, 30 Jun 2008 19:51:02 +0000 (+0000) Subject: missing update on test X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7bae6dc4236ca2fd02da9af845cafa8f58e35c0c;p=dealii-svn.git missing update on test git-svn-id: https://svn.dealii.org/trunk@16408 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/deal.II/project_boundary_rt_01.cc b/tests/deal.II/project_boundary_rt_01.cc index 56be52a028..9f4f2c7143 100644 --- a/tests/deal.II/project_boundary_rt_01.cc +++ b/tests/deal.II/project_boundary_rt_01.cc @@ -57,13 +57,14 @@ void TestFunction::vector_value_list (const std::vector > &points, std::vector > &values) const { - if (degree < 2) + for (unsigned int k=0;k void test_projection (const Triangulation& tr, const FiniteElement& fe) { - deallog << fe.get_name() << std::endl; + deallog << fe.get_name() << std::endl << "Cells: " << tr.n_active_cells() << std::endl; const unsigned int degree = fe.tensor_degree(); @@ -138,9 +139,12 @@ void test_projection (const Triangulation& tr, VectorTools::project_boundary_values(mapping, dof, boundary_map, quadrature, boundary_constraints); + deallog << "Constraints: " << boundary_constraints.size() << std::endl; + // Fill a vector with the projected // boundary values Vector u(dof.n_dofs()); + u = -1.; for (typename std::map::const_iterator i = boundary_constraints.begin(); i != boundary_constraints.end(); ++i) u(i->first) = i->second; @@ -181,6 +185,8 @@ int main() deallog.depth_console (0); deallog.threshold_double(1.e-12); + FE_RaviartThomasNodal<2> rt21(1); + test_hyper_cube(rt21); FE_RaviartThomasNodal<2> rt22(2); test_hyper_cube(rt22); FE_RaviartThomasNodal<2> rt23(3);