[dune-pdelab] Finite Volume Method with slope reconstruction / alpha_volume_post_skeleton

[Markus Blatt]

Hi,

On Tue, May 09, 2017 at 01:54:53PM +0200, Gregor Corbin wrote:
> 
> As we discussed before, it uses the alpha_skeleton method to compute slopes
> at each face, and stores this information. Then in the
> alpha_volume_post_skeleton method, it uses the previously computed slopes.
> 
> But unfortunately, the alpha_volume_post_skeleton will not be executed after
> all the skeleton terms are done, but rather in between. To test this, I
> created a minimal 3x3 YASP-Grid  example and let the methods of the local
> operator print their name and some info about the current cell.
>

The problem might be that there is only one loop over the grid when
assembling the residual. That means for each codim-0 entity the
alpha/lamba methods are applied consecutively.

I guess in your case you would require more passes over the grid as
for an intersection you would need precomputed values of both
neighbors.

You might benefit from hijacking the preAssembly method, too. See my
mail/question from yesterday.

Markus


-- 
Dr. Markus Blatt - HPC-Simulation-Software & Services http://www.dr-blatt.de
Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany
Tel.: +49 (0) 160 97590858
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 819 bytes
Desc: Digital signature
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20170509/33106be6/attachment.sig>