2D P-spline fitting with anisotropic penalties

Usage

pspline_fit_2d(
  data,
  argvals_s,
  argvals_t,
  nbasis_s,
  nbasis_t,
  lambda_s,
  lambda_t,
  order
)