fix getMaxDiagCoeff() in LisMatrix. prepare a diagonal vector by ourself...

a bug fix after #165. Because lis_matrix_get_diagonal() is not available before assembling a matrix, we prepare a full diagonal vector in our side. Any update on diagonal entries is stored both in a matrix and a vector.