Internal Function
It efficiently calculates the working model quantities needed to implement the automatic multiple smoothing parameter estimation procedure by exploiting a result which leads to very fast and stable calculations.
Maintainer: Giampiero Marra giampiero.marra@ucl.ac.uk
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.