Checked that Parallel computation of BivariateSpline takes a bit longer, keeping singlecore version
Now computing sigma0 for both polarisations