[dune-pdelab] Accessing solution vector for MPFA

[Christoph Grüninger]

Hi Shubhangi,
does your MPFA method fulfil "FE and FV schemes requiring at most
face-neighbors" from [1]?
My experience is, that complicated FV schemes exceed this restriction.
You can try to program around this inside the local operator (like your
pointer to the global solution, that should be never needed), but the
resulting stencil generated by PDELab will remain too small.

I still try to answer your question:
The local solution is altered to calculate the finite difference.
Usually, the local solution is changed by 1e-6, 1e-8 or similar.

Bye
Christoph

[1] https://www.dune-project.org/modules/dune-pdelab/

-- 
The method has been developed for use on a high-speed electronic
computer and would be impractical for hand-solution purposes.
                                           [Harlow & Welch 1965]

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 819 bytes
Desc: OpenPGP digital signature
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20160920/1f8ca714/attachment.sig>