ПАРАЛЛЕЛЬНЫЙ АЛГОРИТМ (2,1)-МЕТОДА ПЕРЕМЕННОГО ШАГА > Полезные советы
Тысяча полезных мелочей    

ПАРАЛЛЕЛЬНЫЙ АЛГОРИТМ (2,1)-МЕТОДА ПЕРЕМЕННОГО ШАГА

ПАРАЛЛЕЛЬНЫЙ АЛГОРИТМ (2,1)-МЕТОДА ПЕРЕМЕННОГО ШАГА

Ващенко Г.В. Статья в формате PDF 505 KB

Предложен параллельный алгоритм переменного шага на основе (2,1)-метода. В предлагаемом параллельном алгоритме изменение величины шага построено на основе контроля точности численной схемы.

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

Рассматривается задача Коши для автономной системы обыкновенных дифференциаль ных уравнений первого порядка

 (1)

где y:[t0, tk] → RN, f:[t0, tk]×RN → RN, [t0, tk] - отрезок интегрирования. Для численного решения (1) применим схему (2,1)-метода

 (2)

 

где коэффициенты a, p1 и p2 определяют свойства точности и устойчивости схемы (2), h - шаг интегрирования, fn′ = ∂f(yn)/∂y - матрица Якоби системы (1). Будем считать, что (1) имеет единственное решение. Пусть известны условия для контроля точности вычислений, именно p1 + p2 = 1 и ap1 + 2ap2 = 0,5, . Изменение величины шага основано на оценке локальной ошибки δn. Учитывая соотношение, , новый шаг hnew определяем по формуле hnew = qh, где значение q находится из уравнения q2||δn|| = ε. Если q < 1, то осуществляем повторное вычисление решения с шагом h = hnew. При q > 1 выполняем следующий шаг интегрирования с шагом hnew. Введем функции Par_LU_Decompos(), Par_LU_Solution(), реализующие декомпозицию матрицы Dn и нахождение векторов , . Для контроля точности численной схемы (2) введем функцию accur_control (), для выполнения которой назначим процессор proc(1). Параллельный алгоритм вычисления приближенного решения y(n+1) переменного шага формулируем следующим образом.

Алгоритм. Пусть для численного решения системы (1) используется (2.1)-метод с контролем точности, и известно решение y(n) в точке tn с шагом hn. Тогда для получения значения y(n+1) в точке tn+1 справедлив параллельный алгоритм, в котором на каждом процессоре proc(j) формируется своя j-я часть вектора решения.

Шаг 1. В каждом proc(j), 1 ≤ j ≤ p; (j-1) s + 1 ≤ sj ≤ j⋅s: выполнить recv(, h; 1,..., p), вычислить  и матрицу Якоби Jj, 1 ≤ j ≤ p.

Шаг 2. Сформировать матрицу .

Шаг 3. Разложить матрицу Dn, Dn = Par_LU_Decompos().

Шаг 4. Вычислить ,

.

Шаг 5. Вычислить ,

Шаг 6. В каждом proc(j), 1 ≤ j ≤ p;

(j-1)⋅s +1 ≤ sj ≤ j⋅s:

определить

,

и выполнить .

Шаг 7. В proc(1): выполнить accur_control () и, если необходимо, вывести вектор y(n+1).

Шаг 8. В каждом proc(j), 1 ≤ j » p;

(j-1)⋅s +1 ≤ sj ≤ j⋅s:

вычислить

и выполнить .

Шаг 9. Выполнить следующий шаг интегрирования.

Как показывают теоретические и пpaктические расчеты, выполняемые на кластере ИВМ СО РАН [3] показывают, что основные вычислительные затраты связаны с реализацией
LU-факторизации и последующем решении систем для определения шаговых коэффициентов.

Работа выполнена при финансовой поддержке РФФИ проект №11-01-00106.

Список литературы

  1. Новиков Е.А. Явные методы для жестких систем. - Новосибирск: Наука, 1997.
  2. Хайрер Э., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-алгебраические задачи. - М.: Мир, 1999.
  3. Исаев С.В., Малышев А.В., Шайдуров В.В. Развитие Красноярского центра параллельных вычислений // Вычислительные технологии. - 2006. - №11. - С. 28-33.


