Compute (a[0] + a[1]) (b[0] + b[1]) in extended precision. @param a first term of the multiplication @param b second term of the multiplication @param result placeholder where to put the result
(final double a[], final double b[], final double result[])