output of CouplingOperator functions are now directly sparse matrices
got rid of all FE2D routines that were conserved on git and points directly to the original source
implemented the FE2D case

however, the alpha function (and variable parameters in general) still
have to be introduced