//'
//' @return The filtered CWT, in a matrix of same size (LxD)
// [[Rcpp::export]]
-RcppExport SEXP epclustFilter(SEXP cwt_)
+RcppExport SEXP filterMA(SEXP cwt_)
{
// NOTE: C-style for better efficiency (this must be as fast as possible)
int L = INTEGER(Rf_getAttrib(cwt_, R_DimSymbol))[0],
double* fcwt = REAL(fcwt_); //pointer to the encapsulated vector
// NOTE: unused loop parameter: shifting at the end of the loop is more efficient
- for (int col=D-1; col>=0; col++)
+ for (int col=D-1; col>=0; col--)
{
double v1 = cwt[0]; //first value
double ms = v1 + cwt[1] + cwt[2]; //moving sum at second value