Method readNewData

src/dspcore/fullinputMean.cpp:130–194

Reads new data, computes stats in means_ and n_means_.

// Reads new data, computes stats in means_ and n_means_.
eTickResult cFullinputMean::readNewData()
{
  cVector *vec = reader_->getNextFrame();
  if (vec == NULL) 
    return TICK_SOURCE_NOT_AVAIL;
  if (means_ == NULL) {
    means_ = new cVector(vec->N);
    n_means_arr_ = new long[vec->N];
    if (mean_type_ == MEANTYPE_RQMEAN) {
      for (int i = 0; i < vec->N; i++) {
        means_->data[i] = vec->data[i] * vec->data[i];
      }
    } else if (mean_type_ == MEANTYPE_ABSMEAN) {
      for (int i = 0; i < vec->N; i++) {
        means_->data[i] = (FLOAT_DMEM)fabs(vec->data[i]);
      }
    } else {
      if (mean_type_ == MEANTYPE_AMEAN && exclude_zeros_) {
        for (int i = 0; i < vec->N; i++) {
          if (vec->data[i] != 0.0) {
            means_->data[i] = vec->data[i];
            n_means_arr_[i] = 1;
          }
        }
      } else {
        for (int i = 0; i < vec->N; i++) {
          means_->data[i] = vec->data[i];
        }
      }
    }
    n_means_ = 1;
  } else {
    if (mean_type_ == MEANTYPE_ENORM) {
      for (int i = 0; i < vec->N; i++) {
        if (vec->data[i] > means_->data[i])
          means_->data[i] = vec->data[i];
      }
    } else if (mean_type_ == MEANTYPE_AMEAN) {
      if (exclude_zeros_) {
        for (int i = 0; i < vec->N; i++) {
          if (vec->data[i] != 0.0) {
            means_->data[i] += vec->data[i];
            n_means_arr_[i]++;
          }
        }
      } else {
        for (int i = 0; i < vec->N; i++) {
          means_->data[i] += vec->data[i];
        }
      }
      n_means_++;
    } else if (mean_type_ == MEANTYPE_RQMEAN) {
      for (int i = 0; i < vec->N; i++) {
        means_->data[i] += vec->data[i] * vec->data[i];
      }
      n_means_++;
    } else if (mean_type_ == MEANTYPE_ABSMEAN) {
      for (int i = 0; i < vec->N; i++) {