Faster pitch init

include/algorithms/public/CepstrumF0.hpp

@@ -32,8 +32,11 @@ class CepstrumF0
 
   void init(index size)
   {
-    // avoid allocation of maxSize^2 at constructor
-    mDCT = DCT(size, size);
+    // avoid allocation of maxSize^2 at constructor, but don't reallocate if not needed
+
+    if (mDCT.maxOutputSize() < size || mDCT.maxInputSize() < size)
+      mDCT = DCT(size, size);
+
     mDCT.init(size, size);
 
     mCepstrum = mCepstrumStorage.segment(0, size);

include/algorithms/public/DCT.hpp

@@ -28,6 +28,8 @@ class DCT
   using MatrixXd = Eigen::MatrixXd;
 
   DCT(index maxInputSize, index maxOutputSize)
+  : mMaxInputSize(maxInputSize)
+  , mMaxOutputSize(maxOutputSize)
   {
     mTableStorage = MatrixXd::Zero(maxOutputSize, maxInputSize);
   }
@@ -36,6 +38,11 @@ class DCT
   {
     using namespace std;
     assert(inputSize >= outputSize);
+    assert(inputSize <= mMaxInputSize);
+    assert(outputSize <= mMaxOutputSize);
+    // Exit if already initialised
+    if (mInputSize == inputSize && mOutputSize == outputSize)
+      return;
     mInputSize = inputSize;
     mOutputSize = outputSize;
     mTable = mTableStorage.block(0, 0, mOutputSize, mInputSize);
@@ -61,8 +68,14 @@ class DCT
   {
     output = (mTable * input.matrix()).array();
   }
-  index    mInputSize{40};
-  index    mOutputSize{13};
+
+  index maxInputSize() const { return mMaxInputSize; }
+  index maxOutputSize() const { return mMaxOutputSize; }
+
+  index    mInputSize{0};
+  index    mOutputSize{0};
+  index    mMaxInputSize{0};
+  index    mMaxOutputSize{0};
   MatrixXd mTable;
   MatrixXd mTableStorage;
 };

include/clients/rt/PitchClient.hpp

@@ -150,7 +150,7 @@ class PitchClient : public FluidBaseClient, public AudioIn, public ControlOut
     return { get<kFFT>().winSize(), get<kFFT>().hopSize() }; 
   }
 
-  void  reset()
+  void reset()
   {
     mSTFTBufferedProcess.reset();
     cepstrumF0.init(get<kFFT>().frameSize());