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