+++ /dev/null
-//---------------------------- 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 <base/logstream.h>
-#include <fe/fe_raviart_thomas.h>
-
-#include <fstream>
-#include <string>
-
-#define PRECISION 5
-
-
-
-template<int dim>
-void
-test(const unsigned int degree)
-{
- deallog << "FE_RaviartThomas<" << dim << "> (" << degree << ")"
- << std::endl;
-
- FE_RaviartThomas<dim> fe_rt(degree);
- const FullMatrix<double> & constraints = fe_rt.constraints();
-
- for (unsigned int i=0; i<constraints.m(); ++i)
- {
- for (unsigned int j=0; j<constraints.n(); ++j)
- deallog << constraints(i,j) << ' ';
- deallog << std::endl;
- }
-
- deallog << std::endl;
-}
-
-
-int
-main()
-{
- std::ofstream logfile ("rt_3.output");
- logfile.precision (PRECISION);
- logfile.setf(std::ios::fixed);
- deallog.attach(logfile);
- deallog.depth_console(0);
-
- for (unsigned int degree=0; degree<4; ++degree)
- test<2>(degree);
-
- return 0;
-}
-
-
-
+++ /dev/null
-//---------------------------- 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 <base/logstream.h>
-#include <fe/fe_raviart_thomas.h>
-
-#include <fstream>
-#include <string>
-
-#define PRECISION 5
-
-
-
-template<int dim>
-void
-test(const unsigned int degree)
-{
- deallog << "FE_RaviartThomas<" << dim << "> (" << degree << ")"
- << std::endl;
-
- FE_RaviartThomas<dim> fe_rt(degree);
-
- for (unsigned int c=0; c<GeometryInfo<dim>::children_per_cell; ++c)
- {
- const FullMatrix<double> & m = fe_rt.prolongate(c);
-
- for (unsigned int i=0; i<m.m(); ++i)
- {
- for (unsigned int j=0; j<m.n(); ++j)
- deallog << m(i,j) << ' ';
- deallog << std::endl;
- }
-
- deallog << std::endl;
- }
-}
-
-
-int
-main()
-{
- std::ofstream logfile ("rt_4.output");
- logfile.precision (PRECISION);
- logfile.setf(std::ios::fixed);
- deallog.attach(logfile);
- deallog.depth_console(0);
-
- for (unsigned int degree=0; degree<4; ++degree)
- test<2>(degree);
-// test<3>(degree);
-
- return 0;
-}
-
-
-