]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Correction to BDM degree method 1767/head
authorAlistair Bentley <alistairbntl@gmail.com>
Tue, 20 Oct 2015 19:57:50 +0000 (15:57 -0400)
committerAlistair Bentley <alistairbntl@gmail.com>
Sat, 24 Oct 2015 19:35:40 +0000 (15:35 -0400)
* The documentation says this function should return one less than
    the highest polynomial degree. This is consistent with the
    notation used in the literature on BDM polynomials and is what I
    would expect the function to return.
* However, the function actually returns the degree of the highest
    polynomial degree because it calls PolynomialSpace.degree() (note:
    in this context, PolynomialSpace is the complete polynomial space
    of degree k, next then we add elements from k+1 to complete the
    BDM space) which specifically states it does not return the
    Maxmial degree of the polynomial, but the number of polynomials
    given to the constructor.
* If I understand correctly, since the number of polynomials given
    to the constructor includes constants, this returns degree+1.
* To correct this, one needs to be subtracted from the degree of the
    polynomial_space.
* A typo in the 3D bdm_grad functions has also been corrected
* Made a correction in the compute method to reflect the change in the
    degree method.
* Corrected out in the test polynomials_bdm_01.output file
* Edited some names in the test bdm_1.cc

doc/news/changes.h
include/deal.II/base/polynomials_bdm.h
source/base/polynomials_bdm.cc
tests/base/polynomials_bdm_01.cc
tests/base/polynomials_bdm_01.output
tests/fe/bdm_1.cc

index 54eb17eef03ca036efdfbd089cc9fd8542d80b35..7185d9d8b992fd00dc42fb995a288688516aafd3 100644 (file)
@@ -379,6 +379,11 @@ inconvenience this causes.
 
 <ol>
 
+  <li> Fixed: PolynomialsBDM::degree() now returns the correct value.
+  <br>
+  (Alistair Bentley, 2015/10/24)
+  </li>
+
   <li> New: Triangulation::set_all_manifold_ids_on_boundary(boundary_id, manifold_id) 
   which sets the manifold_id for all parts of the boundary with a given boundary_id.
   <br>
index 7dc44c800c04ea8dcaa3a66cb0f156dc98dc4c7c..ed1b51cd07a682440ff434a08ec66b5e2b096fe2 100644 (file)
@@ -216,7 +216,7 @@ template <int dim>
 inline unsigned int
 PolynomialsBDM<dim>::degree() const
 {
-  return polynomial_space.degree();
+  return polynomial_space.degree()-1;
 }
 
 
