Dear J-P and Wolfgang, Thank you for your answers! As I see the actual problem is to define a step function for material parameters. In FEM simulations mesh is often aligned to the problem`s inner interfaces (like in a "Better mesh" extension of step-6). Is it possible to define for such a case with a discontinuous coefficient using continuous element? E.g. is it valid to take coefficient evaluation out of the q_index cycle and use
current_coefficient = coefficient<dim> (cell->center()) in step-6 code for the case of mesh aligned to the interface of coefficient drop? Best regards, Kostya -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.