added lav_matrix_duplication_ginv() to compute the generalized inverse of the duplication matrix