The problem is generally the matrix corresponding to the fourier mode is full in no fourier cas
so the matrix A associe to “a” is full, it is very bad to solve the problem, but you can solve a problem with bilinear form a , without construction the matrix A only the operator A*u with CG or GMRES algorithm.