From 5186593e9eefd4464bdfaeebfed57949cea77bb0 Mon Sep 17 00:00:00 2001 From: Graham Harper Date: Tue, 3 Jul 2018 21:23:11 +0000 Subject: [PATCH] Add tests for the BR polynomials. --- tests/base/polynomials_br.cc | 181 +++++++++++++++++++ tests/base/polynomials_br.output | 299 +++++++++++++++++++++++++++++++ 2 files changed, 480 insertions(+) create mode 100644 tests/base/polynomials_br.cc create mode 100644 tests/base/polynomials_br.output diff --git a/tests/base/polynomials_br.cc b/tests/base/polynomials_br.cc new file mode 100644 index 0000000000..1350b77592 --- /dev/null +++ b/tests/base/polynomials_br.cc @@ -0,0 +1,181 @@ +// --------------------------------------------------------------------- +// +// HEADER SPACE FOR LATER USE +// +// +// +// +// +// +// +// +// +// +// --------------------------------------------------------------------- + + +// evaluate polynomials_bernardi_raugel on the reference cell +// watch for the first dim*vertices_per_cell shape functions to +// yield delta functions when evaluated on the nodes +// watch for the remaining faces_per_cell shape functions to +// have magnitude 1 on the face centers + +#include +#include +#include +#include + +#include + +#include "../tests.h" + +template +void +check_poly_q(const PolynomialsBernardiRaugel &poly) +{ + std::vector> points; + std::vector> values; + std::vector> grads; + std::vector> grads2; + std::vector> thirds; + std::vector> fourths; + + // set up reference cell vertices + if (dim == 2) + { + for (unsigned int i = 0; i < 2; ++i) + { + for (unsigned int j = 0; j < 2; ++j) + { + Point p; + p[0] = j; + p[1] = i; + points.push_back(p); + } + } + } + else if (dim == 3) + { + for (unsigned int i = 0; i < 2; ++i) + { + for (unsigned int j = 0; j < 2; ++j) + { + for (unsigned int k = 0; k < 2; ++k) + { + Point p; + p[0] = k; + p[1] = j; + p[2] = i; + points.push_back(p); + } + } + } + } + + // loop over vertices + for (unsigned int i = 0; i < points.size(); ++i) + { + values.clear(); + grads.clear(); + grads2.clear(); + thirds.clear(); + fourths.clear(); + values.resize(dim * GeometryInfo::vertices_per_cell + + GeometryInfo::faces_per_cell); + + deallog << "BR1<" << dim << "> point " << i << " (" << points[i][0]; + for (unsigned int d = 1; d < dim; ++d) + deallog << ", " << points[i][d]; + deallog << ")" << std::endl; + + poly.compute(points[i], values, grads, grads2, thirds, fourths); + + // loop through the Q_1^d shape functions + for (unsigned int j = 0; j < dim * GeometryInfo::vertices_per_cell; + ++j) + { + deallog << "BR1<" << dim << "> shape fxn " << j << ": "; + for (unsigned int d = 0; d < dim; ++d) + deallog << '\t' << values[j][d]; + deallog << std::endl; + } + } +} + + + +template +void +check_poly_bubble(const PolynomialsBernardiRaugel &poly) +{ + std::vector> points; + std::vector> values; + std::vector> grads; + std::vector> grads2; + std::vector> thirds; + std::vector> fourths; + + // set up reference cell face midpoints + for (unsigned int i = 0; i < dim; ++i) + { + for (unsigned int j = 0; j < 2; ++j) + { + Point p; + p[0] = 0.5; + p[1] = 0.5; + if (dim == 3) + p[2] = 0.5; + p[i] = j; + points.push_back(p); + } + } + + // loop over vertices + for (unsigned int i = 0; i < points.size(); ++i) + { + values.clear(); + grads.clear(); + grads2.clear(); + thirds.clear(); + fourths.clear(); + values.resize(dim * GeometryInfo::vertices_per_cell + + GeometryInfo::faces_per_cell); + + deallog << "BR1<" << dim << "> point " + << (i + GeometryInfo::vertices_per_cell) << " (" + << points[i][0]; + for (unsigned int d = 1; d < dim; ++d) + deallog << ", " << points[i][d]; + deallog << ")" << std::endl; + + poly.compute(points[i], values, grads, grads2, thirds, fourths); + + // loop through the bubble shape functions + for (unsigned int j = dim * GeometryInfo::vertices_per_cell; + j < dim * GeometryInfo::vertices_per_cell + + GeometryInfo::faces_per_cell; + ++j) + { + deallog << "BR1<" << dim << "> shape fxn " << j << ": "; + for (unsigned int d = 0; d < dim; ++d) + deallog << '\t' << values[j][d]; + deallog << std::endl; + } + } +} + +int +main() +{ + std::ofstream logfile("output"); + deallog << std::setprecision(3); + deallog.attach(logfile); + + PolynomialsBernardiRaugel<2> p_2d(1); + check_poly_q(p_2d); + check_poly_bubble(p_2d); + + PolynomialsBernardiRaugel<3> p_3d(1); + check_poly_q(p_3d); + check_poly_bubble(p_3d); +} diff --git a/tests/base/polynomials_br.output b/tests/base/polynomials_br.output new file mode 100644 index 0000000000..04cd7adc84 --- /dev/null +++ b/tests/base/polynomials_br.output @@ -0,0 +1,299 @@ + +DEAL::BR1<2> point 0 (0.00, 0.00) +DEAL::BR1<2> shape fxn 0: 1.00 0.00 +DEAL::BR1<2> shape fxn 1: 0.00 1.00 +DEAL::BR1<2> shape fxn 2: 0.00 0.00 +DEAL::BR1<2> shape fxn 3: 0.00 0.00 +DEAL::BR1<2> shape fxn 4: 0.00 0.00 +DEAL::BR1<2> shape fxn 5: 0.00 0.00 +DEAL::BR1<2> shape fxn 6: 0.00 0.00 +DEAL::BR1<2> shape fxn 7: 0.00 0.00 +DEAL::BR1<2> point 1 (1.00, 0.00) +DEAL::BR1<2> shape fxn 0: 0.00 0.00 +DEAL::BR1<2> shape fxn 1: 0.00 0.00 +DEAL::BR1<2> shape fxn 2: 1.00 0.00 +DEAL::BR1<2> shape fxn 3: 0.00 1.00 +DEAL::BR1<2> shape fxn 4: 0.00 0.00 +DEAL::BR1<2> shape fxn 5: 0.00 0.00 +DEAL::BR1<2> shape fxn 6: 0.00 0.00 +DEAL::BR1<2> shape fxn 7: 0.00 0.00 +DEAL::BR1<2> point 2 (0.00, 1.00) +DEAL::BR1<2> shape fxn 0: 0.00 0.00 +DEAL::BR1<2> shape fxn 1: 0.00 0.00 +DEAL::BR1<2> shape fxn 2: 0.00 0.00 +DEAL::BR1<2> shape fxn 3: 0.00 0.00 +DEAL::BR1<2> shape fxn 4: 1.00 0.00 +DEAL::BR1<2> shape fxn 5: 0.00 1.00 +DEAL::BR1<2> shape fxn 6: 0.00 0.00 +DEAL::BR1<2> shape fxn 7: 0.00 0.00 +DEAL::BR1<2> point 3 (1.00, 1.00) +DEAL::BR1<2> shape fxn 0: 0.00 0.00 +DEAL::BR1<2> shape fxn 1: 0.00 0.00 +DEAL::BR1<2> shape fxn 2: 0.00 0.00 +DEAL::BR1<2> shape fxn 3: 0.00 0.00 +DEAL::BR1<2> shape fxn 4: 0.00 0.00 +DEAL::BR1<2> shape fxn 5: 0.00 0.00 +DEAL::BR1<2> shape fxn 6: 1.00 0.00 +DEAL::BR1<2> shape fxn 7: 0.00 1.00 +DEAL::BR1<2> point 4 (0.00, 0.500) +DEAL::BR1<2> shape fxn 8: 1.00 0.00 +DEAL::BR1<2> shape fxn 9: 0.00 0.00 +DEAL::BR1<2> shape fxn 10: 0.00 0.00 +DEAL::BR1<2> shape fxn 11: 0.00 0.00 +DEAL::BR1<2> point 5 (1.00, 0.500) +DEAL::BR1<2> shape fxn 8: 0.00 0.00 +DEAL::BR1<2> shape fxn 9: 1.00 0.00 +DEAL::BR1<2> shape fxn 10: 0.00 0.00 +DEAL::BR1<2> shape fxn 11: 0.00 0.00 +DEAL::BR1<2> point 6 (0.500, 0.00) +DEAL::BR1<2> shape fxn 8: 0.00 0.00 +DEAL::BR1<2> shape fxn 9: 0.00 0.00 +DEAL::BR1<2> shape fxn 10: 0.00 1.00 +DEAL::BR1<2> shape fxn 11: 0.00 0.00 +DEAL::BR1<2> point 7 (0.500, 1.00) +DEAL::BR1<2> shape fxn 8: 0.00 0.00 +DEAL::BR1<2> shape fxn 9: 0.00 0.00 +DEAL::BR1<2> shape fxn 10: 0.00 0.00 +DEAL::BR1<2> shape fxn 11: 0.00 1.00 +DEAL::BR1<3> point 0 (0.00, 0.00, 0.00) +DEAL::BR1<3> shape fxn 0: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 3: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 6: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 9: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 12: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 15: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 18: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 21: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 0.00 +DEAL::BR1<3> point 1 (1.00, 0.00, 0.00) +DEAL::BR1<3> shape fxn 0: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 3: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 6: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 9: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 12: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 15: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 18: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 21: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 0.00 +DEAL::BR1<3> point 2 (0.00, 1.00, 0.00) +DEAL::BR1<3> shape fxn 0: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 3: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 6: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 9: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 12: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 15: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 18: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 21: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 0.00 +DEAL::BR1<3> point 3 (1.00, 1.00, 0.00) +DEAL::BR1<3> shape fxn 0: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 3: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 6: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 9: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 12: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 15: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 18: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 21: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 0.00 +DEAL::BR1<3> point 4 (0.00, 0.00, 1.00) +DEAL::BR1<3> shape fxn 0: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 3: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 6: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 9: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 12: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 15: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 18: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 21: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 0.00 +DEAL::BR1<3> point 5 (1.00, 0.00, 1.00) +DEAL::BR1<3> shape fxn 0: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 3: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 6: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 9: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 12: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 15: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 18: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 21: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 0.00 +DEAL::BR1<3> point 6 (0.00, 1.00, 1.00) +DEAL::BR1<3> shape fxn 0: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 3: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 6: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 9: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 12: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 15: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 18: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 21: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 0.00 +DEAL::BR1<3> point 7 (1.00, 1.00, 1.00) +DEAL::BR1<3> shape fxn 0: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 1: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 2: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 3: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 4: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 5: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 6: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 7: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 8: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 9: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 10: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 11: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 12: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 13: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 14: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 15: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 16: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 17: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 18: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 19: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 20: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 21: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 22: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 23: 0.00 0.00 1.00 +DEAL::BR1<3> point 8 (0.00, 0.500, 0.500) +DEAL::BR1<3> shape fxn 24: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 25: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 26: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 27: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 28: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 29: 0.00 0.00 0.00 +DEAL::BR1<3> point 9 (1.00, 0.500, 0.500) +DEAL::BR1<3> shape fxn 24: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 25: 1.00 0.00 0.00 +DEAL::BR1<3> shape fxn 26: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 27: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 28: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 29: 0.00 0.00 0.00 +DEAL::BR1<3> point 10 (0.500, 0.00, 0.500) +DEAL::BR1<3> shape fxn 24: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 25: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 26: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 27: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 28: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 29: 0.00 0.00 0.00 +DEAL::BR1<3> point 11 (0.500, 1.00, 0.500) +DEAL::BR1<3> shape fxn 24: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 25: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 26: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 27: 0.00 1.00 0.00 +DEAL::BR1<3> shape fxn 28: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 29: 0.00 0.00 0.00 +DEAL::BR1<3> point 12 (0.500, 0.500, 0.00) +DEAL::BR1<3> shape fxn 24: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 25: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 26: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 27: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 28: 0.00 0.00 1.00 +DEAL::BR1<3> shape fxn 29: 0.00 0.00 0.00 +DEAL::BR1<3> point 13 (0.500, 0.500, 1.00) +DEAL::BR1<3> shape fxn 24: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 25: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 26: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 27: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 28: 0.00 0.00 0.00 +DEAL::BR1<3> shape fxn 29: 0.00 0.00 1.00 -- 2.39.5