#include <iostream>
#include "tbb/tbb.h"

double IntegralSeq(const double a1, const double b1,
                const double a2, const double b2,
                const double h)
{
    int n1, n2;
    // локальная переменная
    // для подсчета интеграла
    double sum = 0.0f;

    // координата точки сетки по оси x
    // координата точки сетки по оси y
    double x, y;
    // количество точек сетки интегрирования
    // n1 - по координате x
    // n2 - по координате y
    n1 = (int)(b1 - a1) / h;
    n2 = (int)(b2 - a2) / h;

    for(int i = 0; i < n1; ++i)
    {
        for(int j = 0; j < n2; ++j)
        {
            //вычисление координат точки
            x = a1 + i * h + h / 2;
            y = a2 + j * h + h / 2;
            //вычисление интеграла
            sum += ((exp(sin(x * M_PI) * cos(y * M_PI)) + 1) /
                    ((b1 - a1) * (b2 - a2))) * h * h;
        }
    }
    return sum;
}

double IntegralTBB(const double a1, const double b1,
            const double a2, const double b2,
            const double h)
{
    int n1, n2;

    n1 = (int)(b1 - a1) / h;
    n2 = (int)(b2 - a2) / h;

    return  tbb::parallel_reduce(tbb::blocked_range2d<int>(0,n1,0,n2),0.f,[&](const tbb::blocked_range2d<int> & r, float sum)
    {
        for(int i = r.rows().begin(); i != r.rows().end(); ++i) {
            for(int j = r.cols().begin(); j != r.cols().end() ; ++j)
            {
                double x = a1 + i * h + h / 2;
                double y = a2 + j * h + h / 2;
                sum += (exp(sin(x * M_PI) * cos(y * M_PI)) + 1) /
                        ((b1 - a1)*(b2 - a2) * h * h);

            }
        }
        return sum;
    },std::plus<int>());
}

int main()
{
    tbb::tick_count t0 = tbb::tick_count::now();

    double res = IntegralSeq(0.0f, 6.0f, 0.0, 6.0, 0.001);

    tbb::tick_count t1 = tbb::tick_count::now();

    std:: cout  << " result seq:" << res << " time:" << (t1-t0).seconds();
    t0 = tbb::tick_count::now();
    res = IntegralTBB(0.0f, 6.0f, 0.0, 6.0, 0.001);

    t1 = tbb::tick_count::now();

    std:: cout  << " result tbb:" << res << " time:" << (t1-t0).seconds();
    return 0;
}