reduceMatrix function

Apply reduction method from Scholtus (2008)