PROTECT(A2 = allocArray(REALSXP, dimA));
PROTECT(rho = allocArray(REALSXP, dimRho);
PROTECT(pi = allocMatrix(REALSXP, k, L));
PROTECT(A2 = allocArray(REALSXP, dimA));
PROTECT(rho = allocArray(REALSXP, dimRho);
PROTECT(pi = allocMatrix(REALSXP, k, L));