SAP function

Define a two-dimensional penalised tensor-product of marginal B-Spline basis functions.