void create_rp(double *X, int M, int Mp)
{
int j;
for (j = 0; j < M * Mp; j++)
X[j] = rand_normal_r();
}
int initializeVariables__f(double **U, double **Sigma, double **V, double **SNPSd, double **Cov, double **miss, double **mAF, int K, int nSNP, int nIND){
int i, j;
init_random();
*Sigma = calloc(K, sizeof(double));
*U = malloc(sizeof(double)*nSNP*K);
for (i=0; i<K; i++){
for (j=0; j<nSNP; j++){
*(*(U) + i*nSNP + j) = rand_normal(0, 1);
}
}
*SNPSd = malloc(sizeof(double)*nSNP);
*V = calloc(K*nIND, sizeof(double));
*Cov = calloc(nIND*nIND, sizeof(double));
for (i=0; i<K*nIND; i++) *(*(V) + i) = rand_normal_r();
*miss = calloc(nSNP, sizeof(double));
*mAF = calloc(nSNP, sizeof(double));
return 0;
}
