<div dir="ltr">Hi Christoph,<div>Thank you for the prompt response! As a follow-up question, were you able to find a work around when you encountered this limitation? Were there any specific strategies that you looked into?</div><div><br></div><div>In my particular case, I only have a handful of integrals of the form I described above, so I was thinking that it may be possible to add in a block of code (possibly somewhere in the Newton solver?) that looped over the mesh edges after each iteration and added these terms to the residual vector, outside of the calls to alpha_skeleton. Would this be something that (a) could work in theory, and (b) could be implemented without developer-level knowledge of the DUNE code?</div><div><br></div><div>Thank you for your time,</div><div>Bryan</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Sep 10, 2018 at 11:40 PM Christoph Grüninger <<a href="mailto:foss@grueninger.de" target="_blank">foss@grueninger.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><div></div><div>Hi Bryan,</div><div>as far as I know, this is not possible.</div><div><br></div><div><span style="background-color:rgba(255,255,255,0)">> FE and FV schemes requiring at most face-neighbors</span></div><div><span style="background-color:rgba(255,255,255,0)"><br></span></div><div><a href="https://www.dune-project.org/modules/dune-pdelab/" target="_blank">https://www.dune-project.org/modules/dune-pdelab/</a></div><div><br></div><div>The rationale is, that this limitation keeps the interface clean, the number of evaluations low and the sparsity pattern known.</div><div><br></div><div>Bye</div><div>Christoph, who struggled with this limitation already</div><div><br>Am 11.09.2018 um 04:09 schrieb Bryan Doyle <<a href="mailto:btd2@rice.edu" target="_blank">btd2@rice.edu</a>>:<br><br></div><blockquote type="cite"><div><div dir="ltr">Hello,<div>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?</div><div><br></div><div>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</div><div>\int_e u_1*u_2,</div><div>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?</div><div><br></div><div>Any advice or input would be greatly appreciated!</div><div><br></div><div>Thank you for your time,</div><div>Bryan</div></div>
</div></blockquote><blockquote type="cite"><div><span>_______________________________________________</span><br><span>Dune mailing list</span><br><span><a href="mailto:Dune@lists.dune-project.org" target="_blank">Dune@lists.dune-project.org</a></span><br><span><a href="https://lists.dune-project.org/mailman/listinfo/dune" target="_blank">https://lists.dune-project.org/mailman/listinfo/dune</a></span></div></blockquote></div></blockquote></div>