[Dune] Dirichlet conditions on curved boundaries

conf86 at web.de conf86 at web.de
Tue Jun 2 00:03:59 CEST 2015


Thank you for your answer.

> Never directly compare doubles - use something like abs(...-1,0) < eps.
Ok, now it makes sense to me not to do this.

> But you have to use the intersection iterator to find the boundary (use 
> it->boundary())
> and then you can set your data.
So basically, I just need to replace the condition?

 if (it->boundary()) y=1.0;

Sorry, if this might be a stupid question, but I couldn't find any sample code for this.

Best regards.



More information about the Dune mailing list