Обратное быстрое преобразование Фурье

//f-frame, s-frame size, d - direction
double** FFT(double** f, int s, bool d) {
    if (s == 1) return f;
    int sH = s / 2;
 
    double** fOdd = new double*[sH];
    double** fEven = new double*[sH];
    for (int i = 0; i < sH; i++) {
        int j = 2 * i;
        fOdd[i] = f[j];
        fEven[i] = f[j + 1];
    }
 
    double** sOdd = FFT2(fOdd, sH, d);
    double** sEven = FFT2(fEven, sH, d);
 
    double**spectr = new double*[s];
 
    double arg = d ? -DoublePI / s : DoublePI / s;
    double*oBase = new double[2]{ cos(arg),sin(arg) };
    double*o = new double[2]{ 1,0 };
 
 //Mul-произведение, Sum-сумма, Dif-разница
    for (int i = 0; i < sH; i++) {
        double* sO1 = Mul(o, sOdd[i]);
 
        spectr[i] = Sum(sEven[i], sO1);
        spectr[i + sH] = Dif(sEven[i], sO1);
 
        o = Mul(o, oBase);
    }
 
    return spectr;
}