Add method for extrapolation at boundaries
Add a brand new bicubic interpolant (not properly working yet :-)
A little refactoring
CubicInterpolant now can be built in constructor
CubicInterpolant now gets a matrix of nodes instead of row-by-row given vectors of nodes