`scipy.interpolate._fitpack:surfit`

Summary

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

This is the legacy bisplrep-compatible wrapper that handles iopt=0,1,-1 and maintains wrk state.