change data_index= to index=, and better tests
Changed my mind and reverted the addition of index=

Instead I think what I really want is proper support for ancillary
matrices, which provide an elegant solution to the "implicit data"
problem I'm running into in pyrerp, and are also useful for things
like weights= arguments.

Did keep some parts of the changes, though; see changes.rst.