From: kanschat Date: Mon, 20 Sep 2010 16:52:48 +0000 (+0000) Subject: test BDM and RT in 3D X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=e73834e3760e9769356074765ad27f6f2a1cb938;p=dealii-svn.git test BDM and RT in 3D git-svn-id: https://svn.dealii.org/trunk@22088 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/base/bdm.cc b/tests/base/bdm.cc deleted file mode 100644 index 35169a88d8..0000000000 --- a/tests/base/bdm.cc +++ /dev/null @@ -1,88 +0,0 @@ -//----------------------------------------------------------------------------- -// $Id$ -// Version: $Name$ -// -// Copyright (C) 2000, 2001, 2002, 2003, 2004, 2009 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. -// -//----------------------------------------------------------------------------- - -#include "../tests.h" -#include -#include -#include - -#include -#include -#include -#include -#include - -using namespace Polynomials; - -template -void moments (PolynomialsBDM& poly) -{ - std::vector > legendre(2); - for (unsigned int i=0;i(i); - - QGauss<1> qface(2); - - Table<2,double> integrals (poly.n(), poly.n()); - - std::vector > values(poly.n()); - std::vector > grads; - std::vector > grad_grads; - values.resize(poly.n()); - - for (unsigned int face=0;face<2*dim;++face) - for (unsigned int k=0;k p; - switch (face) - { - case 2: - p(1) = 1.; - case 0: - p(0) = x; - break; - case 1: - p(0) = 1.; - case 3: - p(1) = x; - break; - } - std::cerr << p << std::endl; - - poly.compute (p, values, grads, grad_grads); - for (unsigned int i=0;i bdm1(1); - - moments(bdm1); -} diff --git a/tests/base/polynomials_bdm_01.cc b/tests/base/polynomials_bdm_01.cc new file mode 100644 index 0000000000..b8a631a9eb --- /dev/null +++ b/tests/base/polynomials_bdm_01.cc @@ -0,0 +1,70 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2010 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. +// +//---------------------------------------------------------------------- + +// plot PolynomialsBDM on the reference cell + +#include "../tests.h" +#include +#include +#include +#include +#include + +#include +#include +#include + +using namespace std; + +template +void plot(const PolynomialsBDM& poly) +{ + QTrapez<1> base_quadrature; + QIterated quadrature(base_quadrature, poly.degree()+3); + std::vector > values(poly.n()); + std::vector > grads; + std::vector > grads2; + + + for (unsigned int k=0;k' << std::endl; + + deallog << "BDM" << poly.degree() << '<' << dim << '>' + << '\t' << quadrature.point(k); + poly.compute(quadrature.point(k), values, grads, grads2); + + for (unsigned int i=0;i p20(0); + PolynomialsBDM<2> p21(1); + PolynomialsBDM<2> p22(2); + + plot(p20); + plot(p21); + plot(p22); +} diff --git a/tests/base/polynomials_bdm_01/cmp/generic b/tests/base/polynomials_bdm_01/cmp/generic new file mode 100644 index 0000000000..172b5dc977 --- /dev/null +++ b/tests/base/polynomials_bdm_01/cmp/generic @@ -0,0 +1,131 @@ + +DEAL::Name: polynomials_bdm_01.cc 18 polynomials_bdm_01 18446744073709551615 polynomials_bdm_01 +DEAL::File:polynomials_bdm_01 +DEAL::BDM1<2> +DEAL::BDM1<2> 0.00 0.00 1.00 0 0 1.00 0 0 0 0 +DEAL::BDM1<2> 0.250 0.00 1.00 0 0 1.00 0.250 0 -0.250 0 +DEAL::BDM1<2> 0.500 0.00 1.00 0 0 1.00 0.500 0 -0.500 0 +DEAL::BDM1<2> 0.750 0.00 1.00 0 0 1.00 0.750 0 -0.750 0 +DEAL::BDM1<2> 1.00 0.00 1.00 0 0 1.00 1.00 0 -1.00 0 +DEAL::BDM1<2> +DEAL::BDM1<2> 0.00 0.250 1.00 0 0 1.00 0 -0.250 0 0.250 +DEAL::BDM1<2> 0.250 0.250 1.00 0 0 1.00 0.250 -0.250 -0.250 0.250 +DEAL::BDM1<2> 0.500 0.250 1.00 0 0 1.00 0.500 -0.250 -0.500 0.250 +DEAL::BDM1<2> 0.750 0.250 1.00 0 0 1.00 0.750 -0.250 -0.750 0.250 +DEAL::BDM1<2> 1.00 0.250 1.00 0 0 1.00 1.00 -0.250 -1.00 0.250 +DEAL::BDM1<2> +DEAL::BDM1<2> 0.00 0.500 1.00 0 0 1.00 0 -0.500 0 0.500 +DEAL::BDM1<2> 0.250 0.500 1.00 0 0 1.00 0.250 -0.500 -0.250 0.500 +DEAL::BDM1<2> 0.500 0.500 1.00 0 0 1.00 0.500 -0.500 -0.500 0.500 +DEAL::BDM1<2> 0.750 0.500 1.00 0 0 1.00 0.750 -0.500 -0.750 0.500 +DEAL::BDM1<2> 1.00 0.500 1.00 0 0 1.00 1.00 -0.500 -1.00 0.500 +DEAL::BDM1<2> +DEAL::BDM1<2> 0.00 0.750 1.00 0 0 1.00 0 -0.750 0 0.750 +DEAL::BDM1<2> 0.250 0.750 1.00 0 0 1.00 0.250 -0.750 -0.250 0.750 +DEAL::BDM1<2> 0.500 0.750 1.00 0 0 1.00 0.500 -0.750 -0.500 0.750 +DEAL::BDM1<2> 0.750 0.750 1.00 0 0 1.00 0.750 -0.750 -0.750 0.750 +DEAL::BDM1<2> 1.00 0.750 1.00 0 0 1.00 1.00 -0.750 -1.00 0.750 +DEAL::BDM1<2> +DEAL::BDM1<2> 0.00 1.00 1.00 0 0 1.00 0 -1.00 0 1.00 +DEAL::BDM1<2> 0.250 1.00 1.00 0 0 1.00 0.250 -1.00 -0.250 1.00 +DEAL::BDM1<2> 0.500 1.00 1.00 0 0 1.00 0.500 -1.00 -0.500 1.00 +DEAL::BDM1<2> 0.750 1.00 1.00 0 0 1.00 0.750 -1.00 -0.750 1.00 +DEAL::BDM1<2> 1.00 1.00 1.00 0 0 1.00 1.00 -1.00 -1.00 1.00 +DEAL::BDM2<2> +DEAL::BDM2<2> 0.00 0.00 1.00 0 -1.73 0 -1.73 0 0 1.00 0 -1.73 0 -1.73 0 0 0 0 +DEAL::BDM2<2> 0.200 0.00 1.00 0 -1.04 0 -1.73 0 0 1.00 0 -1.04 0 -1.73 0.0400 0 0 0 +DEAL::BDM2<2> 0.400 0.00 1.00 0 -0.346 0 -1.73 0 0 1.00 0 -0.346 0 -1.73 0.160 0 0 0 +DEAL::BDM2<2> 0.600 0.00 1.00 0 0.346 0 -1.73 0 0 1.00 0 0.346 0 -1.73 0.360 0 0 0 +DEAL::BDM2<2> 0.800 0.00 1.00 0 1.04 0 -1.73 0 0 1.00 0 1.04 0 -1.73 0.640 0 0 0 +DEAL::BDM2<2> 1.00 0.00 1.00 0 1.73 0 -1.73 0 0 1.00 0 1.73 0 -1.73 1.00 0 0 0 +DEAL::BDM2<2> +DEAL::BDM2<2> 0.00 0.200 1.00 0 -1.73 0 -1.04 0 0 1.00 0 -1.73 0 -1.04 0 0 0 0.0400 +DEAL::BDM2<2> 0.200 0.200 1.00 0 -1.04 0 -1.04 0 0 1.00 0 -1.04 0 -1.04 0.0400 -0.0800 -0.0800 0.0400 +DEAL::BDM2<2> 0.400 0.200 1.00 0 -0.346 0 -1.04 0 0 1.00 0 -0.346 0 -1.04 0.160 -0.160 -0.160 0.0400 +DEAL::BDM2<2> 0.600 0.200 1.00 0 0.346 0 -1.04 0 0 1.00 0 0.346 0 -1.04 0.360 -0.240 -0.240 0.0400 +DEAL::BDM2<2> 0.800 0.200 1.00 0 1.04 0 -1.04 0 0 1.00 0 1.04 0 -1.04 0.640 -0.320 -0.320 0.0400 +DEAL::BDM2<2> 1.00 0.200 1.00 0 1.73 0 -1.04 0 0 1.00 0 1.73 0 -1.04 1.00 -0.400 -0.400 0.0400 +DEAL::BDM2<2> +DEAL::BDM2<2> 0.00 0.400 1.00 0 -1.73 0 -0.346 0 0 1.00 0 -1.73 0 -0.346 0 0 0 0.160 +DEAL::BDM2<2> 0.200 0.400 1.00 0 -1.04 0 -0.346 0 0 1.00 0 -1.04 0 -0.346 0.0400 -0.160 -0.160 0.160 +DEAL::BDM2<2> 0.400 0.400 1.00 0 -0.346 0 -0.346 0 0 1.00 0 -0.346 0 -0.346 0.160 -0.320 -0.320 0.160 +DEAL::BDM2<2> 0.600 0.400 1.00 0 0.346 0 -0.346 0 0 1.00 0 0.346 0 -0.346 0.360 -0.480 -0.480 0.160 +DEAL::BDM2<2> 0.800 0.400 1.00 0 1.04 0 -0.346 0 0 1.00 0 1.04 0 -0.346 0.640 -0.640 -0.640 0.160 +DEAL::BDM2<2> 1.00 0.400 1.00 0 1.73 0 -0.346 0 0 1.00 0 1.73 0 -0.346 1.00 -0.800 -0.800 0.160 +DEAL::BDM2<2> +DEAL::BDM2<2> 0.00 0.600 1.00 0 -1.73 0 0.346 0 0 1.00 0 -1.73 0 0.346 0 0 0 0.360 +DEAL::BDM2<2> 0.200 0.600 1.00 0 -1.04 0 0.346 0 0 1.00 0 -1.04 0 0.346 0.0400 -0.240 -0.240 0.360 +DEAL::BDM2<2> 0.400 0.600 1.00 0 -0.346 0 0.346 0 0 1.00 0 -0.346 0 0.346 0.160 -0.480 -0.480 0.360 +DEAL::BDM2<2> 0.600 0.600 1.00 0 0.346 0 0.346 0 0 1.00 0 0.346 0 0.346 0.360 -0.720 -0.720 0.360 +DEAL::BDM2<2> 0.800 0.600 1.00 0 1.04 0 0.346 0 0 1.00 0 1.04 0 0.346 0.640 -0.960 -0.960 0.360 +DEAL::BDM2<2> 1.00 0.600 1.00 0 1.73 0 0.346 0 0 1.00 0 1.73 0 0.346 1.00 -1.20 -1.20 0.360 +DEAL::BDM2<2> +DEAL::BDM2<2> 0.00 0.800 1.00 0 -1.73 0 1.04 0 0 1.00 0 -1.73 0 1.04 0 0 0 0.640 +DEAL::BDM2<2> 0.200 0.800 1.00 0 -1.04 0 1.04 0 0 1.00 0 -1.04 0 1.04 0.0400 -0.320 -0.320 0.640 +DEAL::BDM2<2> 0.400 0.800 1.00 0 -0.346 0 1.04 0 0 1.00 0 -0.346 0 1.04 0.160 -0.640 -0.640 0.640 +DEAL::BDM2<2> 0.600 0.800 1.00 0 0.346 0 1.04 0 0 1.00 0 0.346 0 1.04 0.360 -0.960 -0.960 0.640 +DEAL::BDM2<2> 0.800 0.800 1.00 0 1.04 0 1.04 0 0 1.00 0 1.04 0 1.04 0.640 -1.28 -1.28 0.640 +DEAL::BDM2<2> 1.00 0.800 1.00 0 1.73 0 1.04 0 0 1.00 0 1.73 0 1.04 1.00 -1.60 -1.60 0.640 +DEAL::BDM2<2> +DEAL::BDM2<2> 0.00 1.00 1.00 0 -1.73 0 1.73 0 0 1.00 0 -1.73 0 1.73 0 0 0 1.00 +DEAL::BDM2<2> 0.200 1.00 1.00 0 -1.04 0 1.73 0 0 1.00 0 -1.04 0 1.73 0.0400 -0.400 -0.400 1.00 +DEAL::BDM2<2> 0.400 1.00 1.00 0 -0.346 0 1.73 0 0 1.00 0 -0.346 0 1.73 0.160 -0.800 -0.800 1.00 +DEAL::BDM2<2> 0.600 1.00 1.00 0 0.346 0 1.73 0 0 1.00 0 0.346 0 1.73 0.360 -1.20 -1.20 1.00 +DEAL::BDM2<2> 0.800 1.00 1.00 0 1.04 0 1.73 0 0 1.00 0 1.04 0 1.73 0.640 -1.60 -1.60 1.00 +DEAL::BDM2<2> 1.00 1.00 1.00 0 1.73 0 1.73 0 0 1.00 0 1.73 0 1.73 1.00 -2.00 -2.00 1.00 +DEAL::BDM3<2> +DEAL::BDM3<2> 0.00 0.00 1.00 0 -1.73 0 2.24 0 -1.73 0 3.00 0 2.24 0 0 1.00 0 -1.73 0 2.24 0 -1.73 0 3.00 0 2.24 0 0 0 0 +DEAL::BDM3<2> 0.167 0.00 1.00 0 -1.15 0 0.373 0 -1.73 0 2.00 0 2.24 0 0 1.00 0 -1.15 0 0.373 0 -1.73 0 2.00 0 2.24 0.00463 0 0 0 +DEAL::BDM3<2> 0.333 0.00 1.00 0 -0.577 0 -0.745 0 -1.73 0 1.00 0 2.24 0 0 1.00 0 -0.577 0 -0.745 0 -1.73 0 1.00 0 2.24 0.0370 0 0 0 +DEAL::BDM3<2> 0.500 0.00 1.00 0 0 0 -1.12 0 -1.73 0 0 0 2.24 0 0 1.00 0 0 0 -1.12 0 -1.73 0 0 0 2.24 0.125 0 0 0 +DEAL::BDM3<2> 0.667 0.00 1.00 0 0.577 0 -0.745 0 -1.73 0 -1.00 0 2.24 0 0 1.00 0 0.577 0 -0.745 0 -1.73 0 -1.00 0 2.24 0.296 0 0 0 +DEAL::BDM3<2> 0.833 0.00 1.00 0 1.15 0 0.373 0 -1.73 0 -2.00 0 2.24 0 0 1.00 0 1.15 0 0.373 0 -1.73 0 -2.00 0 2.24 0.579 0 0 0 +DEAL::BDM3<2> 1.00 0.00 1.00 0 1.73 0 2.24 0 -1.73 0 -3.00 0 2.24 0 0 1.00 0 1.73 0 2.24 0 -1.73 0 -3.00 0 2.24 1.00 0 0 0 +DEAL::BDM3<2> +DEAL::BDM3<2> 0.00 0.167 1.00 0 -1.73 0 2.24 0 -1.15 0 2.00 0 0.373 0 0 1.00 0 -1.73 0 2.24 0 -1.15 0 2.00 0 0.373 0 0 0 0.00463 +DEAL::BDM3<2> 0.167 0.167 1.00 0 -1.15 0 0.373 0 -1.15 0 1.33 0 0.373 0 0 1.00 0 -1.15 0 0.373 0 -1.15 0 1.33 0 0.373 0.00463 -0.0139 -0.0139 0.00463 +DEAL::BDM3<2> 0.333 0.167 1.00 0 -0.577 0 -0.745 0 -1.15 0 0.667 0 0.373 0 0 1.00 0 -0.577 0 -0.745 0 -1.15 0 0.667 0 0.373 0.0370 -0.0556 -0.0278 0.00463 +DEAL::BDM3<2> 0.500 0.167 1.00 0 0 0 -1.12 0 -1.15 0 0 0 0.373 0 0 1.00 0 0 0 -1.12 0 -1.15 0 0 0 0.373 0.125 -0.125 -0.0417 0.00463 +DEAL::BDM3<2> 0.667 0.167 1.00 0 0.577 0 -0.745 0 -1.15 0 -0.667 0 0.373 0 0 1.00 0 0.577 0 -0.745 0 -1.15 0 -0.667 0 0.373 0.296 -0.222 -0.0556 0.00463 +DEAL::BDM3<2> 0.833 0.167 1.00 0 1.15 0 0.373 0 -1.15 0 -1.33 0 0.373 0 0 1.00 0 1.15 0 0.373 0 -1.15 0 -1.33 0 0.373 0.579 -0.347 -0.0694 0.00463 +DEAL::BDM3<2> 1.00 0.167 1.00 0 1.73 0 2.24 0 -1.15 0 -2.00 0 0.373 0 0 1.00 0 1.73 0 2.24 0 -1.15 0 -2.00 0 0.373 1.00 -0.500 -0.0833 0.00463 +DEAL::BDM3<2> +DEAL::BDM3<2> 0.00 0.333 1.00 0 -1.73 0 2.24 0 -0.577 0 1.00 0 -0.745 0 0 1.00 0 -1.73 0 2.24 0 -0.577 0 1.00 0 -0.745 0 0 0 0.0370 +DEAL::BDM3<2> 0.167 0.333 1.00 0 -1.15 0 0.373 0 -0.577 0 0.667 0 -0.745 0 0 1.00 0 -1.15 0 0.373 0 -0.577 0 0.667 0 -0.745 0.00463 -0.0278 -0.0556 0.0370 +DEAL::BDM3<2> 0.333 0.333 1.00 0 -0.577 0 -0.745 0 -0.577 0 0.333 0 -0.745 0 0 1.00 0 -0.577 0 -0.745 0 -0.577 0 0.333 0 -0.745 0.0370 -0.111 -0.111 0.0370 +DEAL::BDM3<2> 0.500 0.333 1.00 0 0 0 -1.12 0 -0.577 0 0 0 -0.745 0 0 1.00 0 0 0 -1.12 0 -0.577 0 0 0 -0.745 0.125 -0.250 -0.167 0.0370 +DEAL::BDM3<2> 0.667 0.333 1.00 0 0.577 0 -0.745 0 -0.577 0 -0.333 0 -0.745 0 0 1.00 0 0.577 0 -0.745 0 -0.577 0 -0.333 0 -0.745 0.296 -0.444 -0.222 0.0370 +DEAL::BDM3<2> 0.833 0.333 1.00 0 1.15 0 0.373 0 -0.577 0 -0.667 0 -0.745 0 0 1.00 0 1.15 0 0.373 0 -0.577 0 -0.667 0 -0.745 0.579 -0.694 -0.278 0.0370 +DEAL::BDM3<2> 1.00 0.333 1.00 0 1.73 0 2.24 0 -0.577 0 -1.00 0 -0.745 0 0 1.00 0 1.73 0 2.24 0 -0.577 0 -1.00 0 -0.745 1.00 -1.00 -0.333 0.0370 +DEAL::BDM3<2> +DEAL::BDM3<2> 0.00 0.500 1.00 0 -1.73 0 2.24 0 0 0 0 0 -1.12 0 0 1.00 0 -1.73 0 2.24 0 0 0 0 0 -1.12 0 0 0 0.125 +DEAL::BDM3<2> 0.167 0.500 1.00 0 -1.15 0 0.373 0 0 0 0 0 -1.12 0 0 1.00 0 -1.15 0 0.373 0 0 0 0 0 -1.12 0.00463 -0.0417 -0.125 0.125 +DEAL::BDM3<2> 0.333 0.500 1.00 0 -0.577 0 -0.745 0 0 0 0 0 -1.12 0 0 1.00 0 -0.577 0 -0.745 0 0 0 0 0 -1.12 0.0370 -0.167 -0.250 0.125 +DEAL::BDM3<2> 0.500 0.500 1.00 0 0 0 -1.12 0 0 0 0 0 -1.12 0 0 1.00 0 0 0 -1.12 0 0 0 0 0 -1.12 0.125 -0.375 -0.375 0.125 +DEAL::BDM3<2> 0.667 0.500 1.00 0 0.577 0 -0.745 0 0 0 0 0 -1.12 0 0 1.00 0 0.577 0 -0.745 0 0 0 0 0 -1.12 0.296 -0.667 -0.500 0.125 +DEAL::BDM3<2> 0.833 0.500 1.00 0 1.15 0 0.373 0 0 0 0 0 -1.12 0 0 1.00 0 1.15 0 0.373 0 0 0 0 0 -1.12 0.579 -1.04 -0.625 0.125 +DEAL::BDM3<2> 1.00 0.500 1.00 0 1.73 0 2.24 0 0 0 0 0 -1.12 0 0 1.00 0 1.73 0 2.24 0 0 0 0 0 -1.12 1.00 -1.50 -0.750 0.125 +DEAL::BDM3<2> +DEAL::BDM3<2> 0.00 0.667 1.00 0 -1.73 0 2.24 0 0.577 0 -1.00 0 -0.745 0 0 1.00 0 -1.73 0 2.24 0 0.577 0 -1.00 0 -0.745 0 0 0 0.296 +DEAL::BDM3<2> 0.167 0.667 1.00 0 -1.15 0 0.373 0 0.577 0 -0.667 0 -0.745 0 0 1.00 0 -1.15 0 0.373 0 0.577 0 -0.667 0 -0.745 0.00463 -0.0556 -0.222 0.296 +DEAL::BDM3<2> 0.333 0.667 1.00 0 -0.577 0 -0.745 0 0.577 0 -0.333 0 -0.745 0 0 1.00 0 -0.577 0 -0.745 0 0.577 0 -0.333 0 -0.745 0.0370 -0.222 -0.444 0.296 +DEAL::BDM3<2> 0.500 0.667 1.00 0 0 0 -1.12 0 0.577 0 0 0 -0.745 0 0 1.00 0 0 0 -1.12 0 0.577 0 0 0 -0.745 0.125 -0.500 -0.667 0.296 +DEAL::BDM3<2> 0.667 0.667 1.00 0 0.577 0 -0.745 0 0.577 0 0.333 0 -0.745 0 0 1.00 0 0.577 0 -0.745 0 0.577 0 0.333 0 -0.745 0.296 -0.889 -0.889 0.296 +DEAL::BDM3<2> 0.833 0.667 1.00 0 1.15 0 0.373 0 0.577 0 0.667 0 -0.745 0 0 1.00 0 1.15 0 0.373 0 0.577 0 0.667 0 -0.745 0.579 -1.39 -1.11 0.296 +DEAL::BDM3<2> 1.00 0.667 1.00 0 1.73 0 2.24 0 0.577 0 1.00 0 -0.745 0 0 1.00 0 1.73 0 2.24 0 0.577 0 1.00 0 -0.745 1.00 -2.00 -1.33 0.296 +DEAL::BDM3<2> +DEAL::BDM3<2> 0.00 0.833 1.00 0 -1.73 0 2.24 0 1.15 0 -2.00 0 0.373 0 0 1.00 0 -1.73 0 2.24 0 1.15 0 -2.00 0 0.373 0 0 0 0.579 +DEAL::BDM3<2> 0.167 0.833 1.00 0 -1.15 0 0.373 0 1.15 0 -1.33 0 0.373 0 0 1.00 0 -1.15 0 0.373 0 1.15 0 -1.33 0 0.373 0.00463 -0.0694 -0.347 0.579 +DEAL::BDM3<2> 0.333 0.833 1.00 0 -0.577 0 -0.745 0 1.15 0 -0.667 0 0.373 0 0 1.00 0 -0.577 0 -0.745 0 1.15 0 -0.667 0 0.373 0.0370 -0.278 -0.694 0.579 +DEAL::BDM3<2> 0.500 0.833 1.00 0 0 0 -1.12 0 1.15 0 0 0 0.373 0 0 1.00 0 0 0 -1.12 0 1.15 0 0 0 0.373 0.125 -0.625 -1.04 0.579 +DEAL::BDM3<2> 0.667 0.833 1.00 0 0.577 0 -0.745 0 1.15 0 0.667 0 0.373 0 0 1.00 0 0.577 0 -0.745 0 1.15 0 0.667 0 0.373 0.296 -1.11 -1.39 0.579 +DEAL::BDM3<2> 0.833 0.833 1.00 0 1.15 0 0.373 0 1.15 0 1.33 0 0.373 0 0 1.00 0 1.15 0 0.373 0 1.15 0 1.33 0 0.373 0.579 -1.74 -1.74 0.579 +DEAL::BDM3<2> 1.00 0.833 1.00 0 1.73 0 2.24 0 1.15 0 2.00 0 0.373 0 0 1.00 0 1.73 0 2.24 0 1.15 0 2.00 0 0.373 1.00 -2.50 -2.08 0.579 +DEAL::BDM3<2> +DEAL::BDM3<2> 0.00 1.00 1.00 0 -1.73 0 2.24 0 1.73 0 -3.00 0 2.24 0 0 1.00 0 -1.73 0 2.24 0 1.73 0 -3.00 0 2.24 0 0 0 1.00 +DEAL::BDM3<2> 0.167 1.00 1.00 0 -1.15 0 0.373 0 1.73 0 -2.00 0 2.24 0 0 1.00 0 -1.15 0 0.373 0 1.73 0 -2.00 0 2.24 0.00463 -0.0833 -0.500 1.00 +DEAL::BDM3<2> 0.333 1.00 1.00 0 -0.577 0 -0.745 0 1.73 0 -1.00 0 2.24 0 0 1.00 0 -0.577 0 -0.745 0 1.73 0 -1.00 0 2.24 0.0370 -0.333 -1.00 1.00 +DEAL::BDM3<2> 0.500 1.00 1.00 0 0 0 -1.12 0 1.73 0 0 0 2.24 0 0 1.00 0 0 0 -1.12 0 1.73 0 0 0 2.24 0.125 -0.750 -1.50 1.00 +DEAL::BDM3<2> 0.667 1.00 1.00 0 0.577 0 -0.745 0 1.73 0 1.00 0 2.24 0 0 1.00 0 0.577 0 -0.745 0 1.73 0 1.00 0 2.24 0.296 -1.33 -2.00 1.00 +DEAL::BDM3<2> 0.833 1.00 1.00 0 1.15 0 0.373 0 1.73 0 2.00 0 2.24 0 0 1.00 0 1.15 0 0.373 0 1.73 0 2.00 0 2.24 0.579 -2.08 -2.50 1.00 +DEAL::BDM3<2> 1.00 1.00 1.00 0 1.73 0 2.24 0 1.73 0 3.00 0 2.24 0 0 1.00 0 1.73 0 2.24 0 1.73 0 3.00 0 2.24 1.00 -3.00 -3.00 1.00 diff --git a/tests/base/polynomials_rt.cc b/tests/base/polynomials_rt.cc index f5f344dcc0..b0923d10f4 100644 --- a/tests/base/polynomials_rt.cc +++ b/tests/base/polynomials_rt.cc @@ -2,7 +2,7 @@ // $Id$ // Version: $Name$ // -// Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008 by the deal.II authors +// Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008, 2010 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -66,4 +66,12 @@ int main() plot(p20); plot(p21); plot(p22); + + PolynomialsRaviartThomas<3> p30(0); + PolynomialsRaviartThomas<3> p31(1); + PolynomialsRaviartThomas<3> p32(2); + + plot(p30); + plot(p31); + plot(p32); }