gi.attach_triangulation (tr);
std::ifstream in ("../deal.II/grid_in_3d_02/747.ucd");
gi.read (in);
- tr.refine_global (1);
+ //tr.refine_global (1);
}
{
gi.attach_triangulation (tr2);
std::ifstream in ("../deal.II/grid_in_3d_02/747.ucd");
gi.read (in);
- tr2.refine_global (1);
+ //tr2.refine_global (1);
}
Assert (tr.n_active_cells() == tr2.n_active_cells(),
ExcInternalError());
+ deallog << " Number of cells: "
+ << tr.n_active_cells() << ' '
+ << tr2.n_active_cells()
+ << std::endl;
- for (unsigned int i=0; i<1; ++i)
+ for (unsigned int i=0; i<2; ++i)
{
// refine one-tenth of cells randomly
std::vector<bool> flags (tr.n_active_cells(), false);
- for (unsigned int k=0; k<flags.size()/10 + 1; ++k)
+ for (unsigned int k=0; k<flags.size()/30 + 1; ++k)
flags[rand() % flags.size()] = true;
// make sure there's at least one that
// will be refined
tr.execute_coarsening_and_refinement ();
tr2.execute_coarsening_and_refinement ();
- deallog << std::endl;
-
deallog << i << " Number of cells: "
<< tr.n_active_cells() << ' '
<< tr2.n_active_cells()