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