B-spline surface interpolation via 2D de Boor's algorithm

Question

I am implementing the de Boor\'s algorithm for B-spline fitting for 2D surfaces.

1D case (Correct)

I modified the 1D case from Wikipedia, test it against sci