Распараллеливание вычисления интеграла используя редукции (OpenMP) - OpenMP - Ответ 12351569

double integralReduction(double a, double b, double steps)
{
    int n = 3;
    double result = 0;
    double function = 0;
    double dx = (b - a) / steps;
    double startTime = 0, endTime = 0, time;
    omp_set_num_threads(n);
    double i;
#pragma omp parallel
    {
        startTime = omp_get_wtime();
#pragma omp for reduction (+:function)
            for (i = a; i <b; i += dx)
            {
                function += log(i) - (3 * sin(3 * i));
            }
            endTime = omp_get_wtime();
            printf("Time work: %f\nNumber potok: %d\n", endTime - startTime, omp_get_thread_num());
            result = function * dx;
    }
    return result;
}

double integralReduction(double a, double b, double steps)
{
    int n = 3;
    double result = 0;
    double function = 0;
    double dx = (b - a) / steps;
    double x = a;
    double startTime = 0, endTime = 0;
    omp_set_num_threads(n);
    int i;
#pragma omp parallel
    {
        startTime = omp_get_wtime();
#pragma omp for reduction (+:result,x) private (i)
            for (i = 0; i <steps; i++)
            {
                x += dx;
                function += log(x) - (3 * sin(3 * x));
                result += function * dx;
            }
            endTime = omp_get_wtime();
            printf("Time work: %f\nNumber potok: %d\n", endTime - startTime, omp_get_thread_num());
    }
    return result;
}