[Dune] Using PDELab's alpha_skeleton with larger stencil?

Christoph GrĂ¼ninger foss at grueninger.de
Tue Sep 11 06:40:22 CEST 2018


Hi Bryan,
as far as I know, this is not possible.

> FE and FV schemes requiring at most face-neighbors

https://www.dune-project.org/modules/dune-pdelab/

The rationale is, that this limitation keeps the interface clean, the number of evaluations low and the sparsity pattern known.

Bye
Christoph, who struggled with this limitation already

> Am 11.09.2018 um 04:09 schrieb Bryan Doyle <btd2 at rice.edu>:
> 
> Hello,
> I am using DUNE's PDELab, and noticed that in the alpha_skeleton calls, only 2 coefficient vectors are available: x_s and x_n. This implies that only the unknowns evaluated on the self and neighbor cells are available for use in alpha_skeleton's calculations. Suppose my weak form had an integral over an edge which depended on the unknowns evaluated on cells other than the self and neighbor elements; would it be possible to implement this using PDELab?
> 
> As a more concrete example, assume I have a uniform square grid, with elements E_1 and E_2 which do not share an edge. Suppose I wish to calculate
> \int_e u_1*u_2,
> where u_1 and u_2 are the unknowns evaluated in the centers of elements E_1 and E_2, respectively, and e is some edge of E_1. Is calculating this integral possible using PDELab's alpha_skeleton function, or is there another way to do so?
> 
> Any advice or input would be greatly appreciated!
> 
> Thank you for your time,
> Bryan
> _______________________________________________
> Dune mailing list
> Dune at lists.dune-project.org
> https://lists.dune-project.org/mailman/listinfo/dune
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20180911/f73f7489/attachment.htm>


More information about the Dune mailing list