working.comp1: Internal Function
Description
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.