Wrote a unified version of the vectorization test to be used as a basis for the tests of @jwillis, @rgb and @alepper
The test also contains a brute-force calculation to check that any approximation made even in the stock DOPAIR1() function is correct.