Method processVector

src/lldcore/spectral.cpp:586–1555

a derived class should override this method, in order to implement the actual processing

// a derived class should override this method, in order to implement the actual processing
int cSpectral::processVector(const FLOAT_DMEM *src, FLOAT_DMEM *dst, long Nsrc, long Ndst, int idxi) // idxi=input field index
{
  long n = 0;
  long i = 0;
  long j = 0;

  double N = (double) ((Nsrc-1)*2);  // assumes FFT magnitude input array!!
  double F0 = 1.0/fsSec;

  if (frq == NULL) { // input block's frequency axis
    // TODO: use frqSet flag to avoid doing the code below on every tick if info==NULL
    const FrameMetaInfo * fmeta = reader_->getFrameMetaInfo();
    if ((fmeta != NULL)&&(idxi < fmeta->N)) {
      // TODO: check input block's type (in fmeta)
      // fmeta->field[idxi].dataType & SPECTRAL == 1 ?
      nScale = fmeta->field[idxi].infoSize / sizeof(double);
      frq = (double *)(fmeta->field[idxi].info);
      if (nScale > 0) {
        // old frequency axis computation mode:
        //nScale=0;
        if (nScale != Nsrc) {
          SMILE_IWRN(2,"number of frequency axis points (from info struct) [%i] is not equal to Nsrc [%i] ! Field index: %i (check the processArrayFields option).",nScale,Nsrc,idxi);
          nScale = MIN(nScale,Nsrc);
        }
      }
    }
  }

  if (frqScale == -1) {
      cVectorMeta *mdata = writer_->getLevelMetaDataPtr();
      if (mdata != NULL && mdata->ID == 1001 /* SCALED_SPEC */) {
        frqScale = (int)mdata->fData[6];
        frqScaleParam = (double)mdata->fData[7];
      } else {
        frqScale = SPECTSCALE_LINEAR;
        frqScaleParam = 0.0;
      }
  }

  if (specRangeLowerBin == -1) {
    if (specRangeLower == specRangeUpper && specRangeUpper == 0) {
      specRangeLowerBin = 1;
      specRangeUpperBin = Nsrc - 1;
    } else {
      for (i = 0; i < Nsrc; i++) {
        if ((double)specRangeLower >= frq[i]) {
          specRangeLowerBin = i;
        }
        if ((double)specRangeUpper > frq[i]) {
          specRangeUpperBin = i;
        }
      }
      if (specRangeUpperBin == -1 || specRangeUpperBin >= Nsrc) {
        specRangeUpperBin = Nsrc - 1;
      }
      if (specRangeLowerBin < 0) {
        specRangeLowerBin = 0;
      }