Fixed the fact that the C testing scheme was not actually comparing all the columns of a matrix, merely the first column!  Also, collected more matrix comparisons and made them all call the same function.  Similarly in the python tests, grouped the tests for matrices.

Added one new test based on a model Herbert set up.