Fast smoothing parameter separation in multidimensional generalized P-splines: the SAP algorithm | Publicación