Fix band matrix for addentries!

This changes over to a new data structure to represent Banded Matrices that is a lot faster for multiplication and a lot less hacked together.

It means up to 2x faster \ just due to TimesOperator being faster

The next step in this fork is to rewrite AlmostBandedOperator and givensreduce, which will hopefully lead to another 2x speed up (or more).