added matrix base and matrix-vector product, now ready for generic
multiplication of index/position vectors
added test of img iterator
rethrow same exception, just in case in matrix-vector product
going to switch to single const and iter types rather than one for each
type
unified iterators, think this is more intuitive
fix to writer, removed a little debugging