Internal Function
It provides an overall penalty matrix in a format suitable for estimation conditional on smoothing parameters.
Maintainer: Giampiero Marra giampiero.marra@ucl.ac.uk