+++ /dev/null
-//-----------------------------------------------------------------------------
-// $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);
-}
--- /dev/null
+//----------------------------------------------------------------------
+// $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);
+}
--- /dev/null
+
+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