From: bangerth Date: Sun, 19 Sep 2010 17:04:14 +0000 (+0000) Subject: Remove a test that can't have worked for many years (close to 10). X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=201b4681552d6f92437ed7bb250d604b9fcafea2;p=dealii-svn.git Remove a test that can't have worked for many years (close to 10). git-svn-id: https://svn.dealii.org/trunk@22065 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/deal.II/fe_tables.cc b/tests/deal.II/fe_tables.cc deleted file mode 100644 index 810514c973..0000000000 --- a/tests/deal.II/fe_tables.cc +++ /dev/null @@ -1,167 +0,0 @@ -//---------------------------- fe_tables.cc --------------------------- -// $Id$ -// Version: $Name$ -// -// Copyright (C) 1998, 1999, 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. -// -//---------------------------- fe_tables.cc --------------------------- - - -#include "../tests.h" -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include - -#include - -std::ofstream logfile("fe_tables/output"); - -#define TEST_ELEMENT(e) { deallog.push(#e); e el;\ - print_fe_statistics(el); deallog.pop(); deallog << std::endl; } -#define TEST_MULTIPLE(e,n,d) { deallog.push(#e "x" #n); e eb; FESystem el(eb,n); \ - print_fe_statistics(el); deallog.pop(); deallog << std::endl; } -#define TEST_MIXED2(e1,n1,e2,n2,d) { deallog.push(#e1 "x" #n1 "-" #e2 "x" #n2);\ - e1 eb1; e2 eb2; FESystem el(eb1,n1,eb2,n2);\ - print_fe_statistics(el); deallog.pop(); deallog << std::endl; } -#define TEST_MATRIX(e1,e2) { deallog.push( #e1 " onto " #e2); e1 el1; e2 el2;\ - print_fe_matrices(el1,el2); deallog.pop(); deallog << std::endl; } - -template -inline void -print_fe_statistics(const FiniteElement& fe) -{ - Triangulation tr; - GridGenerator::hyper_cube(tr,-1,1); - DoFHandler dof(tr); - dof.distribute_dofs(fe); - StraightBoundary boundary; - typename DoFHandler::active_cell_iterator cell = dof.begin_active(); - typename DoFHandler::active_face_iterator face = dof.begin_active_face(); - - const std::vector > unit_points = fe.get_unit_support_points(); - const std::vector > support_points = fe.get_support_points(); - const std::vector > face_support_points = fe.get_face_support_points(); - - deallog << "dofs_per_cell" << " " << fe.dofs_per_cell; - deallog << ": vertex" << " " << fe.dofs_per_vertex; - deallog << " line" << " " << fe.dofs_per_line; - deallog << " quad" << " " < p = fe.system_to_component_index(i); - deallog << "Index " << i << " (" - << p.first << "," << p.second << ") -> " - << fe.component_to_system_index(p.first, p.second) - << " support " << support_points[i] << " unit: " << unit_points[i] - << std::endl; - } - for (unsigned i=0;i p = fe.face_system_to_component_index(i); - deallog << "FaceIndex " << i << " (" - << p.first << "," << p.second << ") -> " - << fe.face_component_to_system_index(p.first, p.second) - << " support " << face_support_points[i] - << std::endl; - } - deallog.pop(); -} - -template -inline void -print_fe_matrices(const FiniteElement& high, - const FiniteElement& low) -{ - FullMatrix interpolation(low.dofs_per_cell, high.dofs_per_cell); - FETools::get_interpolation_matrix(high, low, interpolation); - deallog << "Interpolation" << std::endl; - interpolation.print(logfile); -} - - -int main() -{ - deallog.attach(logfile); - deallog.depth_console(0); - deallog.threshold_double(1.e-10); - - deallog << std::setprecision(4); - deallog << std::fixed; - - deallog.push("GeometryInfo"); - - deallog.push("1D"); - deallog << " vertices: " << GeometryInfo<1>::vertices_per_cell - << " lines: " << GeometryInfo<1>::lines_per_cell - << " quads: " << GeometryInfo<1>::quads_per_cell - << " hexes: " << GeometryInfo<1>::hexes_per_cell - << std::endl; - deallog.pop(); - - deallog.push("2D"); - deallog << " vertices: " << GeometryInfo<2>::vertices_per_cell - << " lines: " << GeometryInfo<2>::lines_per_cell - << " quads: " << GeometryInfo<2>::quads_per_cell - << " hexes: " << GeometryInfo<2>::hexes_per_cell - << std::endl; - deallog.pop(); - - deallog.push("3D"); - deallog << " vertices: " << GeometryInfo<3>::vertices_per_cell - << " lines: " << GeometryInfo<3>::lines_per_cell - << " quads: " << GeometryInfo<3>::quads_per_cell - << " hexes: " << GeometryInfo<3>::hexes_per_cell - << std::endl; - deallog.pop(); - - deallog.push("4D"); - deallog << " vertices: " << GeometryInfo<4>::vertices_per_cell - << " lines: " << GeometryInfo<4>::lines_per_cell - << " quads: " << GeometryInfo<4>::quads_per_cell - << " hexes: " << GeometryInfo<4>::hexes_per_cell - << std::endl; - deallog.pop(); - - deallog.pop(); - - TEST_ELEMENT(FEDG_Q0<2>); - TEST_ELEMENT(FEDG_Q1<2>); - - TEST_ELEMENT(FEQ1<2>); - TEST_ELEMENT(FEQ2<2>); - TEST_ELEMENT(FEQ3<2>); - TEST_ELEMENT(FEQ4<2>); - - TEST_MULTIPLE(FEQ1<2>,3,2); - TEST_MULTIPLE(FEQ2<2>,3,2); - TEST_MULTIPLE(FEQ3<2>,3,2); - - TEST_MIXED2(FEQ1<2>,1,FEDG_Q0<2>,1,2); - TEST_MIXED2(FEQ2<2>,3,FEQ1<2>,1,2); - TEST_MIXED2(FEQ3<2>,3,FEQ2<2>,2,2); - - deallog.push("Matrices"); - TEST_MATRIX(FEQ2<2>, FEQ1<2>); - TEST_MATRIX(FEQ3<2>, FEQ2<2>); - TEST_MATRIX(FEQ4<2>, FEQ3<2>); - deallog.pop(); -}