Ulrich, > The error occurs in the debuging mode in the case of non-global (but > without hanging nodes) refinement with cut_y using > cell->set_refine_flag( RefinementCase<2>::cut_y ); > and > m_triangulation.execute_coarsening_and_refinement (); > in combination with the use of Raviart-Thomas elements.
I'm afraid the Raviart-Thomas element has never been converted to support anisotropic refinement :-( > The error does not occur without the debuging mode. In this case the > program runs and produces at least plausible results. The check isn't run in non-debug mode, so the program would not abort. I would at least doubt the correctness of results in that case, though it is entirely possible that the check is actually in error -- it was not written with anisotropic refinement in mind... I guess we'd need a volunteer to try to understand what would have to be done to make this work. > The error does not occur in the case of the deal.II tutorial Step-3 even > for non-global refinement with cut_y. A probable reason is that this > step does not use Raviart-Thomas elements. Yes. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
