fix numbering in local interpolant ( u^1_x, u^2_x, ..., u^N_x, u^1_y u^2 ... u^N_y )
add operator<< for Dof
Merge branch 'feature/interpolate' of https://github.com/feelpp/feelpp into feature/interpolate