PROTECT(dimPhi = allocVector(INTSXP, 3));
int* pDimPhi = INTEGER(dimPhi);
pDimPhi[0] = p; pDimPhi[1] = m; pDimPhi[2] = k;
PROTECT(phi = allocArray(REALSXP, dimPhi));
PROTECT(LLF = allocVector(REALSXP, 1));
PROTECT(dimPhi = allocVector(INTSXP, 3));
int* pDimPhi = INTEGER(dimPhi);
pDimPhi[0] = p; pDimPhi[1] = m; pDimPhi[2] = k;
PROTECT(phi = allocArray(REALSXP, dimPhi));
PROTECT(LLF = allocVector(REALSXP, 1));