Some flux schemes need a term $$\int_\Gamma (\partial_n \tilde u) f_2~dS$$ This is complicated, since nodes on the interior have nonzero $\frac{\partial \tilde {u}}{\partial n}$. ## Tasks - [ ] create a subroutine in the integration algorithms for this computation - [ ] write a test for it ## Acceptance Criteria - [ ] tests pass