From: Wolfgang Bangerth Date: Wed, 9 Feb 2011 14:15:50 +0000 (+0000) Subject: Fix a problem with refining surfaces. X-Git-Tag: v8.0.0~4390 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=3f074f3655db31cd7f2ac121c4fa56e447ae2876;p=dealii.git Fix a problem with refining surfaces. git-svn-id: https://svn.dealii.org/trunk@23308 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/grid/tria.cc b/deal.II/source/grid/tria.cc index c780ee03f4..775cbf2382 100644 --- a/deal.II/source/grid/tria.cc +++ b/deal.II/source/grid/tria.cc @@ -4660,11 +4660,31 @@ namespace internal if (line->has_children() == false) { line->set_user_flag (); -//TODO[WB]: we overwrite the user_index here because we later on need to find -// out which boundary object we have to ask to refine this line. find a better -// way to do that +//TODO[WB]: we overwrite the user_index here because we later on need +// to find out which boundary object we have to ask to refine this +// line. we can't use the boundary_indicator field because that can +// only be used for lines at the boundary of the domain, but we also +// need a domain description for interior lines in the codim-1 case if (spacedim > dim) - line->set_user_index(cell->material_id()); + { + if (line->at_boundary()) + // if + // possible + // honor + // boundary + // indicator + line->set_user_index(line->boundary_indicator()); + else + // otherwise + // take + // manifold + // description + // from + // the + // adjacent + // cell + line->set_user_index(cell->material_id()); + } } } } diff --git a/tests/codim_one/boundary_indicator_01.cc b/tests/codim_one/boundary_indicator_01.cc new file mode 100644 index 0000000000..207970fbeb --- /dev/null +++ b/tests/codim_one/boundary_indicator_01.cc @@ -0,0 +1,116 @@ +//---------------------------- boundary_indicator_01.cc --------------------------- +// $Id: bem.cc 22693 2010-11-11 20:11:47Z kanschat $ +// Version: $Name$ +// +// Copyright (C) 2010, 2011 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. +// +//---------------------------- boundary_indicator_01.cc --------------------------- + + +// for surfaces, we need some sort of mapping also for interior cells +// and faces, but at the time of writing this, we can only set +// boundary_indicators of faces and edges that are truly at the +// boundary of the domain. to work around this, we use the material_id +// field that one can set on cells, and copy it also to the adjacent +// faces. initially, however, we also copied the material_id to the +// boundary_indicator of adjacent faces that truly were at the +// boundary of the domain, and which might have had something +// purposefully set already +// +// this test verifies that this is now fixed. + +#include "../tests.h" + +#include +#include +#include +#include +#include + +using namespace std; + + + +template +void save_mesh(const Triangulation& tria) +{ + GridOut grid_out; + grid_out.write_gnuplot (tria, deallog.get_file_stream()); +} + + +int main () +{ + ofstream logfile("boundary_indicator_01/output"); + deallog.attach(logfile); + deallog.depth_console(0); + + // Extract the boundary of 3/4 of a sphere + { + const int dim = 3; + deallog << "Testing hyper_cube in dim: " << dim << "..."<< endl; + + const HyperBallBoundary boundary_description; + Triangulation volume_mesh; + GridGenerator::hyper_ball(volume_mesh); + volume_mesh.set_boundary (0, boundary_description); + + // exclude one of the 6 faces + // from the surface mesh + // extraction + for (typename Triangulation::active_cell_iterator + cell = volume_mesh.begin_active(); + cell != volume_mesh.end(); ++cell) + { + bool done = false; + for (unsigned int f=0; f::faces_per_cell; ++f) + if (cell->at_boundary(f)) + { + cell->face(f)->set_boundary_indicator(1); + done = true; + break; + } + if (done) + break; + } + + const HyperBallBoundary surface_description; + Triangulation boundary_mesh; + boundary_mesh.set_boundary (0, surface_description); + + // now extract a mesh of the 5 + // surface faces + std::set boundary_indicators; + boundary_indicators.insert (0); + GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh, + boundary_indicators); + deallog << volume_mesh.n_active_cells () << std::endl; + deallog << boundary_mesh.n_active_cells () << std::endl; + + // at this point, all cells and + // edges of the surface mesh + // should have boundary indicator + // 0. set those at the boundary + // of the mesh to 1 to force + // straight line refinement, then + // refine + for (typename Triangulation::active_cell_iterator + cell = boundary_mesh.begin_active(); + cell != boundary_mesh.end(); ++cell) + for (unsigned int f=0; f::faces_per_cell; ++f) + if (cell->at_boundary(f)) + cell->face(f)->set_boundary_indicator(1); + + boundary_mesh.refine_global (2); + + save_mesh(boundary_mesh); + } + + + return 0; +} diff --git a/tests/codim_one/boundary_indicator_01/cmp/generic b/tests/codim_one/boundary_indicator_01/cmp/generic new file mode 100644 index 0000000000..29384a7a22 --- /dev/null +++ b/tests/codim_one/boundary_indicator_01/cmp/generic @@ -0,0 +1,564 @@ + +DEAL::Testing hyper_cube in dim: 3... +DEAL::7 +DEAL::5 +0.577350 -0.577350 -0.577350 2 0 +0.673887 -0.673887 -0.302905 2 0 +0.862418 -0.379225 -0.335298 2 0 +0.577350 -0.288675 -0.577350 2 0 +0.577350 -0.577350 -0.577350 2 0 + + +0.673887 -0.673887 -0.302905 2 0 +0.707107 -0.707107 0.00000 2 0 +0.923762 -0.382710 0.0140184 2 0 +0.862418 -0.379225 -0.335298 2 0 +0.673887 -0.673887 -0.302905 2 0 + + +0.577350 -0.288675 -0.577350 2 0 +0.862418 -0.379225 -0.335298 2 0 +0.943953 0.00000 -0.330081 2 0 +0.577350 0.00000 -0.577350 2 0 +0.577350 -0.288675 -0.577350 2 0 + + +0.862418 -0.379225 -0.335298 2 0 +0.923762 -0.382710 0.0140184 2 0 +0.999665 0.00000 0.0259008 2 0 +0.943953 0.00000 -0.330081 2 0 +0.862418 -0.379225 -0.335298 2 0 + + +0.707107 -0.707107 0.00000 2 0 +0.673887 -0.673887 0.302905 2 0 +0.853386 -0.365005 0.372160 2 0 +0.923762 -0.382710 0.0140184 2 0 +0.707107 -0.707107 0.00000 2 0 + + +0.673887 -0.673887 0.302905 2 0 +0.577350 -0.577350 0.577350 2 0 +0.673887 -0.302905 0.673887 2 0 +0.853386 -0.365005 0.372160 2 0 +0.673887 -0.673887 0.302905 2 0 + + +0.923762 -0.382710 0.0140184 2 0 +0.853386 -0.365005 0.372160 2 0 +0.918846 0.00000 0.394617 2 0 +0.999665 0.00000 0.0259008 2 0 +0.923762 -0.382710 0.0140184 2 0 + + +0.853386 -0.365005 0.372160 2 0 +0.673887 -0.302905 0.673887 2 0 +0.707107 0.00000 0.707107 2 0 +0.918846 0.00000 0.394617 2 0 +0.853386 -0.365005 0.372160 2 0 + + +0.577350 0.00000 -0.577350 2 0 +0.943953 0.00000 -0.330081 2 0 +0.862418 0.379225 -0.335298 2 0 +0.577350 0.288675 -0.577350 2 0 +0.577350 0.00000 -0.577350 2 0 + + +0.943953 0.00000 -0.330081 2 0 +0.999665 0.00000 0.0259008 2 0 +0.923762 0.382710 0.0140184 2 0 +0.862418 0.379225 -0.335298 2 0 +0.943953 0.00000 -0.330081 2 0 + + +0.577350 0.288675 -0.577350 2 0 +0.862418 0.379225 -0.335298 2 0 +0.673887 0.673887 -0.302905 2 0 +0.577350 0.577350 -0.577350 2 0 +0.577350 0.288675 -0.577350 2 0 + + +0.862418 0.379225 -0.335298 2 0 +0.923762 0.382710 0.0140184 2 0 +0.707107 0.707107 0.00000 2 0 +0.673887 0.673887 -0.302905 2 0 +0.862418 0.379225 -0.335298 2 0 + + +0.999665 0.00000 0.0259008 2 0 +0.918846 0.00000 0.394617 2 0 +0.853386 0.365005 0.372160 2 0 +0.923762 0.382710 0.0140184 2 0 +0.999665 0.00000 0.0259008 2 0 + + +0.918846 0.00000 0.394617 2 0 +0.707107 0.00000 0.707107 2 0 +0.673887 0.302905 0.673887 2 0 +0.853386 0.365005 0.372160 2 0 +0.918846 0.00000 0.394617 2 0 + + +0.923762 0.382710 0.0140184 2 0 +0.853386 0.365005 0.372160 2 0 +0.673887 0.673887 0.302905 2 0 +0.707107 0.707107 0.00000 2 0 +0.923762 0.382710 0.0140184 2 0 + + +0.853386 0.365005 0.372160 2 0 +0.673887 0.302905 0.673887 2 0 +0.577350 0.577350 0.577350 2 0 +0.673887 0.673887 0.302905 2 0 +0.853386 0.365005 0.372160 2 0 + + +-0.577350 -0.577350 0.577350 2 0 +-0.673887 -0.302905 0.673887 2 0 +-0.365731 -0.365731 0.855851 2 0 +-0.302905 -0.673887 0.673887 2 0 +-0.577350 -0.577350 0.577350 2 0 + + +-0.673887 -0.302905 0.673887 2 0 +-0.707107 0.00000 0.707107 2 0 +-0.382683 0.00000 0.923880 2 0 +-0.365731 -0.365731 0.855851 2 0 +-0.673887 -0.302905 0.673887 2 0 + + +-0.302905 -0.673887 0.673887 2 0 +-0.365731 -0.365731 0.855851 2 0 +0.00000 -0.382683 0.923880 2 0 +0.00000 -0.707107 0.707107 2 0 +-0.302905 -0.673887 0.673887 2 0 + + +-0.365731 -0.365731 0.855851 2 0 +-0.382683 0.00000 0.923880 2 0 +0.00000 0.00000 1.00000 2 0 +0.00000 -0.382683 0.923880 2 0 +-0.365731 -0.365731 0.855851 2 0 + + +-0.707107 0.00000 0.707107 2 0 +-0.673887 0.302905 0.673887 2 0 +-0.365731 0.365731 0.855851 2 0 +-0.382683 0.00000 0.923880 2 0 +-0.707107 0.00000 0.707107 2 0 + + +-0.673887 0.302905 0.673887 2 0 +-0.577350 0.577350 0.577350 2 0 +-0.302905 0.673887 0.673887 2 0 +-0.365731 0.365731 0.855851 2 0 +-0.673887 0.302905 0.673887 2 0 + + +-0.382683 0.00000 0.923880 2 0 +-0.365731 0.365731 0.855851 2 0 +0.00000 0.382683 0.923880 2 0 +0.00000 0.00000 1.00000 2 0 +-0.382683 0.00000 0.923880 2 0 + + +-0.365731 0.365731 0.855851 2 0 +-0.302905 0.673887 0.673887 2 0 +0.00000 0.707107 0.707107 2 0 +0.00000 0.382683 0.923880 2 0 +-0.365731 0.365731 0.855851 2 0 + + +0.00000 -0.707107 0.707107 2 0 +0.00000 -0.382683 0.923880 2 0 +0.365731 -0.365731 0.855851 2 0 +0.302905 -0.673887 0.673887 2 0 +0.00000 -0.707107 0.707107 2 0 + + +0.00000 -0.382683 0.923880 2 0 +0.00000 0.00000 1.00000 2 0 +0.382683 0.00000 0.923880 2 0 +0.365731 -0.365731 0.855851 2 0 +0.00000 -0.382683 0.923880 2 0 + + +0.302905 -0.673887 0.673887 2 0 +0.365731 -0.365731 0.855851 2 0 +0.673887 -0.302905 0.673887 2 0 +0.577350 -0.577350 0.577350 2 0 +0.302905 -0.673887 0.673887 2 0 + + +0.365731 -0.365731 0.855851 2 0 +0.382683 0.00000 0.923880 2 0 +0.707107 0.00000 0.707107 2 0 +0.673887 -0.302905 0.673887 2 0 +0.365731 -0.365731 0.855851 2 0 + + +0.00000 0.00000 1.00000 2 0 +0.00000 0.382683 0.923880 2 0 +0.365731 0.365731 0.855851 2 0 +0.382683 0.00000 0.923880 2 0 +0.00000 0.00000 1.00000 2 0 + + +0.00000 0.382683 0.923880 2 0 +0.00000 0.707107 0.707107 2 0 +0.302905 0.673887 0.673887 2 0 +0.365731 0.365731 0.855851 2 0 +0.00000 0.382683 0.923880 2 0 + + +0.382683 0.00000 0.923880 2 0 +0.365731 0.365731 0.855851 2 0 +0.673887 0.302905 0.673887 2 0 +0.707107 0.00000 0.707107 2 0 +0.382683 0.00000 0.923880 2 0 + + +0.365731 0.365731 0.855851 2 0 +0.302905 0.673887 0.673887 2 0 +0.577350 0.577350 0.577350 2 0 +0.673887 0.302905 0.673887 2 0 +0.365731 0.365731 0.855851 2 0 + + +-0.577350 -0.577350 -0.577350 2 0 +-0.577350 -0.288675 -0.577350 2 0 +-0.862418 -0.379225 -0.335298 2 0 +-0.673887 -0.673887 -0.302905 2 0 +-0.577350 -0.577350 -0.577350 2 0 + + +-0.577350 -0.288675 -0.577350 2 0 +-0.577350 0.00000 -0.577350 2 0 +-0.943953 0.00000 -0.330081 2 0 +-0.862418 -0.379225 -0.335298 2 0 +-0.577350 -0.288675 -0.577350 2 0 + + +-0.673887 -0.673887 -0.302905 2 0 +-0.862418 -0.379225 -0.335298 2 0 +-0.923762 -0.382710 0.0140184 2 0 +-0.707107 -0.707107 0.00000 2 0 +-0.673887 -0.673887 -0.302905 2 0 + + +-0.862418 -0.379225 -0.335298 2 0 +-0.943953 0.00000 -0.330081 2 0 +-0.999665 0.00000 0.0259008 2 0 +-0.923762 -0.382710 0.0140184 2 0 +-0.862418 -0.379225 -0.335298 2 0 + + +-0.577350 0.00000 -0.577350 2 0 +-0.577350 0.288675 -0.577350 2 0 +-0.862418 0.379225 -0.335298 2 0 +-0.943953 0.00000 -0.330081 2 0 +-0.577350 0.00000 -0.577350 2 0 + + +-0.577350 0.288675 -0.577350 2 0 +-0.577350 0.577350 -0.577350 2 0 +-0.673887 0.673887 -0.302905 2 0 +-0.862418 0.379225 -0.335298 2 0 +-0.577350 0.288675 -0.577350 2 0 + + +-0.943953 0.00000 -0.330081 2 0 +-0.862418 0.379225 -0.335298 2 0 +-0.923762 0.382710 0.0140184 2 0 +-0.999665 0.00000 0.0259008 2 0 +-0.943953 0.00000 -0.330081 2 0 + + +-0.862418 0.379225 -0.335298 2 0 +-0.673887 0.673887 -0.302905 2 0 +-0.707107 0.707107 0.00000 2 0 +-0.923762 0.382710 0.0140184 2 0 +-0.862418 0.379225 -0.335298 2 0 + + +-0.707107 -0.707107 0.00000 2 0 +-0.923762 -0.382710 0.0140184 2 0 +-0.853386 -0.365005 0.372160 2 0 +-0.673887 -0.673887 0.302905 2 0 +-0.707107 -0.707107 0.00000 2 0 + + +-0.923762 -0.382710 0.0140184 2 0 +-0.999665 0.00000 0.0259008 2 0 +-0.918846 0.00000 0.394617 2 0 +-0.853386 -0.365005 0.372160 2 0 +-0.923762 -0.382710 0.0140184 2 0 + + +-0.673887 -0.673887 0.302905 2 0 +-0.853386 -0.365005 0.372160 2 0 +-0.673887 -0.302905 0.673887 2 0 +-0.577350 -0.577350 0.577350 2 0 +-0.673887 -0.673887 0.302905 2 0 + + +-0.853386 -0.365005 0.372160 2 0 +-0.918846 0.00000 0.394617 2 0 +-0.707107 0.00000 0.707107 2 0 +-0.673887 -0.302905 0.673887 2 0 +-0.853386 -0.365005 0.372160 2 0 + + +-0.999665 0.00000 0.0259008 2 0 +-0.923762 0.382710 0.0140184 2 0 +-0.853386 0.365005 0.372160 2 0 +-0.918846 0.00000 0.394617 2 0 +-0.999665 0.00000 0.0259008 2 0 + + +-0.923762 0.382710 0.0140184 2 0 +-0.707107 0.707107 0.00000 2 0 +-0.673887 0.673887 0.302905 2 0 +-0.853386 0.365005 0.372160 2 0 +-0.923762 0.382710 0.0140184 2 0 + + +-0.918846 0.00000 0.394617 2 0 +-0.853386 0.365005 0.372160 2 0 +-0.673887 0.302905 0.673887 2 0 +-0.707107 0.00000 0.707107 2 0 +-0.918846 0.00000 0.394617 2 0 + + +-0.853386 0.365005 0.372160 2 0 +-0.673887 0.673887 0.302905 2 0 +-0.577350 0.577350 0.577350 2 0 +-0.673887 0.302905 0.673887 2 0 +-0.853386 0.365005 0.372160 2 0 + + +-0.577350 -0.577350 -0.577350 2 0 +-0.673887 -0.673887 -0.302905 2 0 +-0.379225 -0.862418 -0.335298 2 0 +-0.288675 -0.577350 -0.577350 2 0 +-0.577350 -0.577350 -0.577350 2 0 + + +-0.673887 -0.673887 -0.302905 2 0 +-0.707107 -0.707107 0.00000 2 0 +-0.382710 -0.923762 0.0140184 2 0 +-0.379225 -0.862418 -0.335298 2 0 +-0.673887 -0.673887 -0.302905 2 0 + + +-0.288675 -0.577350 -0.577350 2 0 +-0.379225 -0.862418 -0.335298 2 0 +0.00000 -0.943953 -0.330081 2 0 +0.00000 -0.577350 -0.577350 2 0 +-0.288675 -0.577350 -0.577350 2 0 + + +-0.379225 -0.862418 -0.335298 2 0 +-0.382710 -0.923762 0.0140184 2 0 +0.00000 -0.999665 0.0259008 2 0 +0.00000 -0.943953 -0.330081 2 0 +-0.379225 -0.862418 -0.335298 2 0 + + +-0.707107 -0.707107 0.00000 2 0 +-0.673887 -0.673887 0.302905 2 0 +-0.365005 -0.853386 0.372160 2 0 +-0.382710 -0.923762 0.0140184 2 0 +-0.707107 -0.707107 0.00000 2 0 + + +-0.673887 -0.673887 0.302905 2 0 +-0.577350 -0.577350 0.577350 2 0 +-0.302905 -0.673887 0.673887 2 0 +-0.365005 -0.853386 0.372160 2 0 +-0.673887 -0.673887 0.302905 2 0 + + +-0.382710 -0.923762 0.0140184 2 0 +-0.365005 -0.853386 0.372160 2 0 +0.00000 -0.918846 0.394617 2 0 +0.00000 -0.999665 0.0259008 2 0 +-0.382710 -0.923762 0.0140184 2 0 + + +-0.365005 -0.853386 0.372160 2 0 +-0.302905 -0.673887 0.673887 2 0 +0.00000 -0.707107 0.707107 2 0 +0.00000 -0.918846 0.394617 2 0 +-0.365005 -0.853386 0.372160 2 0 + + +0.00000 -0.577350 -0.577350 2 0 +0.00000 -0.943953 -0.330081 2 0 +0.379225 -0.862418 -0.335298 2 0 +0.288675 -0.577350 -0.577350 2 0 +0.00000 -0.577350 -0.577350 2 0 + + +0.00000 -0.943953 -0.330081 2 0 +0.00000 -0.999665 0.0259008 2 0 +0.382710 -0.923762 0.0140184 2 0 +0.379225 -0.862418 -0.335298 2 0 +0.00000 -0.943953 -0.330081 2 0 + + +0.288675 -0.577350 -0.577350 2 0 +0.379225 -0.862418 -0.335298 2 0 +0.673887 -0.673887 -0.302905 2 0 +0.577350 -0.577350 -0.577350 2 0 +0.288675 -0.577350 -0.577350 2 0 + + +0.379225 -0.862418 -0.335298 2 0 +0.382710 -0.923762 0.0140184 2 0 +0.707107 -0.707107 0.00000 2 0 +0.673887 -0.673887 -0.302905 2 0 +0.379225 -0.862418 -0.335298 2 0 + + +0.00000 -0.999665 0.0259008 2 0 +0.00000 -0.918846 0.394617 2 0 +0.365005 -0.853386 0.372160 2 0 +0.382710 -0.923762 0.0140184 2 0 +0.00000 -0.999665 0.0259008 2 0 + + +0.00000 -0.918846 0.394617 2 0 +0.00000 -0.707107 0.707107 2 0 +0.302905 -0.673887 0.673887 2 0 +0.365005 -0.853386 0.372160 2 0 +0.00000 -0.918846 0.394617 2 0 + + +0.382710 -0.923762 0.0140184 2 0 +0.365005 -0.853386 0.372160 2 0 +0.673887 -0.673887 0.302905 2 0 +0.707107 -0.707107 0.00000 2 0 +0.382710 -0.923762 0.0140184 2 0 + + +0.365005 -0.853386 0.372160 2 0 +0.302905 -0.673887 0.673887 2 0 +0.577350 -0.577350 0.577350 2 0 +0.673887 -0.673887 0.302905 2 0 +0.365005 -0.853386 0.372160 2 0 + + +-0.577350 0.577350 -0.577350 2 0 +-0.288675 0.577350 -0.577350 2 0 +-0.379225 0.862418 -0.335298 2 0 +-0.673887 0.673887 -0.302905 2 0 +-0.577350 0.577350 -0.577350 2 0 + + +-0.288675 0.577350 -0.577350 2 0 +0.00000 0.577350 -0.577350 2 0 +0.00000 0.943953 -0.330081 2 0 +-0.379225 0.862418 -0.335298 2 0 +-0.288675 0.577350 -0.577350 2 0 + + +-0.673887 0.673887 -0.302905 2 0 +-0.379225 0.862418 -0.335298 2 0 +-0.382710 0.923762 0.0140184 2 0 +-0.707107 0.707107 0.00000 2 0 +-0.673887 0.673887 -0.302905 2 0 + + +-0.379225 0.862418 -0.335298 2 0 +0.00000 0.943953 -0.330081 2 0 +0.00000 0.999665 0.0259008 2 0 +-0.382710 0.923762 0.0140184 2 0 +-0.379225 0.862418 -0.335298 2 0 + + +0.00000 0.577350 -0.577350 2 0 +0.288675 0.577350 -0.577350 2 0 +0.379225 0.862418 -0.335298 2 0 +0.00000 0.943953 -0.330081 2 0 +0.00000 0.577350 -0.577350 2 0 + + +0.288675 0.577350 -0.577350 2 0 +0.577350 0.577350 -0.577350 2 0 +0.673887 0.673887 -0.302905 2 0 +0.379225 0.862418 -0.335298 2 0 +0.288675 0.577350 -0.577350 2 0 + + +0.00000 0.943953 -0.330081 2 0 +0.379225 0.862418 -0.335298 2 0 +0.382710 0.923762 0.0140184 2 0 +0.00000 0.999665 0.0259008 2 0 +0.00000 0.943953 -0.330081 2 0 + + +0.379225 0.862418 -0.335298 2 0 +0.673887 0.673887 -0.302905 2 0 +0.707107 0.707107 0.00000 2 0 +0.382710 0.923762 0.0140184 2 0 +0.379225 0.862418 -0.335298 2 0 + + +-0.707107 0.707107 0.00000 2 0 +-0.382710 0.923762 0.0140184 2 0 +-0.365005 0.853386 0.372160 2 0 +-0.673887 0.673887 0.302905 2 0 +-0.707107 0.707107 0.00000 2 0 + + +-0.382710 0.923762 0.0140184 2 0 +0.00000 0.999665 0.0259008 2 0 +0.00000 0.918846 0.394617 2 0 +-0.365005 0.853386 0.372160 2 0 +-0.382710 0.923762 0.0140184 2 0 + + +-0.673887 0.673887 0.302905 2 0 +-0.365005 0.853386 0.372160 2 0 +-0.302905 0.673887 0.673887 2 0 +-0.577350 0.577350 0.577350 2 0 +-0.673887 0.673887 0.302905 2 0 + + +-0.365005 0.853386 0.372160 2 0 +0.00000 0.918846 0.394617 2 0 +0.00000 0.707107 0.707107 2 0 +-0.302905 0.673887 0.673887 2 0 +-0.365005 0.853386 0.372160 2 0 + + +0.00000 0.999665 0.0259008 2 0 +0.382710 0.923762 0.0140184 2 0 +0.365005 0.853386 0.372160 2 0 +0.00000 0.918846 0.394617 2 0 +0.00000 0.999665 0.0259008 2 0 + + +0.382710 0.923762 0.0140184 2 0 +0.707107 0.707107 0.00000 2 0 +0.673887 0.673887 0.302905 2 0 +0.365005 0.853386 0.372160 2 0 +0.382710 0.923762 0.0140184 2 0 + + +0.00000 0.918846 0.394617 2 0 +0.365005 0.853386 0.372160 2 0 +0.302905 0.673887 0.673887 2 0 +0.00000 0.707107 0.707107 2 0 +0.00000 0.918846 0.394617 2 0 + + +0.365005 0.853386 0.372160 2 0 +0.673887 0.673887 0.302905 2 0 +0.577350 0.577350 0.577350 2 0 +0.302905 0.673887 0.673887 2 0 +0.365005 0.853386 0.372160 2 0 + +