[AMP] Распараллеливание вложенных циклов

for (int k = 0; k <= Nfi; k++)
    {
        TmpPsi[k] = 0.0;
        for (int n = 0; n <= Nx; n++)
        {
            complex<double> S = 0.0;
            for (int nu = 0; nu <= Nx + 2; nu++)
                for (int j = -Nfi; j <= Nfi; j++)
                    S += B2_Func(n, -k, nu, j)*Vfi_(nu, j);
            for (int nu = 0; nu <= Nx + 2; nu++)
                for (int m = 0; m <= Nx + 2; m++)
                {
                    double AA0 = A0_Func(n, nu, m);
                    double AA1 = AA1_Func(n, nu, m);
                    double AA2 = AA2_Func(n, nu, m);
                    double AA3 = AA3_Func(n, nu, m);
                    double AA4 = AA4_Func(n, nu, m);
                    double AA5 = AA5_Func(n, nu, m);
                    double AA6 = AA6_Func(n, nu, m);
                    double AA7 = AA7_Func(n, nu, m);
 
                    for (int j = -Nfi; j <= Nfi; j++)
                        for (int l = -Nfi; l <= Nfi; l++)
                            {
                                complex<double> B3 = AA1*Jz_(j + l - k);
                                if ((j + l) == k)
                                    B3 += complex<double>(0, M_PI*(2 * j*AA1 - l*AA0));
                                if (abs(j + l - k) <= 2)
                                    B3 += complex<double>(0, beta / 4)*(AA4*l + AA5*(l - k) + AA2*(k - j) - (k*AA3) / 2)*J_(j + l - k, 1);
                                if (abs(j + l - k) <= 4)
                                    B3 += complex<double>(0, beta*beta / 16)*(AA6*(j - l - k) + (AA7*(j + l + k)) / 4)*J_(j + l - k, 2);
                                S += B3*Vfi_(nu, j)*Vfi_(m, l);
                }
            }
            TmpPsi[k](n) = S / M_PI;
        }
        if (k == 0)
            TmpPsi[k](Nx + 1) = -wt;
        if (k == 1)
            TmpPsi[k](Nx + 2) = (beta / 2)*conj(yt);
        Blas_Mat_Vec_Mult(A_Inv, TmpPsi[k], Psi[k]);
    }