Remove copy from add_to_diagonal
The recently-added function for FieldMatrix, add_to_diagonal( double ), should simply augment the diagonal of the FieldMatrix from which it is called, rather than returning a new FieldMatrix.
This is on the DenseLinAlgMatlabLikeOperators branch, which has not been merged to master yet.