- int L = cwt.nrow(),
- D = cwt.ncol();
- NumericMatrix fcwt(L, D); //fill with 0... TODO: back to SEXP C-style?
- double *cwt_c = cwt.begin(),
- *fcwt_c = fcwt.begin();
+ int L = INTEGER(Rf_getAttrib(cwt_, R_DimSymbol))[0],
+ D = INTEGER(Rf_getAttrib(cwt_, R_DimSymbol))[1];
+ double *cwt = REAL(cwt_);
+ SEXP fcwt_;
+ PROTECT(fcwt_ = Rf_allocMatrix(REALSXP, L, D));
+ double* fcwt = REAL(fcwt_); //(double*)malloc(L*D*sizeof(double));