Method Process

Source/FFTtoAdditive.cpp:86–146

}

void FFTtoAdditive::Process(double time)
{
   PROFILER(FFTtoAdditive);

   IAudioReceiver* target = GetTarget();

   if (target == nullptr || !mEnabled)
      return;

   ComputeSliders(0);
   SyncBuffers();

   float inputPreampSq = mInputPreamp * mInputPreamp;
   float volSq = mVolume * mVolume;

   int bufferSize = GetBuffer()->BufferSize();

   mRollingInputBuffer.WriteChunk(GetBuffer()->GetChannel(0), bufferSize, 0);

   //copy rolling input buffer into working buffer and window it
   mRollingInputBuffer.ReadChunk(mFFTData.mTimeDomain, fftWindowSize, 0, 0);
   Mult(mFFTData.mTimeDomain, mWindower, fftWindowSize);
   Mult(mFFTData.mTimeDomain, inputPreampSq, fftWindowSize);

   mFFT.Forward(mFFTData.mTimeDomain,
                mFFTData.mRealValues,
                mFFTData.mImaginaryValues);

   for (int i = 0; i < fftFreqDomainSize; ++i)
   {
      float real = mFFTData.mRealValues[i];
      float imag = mFFTData.mImaginaryValues[i];

      //cartesian to polar
      float amp = 2. * sqrtf(real * real + imag * imag);
      float phase = atan2(imag, real);

      mFFTData.mRealValues[i] = amp / (fftWindowSize / 2);
      mFFTData.mImaginaryValues[i] = phase;
   }

   float* out = target->GetBuffer()->GetChannel(0);
   for (int i = 0; i < bufferSize; ++i)
   {
      float write = 0;
      for (int j = 1; j < numPartials; ++j)
      {
         float phase = ((mFFTData.mImaginaryValues[j + 1] + i * mPhaseInc[j]) / FTWO_PI) * 512;
         float sample = SinSample(phase) * mFFTData.mRealValues[j + 1] * volSq * .4f;
         write += sample;
      }

      GetVizBuffer()->Write(write, 0);

      out[i] += write;

      time += gInvSampleRateMs;
   }