TPC: fix initialization of FFTW arrays

Author: ehellbar

Detectors/TPC/calibration/include/TPCCalibration/IDCFourierTransform.h

@@ -43,19 +43,13 @@ class IDCFourierTransform : public IDCFourierTransformBase<Type>
   /// \param timeFrames number of time frames which will be stored
   /// \param nFourierCoefficientsStore number of courier coefficients (real+imag) which will be stored (the maximum can be 'rangeIDC + 2', should be an even number when using naive FT). If less than maximum is setn the inverse fourier transform will not work.
   template <bool IsEnabled = true, typename std::enable_if<(IsEnabled && (std::is_same<Type, IDCFourierTransformBaseAggregator>::value)), int>::type = 0>
-  IDCFourierTransform(const unsigned int rangeIDC = 200, const unsigned int nFourierCoefficientsStore = 200 + 2) : IDCFourierTransformAggregator(rangeIDC), mFourierCoefficients{1, nFourierCoefficientsStore}, mVal1DIDCs(sNThreads), mCoefficients(sNThreads)
-  {
-    initFFTW3Members();
-  };
+  IDCFourierTransform(const unsigned int rangeIDC = 200, const unsigned int nFourierCoefficientsStore = 200 + 2) : IDCFourierTransformAggregator(rangeIDC), mFourierCoefficients{1, nFourierCoefficientsStore}, mVal1DIDCs(sNThreads), mCoefficients(sNThreads) {}
 
   /// constructor for  EPN type
   /// \param rangeIDC number of IDCs for each interval which will be used to calculate the fourier coefficients
   /// \param nFourierCoefficientsStore number of courier coefficients (real+imag) which will be stored (the maximum can be 'rangeIDC + 2', should be an even number when using naive FT). If less than maximum is setn the inverse fourier transform will not work.
   template <bool IsEnabled = true, typename std::enable_if<(IsEnabled && (std::is_same<Type, IDCFourierTransformBaseEPN>::value)), int>::type = 0>
-  IDCFourierTransform(const unsigned int rangeIDC = 200, const unsigned int nFourierCoefficientsStore = 200 + 2) : IDCFourierTransformEPN(rangeIDC), mFourierCoefficients{1, nFourierCoefficientsStore}, mVal1DIDCs(sNThreads), mCoefficients(sNThreads)
-  {
-    initFFTW3Members();
-  };
+  IDCFourierTransform(const unsigned int rangeIDC = 200, const unsigned int nFourierCoefficientsStore = 200 + 2) : IDCFourierTransformEPN(rangeIDC), mFourierCoefficients{1, nFourierCoefficientsStore}, mVal1DIDCs(sNThreads), mCoefficients(sNThreads) {}
 
   // Destructor
   ~IDCFourierTransform();

Detectors/TPC/calibration/src/IDCFourierTransform.cxx

@@ -35,6 +35,8 @@ o2::tpc::IDCFourierTransform<Type>::~IDCFourierTransform()
 template <class Type>
 void o2::tpc::IDCFourierTransform<Type>::initFFTW3Members()
 {
+  mVal1DIDCs.resize(sNThreads);
+  mCoefficients.resize(sNThreads);
   for (int thread = 0; thread < sNThreads; ++thread) {
     mVal1DIDCs[thread] = fftwf_alloc_real(this->mRangeIDC);
     mCoefficients[thread] = fftwf_alloc_complex(getNMaxCoefficients());

Detectors/TPC/workflow/include/TPCWorkflow/TPCFourierTransformAggregatorSpec.h

@@ -66,6 +66,7 @@ class TPCFourierTransformAggregatorSpec : public o2::framework::Task
     mEnableFFTCCDB = ic.options().get<bool>("enable-fft-CCDB");
     int nthreads = ic.options().get<int>("nthreads");
     TPCFourierTransformAggregatorSpec::IDCFType::setNThreads(nthreads);
+    initFFTW3Members();
     resizeBuffer(mInputLanes);
   }