<div dir="ltr"><div class="gmail_quote">Hi all,<br><div dir="ltr"><div>I really need help with defining my domain in dune. </div><div><br></div><div>I have a problem similar to example 5 of dune pdelab howto (a two-component diffusion-reaction problem) but with different domain. </div>

<div>In howto problem the domain is (0,2)^2 but mine is </div><div><div>{(x; y) : ((x-0.5)/0.25)^2+((y-0.5)/0.5)^2<1}</div><div><br></div></div><div>Also, in howto problem the domain is defined with: </div><div><div>Dune::FieldVector<double,2> L(2.0);</div>

<div>      Dune::FieldVector<int,2> N(1);</div><div>      Dune::FieldVector<bool,2> periodic(false);</div><div>      int overlap=0;</div><div>      Dune::YaspGrid<2> grid(L,N,periodic,overlap);</div><div>

      grid.globalRefine(level);</div><div>      typedef Dune::YaspGrid<2>::LeafGridView GV;</div><div>      const GV& gv=grid.leafView();</div></div><div><br></div><div>Now I shoud replace that YaspGrid (0,2)^2 with something of myown but I am unsure how.</div>

<div><br></div><div>Thank you in advance.</div><div><br></div><div>Best regards,</div><div>Mario</div></div>
</div><br></div>