From: Guido Kanschat Date: Wed, 9 Mar 2005 03:21:20 +0000 (+0000) Subject: embedding matrices X-Git-Tag: v8.0.0~14474 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=26b47257924e25674c1f5452c0fb651bb381d1a0;p=dealii.git embedding matrices git-svn-id: https://svn.dealii.org/trunk@10054 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/fe/fe_tools.cc b/tests/fe/fe_tools.cc index 639dcba9bc..db7d4071f9 100644 --- a/tests/fe/fe_tools.cc +++ b/tests/fe/fe_tools.cc @@ -2,7 +2,7 @@ // fe_tools.cc,v 1.1 2003/11/28 15:03:26 guido Exp // Version: // -// Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004 by the deal.II authors +// Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -15,6 +15,7 @@ #include "../tests.h" #include +#include #include #include @@ -28,6 +29,49 @@ #include +template +void print_formatted (const FullMatrix &A, + const unsigned int precision, + const unsigned int width) +{ + for (unsigned int i=0; i +void test_embedding (const FiniteElement& fe) +{ + const unsigned int n = fe.dofs_per_cell; + const unsigned int nc= GeometryInfo::children_per_cell; + + FullMatrix P[nc]; + for (unsigned int i=0;i void test_projection (const FiniteElement& fe1, const FiniteElement& fe2, @@ -61,6 +105,14 @@ void test_projection (std::ostream& out) FE_DGP p3(3); FE_DGP p4(4); + test_embedding(q0); + test_embedding(q1); + test_embedding(q2); + test_embedding(q3); + test_embedding(p1); + test_embedding(p2); + test_embedding(p3); + test_projection(p1,p0, out); test_projection(p0,p1, out); test_projection(p2,p1, out);