Quad4 finite element 

This PR is extracted from #100. Interpolation and extrapolation functionalities in #100 are removed as they are not used at the moment.

The codes are still not ready for review and waiting #203 and #205.