Dear Wolfgang, > In other words, there are no NaNs for me. What do you get? And what > version of > deal.II are you using? > 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; the half with material_id=0 has not included in the function_map and consequently, it leads to the -nan values. P.S. for the dst vector elements, I considered constant value 1 and 2 for cells with material_id=1 and 0 respectively. Thanks and regards, -- 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.
