]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
test BDM and RT in 3D
authorkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 20 Sep 2010 16:52:48 +0000 (16:52 +0000)
committerkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 20 Sep 2010 16:52:48 +0000 (16:52 +0000)
git-svn-id: https://svn.dealii.org/trunk@22088 0785d39b-7218-0410-832d-ea1e28bc413d

tests/base/bdm.cc [deleted file]
tests/base/polynomials_bdm_01.cc [new file with mode: 0644]
tests/base/polynomials_bdm_01/cmp/generic [new file with mode: 0644]
tests/base/polynomials_rt.cc

diff --git a/tests/base/bdm.cc b/tests/base/bdm.cc
deleted file mode 100644 (file)
index 35169a8..0000000
+++ /dev/null
@@ -1,88 +0,0 @@
-//-----------------------------------------------------------------------------
-//    $Id$
-//    Version: $Name$ 
-//
-//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2009 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 "../tests.h"
-#include <iomanip>
-#include <fstream>
-#include <cmath>
-
-#include <base/logstream.h>
-#include <base/table.h>
-#include <base/polynomial.h>
-#include <base/polynomials_bdm.h>
-#include <base/quadrature_lib.h>
-
-using namespace Polynomials;
-
-template <int dim>
-void moments (PolynomialsBDM<dim>& poly)
-{
-  std::vector<Polynomial<double> > legendre(2);
-  for (unsigned int i=0;i<legendre.size();++i)
-    legendre[i] = Monomial<double>(i);
-
-  QGauss<1> qface(2);
-
-  Table<2,double> integrals (poly.n(), poly.n());
-
-  std::vector<Tensor<1,dim> > values(poly.n());
-  std::vector<Tensor<2,dim> > grads;
-  std::vector<Tensor<3,dim> > grad_grads;
-  values.resize(poly.n());
-
-  for (unsigned int face=0;face<2*dim;++face)
-    for (unsigned int k=0;k<qface.n_quadrature_points;++k)
-      {
-       const double w = qface.weight(k);
-       const double x = qface.point(k)(0);
-       Point<dim> p;
-       switch (face)
-         {
-           case 2:
-             p(1) = 1.;
-           case 0:
-             p(0) = x;
-             break;
-           case 1:
-             p(0) = 1.;
-           case 3:
-             p(1) = x;
-             break;          
-         }
-       std::cerr << p << std::endl;
-       
-       poly.compute (p, values, grads, grad_grads);
-       for (unsigned int i=0;i<poly.n();++i)
-         {
-           integrals(2*face,i) += w * values[i][1-face%2] ;
-           integrals(2*face+1,i) += w * values[i][1-face%2] * legendre[1].value(x);
-         }
-      }
-
-  std::cout.setf(std::ios::fixed);
-  for (unsigned int i=0;i<integrals.n_rows();++i)
-    {
-      for (unsigned int j=0;j<poly.n();++j)
-       std::cout << std::setw(13) << integrals(i,j);
-      std::cout << std::endl;
-    }
-}
-
-
-
-int main()
-{
-  PolynomialsBDM<2> bdm1(1);
-
-  moments(bdm1);
-}
diff --git a/tests/base/polynomials_bdm_01.cc b/tests/base/polynomials_bdm_01.cc
new file mode 100644 (file)
index 0000000..b8a631a
--- /dev/null
@@ -0,0 +1,70 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2010 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.
+//
+//----------------------------------------------------------------------
+
+// plot PolynomialsBDM on the reference cell
+
+#include "../tests.h"
+#include <base/logstream.h>
+#include <base/job_identifier.h>
+#include <base/tensor.h>
+#include <base/polynomials_bdm.h>
+#include <base/quadrature_lib.h>
+
+#include <vector>
+#include <iomanip>
+#include <fstream>
+
+using namespace std;
+
+template <int dim>
+void plot(const PolynomialsBDM<dim>& poly)
+{
+  QTrapez<1> base_quadrature;
+  QIterated<dim> quadrature(base_quadrature, poly.degree()+3);
+  std::vector<Tensor<1,dim> > values(poly.n());
+  std::vector<Tensor<2,dim> > grads;
+  std::vector<Tensor<3,dim> > grads2;
+
+  
+  for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+    {
+      if (k%(poly.degree()+4) == 0)
+       deallog << "BDM" << poly.degree() << '<' << dim << '>' << std::endl;
+      
+      deallog << "BDM" << poly.degree() << '<' << dim << '>'
+             << '\t' << quadrature.point(k);
+      poly.compute(quadrature.point(k), values, grads, grads2);
+      
+      for (unsigned int i=0;i<poly.n();++i)
+       for (unsigned int d=0;d<dim;++d)
+         deallog << '\t' << values[i][d];
+      deallog << std::endl;
+    }
+}
+
+int main()
+{
+  std::ofstream logfile(JobIdentifier::base_name(__FILE__).c_str());
+  deallog << std::setprecision(3);
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+  
+  PolynomialsBDM<2> p20(0);
+  PolynomialsBDM<2> p21(1);
+  PolynomialsBDM<2> p22(2);
+  
+  plot(p20);
+  plot(p21);
+  plot(p22);
+}
diff --git a/tests/base/polynomials_bdm_01/cmp/generic b/tests/base/polynomials_bdm_01/cmp/generic
new file mode 100644 (file)
index 0000000..172b5dc
--- /dev/null
@@ -0,0 +1,131 @@
+
+DEAL::Name: polynomials_bdm_01.cc 18 polynomials_bdm_01 18446744073709551615 polynomials_bdm_01
+DEAL::File:polynomials_bdm_01
+DEAL::BDM1<2>
+DEAL::BDM1<2>  0.00 0.00       1.00    0       0       1.00    0       0       0       0
+DEAL::BDM1<2>  0.250 0.00      1.00    0       0       1.00    0.250   0       -0.250  0
+DEAL::BDM1<2>  0.500 0.00      1.00    0       0       1.00    0.500   0       -0.500  0
+DEAL::BDM1<2>  0.750 0.00      1.00    0       0       1.00    0.750   0       -0.750  0
+DEAL::BDM1<2>  1.00 0.00       1.00    0       0       1.00    1.00    0       -1.00   0
+DEAL::BDM1<2>
+DEAL::BDM1<2>  0.00 0.250      1.00    0       0       1.00    0       -0.250  0       0.250
+DEAL::BDM1<2>  0.250 0.250     1.00    0       0       1.00    0.250   -0.250  -0.250  0.250
+DEAL::BDM1<2>  0.500 0.250     1.00    0       0       1.00    0.500   -0.250  -0.500  0.250
+DEAL::BDM1<2>  0.750 0.250     1.00    0       0       1.00    0.750   -0.250  -0.750  0.250
+DEAL::BDM1<2>  1.00 0.250      1.00    0       0       1.00    1.00    -0.250  -1.00   0.250
+DEAL::BDM1<2>
+DEAL::BDM1<2>  0.00 0.500      1.00    0       0       1.00    0       -0.500  0       0.500
+DEAL::BDM1<2>  0.250 0.500     1.00    0       0       1.00    0.250   -0.500  -0.250  0.500
+DEAL::BDM1<2>  0.500 0.500     1.00    0       0       1.00    0.500   -0.500  -0.500  0.500
+DEAL::BDM1<2>  0.750 0.500     1.00    0       0       1.00    0.750   -0.500  -0.750  0.500
+DEAL::BDM1<2>  1.00 0.500      1.00    0       0       1.00    1.00    -0.500  -1.00   0.500
+DEAL::BDM1<2>
+DEAL::BDM1<2>  0.00 0.750      1.00    0       0       1.00    0       -0.750  0       0.750
+DEAL::BDM1<2>  0.250 0.750     1.00    0       0       1.00    0.250   -0.750  -0.250  0.750
+DEAL::BDM1<2>  0.500 0.750     1.00    0       0       1.00    0.500   -0.750  -0.500  0.750
+DEAL::BDM1<2>  0.750 0.750     1.00    0       0       1.00    0.750   -0.750  -0.750  0.750
+DEAL::BDM1<2>  1.00 0.750      1.00    0       0       1.00    1.00    -0.750  -1.00   0.750
+DEAL::BDM1<2>
+DEAL::BDM1<2>  0.00 1.00       1.00    0       0       1.00    0       -1.00   0       1.00
+DEAL::BDM1<2>  0.250 1.00      1.00    0       0       1.00    0.250   -1.00   -0.250  1.00
+DEAL::BDM1<2>  0.500 1.00      1.00    0       0       1.00    0.500   -1.00   -0.500  1.00
+DEAL::BDM1<2>  0.750 1.00      1.00    0       0       1.00    0.750   -1.00   -0.750  1.00
+DEAL::BDM1<2>  1.00 1.00       1.00    0       0       1.00    1.00    -1.00   -1.00   1.00
+DEAL::BDM2<2>
+DEAL::BDM2<2>  0.00 0.00       1.00    0       -1.73   0       -1.73   0       0       1.00    0       -1.73   0       -1.73   0       0       0       0
+DEAL::BDM2<2>  0.200 0.00      1.00    0       -1.04   0       -1.73   0       0       1.00    0       -1.04   0       -1.73   0.0400  0       0       0
+DEAL::BDM2<2>  0.400 0.00      1.00    0       -0.346  0       -1.73   0       0       1.00    0       -0.346  0       -1.73   0.160   0       0       0
+DEAL::BDM2<2>  0.600 0.00      1.00    0       0.346   0       -1.73   0       0       1.00    0       0.346   0       -1.73   0.360   0       0       0
+DEAL::BDM2<2>  0.800 0.00      1.00    0       1.04    0       -1.73   0       0       1.00    0       1.04    0       -1.73   0.640   0       0       0
+DEAL::BDM2<2>  1.00 0.00       1.00    0       1.73    0       -1.73   0       0       1.00    0       1.73    0       -1.73   1.00    0       0       0
+DEAL::BDM2<2>
+DEAL::BDM2<2>  0.00 0.200      1.00    0       -1.73   0       -1.04   0       0       1.00    0       -1.73   0       -1.04   0       0       0       0.0400
+DEAL::BDM2<2>  0.200 0.200     1.00    0       -1.04   0       -1.04   0       0       1.00    0       -1.04   0       -1.04   0.0400  -0.0800 -0.0800 0.0400
+DEAL::BDM2<2>  0.400 0.200     1.00    0       -0.346  0       -1.04   0       0       1.00    0       -0.346  0       -1.04   0.160   -0.160  -0.160  0.0400
+DEAL::BDM2<2>  0.600 0.200     1.00    0       0.346   0       -1.04   0       0       1.00    0       0.346   0       -1.04   0.360   -0.240  -0.240  0.0400
+DEAL::BDM2<2>  0.800 0.200     1.00    0       1.04    0       -1.04   0       0       1.00    0       1.04    0       -1.04   0.640   -0.320  -0.320  0.0400
+DEAL::BDM2<2>  1.00 0.200      1.00    0       1.73    0       -1.04   0       0       1.00    0       1.73    0       -1.04   1.00    -0.400  -0.400  0.0400
+DEAL::BDM2<2>
+DEAL::BDM2<2>  0.00 0.400      1.00    0       -1.73   0       -0.346  0       0       1.00    0       -1.73   0       -0.346  0       0       0       0.160
+DEAL::BDM2<2>  0.200 0.400     1.00    0       -1.04   0       -0.346  0       0       1.00    0       -1.04   0       -0.346  0.0400  -0.160  -0.160  0.160
+DEAL::BDM2<2>  0.400 0.400     1.00    0       -0.346  0       -0.346  0       0       1.00    0       -0.346  0       -0.346  0.160   -0.320  -0.320  0.160
+DEAL::BDM2<2>  0.600 0.400     1.00    0       0.346   0       -0.346  0       0       1.00    0       0.346   0       -0.346  0.360   -0.480  -0.480  0.160
+DEAL::BDM2<2>  0.800 0.400     1.00    0       1.04    0       -0.346  0       0       1.00    0       1.04    0       -0.346  0.640   -0.640  -0.640  0.160
+DEAL::BDM2<2>  1.00 0.400      1.00    0       1.73    0       -0.346  0       0       1.00    0       1.73    0       -0.346  1.00    -0.800  -0.800  0.160
+DEAL::BDM2<2>
+DEAL::BDM2<2>  0.00 0.600      1.00    0       -1.73   0       0.346   0       0       1.00    0       -1.73   0       0.346   0       0       0       0.360
+DEAL::BDM2<2>  0.200 0.600     1.00    0       -1.04   0       0.346   0       0       1.00    0       -1.04   0       0.346   0.0400  -0.240  -0.240  0.360
+DEAL::BDM2<2>  0.400 0.600     1.00    0       -0.346  0       0.346   0       0       1.00    0       -0.346  0       0.346   0.160   -0.480  -0.480  0.360
+DEAL::BDM2<2>  0.600 0.600     1.00    0       0.346   0       0.346   0       0       1.00    0       0.346   0       0.346   0.360   -0.720  -0.720  0.360
+DEAL::BDM2<2>  0.800 0.600     1.00    0       1.04    0       0.346   0       0       1.00    0       1.04    0       0.346   0.640   -0.960  -0.960  0.360
+DEAL::BDM2<2>  1.00 0.600      1.00    0       1.73    0       0.346   0       0       1.00    0       1.73    0       0.346   1.00    -1.20   -1.20   0.360
+DEAL::BDM2<2>
+DEAL::BDM2<2>  0.00 0.800      1.00    0       -1.73   0       1.04    0       0       1.00    0       -1.73   0       1.04    0       0       0       0.640
+DEAL::BDM2<2>  0.200 0.800     1.00    0       -1.04   0       1.04    0       0       1.00    0       -1.04   0       1.04    0.0400  -0.320  -0.320  0.640
+DEAL::BDM2<2>  0.400 0.800     1.00    0       -0.346  0       1.04    0       0       1.00    0       -0.346  0       1.04    0.160   -0.640  -0.640  0.640
+DEAL::BDM2<2>  0.600 0.800     1.00    0       0.346   0       1.04    0       0       1.00    0       0.346   0       1.04    0.360   -0.960  -0.960  0.640
+DEAL::BDM2<2>  0.800 0.800     1.00    0       1.04    0       1.04    0       0       1.00    0       1.04    0       1.04    0.640   -1.28   -1.28   0.640
+DEAL::BDM2<2>  1.00 0.800      1.00    0       1.73    0       1.04    0       0       1.00    0       1.73    0       1.04    1.00    -1.60   -1.60   0.640
+DEAL::BDM2<2>
+DEAL::BDM2<2>  0.00 1.00       1.00    0       -1.73   0       1.73    0       0       1.00    0       -1.73   0       1.73    0       0       0       1.00
+DEAL::BDM2<2>  0.200 1.00      1.00    0       -1.04   0       1.73    0       0       1.00    0       -1.04   0       1.73    0.0400  -0.400  -0.400  1.00
+DEAL::BDM2<2>  0.400 1.00      1.00    0       -0.346  0       1.73    0       0       1.00    0       -0.346  0       1.73    0.160   -0.800  -0.800  1.00
+DEAL::BDM2<2>  0.600 1.00      1.00    0       0.346   0       1.73    0       0       1.00    0       0.346   0       1.73    0.360   -1.20   -1.20   1.00
+DEAL::BDM2<2>  0.800 1.00      1.00    0       1.04    0       1.73    0       0       1.00    0       1.04    0       1.73    0.640   -1.60   -1.60   1.00
+DEAL::BDM2<2>  1.00 1.00       1.00    0       1.73    0       1.73    0       0       1.00    0       1.73    0       1.73    1.00    -2.00   -2.00   1.00
+DEAL::BDM3<2>
+DEAL::BDM3<2>  0.00 0.00       1.00    0       -1.73   0       2.24    0       -1.73   0       3.00    0       2.24    0       0       1.00    0       -1.73   0       2.24    0       -1.73   0       3.00    0       2.24    0       0       0       0
+DEAL::BDM3<2>  0.167 0.00      1.00    0       -1.15   0       0.373   0       -1.73   0       2.00    0       2.24    0       0       1.00    0       -1.15   0       0.373   0       -1.73   0       2.00    0       2.24    0.00463 0       0       0
+DEAL::BDM3<2>  0.333 0.00      1.00    0       -0.577  0       -0.745  0       -1.73   0       1.00    0       2.24    0       0       1.00    0       -0.577  0       -0.745  0       -1.73   0       1.00    0       2.24    0.0370  0       0       0
+DEAL::BDM3<2>  0.500 0.00      1.00    0       0       0       -1.12   0       -1.73   0       0       0       2.24    0       0       1.00    0       0       0       -1.12   0       -1.73   0       0       0       2.24    0.125   0       0       0
+DEAL::BDM3<2>  0.667 0.00      1.00    0       0.577   0       -0.745  0       -1.73   0       -1.00   0       2.24    0       0       1.00    0       0.577   0       -0.745  0       -1.73   0       -1.00   0       2.24    0.296   0       0       0
+DEAL::BDM3<2>  0.833 0.00      1.00    0       1.15    0       0.373   0       -1.73   0       -2.00   0       2.24    0       0       1.00    0       1.15    0       0.373   0       -1.73   0       -2.00   0       2.24    0.579   0       0       0
+DEAL::BDM3<2>  1.00 0.00       1.00    0       1.73    0       2.24    0       -1.73   0       -3.00   0       2.24    0       0       1.00    0       1.73    0       2.24    0       -1.73   0       -3.00   0       2.24    1.00    0       0       0
+DEAL::BDM3<2>
+DEAL::BDM3<2>  0.00 0.167      1.00    0       -1.73   0       2.24    0       -1.15   0       2.00    0       0.373   0       0       1.00    0       -1.73   0       2.24    0       -1.15   0       2.00    0       0.373   0       0       0       0.00463
+DEAL::BDM3<2>  0.167 0.167     1.00    0       -1.15   0       0.373   0       -1.15   0       1.33    0       0.373   0       0       1.00    0       -1.15   0       0.373   0       -1.15   0       1.33    0       0.373   0.00463 -0.0139 -0.0139 0.00463
+DEAL::BDM3<2>  0.333 0.167     1.00    0       -0.577  0       -0.745  0       -1.15   0       0.667   0       0.373   0       0       1.00    0       -0.577  0       -0.745  0       -1.15   0       0.667   0       0.373   0.0370  -0.0556 -0.0278 0.00463
+DEAL::BDM3<2>  0.500 0.167     1.00    0       0       0       -1.12   0       -1.15   0       0       0       0.373   0       0       1.00    0       0       0       -1.12   0       -1.15   0       0       0       0.373   0.125   -0.125  -0.0417 0.00463
+DEAL::BDM3<2>  0.667 0.167     1.00    0       0.577   0       -0.745  0       -1.15   0       -0.667  0       0.373   0       0       1.00    0       0.577   0       -0.745  0       -1.15   0       -0.667  0       0.373   0.296   -0.222  -0.0556 0.00463
+DEAL::BDM3<2>  0.833 0.167     1.00    0       1.15    0       0.373   0       -1.15   0       -1.33   0       0.373   0       0       1.00    0       1.15    0       0.373   0       -1.15   0       -1.33   0       0.373   0.579   -0.347  -0.0694 0.00463
+DEAL::BDM3<2>  1.00 0.167      1.00    0       1.73    0       2.24    0       -1.15   0       -2.00   0       0.373   0       0       1.00    0       1.73    0       2.24    0       -1.15   0       -2.00   0       0.373   1.00    -0.500  -0.0833 0.00463
+DEAL::BDM3<2>
+DEAL::BDM3<2>  0.00 0.333      1.00    0       -1.73   0       2.24    0       -0.577  0       1.00    0       -0.745  0       0       1.00    0       -1.73   0       2.24    0       -0.577  0       1.00    0       -0.745  0       0       0       0.0370
+DEAL::BDM3<2>  0.167 0.333     1.00    0       -1.15   0       0.373   0       -0.577  0       0.667   0       -0.745  0       0       1.00    0       -1.15   0       0.373   0       -0.577  0       0.667   0       -0.745  0.00463 -0.0278 -0.0556 0.0370
+DEAL::BDM3<2>  0.333 0.333     1.00    0       -0.577  0       -0.745  0       -0.577  0       0.333   0       -0.745  0       0       1.00    0       -0.577  0       -0.745  0       -0.577  0       0.333   0       -0.745  0.0370  -0.111  -0.111  0.0370
+DEAL::BDM3<2>  0.500 0.333     1.00    0       0       0       -1.12   0       -0.577  0       0       0       -0.745  0       0       1.00    0       0       0       -1.12   0       -0.577  0       0       0       -0.745  0.125   -0.250  -0.167  0.0370
+DEAL::BDM3<2>  0.667 0.333     1.00    0       0.577   0       -0.745  0       -0.577  0       -0.333  0       -0.745  0       0       1.00    0       0.577   0       -0.745  0       -0.577  0       -0.333  0       -0.745  0.296   -0.444  -0.222  0.0370
+DEAL::BDM3<2>  0.833 0.333     1.00    0       1.15    0       0.373   0       -0.577  0       -0.667  0       -0.745  0       0       1.00    0       1.15    0       0.373   0       -0.577  0       -0.667  0       -0.745  0.579   -0.694  -0.278  0.0370
+DEAL::BDM3<2>  1.00 0.333      1.00    0       1.73    0       2.24    0       -0.577  0       -1.00   0       -0.745  0       0       1.00    0       1.73    0       2.24    0       -0.577  0       -1.00   0       -0.745  1.00    -1.00   -0.333  0.0370
+DEAL::BDM3<2>
+DEAL::BDM3<2>  0.00 0.500      1.00    0       -1.73   0       2.24    0       0       0       0       0       -1.12   0       0       1.00    0       -1.73   0       2.24    0       0       0       0       0       -1.12   0       0       0       0.125
+DEAL::BDM3<2>  0.167 0.500     1.00    0       -1.15   0       0.373   0       0       0       0       0       -1.12   0       0       1.00    0       -1.15   0       0.373   0       0       0       0       0       -1.12   0.00463 -0.0417 -0.125  0.125
+DEAL::BDM3<2>  0.333 0.500     1.00    0       -0.577  0       -0.745  0       0       0       0       0       -1.12   0       0       1.00    0       -0.577  0       -0.745  0       0       0       0       0       -1.12   0.0370  -0.167  -0.250  0.125
+DEAL::BDM3<2>  0.500 0.500     1.00    0       0       0       -1.12   0       0       0       0       0       -1.12   0       0       1.00    0       0       0       -1.12   0       0       0       0       0       -1.12   0.125   -0.375  -0.375  0.125
+DEAL::BDM3<2>  0.667 0.500     1.00    0       0.577   0       -0.745  0       0       0       0       0       -1.12   0       0       1.00    0       0.577   0       -0.745  0       0       0       0       0       -1.12   0.296   -0.667  -0.500  0.125
+DEAL::BDM3<2>  0.833 0.500     1.00    0       1.15    0       0.373   0       0       0       0       0       -1.12   0       0       1.00    0       1.15    0       0.373   0       0       0       0       0       -1.12   0.579   -1.04   -0.625  0.125
+DEAL::BDM3<2>  1.00 0.500      1.00    0       1.73    0       2.24    0       0       0       0       0       -1.12   0       0       1.00    0       1.73    0       2.24    0       0       0       0       0       -1.12   1.00    -1.50   -0.750  0.125
+DEAL::BDM3<2>
+DEAL::BDM3<2>  0.00 0.667      1.00    0       -1.73   0       2.24    0       0.577   0       -1.00   0       -0.745  0       0       1.00    0       -1.73   0       2.24    0       0.577   0       -1.00   0       -0.745  0       0       0       0.296
+DEAL::BDM3<2>  0.167 0.667     1.00    0       -1.15   0       0.373   0       0.577   0       -0.667  0       -0.745  0       0       1.00    0       -1.15   0       0.373   0       0.577   0       -0.667  0       -0.745  0.00463 -0.0556 -0.222  0.296
+DEAL::BDM3<2>  0.333 0.667     1.00    0       -0.577  0       -0.745  0       0.577   0       -0.333  0       -0.745  0       0       1.00    0       -0.577  0       -0.745  0       0.577   0       -0.333  0       -0.745  0.0370  -0.222  -0.444  0.296
+DEAL::BDM3<2>  0.500 0.667     1.00    0       0       0       -1.12   0       0.577   0       0       0       -0.745  0       0       1.00    0       0       0       -1.12   0       0.577   0       0       0       -0.745  0.125   -0.500  -0.667  0.296
+DEAL::BDM3<2>  0.667 0.667     1.00    0       0.577   0       -0.745  0       0.577   0       0.333   0       -0.745  0       0       1.00    0       0.577   0       -0.745  0       0.577   0       0.333   0       -0.745  0.296   -0.889  -0.889  0.296
+DEAL::BDM3<2>  0.833 0.667     1.00    0       1.15    0       0.373   0       0.577   0       0.667   0       -0.745  0       0       1.00    0       1.15    0       0.373   0       0.577   0       0.667   0       -0.745  0.579   -1.39   -1.11   0.296
+DEAL::BDM3<2>  1.00 0.667      1.00    0       1.73    0       2.24    0       0.577   0       1.00    0       -0.745  0       0       1.00    0       1.73    0       2.24    0       0.577   0       1.00    0       -0.745  1.00    -2.00   -1.33   0.296
+DEAL::BDM3<2>
+DEAL::BDM3<2>  0.00 0.833      1.00    0       -1.73   0       2.24    0       1.15    0       -2.00   0       0.373   0       0       1.00    0       -1.73   0       2.24    0       1.15    0       -2.00   0       0.373   0       0       0       0.579
+DEAL::BDM3<2>  0.167 0.833     1.00    0       -1.15   0       0.373   0       1.15    0       -1.33   0       0.373   0       0       1.00    0       -1.15   0       0.373   0       1.15    0       -1.33   0       0.373   0.00463 -0.0694 -0.347  0.579
+DEAL::BDM3<2>  0.333 0.833     1.00    0       -0.577  0       -0.745  0       1.15    0       -0.667  0       0.373   0       0       1.00    0       -0.577  0       -0.745  0       1.15    0       -0.667  0       0.373   0.0370  -0.278  -0.694  0.579
+DEAL::BDM3<2>  0.500 0.833     1.00    0       0       0       -1.12   0       1.15    0       0       0       0.373   0       0       1.00    0       0       0       -1.12   0       1.15    0       0       0       0.373   0.125   -0.625  -1.04   0.579
+DEAL::BDM3<2>  0.667 0.833     1.00    0       0.577   0       -0.745  0       1.15    0       0.667   0       0.373   0       0       1.00    0       0.577   0       -0.745  0       1.15    0       0.667   0       0.373   0.296   -1.11   -1.39   0.579
+DEAL::BDM3<2>  0.833 0.833     1.00    0       1.15    0       0.373   0       1.15    0       1.33    0       0.373   0       0       1.00    0       1.15    0       0.373   0       1.15    0       1.33    0       0.373   0.579   -1.74   -1.74   0.579
+DEAL::BDM3<2>  1.00 0.833      1.00    0       1.73    0       2.24    0       1.15    0       2.00    0       0.373   0       0       1.00    0       1.73    0       2.24    0       1.15    0       2.00    0       0.373   1.00    -2.50   -2.08   0.579
+DEAL::BDM3<2>
+DEAL::BDM3<2>  0.00 1.00       1.00    0       -1.73   0       2.24    0       1.73    0       -3.00   0       2.24    0       0       1.00    0       -1.73   0       2.24    0       1.73    0       -3.00   0       2.24    0       0       0       1.00
+DEAL::BDM3<2>  0.167 1.00      1.00    0       -1.15   0       0.373   0       1.73    0       -2.00   0       2.24    0       0       1.00    0       -1.15   0       0.373   0       1.73    0       -2.00   0       2.24    0.00463 -0.0833 -0.500  1.00
+DEAL::BDM3<2>  0.333 1.00      1.00    0       -0.577  0       -0.745  0       1.73    0       -1.00   0       2.24    0       0       1.00    0       -0.577  0       -0.745  0       1.73    0       -1.00   0       2.24    0.0370  -0.333  -1.00   1.00
+DEAL::BDM3<2>  0.500 1.00      1.00    0       0       0       -1.12   0       1.73    0       0       0       2.24    0       0       1.00    0       0       0       -1.12   0       1.73    0       0       0       2.24    0.125   -0.750  -1.50   1.00
+DEAL::BDM3<2>  0.667 1.00      1.00    0       0.577   0       -0.745  0       1.73    0       1.00    0       2.24    0       0       1.00    0       0.577   0       -0.745  0       1.73    0       1.00    0       2.24    0.296   -1.33   -2.00   1.00
+DEAL::BDM3<2>  0.833 1.00      1.00    0       1.15    0       0.373   0       1.73    0       2.00    0       2.24    0       0       1.00    0       1.15    0       0.373   0       1.73    0       2.00    0       2.24    0.579   -2.08   -2.50   1.00
+DEAL::BDM3<2>  1.00 1.00       1.00    0       1.73    0       2.24    0       1.73    0       3.00    0       2.24    0       0       1.00    0       1.73    0       2.24    0       1.73    0       3.00    0       2.24    1.00    -3.00   -3.00   1.00
index f5f344dcc0e67116657809f455caabe725f5996a..b0923d10f48ab0de460d23941379cdb802d8b19a 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -66,4 +66,12 @@ int main()
   plot(p20);
   plot(p21);
   plot(p22);
+
+  PolynomialsRaviartThomas<3> p30(0);
+  PolynomialsRaviartThomas<3> p31(1);
+  PolynomialsRaviartThomas<3> p32(2);
+
+  plot(p30);
+  plot(p31);
+  plot(p32);
 }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.