Dear Deal.II community, I'm getting trouble in changing the values of my solution variable Phi. Firstly, I create a R20 circle whose center is located at the center of the square domain by using VectorTools::interpolate(dof_handler, user_defined_function, solution_Phi). After creating this circle and several program loops, the radius of this circle grows up to about 40, as shown below, solution_Phi inside the circle is 1(denoted by red color) while outside the circle is -1(denoted by blue color). Here comes my problem, are there any practical ways to assign those solution_Phi who is located at the top half of the square domain(y>=60) to -1 while others(y<60) keep unchanged? [image: 1.jpeg]
Any suggestions or help will be appreciated! Thank you very much. Best, Toddy -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/ac2bfa5a-fac7-4629-a4c3-529765cc7db8n%40googlegroups.com.
