NumInt by gauss point. As many points as you like.
Fixed bug in MPI communication.
Only plot every i'th timestep.
Added all B3.c .. B7.c to library, And the same for mass matrix.