calculateConstructVCV function

Internal: Calculate construct variance-covariance matrix