Hello Christian,
could you please provide the vertex and cell data? I think there the
error should lie.
Best Regards,
Markus
Am 27.09.10 11:36, schrieb Christian May:
Dear Wolfgang,
thanks for your answer. The way I create the triangulation is the following:
///////////////////// code snippet
std::vector< Point<dim> > vertices;
std::vector< CellData<dim> > cells;
// code which reads in vertices and cells from a file
triangulation.create_triangulation(vertices,cells,SubCellData());
for(typename DoFHandler<dim>::active_cell_iterator
cell=dof_handler.begin_active();cell!=dof_handler.end();++cell) {
for(int f=0;f<GeometryInfo<dim>::faces_per_cell;++f) {
const TriaIterator<dealii::DoFAccessor<dim-1,
dealii::DoFHandler<dim, dim> > >
currentface=cell->face(f);
const int indicator=currentface->boundary_indicator();
std::cout<<"\t"<<indicator;
}
std::cout<< std::endl;
}
///////////////////// end code snippet
The grid it creates looks reasonable and has the outer shape it is
supposed to have. However, if I run this code, which shows a line for
each cell with boundary indicators of all faces,
I mostly get lines like the following:
0 0 255 255 255 255
In particular, if I do a grep -v 0, I realize that there is no element
without a 0 assigned face. In other words, according to this
triangulation all cells have a face on the outer boundary, which is
clearly not the way it should be.
I am aware that this code snippet is not sufficient to be executed. If
the problem lies somewhere in the way the vertex and cell data is
supplied, I can also provide that if you tell me the preferred format.
Many thanks
Christian
On Sun, Sep 26, 2010 at 3:10 AM, Wolfgang Bangerth
<[email protected]> wrote:
Christian,
I don't recall if anyone had already responded to this:
I have a list of vertices and cells in a file which I read into a vector
of Point and a vector of CellData. I then use these two parameters in
the invocation of create_triangulation, along with empty SubCellData.
My problem is now that imposing Dirichlet boundary conditions goes wrong.
Please excuse my ignorance as a new users of deal: I assumed that
create_triangulation would determine the boundary faces itself.
That's exactly what should happen. But you only say that it doesn't work, you
don't describe what actually happens. Can you elaborate?
W.
--
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
