From fed644818be33a9647689130ff5a73178cc90335 Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Tue, 25 Oct 2005 07:13:48 +0000 Subject: [PATCH] Plot Raviart-Thomas polynomials git-svn-id: https://svn.dealii.org/trunk@11655 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/base/polynomials_rt.cc | 67 ++++++++++++++++++++++++++++++++++++ 1 file changed, 67 insertions(+) create mode 100644 tests/base/polynomials_rt.cc diff --git a/tests/base/polynomials_rt.cc b/tests/base/polynomials_rt.cc new file mode 100644 index 0000000000..cbecc50c38 --- /dev/null +++ b/tests/base/polynomials_rt.cc @@ -0,0 +1,67 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 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 +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- + +// plot PolynomialsRaviartThomas on the reference cell + +#include +#include +#include +#include + +#include +#include + +using namespace std; + +template +void plot(const PolynomialsRaviartThomas& 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 << "RT" << poly.degree() << '<' << dim << '>' + << '\t' << quadrature.point(k); + poly.compute(quadrature.point(k), values, grads, grads2); + + for (unsigned int i=0;i p20(0); + PolynomialsRaviartThomas<2> p21(1); + PolynomialsRaviartThomas<2> p22(2); + + plot(p20); + plot(p21); + plot(p22); +} -- 2.39.5