From: Guido Kanschat Date: Thu, 29 Dec 2005 09:33:55 +0000 (+0000) Subject: cylinder_shell fixed X-Git-Tag: v8.0.0~12726 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d7486d06c1c43f4afaeaa369d37d791d6e633df7;p=dealii.git cylinder_shell fixed git-svn-id: https://svn.dealii.org/trunk@11933 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/grid/grid_generator.cc b/deal.II/deal.II/source/grid/grid_generator.cc index a976647ff0..336fd42e4c 100644 --- a/deal.II/deal.II/source/grid/grid_generator.cc +++ b/deal.II/deal.II/source/grid/grid_generator.cc @@ -1384,20 +1384,20 @@ void GridGenerator::cylinder_shell (Triangulation &tria, // first N ones are on the // outer one, and all are // numbered counter-clockwise - std::vector > vertices_2d(2*N_r); + std::vector > vertices_2d(2*N_r); for (unsigned int i=0; i( std::cos(2*pi * i/N_r), + vertices_2d[i] = Point<2>( std::cos(2*pi * i/N_r), std::sin(2*pi * i/N_r)) * outer_radius; vertices_2d[i+N_r] = vertices_2d[i] * (inner_radius/outer_radius); }; - std::vector > vertices_3d; + std::vector > vertices_3d; vertices_3d.reserve (2*N_r*(N_z+1)); for (unsigned int j=0; j<=N_z; ++j) for (unsigned int i=0; i<2*N_r; ++i) { - const Point v (vertices_2d[i][0], + const Point<3> v (vertices_2d[i][0], vertices_2d[i][1], j*length/N_z); vertices_3d.push_back (v);