[Dune] Dirichlet conditions on curved boundaries

[conf86 at web.de]

I found the code, but I have the impression, I can't figure out how to use it with mine.
I am currently using a modified version of example02_bcextension.hh from the course-examples in dune-pdelab-howto. Do I need to fully rewrite bcextension.hh or is there a possibility to integrate IntersectionIterator into it?

> Gesendet: Dienstag, 02. Juni 2015 um 00:08 Uhr
> Von: "Andreas Dedner" <a.s.dedner at warwick.ac.uk>
> An: conf86 at web.de
> Betreff: Re: [Dune] Dirichlet conditions on curved boundaries
>
> Well not quite - you need to know which dofs lie on that intersection 
> (an edge of the boundary in
> 2d or a face in 3d) and then you set those dofs.
> In the dune-grid-howto there is a finiteelement.cc example which should 
> give you an idea.
> Andreas
> 
> On 01/06/15 23:03, conf86 at web.de wrote:
> > 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.
> >
> > _______________________________________________
> > Dune mailing list
> > Dune at dune-project.org
> > http://lists.dune-project.org/mailman/listinfo/dune
> 
>