Оптимизация кода

while (!arr[0]) {
/* Здесь я всеми силами пытаюсь оптимизировать данный цикл вычислений, он очень тормозит программу */
        #pragma omp parallel for
        for (int i = 0; i < 898; i++)
        {
        chast[i] = (arr[6] * sym[7 * i] + arr[5] * sym[7 * i + 1] + arr[4] * sym[7 * i + 2] + arr[3] * sym[7 * i + 3] + arr[2] * sym[7 * i + 4] + arr[1] * sym[7 * i + 5] + arr[0] * sym[7 * i + 6]) % 31;
        }
 
        double mass[31];
        memset(mass, 0, sizeof(mass));
        double q[31];
        for (int i = 0; i < 898; i++)
        {
                mass[chast[i]]++;
        }
        double khi_2 = 1;
        for (int i = 0; i < 31; i++)
        {
            q[i] = mass[i] / 898;
        }
        for (int i = 0; i < 31; i++)
        {
            khi_2 += q[i] * q[i];
        }
        if (khi_2 < 1.0)
        {
            cout << "khi square: " << khi_2 << '\n';
            for (int i = M - 1; i >= 0; i--)
                printf("%d ", arr[i]);
            printf("\n");
        }
        arr[M-1] ++;
        int j = M-1;
        while (arr[j] == N) {
            arr[j] = 0;
            arr[--j]++;
        }
 
    }

        for (int i = 0; i < 31; i++)
        {
            q[i] = ((double)mass[i]) / 898.0;
            khi_2 += q[i] * q[i];
        }