[Dune] Dirichlet conditions on curved boundaries

Andreas Dedner a.s.dedner at warwick.ac.uk
Mon Jun 1 22:10:18 CEST 2015


Hi.

On 01/06/15 20:46, conf86 at web.de wrote:
> if (x[0]*x[0]+x[1]*x[1]==1) y = 1.0;
Never directly compare doubles - use something like abs(...-1,0) < eps.
But you have to use the intersection iterator to find the boundary (use 
it->boundary())
and then you can set your data.
Andreas




More information about the Dune mailing list