РОЛЬ ОХОТОВЕДЕНИЯ В СОХРАНЕНИИ БИОРАЗНООБРАЗИЯ

РОЛЬ ОХОТОВЕДЕНИЯ В СОХРАНЕНИИ БИОРАЗНООБРАЗИЯ Статья в формате PDF 125 KB...

03 07 2026 11:13:34

ДИСТАНЦИОННОЕ УПРАВЛЕНИЕ ЗДОРОВЬЕМ ЧЕЛОВЕКА С ПОМОЩЬЮ КВАНТОВО-ВОЛНОВЫХ ФИЗИОЛОГИЧЕСКИХ ТЕХНОЛОГИЙ (КВАНТОВО-ВОЛНОВАЯ ФИЗИОЛОГИЯ)

ДИСТАНЦИОННОЕ УПРАВЛЕНИЕ ЗДОРОВЬЕМ ЧЕЛОВЕКА С ПОМОЩЬЮ КВАНТОВО-ВОЛНОВЫХ ФИЗИОЛОГИЧЕСКИХ ТЕХНОЛОГИЙ (КВАНТОВО-ВОЛНОВАЯ ФИЗИОЛОГИЯ) В работе рассматриваются вопросы дистанционного управления здоровьем человека с помощью квантово-волновых нейроинформационных технологий – электроакустических импульсов, скопированных у адаптированной к гипоксии нервной клетке. Приведены данные, cсвидетельствующие о нормализующем действии моделей нейроинформационных сигналов на концентрацию СО2 в крови. В результате этого просвет кровеносных сосудов расширяется, в клетках восстанавливается режим нормоксии – основного фактора здоровья человека. ...

27 06 2026 22:13:56

РАЗВИТИЕ ЛЕСНОГО КОМПЛЕКСА В&#8239;РОССИИ

РАЗВИТИЕ ЛЕСНОГО КОМПЛЕКСА В&#8239;РОССИИ Статья в формате PDF 314 KB...

25 06 2026 2:26:47

ПРОБЛЕМЫ ИДЕНТИФИКАЦИИ МОДЕЛЕЙ РАСПРЕДЕЛЕНИЯ СЛУЧАЙНЫХ ВЕЛИЧИН С ПРИМЕНЕНИЕМ СОВРЕМЕННОГО ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ

ПРОБЛЕМЫ ИДЕНТИФИКАЦИИ МОДЕЛЕЙ РАСПРЕДЕЛЕНИЯ СЛУЧАЙНЫХ ВЕЛИЧИН С ПРИМЕНЕНИЕМ СОВРЕМЕННОГО ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ Рассмотрены некоторые проблемы идентификации моделей распределения данных, при использовании современного математического аппарата для решения этой задачи. Показано, что использование методов нелинейной оптимизации для идентификации моделей приводит к улучшению результатов идентификации, но одновременно, изменяет формальную постановку задачи. Выделено три группы проблем, связанных с выбором критериев согласия, их критических значений и проверкой адекватности получаемых моделей. Проанализированы возможные подходы к решению этих проблем. ...

22 06 2026 12:31:19

ИСПОЛЬЗОВАНИЕ ТЕХНОЛОГИИ МУЛЬТИМЕДИА В ОБУЧЕНИИ

ИСПОЛЬЗОВАНИЕ ТЕХНОЛОГИИ МУЛЬТИМЕДИА В ОБУЧЕНИИ Статья в формате PDF 246 KB...

21 06 2026 2:53:44

КОЛЕБАНИЯ УПРУГО-ПЛАСТИЧЕСКОЙ СРЕДЫ

КОЛЕБАНИЯ УПРУГО-ПЛАСТИЧЕСКОЙ СРЕДЫ Статья в формате PDF 147 KB...

16 06 2026 4:46:58

Концепт «время» в русской и китайской культурах

Концепт «время» в русской и китайской культурах Статья в формате PDF 319 KB...

