From: Guido Kanschat Date: Wed, 19 Oct 2005 04:56:45 +0000 (+0000) Subject: tests for RTN and RT X-Git-Tag: v8.0.0~12981 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=36be9af56355dc7f404cf267d6294360574ed86c;p=dealii.git tests for RTN and RT git-svn-id: https://svn.dealii.org/trunk@11622 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/fe/Makefile b/tests/fe/Makefile index 1ebe5c3f2f..5c6c3b1d28 100644 --- a/tests/fe/Makefile +++ b/tests/fe/Makefile @@ -30,7 +30,7 @@ tests = fe_data_test traits fe_tools fe_tools_test mapping \ transfer internals \ dgq_1 \ q_1 q_2 q_3 q_4 interpolate_q1 \ - rt_1 rt_2 rt_3 rt_4 rt_5 rt_6 rtdiff rtn_1 interpolate_rtn \ + rt_1 rt_2 rt_3 rt_4 rt_5 rt_6 rtdiff rtn_1 interpolate_rt interpolate_rtn \ dgp_monomial_1 dgp_monomial_2 \ system_1 interpolate_system \ non_primitive_1 non_primitive_2 \ diff --git a/tests/fe/interpolate_rt.cc b/tests/fe/interpolate_rt.cc new file mode 100644 index 0000000000..4d61217071 --- /dev/null +++ b/tests/fe/interpolate_rt.cc @@ -0,0 +1,73 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 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. +// +//---------------------------------------------------------------------- + +#include "interpolate_common.h" +#include +#include + +#include + +#include + +// FE_RaviartThomas::interpolate(...) + +template +void check1(const Function& f, + const unsigned int degree) +{ + FE_RaviartThomas fe(degree); + deallog << fe.get_name() << ' '; + deallog << fe.get_generalized_support_points().size() << ' '; + + std::vector dofs(fe.dofs_per_cell); + + std::vector > + values(dim, std::vector(fe.get_generalized_support_points().size())); + std::vector > + vectors(fe.get_generalized_support_points().size(), + Vector(dim)); + f.vector_value_list(fe.get_generalized_support_points(), vectors); + + for (unsigned int c=0;c w1; +// check1(w1,1,2); +// check1(w1,2,2); +// check1(w1,3,2); + Q1WedgeFunction<2,1,2> w21; + check1(w21,1); + check1(w21,2); + Q1WedgeFunction<2,2,2> w22; + check1(w22,2); + Q1WedgeFunction<2,3,2> w23; + check1(w23,3); +// Q1WedgeFunction<3,1,3> w3; +// check1(w3,1); +// check1(w3,2); +} diff --git a/tests/fe/interpolate_rtn.cc b/tests/fe/interpolate_rtn.cc index f22e6e5265..f18bf86caa 100644 --- a/tests/fe/interpolate_rtn.cc +++ b/tests/fe/interpolate_rtn.cc @@ -54,18 +54,20 @@ int main() std::ofstream logfile ("interpolate_rtn.output"); deallog.attach(logfile); deallog.depth_console(0); - deallog.threshold_double(1.e-15); + deallog.threshold_double(1.e-13); -// Q1WedgeFunction<1,1,2> w1; -// check1(w1,1,2); -// check1(w1,2,2); -// check1(w1,3,2); +// Q1WedgeFunction<1,1,2> w1; +// check1(w1,1,2); +// check1(w1,2,2); +// check1(w1,3,2); Q1WedgeFunction<2,1,2> w2; check1(w2,1); check1(w2,2); Q1WedgeFunction<2,2,2> w22; check1(w22,2); -// Q1WedgeFunction<3,1,3> w3; -// check1(w3,1); -// check1(w3,2); + Q1WedgeFunction<2,3,2> w23; + check1(w23,3); + Q1WedgeFunction<3,1,3> w3; + check1(w3,1); + check1(w3,2); }