Takes the average matrix Q and computes the jumps
giving matrix deltaQ as output. The entry delta_q[i, j] of the matrix contains the jump of the i-th unknown in the system at the interface j.
Release Notes: