So far we used Qk^d for representing the physical flux. However, we know
that the physical flux is basically the gradient of q Qk function.
To consider this, we add a finite element representing gradients of Q1
functions and use these for interpolating the physical flux. This
greatly reduces the number of basis functions and thus the time and
memory consumption for MEG with physical flux