From: kayser-herold Date: Sat, 29 Jul 2006 22:02:39 +0000 (+0000) Subject: Another test that presently produces wrong results X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=869ae5b3a7f9f2b26436f2d7d33cda3991136ec3;p=dealii-svn.git Another test that presently produces wrong results git-svn-id: https://svn.dealii.org/trunk@13510 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/fe/rt_13.cc b/tests/fe/rt_13.cc new file mode 100644 index 0000000000..58d43c8ce5 --- /dev/null +++ b/tests/fe/rt_13.cc @@ -0,0 +1,98 @@ +//---------------------------- rt_13.cc --------------------------- +// rt_13.cc,v 1.3 2003/06/09 16:00:38 wolf Exp +// Version: +// +// Copyright (C) 2003, 2005, 2006 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_13.cc --------------------------- + +// Like rt_11, but for gradients (just as rt_12 is to rt_10) + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include +#include + +#define PRECISION 2 + + +std::ofstream logfile ("rt_13/output"); + +template +void +test (const unsigned int degree) +{ + FE_RaviartThomas fe_rt(degree); + + deallog << "Degree=" << degree + << std::endl; + + for (double h=1; h>1./128; h/=2) + { + deallog << " h=" << h + << std::endl; + + Triangulation tr; + GridGenerator::hyper_cube(tr, 0., h); + + DoFHandler dof(tr); + dof.distribute_dofs(fe_rt); + + QTrapez quadrature; + + FEFaceValues fe_values (fe_rt, quadrature, update_gradients); + fe_values.reinit (dof.begin_active(), 0); + for (unsigned int q=0; q(i); + + return 0; +} + + +