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.

Reply via email to