Удк 519. 6 Механизмы управления длиной шага в параллельных алгоритмах решения дифференциальных уравнений

Вид материалаДокументы
Подобный материал:
УДК 519.6

МЕХАНИЗМЫ УПРАВЛЕНИЯ ДЛИНОЙ ШАГА В ПАРАЛЛЕЛЬНЫХ АЛГОРИТМАХ РЕШЕНИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Дмитриева О.А., Куприй Я.А.,

Донецкий национальный технический университет


В данной работе рассматривается применение переменного шага к алгоритмам решения уравнений блочными предиктор – корректорными методы Рунге-Кутты.

Одними из параллельных методов решения дифференциальных уравнений являются блочные методы Рунге-Кутты. В методах этой категории, оси решений разделены на последовательности смежных блоков, каждый из которых содержит фиксированной число точек. На каждом шаге процедуры генерируется решение в каждой из точек внутри блока [1].

Рассмотрим обыкновенное дифференциальное уравнение вида (1):


 (1)


где f – некоторая функция, связывающая независимую переменную t, искомую функцию y(t) и ее производные до m-го порядка включительно. Пусть ut – приближенное решение задачи (1) в момент времени t.

S-этапный блочный метод Рунге-Кутты задается формулами (2) [2]:

(2)


где k – число аппроксимаций решения;

Ai, Aij, B0 и Bj – матрицы коэффициентов, размерностью r x k, r x r, k x k и k x r;

ek – k-й единичный вектор;

– k-мерный блочный вектор, компонентами которого являются приближения значений точного решения y(tn+cjh). При этом, ck = 1, другие значения cj – могут быть положены равными любому действительному значению.

Корректор генерируется с помощью формулы (3) [3]:


(3)


Для контроля шага используется оценка погрешности (4)


, (4)


где m – число итераций.


Использование переменного шага является одним из возможных методов повышения точности численного решения дифференциальных уравнений, но, с другой стороны, увеличивает объем вычислений. В приведенном алгоритме для управления длиной шага осуществлять дополнительные итерации не требуется, поскольку компонента Y(i) является элементом итерационного процесса.


Литература

1. Liming Yang. Parallel Block Predictor-Corrector Methods for Numerical Solution of ODE’s // A thesis of the School of Graduate Studies and Research of the University of Ottawa in partial fulfillment of the requirements of Master of science in Systems Science. Ottawa, Ontario, Canada, 1993.

2. P.J. van der Houwen, B.P.Sommeijer. Block Runge-Kutta methods on parallel computers // CWI Centre for Mathematics and Computer Science, Report NM-R8906, March, 1989.

3. P.J. van der Houwen. Block Runge-Kutta methods// CWI Centre for Mathematics and Computer Science, Report NM-R8913, June, 1989.