Method dataProcessorCustomFinalise

src/dsp/specScale.cpp:248–321

}

int cSpecScale::dataProcessorCustomFinalise()
{
  int i;
  int ret = cVectorProcessor::dataProcessorCustomFinalise();
  if (!ret) return ret;

  // check maxF < sampleFreq.
  double samplF = deltaF * (double)(nMag - 1); // sampling frequency
  if ((maxF <= minF)||(maxF > samplF)) {
    maxF = samplF; 
  }
  
  // transform from source to target scale:

  fmin_t = smileDsp_specScaleTransfFwd(minF,scale,param); 
  fmax_t = smileDsp_specScaleTransfFwd(maxF,scale,param);

  // target delta f
  deltaF_t = (fmax_t - fmin_t) / (nPointsTarget - 1);

  // calculate the target frequencies of the linear scale fft input points
  if (f_t == NULL) f_t = (double*)malloc(sizeof(double)*nMag);

  if (scale == SPECTSCALE_LOG) {
    for (i=1; i < nMag; i++) {
      f_t[i] = smileDsp_specScaleTransfFwd( (double)i * (double)deltaF , scale, param );
    }
    f_t[0] = 2.0 * f_t[1] - f_t[2]; // heuristic for the 0th frequency (only valid for log2 ??)
  } else { // generic transform:
    for (i=0; i < nMag; i++) {
      f_t[i] = smileDsp_specScaleTransfFwd( (double)i * (double)deltaF , scale, param );
    }
  }

  splineCache = std::unique_ptr<sSmileMath_splineCache>(new sSmileMath_splineCache());
  smileMath_cspline_init(f_t, nMag, splineCache.get());

  double *x = (double*)malloc(sizeof(double)*nPointsTarget);  
  for (i=0; i < nPointsTarget; i++) {
    x[i] = fmin_t + (double)i * deltaF_t;
  }
  splintCache = std::unique_ptr<sSmileMath_splintCache>(new sSmileMath_splintCache());
  if (!smileMath_csplint_init(f_t, nMag, x, nPointsTarget, splintCache.get())) {
    SMILE_IERR(1,"smileMath_csplint_init failed. Output of this component will be invalid!");
    splintCache = NULL;
  }
  free(x);
  
  double nOctaves = log(maxF / minF)/log(2.0);
  double nPointsPerOctave = nPointsTarget / nOctaves; // this is valid for log(2.0) scale only...
  if (auditoryWeighting) {
    /* auditory weighting function (octave scale only...)*/
    double atan_s = nPointsPerOctave * smileMath_log2(65.0 / 50.0) - 1.0;
    audw = (double*)malloc(sizeof(double)*nPointsTarget);
    for (i=0; i < nPointsTarget; i++) {
      audw[i] = 0.5 + atan (3.0 * (i + 1 - atan_s) / nPointsPerOctave) / M_PI;
    }
  }