38513-LW_Integral_ppt

Параллельная версия. Распараллеливание

оптимизированного алгоритма

 Далее распараллельте основные вычислительные циклы и освободите выделенную память

//вычисление интеграла sum = 0.0;

#pragma omp parallel for private (x, y, j) / reduction(+: sum)

for(i = 0; i < npi; i++)

{

for(j = 0; j < npi; j++)

{

//вычисление интеграла

sum += (exp(sinx[i] * cosy[j]) /

((b1 - a1) * (b2 - a2))) * h * h;

}

}

*res = sum * (n1 / npi) * (n2 / npi) + 1; delete [] sinx;

delete [] cosy;

}

Параллельная версия. Распараллеливание

оптимизированного алгоритма

Пересоберите получившийся код (команда Build→Rebuild Solution) и запустите его на выполнение.

Параллельная версия. Распараллеливание

оптимизированного алгоритма

Сравнение времени численного интегрирования последовательной и параллельной реализаций

Задания для самостоятельной работы…

Реализуйте алгоритмы, использующие буфер, на случай, когда размер буфера не кратен количеству узлов сетки интегрирования.

Предложите дальнейшие алгоритмические методы сокращения времени вычисления последней версии реализации алгоритма численного интегрирования.

Рассмотреть возможность использования Intel MKL для векторного вычисления математических функций (функции VML).

Рассмотреть другие квадратурные формулы.

Использованные источники информации

Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. – Бином. Лаборатория знаний, 2008. – 640c.

Гергель В.П. Теория и практика параллельных вычислений. – Бином. Лаборатория знаний, 2007. – 424c.

Авторский коллектив

Козинов Евгений Александрович ассистент кафедры

Математического обеспечения ЭВМ факультета ВМК ННГУ.

Evgeniy.Kozinov@gmail.com

Сысоев Александр Владимирович,

ассистент кафедры Математического обеспечения ЭВМ факультета ВМК ННГУ. sysoyev@vmk.unn.ru