From 6217205240ad67381169f26b2c10764cabb9a0ec Mon Sep 17 00:00:00 2001 From: wolf Date: Wed, 23 Apr 2003 14:45:56 +0000 Subject: [PATCH] Remove from main branch to avoid problem with building Makefile.depend git-svn-id: https://svn.dealii.org/trunk@7452 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/fe/rt_1.cc | 76 ---------------------- tests/fe/rt_2.cc | 161 ----------------------------------------------- 2 files changed, 237 deletions(-) delete mode 100644 tests/fe/rt_1.cc delete mode 100644 tests/fe/rt_2.cc diff --git a/tests/fe/rt_1.cc b/tests/fe/rt_1.cc deleted file mode 100644 index 1911800c15..0000000000 --- a/tests/fe/rt_1.cc +++ /dev/null @@ -1,76 +0,0 @@ -// $Id$ -// (c) Wolfgang Bangerth -// -// Show the shape functions of the Raviart-Thomas element on the unit cell - -#include -#include - -#include -#include -#include - -#define PRECISION 2 - - - -template -inline void -plot_shape_functions(const unsigned int degree) -{ - FE_RaviartThomas fe_rt(degree); - - const unsigned int div=2; - for (unsigned int mz=0;mz<=((dim>2) ? div : 0) ;++mz) - for (unsigned int my=0;my<=((dim>1) ? div : 0) ;++my) - for (unsigned int mx=0;mx<=div;++mx) - { - const Point p = (dim==2 ? - Point(1.*mx/div, 1.*my/div) : - Point(1.*mx/div, 1.*my/div, 1.*mz/div)); - deallog << "q_point = " << p << std::endl; - - for (unsigned int i=0;i(degree); -// plot_shape_functions<3>(degree); - - return 0; -} - - - diff --git a/tests/fe/rt_2.cc b/tests/fe/rt_2.cc deleted file mode 100644 index 1cb6e336db..0000000000 --- a/tests/fe/rt_2.cc +++ /dev/null @@ -1,161 +0,0 @@ -// $Id$ -// (c) Wolfgang Bangerth -// -// Show the shape functions of the Raviart-Thomas element on a grid -// with only one cell. This cell is rotated, stretched, scaled, etc, -// and on each of these cells each time we evaluate the shape -// functions. - -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include - -#include -#include -#include - -#define PRECISION 2 - - -Point<2> stretch_coordinates (const Point<2> p) -{ - return Point<2>(2*p(0), p(1)); -} - - - -Point<2> tilt_coordinates (const Point<2> p) -{ - return Point<2>(p(0)+p(1), p(1)); -} - - -void -transform_grid (Triangulation<2> &tria, - const unsigned int transform) -{ - switch (transform) - { - // first round: take - // original grid - case 0: - break; - - // second round: rotate - // triangulation - case 1: - GridTools::rotate (3.14159265358/2, tria); - break; - - // third round: inflate - // by a factor of 2 - case 2: - GridTools::scale (2, tria); - break; - - // third round: scale - // back, rotate back, - // stretch - case 3: - GridTools::scale (.5, tria); - GridTools::rotate (-3.14159265358/2, tria); - GridTools::transform (&stretch_coordinates, tria); - - break; - - default: - Assert (false, ExcNotImplemented()); - }; -} - - - - -template -inline void -plot_shape_functions() -{ - FE_RaviartThomas fe_rt(0); - Triangulation tr; - GridGenerator::hyper_cube(tr, 0., 1.); - - // check the following with a - // number of transformed - // triangulations - for (unsigned int transform=0; transform<4; ++transform) - { - deallog << "GRID TRANSFORMATION " << transform << std::endl; - - transform_grid (tr, transform); - - DoFHandler dof(tr); - typename DoFHandler::cell_iterator c = dof.begin(); - dof.distribute_dofs(fe_rt); - - QTrapez<1> q_trapez; - const unsigned int div=2; - QIterated q(q_trapez, div); - FEValues fe(fe_rt, q, update_values|update_gradients|update_q_points); - fe.reinit(c); - - unsigned int q_point=0; - for (unsigned int mz=0;mz<=((dim>2) ? div : 0) ;++mz) - for (unsigned int my=0;my<=((dim>1) ? div : 0) ;++my) - for (unsigned int mx=0;mx<=div;++mx) - { - deallog << "q_point(" << q_point << ")=" << fe.quadrature_point(q_point) - << std::endl; - - for (unsigned int i=0;i(); -// plot_shape_functions<3>(); - - return 0; -} - - - -- 2.39.5