pwUserFunctionFitSum

 residuals = pwUserFunctionFitSum(residuals, couples)

 User function to fit the sum of the observables over all couples
 to the data of couple 1.