Integration on bulk and virtual interface terms are performed by a
matrix vector product. Add 2 functions in fullMatrix.h and cpp to
realize this operation