What's Inside A variance component estimation method for sparse matrix applications | tbr*a