`scipy.interpolate._fitpack:surfit_lsq`

Summary

[c,fp,ier] = surfit_lsq(x,y,z,w,xb,xe,yb,ye,kx,ky,s,nxest,nyest,nmax,eps,tx,ty)

This LSQ-only wrapper modifies tx and ty in-place (both must be writable float64 ndarrays of length nmax).