Laplacian matrices joint diagonalization