Sparse Fourier coefficients

Only consider non-zero Fourier coefficients in the computation when assembling the matrix. 
This would allow to cheaply explore high frequency perturbations to an equilibrium, or global perturbations (modes that are not multiples of NFP) 