- PROTECT(A1 = allocArray(REALSXP, dimA));
- PROTECT(A2 = allocArray(REALSXP, dimA));
- PROTECT(rho = allocArray(REALSXP, dimRho);
+ PROTECT(A1 = allocArray(INTSXP, dimA));
+ PROTECT(A2 = allocArray(INTSXP, dimA));
+ PROTECT(dimRho = allocVector(INTSXP, 4));
+ int* pDimRho = INTEGER(dimRho);
+ pDimRho[0] = m; pDimRho[1] = m; pDimRho[2] = k; pDimRho[3] = L;
+ PROTECT(rho = allocArray(REALSXP, dimRho));