11 06 2026 11:11:13

HLA-АССОЦИАЦИИ ПРИ ТЯЖЕЛОМ ТЕЧЕНИИ АТОПИЧЕСКОГО ДЕРМАТИТА И АТОПИЧЕСКОЙ БРОНХИАЛЬНОЙ АСТМЫ У ДЕТЕЙ

HLA-АССОЦИАЦИИ ПРИ ТЯЖЕЛОМ ТЕЧЕНИИ АТОПИЧЕСКОГО ДЕРМАТИТА И АТОПИЧЕСКОЙ БРОНХИАЛЬНОЙ АСТМЫ У ДЕТЕЙ У детей установлено существование ассоциативной связи тяжёлого течения атопического дерматита и атопической бронхиальной астмы с иммуногенетическими параметрами. ...

09 06 2026 18:28:23

ЭКОНОМИЧЕСКАЯ ОЦЕНКА СОРТОВ И ЭЛИТ ЯБЛОНИ В АГРОЛАНДШАФТАХ В РАЗНЫХ ЭКОЛОГИЧЕСКИХ ЗОНАХ СЕВЕРНОГО КАВКАЗА

ЭКОНОМИЧЕСКАЯ ОЦЕНКА СОРТОВ И ЭЛИТ ЯБЛОНИ В АГРОЛАНДШАФТАХ В РАЗНЫХ ЭКОЛОГИЧЕСКИХ ЗОНАХ СЕВЕРНОГО КАВКАЗА При оценке экономической эффективности культивирования сортов яблони в разных экологических условиях Северного Кавказа важными показателями являются продуктивность и качество продукции. Значительный интерес в этом отношении представляют сорта селекции Северо-Кавказского НИИ горного и предгорного садоводства: Фестиваль гор, Долинское, Златогор, Лескенское и другие, которые являются аборигенными и конкурентоспособными с культивируемыми в регионе интродуцированными зарубежными и отечественными сортами яблони. ...

06 06 2026 20:34:31

Опыт применения ГИС технологий и&#8239;ДДЗЗ для оценки влияния антропогенного фактора на растительность отдельных территорий Якутии

Опыт применения ГИС технологий и&#8239;ДДЗЗ для оценки влияния антропогенного фактора на растительность отдельных территорий Якутии В работе приведены результаты применения ГИС технологий в различном масштабе для анализа структуры растительности и влияние антропогенной нагрузки на параметры растительного покрова регионов, в целом, и отдельных сообществ Якутии, в частности. Примененные подходы могут быть использованы в различном масштабе для анализа степени антропогенного пресса территорий и анализа растительности. ...

05 06 2026 9:11:45

ВИКТОР СТЕПАНОВИЧ ДМИТРИЕВ

ВИКТОР СТЕПАНОВИЧ ДМИТРИЕВ Статья в формате PDF 317 KB...

02 06 2026 13:38:40

БИОТЕХНИЧЕСКИЙ ЗАКОН И ЧИСЛЕННОСТЬ НАБЛЮДЕНИЙ

БИОТЕХНИЧЕСКИЙ ЗАКОН И ЧИСЛЕННОСТЬ НАБЛЮДЕНИЙ Статья в формате PDF 390 KB...

27 05 2026 4:26:46

ОПРЕДЕЛЕНИЕ ФАКТОРОВ, ВЛИЯЮЩИХ НА КАЧЕСТВО ИЗМЕРЕНИЙ ПРИБОРА МАЭС

ОПРЕДЕЛЕНИЕ ФАКТОРОВ, ВЛИЯЮЩИХ НА КАЧЕСТВО ИЗМЕРЕНИЙ ПРИБОРА МАЭС Существующие методы атомной эмиссионной спектроскопии для исследования состава металлов и сплавов используются во всех отраслях машиностроения. По мнению авторов, современные методы уже не обеспечивают необходимых точностей измерений. В данной работе авторами проведены исследования влияния внешних факторов на точность измерений прибора атомно-эмиссионной спектроскопии. ...

25 05 2026 20:49:54

Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::