Hi, Wolfgang
>> I found one problem with inp file - in line "1 0 quad 1 2 3 4". I
>> guess the second number should be nonzero here. I have try to set 1 and
>> get better results. Now I can see velocity changes around the
>> input-output bounds.
> Is this an input file you created yourself, or are you saying that there
> is a mistake in the input files we distribute as part of step-35?
I created this input file by myself. It's my fault I used only following
code for set bounds:
Triangulation<2>::active_cell_iterator
cell = triangulation.begin_active(),
endc = triangulation.end();
for (; cell != endc; ++cell)
{
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
{
if (cell->face(f)->at_boundary())
{
cell->face(f)->set_boundary_id(2);
...
Now I think I need to use Triangulation<2>::::active_quad_iterator to set
boundary id to every quad.
Thank you!
E.P.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
