added lav_matrix_bdiag() function
further tricks in lav_matrix.R