index a82d87aafb64cda216bf9a0a9fa9e3c229e6c3f2..3653da62f1e30cd97245c0955b0bc77e9732f1bf 100644 (file)
@@ -188,7 +188,7 @@ PolynomialsBDM<dim>::compute (const Point<dim>            &unit_point,
               // p(t) = t^(i+1)
               monomials[i+1].value(unit_point(d), monovali[d]);
               // q(t) = t^(k-i)
-              monomials[degree()-i-1].value(unit_point(d), monovalk[d]);
+              monomials[degree()-i].value(unit_point(d), monovalk[d]);
             }
           if (values.size() != 0)
             {
@@ -218,9 +218,9 @@ PolynomialsBDM<dim>::compute (const Point<dim>            &unit_point,
               grads[start][1][0] = 0.;
               grads[start][1][1] = -monovali[1][1] * monovalk[2][0];
               grads[start][1][2] = -monovali[1][0] * monovalk[2][1];
-              grads[start+2][2][0] = 0.;
-              grads[start+2][2][1] = 0.;
-              grads[start+2][2][2] = 0.;
+              grads[start][2][0] = 0.;
+              grads[start][2][1] = 0.;
+              grads[start][2][2] = 0.;
 
               grads[start+1][1][1] = monovali[2][1] * monovalk[0][0];
               grads[start+1][1][2] = unit_point(1) * monovali[2][2] * monovalk[0][0];
index 10ba6c4bef86b0a1314047e1008fec00b78ac0ef..eeb3590ffadf05cab6843c7fd482339761db5927 100644 (file)
@@ -33,17 +33,16 @@ template <int dim>
 void plot(const PolynomialsBDM<dim> &poly)
 {
   QTrapez<1> base_quadrature;
-  QIterated<dim> quadrature(base_quadrature, poly.degree()+3);
+  QIterated<dim> quadrature(base_quadrature, poly.degree()+4);
   std::vector<Tensor<1,dim> > values(poly.n());
   std::vector<Tensor<2,dim> > grads;
   std::vector<Tensor<3,dim> > grads2;
   std::vector<Tensor<4,dim> > thirds;
   std::vector<Tensor<5,dim> > fourths;
 
-
   for (unsigned int k=0; k<quadrature.size(); ++k)
     {
-      if (k%(poly.degree()+4) == 0)
+      if (k%(poly.degree()+5) == 0)
         deallog << "BDM" << poly.degree() << '<' << dim << '>' << std::endl;
 
       deallog << "BDM" << poly.degree() << '<' << dim << '>'
index d7a51a8eacd974e6d49f8ddf56d0b7f57d60a599..bb02e48278d76654441feb3c588632db8b798349 100644 (file)
 
+DEAL::BDM0<2>
+DEAL::BDM0<2>  0.00 0.00       1.00    0       0       1.00    0       0       0       0
+DEAL::BDM0<2>  0.250 0.00      1.00    0       0       1.00    0.250   0       0.250   0
+DEAL::BDM0<2>  0.500 0.00      1.00    0       0       1.00    0.500   0       0.500   0
+DEAL::BDM0<2>  0.750 0.00      1.00    0       0       1.00    0.750   0       0.750   0
+DEAL::BDM0<2>  1.00 0.00       1.00    0       0       1.00    1.00    0       1.00    0
+DEAL::BDM0<2>
+DEAL::BDM0<2>  0.00 0.250      1.00    0       0       1.00    0       -0.250  0       -0.250
+DEAL::BDM0<2>  0.250 0.250     1.00    0       0       1.00    0.250   -0.250  0.250   -0.250
+DEAL::BDM0<2>  0.500 0.250     1.00    0       0       1.00    0.500   -0.250  0.500   -0.250
+DEAL::BDM0<2>  0.750 0.250     1.00    0       0       1.00    0.750   -0.250  0.750   -0.250
+DEAL::BDM0<2>  1.00 0.250      1.00    0       0       1.00    1.00    -0.250  1.00    -0.250
+DEAL::BDM0<2>
+DEAL::BDM0<2>  0.00 0.500      1.00    0       0       1.00    0       -0.500  0       -0.500
+DEAL::BDM0<2>  0.250 0.500     1.00    0       0       1.00    0.250   -0.500  0.250   -0.500
+DEAL::BDM0<2>  0.500 0.500     1.00    0       0       1.00    0.500   -0.500  0.500   -0.500
+DEAL::BDM0<2>  0.750 0.500     1.00    0       0       1.00    0.750   -0.500  0.750   -0.500
+DEAL::BDM0<2>  1.00 0.500      1.00    0       0       1.00    1.00    -0.500  1.00    -0.500
+DEAL::BDM0<2>
+DEAL::BDM0<2>  0.00 0.750      1.00    0       0       1.00    0       -0.750  0       -0.750
+DEAL::BDM0<2>  0.250 0.750     1.00    0       0       1.00    0.250   -0.750  0.250   -0.750
+DEAL::BDM0<2>  0.500 0.750     1.00    0       0       1.00    0.500   -0.750  0.500   -0.750
+DEAL::BDM0<2>  0.750 0.750     1.00    0       0       1.00    0.750   -0.750  0.750   -0.750
+DEAL::BDM0<2>  1.00 0.750      1.00    0       0       1.00    1.00    -0.750  1.00    -0.750
+DEAL::BDM0<2>
+DEAL::BDM0<2>  0.00 1.00       1.00    0       0       1.00    0       -1.00   0       -1.00
+DEAL::BDM0<2>  0.250 1.00      1.00    0       0       1.00    0.250   -1.00   0.250   -1.00
+DEAL::BDM0<2>  0.500 1.00      1.00    0       0       1.00    0.500   -1.00   0.500   -1.00
+DEAL::BDM0<2>  0.750 1.00      1.00    0       0       1.00    0.750   -1.00   0.750   -1.00
+DEAL::BDM0<2>  1.00 1.00       1.00    0       0       1.00    1.00    -1.00   1.00    -1.00
 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>  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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::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>  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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::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>  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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::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>  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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<2>
+DEAL::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM1<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::BDM2<2>
+DEAL::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::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>  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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::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>  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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::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>  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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::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>  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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::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>  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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::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
+DEAL::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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::BDM2<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 e37cda7f4e8eaa8a73869e9dfe589cf675f71823..848c5f511a7f50fec51de556c058ac3e8c4496c0 100644 (file)
@@ -35,8 +35,8 @@ template<int dim>
 inline void
 plot_shape_functions(const unsigned int degree)
 {
-  FE_BDM<dim> fe_rt(degree);
-  deallog.push(fe_rt.get_name());
+  FE_BDM<dim> fe_bdm(degree);
+  deallog.push(fe_bdm.get_name());
 
   const unsigned int div=2;
   for (unsigned int mz=0; mz<=((dim>2) ? div : 0) ; ++mz)
@@ -52,22 +52,22 @@ plot_shape_functions(const unsigned int degree)
             // values contain
             // quadrature point and one
             // vector of dim entries
-            // for each chape function
+            // for each shape function
             deallog << "value " << p;
-            for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
+            for (unsigned int i=0; i<fe_bdm.dofs_per_cell; ++i)
               {
                 for (unsigned int c=0; c<dim; ++c)
-                  deallog << " " << fe_rt.shape_value_component(i,p,c);
+                  deallog << " " << fe_bdm.shape_value_component(i,p,c);
                 deallog << "  ";
               }
             deallog << std::endl << "grad "  << p;
-            for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
+            for (unsigned int i=0; i<fe_bdm.dofs_per_cell; ++i)
               {
                 for (unsigned int c=0; c<dim; ++c)
                   {
                     deallog << ' ';
                     for (unsigned int d=0; d<dim; ++d)
-                      deallog << ' ' << fe_rt.shape_grad_component(i,p,c)[d];
+                      deallog << ' ' << fe_bdm.shape_grad_component(i,p,c)[d];
                   }
               }
             deallog << std::endl;

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.