Compute the variance-covariance matrix of the incomplete score vector involved in Louis' identity for the observed information matrix
invmatgen(G,H,listE)
The inverse of matrix G + H computed recursively thanks to matrices listed in listE
Primary matrix for the sum decomposition of $G+H$
Secondary matrix for the sum decomposition of $G+H$
Auxiliary matrices that sum up to H
Miller K. (1981). On the inverse of the sum of matrices,
Mathematics Magazine, 54, 67--72
fastCUB