From 9c3ea22ff77ca6c05d018e657bc3bcc072a82661 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Fri, 2 May 2003 23:12:52 +0000 Subject: [PATCH] . git-svn-id: https://svn.dealii.org/trunk@7557 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/fe/rt_3.cc | 64 ++++++++++++++++++++++++++++++++++++++++++++ tests/fe/rt_4.cc | 69 ++++++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 133 insertions(+) create mode 100644 tests/fe/rt_3.cc create mode 100644 tests/fe/rt_4.cc diff --git a/tests/fe/rt_3.cc b/tests/fe/rt_3.cc new file mode 100644 index 0000000000..98676319ce --- /dev/null +++ b/tests/fe/rt_3.cc @@ -0,0 +1,64 @@ +//---------------------------- rt_3.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2003 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. +// +//---------------------------- rt_3.cc --------------------------- + + +// Just output the constraint matrices of the RT element + +#include +#include + +#include +#include + +#define PRECISION 5 + + + +template +void +test(const unsigned int degree) +{ + deallog << "FE_RaviartThomas<" << dim << "> (" << degree << ")" + << std::endl; + + FE_RaviartThomas fe_rt(degree); + const FullMatrix & constraints = fe_rt.constraints(); + + for (unsigned int i=0; i(degree); + + return 0; +} + + + diff --git a/tests/fe/rt_4.cc b/tests/fe/rt_4.cc new file mode 100644 index 0000000000..99fab02035 --- /dev/null +++ b/tests/fe/rt_4.cc @@ -0,0 +1,69 @@ +//---------------------------- rt_4.cc --------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2003 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. +// +//---------------------------- rt_4.cc --------------------------- + + +// Just output the prolongation matrices of the RT element + +#include +#include + +#include +#include + +#define PRECISION 5 + + + +template +void +test(const unsigned int degree) +{ + deallog << "FE_RaviartThomas<" << dim << "> (" << degree << ")" + << std::endl; + + FE_RaviartThomas fe_rt(degree); + + for (unsigned int c=0; c::children_per_cell; ++c) + { + const FullMatrix & m = fe_rt.prolongate(c); + + for (unsigned int i=0; i(degree); +// test<3>(degree); + + return 0; +} + + + -- 2.39.5