void load_sigma2(const fmat &V, fvec &sigma2){ int nsnp = V.size(); sigma2 = fvec (nsnp, .0f); for(int i = 0; i < nsnp; ++i){ sigma2[i] = V[i][i]; } }