add documentation on the effect of ConstraintMatrix::distribute_local_to...

Documentation update highlighting the difference between the two ConstraintMatrix::distribute_local_to_global() functions for matrices and their effect on the spectrum of a matrix
The documentation of diag_element() says that you can't call this function for non-square matrices.

What's a "source problem"? A regular linear system of the form $Ax=b$?