[dune-pdelab] Finite Volume Method with slope reconstruction

[Christian Engwer]

Hi Gregor,

> If I use a grid operator for this, how do I get the solution on neighboring
> cells?  I suppose I could use the alpha_skeleton method to communicate this
> information and then doing the reconstruction in the
> alpha_volume_post_skeleton method. But this would mean a lot of storage and
> copying overhead.

I don't understand why you think you have to communicate information?!
You use alpha_skeleton to compute the slope w.r.t. each neighbor. You
can store this information in a stateful local operator and reuse the
information on the same cell in the post-skeleton.

> Also, the reconstruction has to be done before each stage of the time
> stepping. So would I have to call the reconstruction operator from the
> preStage() method of my original local operator?

Yes.

> Another question related to this: Are the intersections in a YASPGrid
> ordered in a consistent way? I.e. has intersections(grid_view,cell)[0]
> always the same normal?

Yes.

Ciao
Christian