Hamed,
> Thank you very much for your help. In fact, what you got is for the case that > the material_id for the entire domain (both half of cubes) is included in > function_map. My domain is a cube divided into two equal part at the plane > x=0.5. One half has material_id=0 and the other half has material_id=1. So if > you comment the following lines in Run: > > const ConstantFunction<dim> constant_function_2(2.0); > > function_map[0] = &constant_function_2; Ah, I see. I also see that you had left a comment there -- though, in the future, make it easy on people by making the default case the one that shows the problem. I've opened a github issue for this here: https://github.com/dealii/dealii/issues/7522 I don't have the time soon to debug this, but with a good testcase like yours it shouldn't be very difficult to do. Have you looked at the implementation of the function to see whether you can find yourself what the problem is? We'd of course be happy to help you find your way through the code if you run into problems! Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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.
