If i’m correct, jump() and average() functions give the jump and mean value of a given function at an edge for discontinuous Galerkin schemes. But how to retrieve these fundamental quantities when working with a vector of functions ? I tried something like jump([u1,u2]) but this doesn’t seem to work. Any idea ?