Method process

src/functionals/functionalModulation.cpp:487–562

}

long cFunctionalModulation::process(FLOAT_DMEM *in,
    FLOAT_DMEM *inSorted, FLOAT_DMEM *out, long Nin, long Nout)
{
  if (mappedSpec_ == NULL ||
      (stftWinSizeFrames_ == 0 && mappedSpec_->getNin() != Nin)) {
    // we need to init
    FLOAT_DMEM T = (FLOAT_DMEM)getInputPeriod();
    if (T == 0.0) {
      SMILE_IERR(1, "Cannot compute modulation spectrum when input period is unknown (asynchronous input level!). T = 0.0");
      return 0;
    }
    if (stftWinSizeFrames_ == 0 && T > 0) {
      stftWinSizeFrames_ = stftWinSizeSec_ / T;
      stftWinStepFrames_ = stftWinStepSec_ / T;
    }
    long N = Nin;
    if (stftWinSizeFrames_ > 0) {
      N = stftWinSizeFrames_;
      //fprintf(stderr, "XXX stftWinSize: %ld\n", stftWinSizeFrames_);
    }
    //fprintf(stderr, "XXX allocate with N = %ld\n", N);
    mappedSpec_ = new cSmileUtilMappedMagnitudeSpectrum(N, modSpecNumBins_,
        winFuncId_, modSpecMinFreq_, modSpecMaxFreq_, T);
  }
  if (avgModSpec_ == NULL) {
    avgModSpec_ = (FLOAT_DMEM*)malloc(sizeof(FLOAT_DMEM) * modSpecNumBins_);
  }
  if (removeNonZeroMean_) {
    if (inNormN_ != Nin) {
      if (inNorm_ != NULL) {
        free(inNorm_);
      }
      inNorm_ = NULL;
    }
    if (inNorm_ == NULL) {
      inNormN_ = Nin;
      inNorm_ = (FLOAT_DMEM*)malloc(sizeof(FLOAT_DMEM) * Nin);
    }
    FLOAT_DMEM mean = 0.0;
    long nMean = 0;
    for (long i = 0; i < Nin; i++) {
      if (in[i] != 0.0) {
        mean += in[i];
        nMean++;
      }
    }
    if (nMean > 0) {
      mean /= (FLOAT_DMEM)nMean;
    }
    for (long i = 0; i < Nin; i++) {
      if (in[i] != 0.0) {
        inNorm_[i] = in[i] - mean;
      } else {
        inNorm_[i] = 0.0;
      }
    }
    in = inNorm_;
  }