Код вызвал срабатывание точки останова

// AVal - массив анализируемых данных, Nvl - длина массива должна быть кратна степени 2.
// FTvl - массив полученных значений, Nft - длина массива должна быть равна Nvl.
 
const double TwoPi = 6.283185307179586;
#include <cmath> 
#include <iostream>
#include <iostream>
#include <fstream>
using namespace std;
 
void FFT(double *AVal, double *FTvl, int Nvl, int Nft) {
    int i, j, n, m, Mmax, Istp;
    double Tmpr, Tmpi, Wtmp, Theta;
    double Wpr, Wpi, Wr, Wi;
    double *Tmvl;
 
    n = Nvl * 2; Tmvl = new double[n];
 
    for (i = 0; i < Nvl; i++) {
        j = i * 2; Tmvl[j] = 0; Tmvl[j + 1] = AVal[i];
    }
 
    i = 1; j = 1;
    while (i < n) {
        if (j > i) {
            Tmpr = Tmvl[i]; Tmvl[i] = Tmvl[j]; Tmvl[j] = Tmpr;
            Tmpr = Tmvl[i + 1]; Tmvl[i + 1] = Tmvl[j + 1]; Tmvl[j + 1] = Tmpr;
        }
        i = i + 2; m = Nvl;
        while ((m >= 2) && (j > m)) {
            j = j - m; m = m >> 1;
        }
        j = j + m;
    }
 
    Mmax = 2;
    while (n > Mmax) {
        Theta = -TwoPi / Mmax; Wpi = sin(Theta);
        Wtmp = sin(Theta / 2); Wpr = Wtmp * Wtmp * 2;
        Istp = Mmax * 2; Wr = 1; Wi = 0; m = 1;
 
        while (m < Mmax) {
            i = m; m = m + 2; Tmpr = Wr; Tmpi = Wi;
            Wr = Wr - Tmpr * Wpr - Tmpi * Wpi;
            Wi = Wi + Tmpr * Wpi - Tmpi * Wpr;
 
            while (i < n) {
                j = i + Mmax;
                Tmpr = Wr * Tmvl[j] - Wi * Tmvl[j - 1];
                Tmpi = Wi * Tmvl[j] + Wr * Tmvl[j - 1];
 
                Tmvl[j] = Tmvl[i] - Tmpr; Tmvl[j - 1] = Tmvl[i - 1] - Tmpi;
                Tmvl[i] = Tmvl[i] + Tmpr; Tmvl[i - 1] = Tmvl[i - 1] + Tmpi;
                i = i + Istp;
            }
        }
 
        Mmax = Istp;
    }
 
    for (i = 0; i < Nft; i++) {
        j = i * 2; FTvl[Nft - i - 1] = sqrt((Tmvl[j])*(Tmvl[j]) + (Tmvl[j + 1])*(Tmvl[j + 1]));
    }
 
    delete[] Tmvl;
}
int main()
{
    int Nvl;
    ofstream f;
    cin >> Nvl;
    double *AVal;
    AVal = new double[Nvl];
    for (int i = 0; i < Nvl; ++i) {
        int elem;
        cin >> elem;
        AVal[i] = elem;
    }
    int Nft = Nvl;
    double *FTvl;
    FTvl = new double[Nft];
    FFT(AVal, FTvl, Nvl, Nft);
    f.open("FF.txt");
    for (int i = 0; i < Nft; ++i)
        f << FTvl[i] << " ";
    f.close();
    cout << endl;
}