[Dune-devel] help with defining a ellipse grid in dune

[Mario Rogulj]

Hi all,
I really need help with defining my domain in dune.

I have a problem similar to example 5 of dune pdelab howto (a two-component
diffusion-reaction problem) but with different domain.
In howto problem the domain is (0,2)^2 but mine is
{(x; y) : ((x-0.5)/0.25)^2+((y-0.5)/0.5)^2<1}

Also, in howto problem the domain is defined with:
Dune::FieldVector<double,2> L(2.0);
      Dune::FieldVector<int,2> N(1);
      Dune::FieldVector<bool,2> periodic(false);
      int overlap=0;
      Dune::YaspGrid<2> grid(L,N,periodic,overlap);
      grid.globalRefine(level);
      typedef Dune::YaspGrid<2>::LeafGridView GV;
      const GV& gv=grid.leafView();

Now I shoud replace that YaspGrid (0,2)^2 with something of myown but I am
unsure how.

Thank you in advance.

Best regards,
Mario
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune-devel/attachments/20130909/153515a5/attachment.htm>