Honor this preference in operator=. - add several methods to the API - rework API for diagonal matrices - add benchmarking code