ПАРАЛЛЕЛЬНЫЕ ЯВНЫЕ ОДНОШАГОВЫЕ МЕТОДЫ ДЛЯ ЧИСЛЕННОГО РЕШЕНИЯ ЖЕСТКИХ СИСТЕМ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ > Полезные советы
Тысяча полезных мелочей    

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

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

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

Предложены параллельные явные одношаговые методы первого, второго порядков, обеспечивающие возможность с минимальными вычислительными затратами интегрировать жесткие системы обыкновенных дифференциальных уравнений. В предлагаемых параллельных алгоритмах изменение величины шага построены на основе контроля точности и устойчивости численной схемы, а в неравенстве для контроля точности применяется оценка локальной ошибки метода.

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

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

y′ = f(y), y(t0) = y0, t0 ≤ t ≤ tk (1)

где y: [t0, tk] →RN, f: [t0, tk]× RN →RN, [t0, tk] -отрезок интегрирования. В предположении существования и единственности решения задачи (1) параллельная схема метода первого порядка с контролем точности для численного решения (1) в вычислительной системе из p процессоров, N > p и s = N/p, если N кратно p, или s = [N/p] + g, в противном случае, записывается в виде [1]

 (2)

где yjs(n) ∈ Comp(j), || δn ||= 0.5h || fn + 1 - fn || ≤ ε, 1 ≤ j ≤ p, (j-1)⋅s + 1 ≤ js ≤ j⋅ s, ||⋅|| - некоторая норма в RN , || δn ||- норма вектора локальной погрешности, fn + 1 и fn - значения правой части системы (1) соответственно в точках t n+1 и tn, ε требуемая точность. Параллельная схема второго порядка для численного решения (1) имеет вид

 (3)

Неравенство для оценки устойчивости h | λmax |≤ D, где | λmax | -наибольшее собственное число якобиана, D - размер области устойчивости (для схемы (3) он равен 2). Выбор величины шага hn для схемы (2) определяется по формуле hn = qhn/1.1, где q = (ε /|| δn ||)1/2, а для схемы (3) по формуле hn = max(hn, qhn)/1.1, где q = (D / hnmax | )1/2 .

Укрупненная схема параллельных алгоритмов предложенных вычислительных схем (2), (3) состоит в следующем. Компоненты yjs(n) распределяются по p процессорам согласно блочной схеме распределения по s компонентов в каждом. Каждая задача Uj выполняется на proc(j), Uj ∈ proc(j). Proc(1) определяет значение шага hn и передает всем proc(j), используя коммуникационную операцию one-to- all. В каждом proc(j) вычисляются yjs(n), т.е. решается задача Uj, вычисляется значение локальной нормы || δn ||j и выполняется операция all-to-all. Для вычисления значений элементов fjs(y(n) ) вектора правой части разpaбатывается отдель ная функция. Таким образом, общая схема параллельного алгоритма сводится к линейной форме и обеспечивается возможность анализа и оценки его эффективности алгоритма.

Алгоритмы реализованы в виде отдельных функций языка С и включены в комплекс программ, предназначенных для численного моделирования процессов, описываемых жесткими системами на многопроцессорных вычислительных системах кластерной архитектуры. Коммуникационные операции реализованы функциями библиотеки MPI.

Расчеты, выполняемые на 99-процессорном кластере ИВМ СО РАН [4] показали, что параллельные схемы (2), (3) применяться в случаях, когда расчеты требуется проводить с невысокой точностью - порядка 1 % и ниже.

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

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


Концепт «удача» в русских и китайских песнях

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

11 06 2026 14:19:17

ЭКОНОМИЧЕСКАЯ БЕЗОПАСНОСТЬ

ЭКОНОМИЧЕСКАЯ БЕЗОПАСНОСТЬ Статья в формате PDF 100 KB...

10 06 2026 5:58:24

ВИРТУАЛЬНЫЕ ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЕ ПРИБОРЫ

ВИРТУАЛЬНЫЕ ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЕ ПРИБОРЫ Статья в формате PDF 265 KB...

31 05 2026 18:38:30

ИСПОЛЬЗОВАНИЕ МОЛЕКУЛЯРНЫХ ТЕСТОВ В РАЗВЕДЕНИИ ЖИВОТНЫХ

ИСПОЛЬЗОВАНИЕ МОЛЕКУЛЯРНЫХ ТЕСТОВ В РАЗВЕДЕНИИ ЖИВОТНЫХ Изучена активность оксидоредуктаз в митохондриях различных органов свиней трех линий породы СМ-1 новосибирской селекции. Исследована активность цитохромоксидазы, сукцинатдегидрогеназы в митохондриях, супернатанте печении и сердца животных. Анализ всех экспериментальных групп показал, что по изменению ферментативной активности митохондрий лучшими являются свиньи линий Светлого и Совета. Энергию клетке поставляют митохондоии. В состав митохондрий входят цитохромы, в частности, цитохром аа3(цитохромоксидаза), сукцинатдегидрогеназа. Во внутренней митохондриальной мембране приблизительно четвертую часть от общего белка составляют ферменты, которые принимают участие в трaнcпорте электронов и тканевом дыхании: флавопротеиды, цитохромы и ферменты, участвующие в синтезе макроэргов. Остальная часть общего белка внутренней мембраны митохондрий выполняет структурные функции вместе с входящими в ее состав липидами [1]. ...

29 05 2026 11:35:29

ДАШКЕВИЧ ЮРИЙ МИХАЙЛОВИЧ

ДАШКЕВИЧ ЮРИЙ МИХАЙЛОВИЧ Статья в формате PDF 64 KB...

21 05 2026 8:51:15

СТОИТ ЛИ ИЗУЧАТЬ ГЕОМЕТРИЮ С ПЯТОГО КЛАССА?

СТОИТ ЛИ ИЗУЧАТЬ ГЕОМЕТРИЮ С ПЯТОГО КЛАССА? Статья в формате PDF 250 KB...

16 05 2026 7:52:55

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

ОСОБЕННОСТИ АНАТОМИЧЕСКОЙ СТРУКТУРЫ РОГОЗА УЗКОЛИСТНОГО В УСЛОВИЯХ НАГРУЗКИ ПО СВИНЦУ Установлено влияние уксуснокислого свинца (2,5∙10–1 мг/л) на анатомическое строение почвенных и водных корней рогоза узколистного (Typha angustifolia L.). Происходит адаптационное перераспределение активности разрушения паренхимных клеток и образования воздухоносных полостей с водных корней, непосредственно контактирующих с растворенной в воде солью, на почвенные. Объем воздухоносных полостей специфичен периоду вегетации растений и возрасту корней. ...

12 05 2026 6:11:50

ФАКТОРНЫЙ АНАЛИЗ ЗАГРЯЗНЕНИЯ РОДНИКОВ

ФАКТОРНЫЙ АНАЛИЗ ЗАГРЯЗНЕНИЯ РОДНИКОВ Получены закономерности взаимного влияния концентрации по 22 видам загрязнения семи родников, отобранных для исследования моделированием взаимосвязей между факторами. Дана полная корреляционная матрица монарных (на основе рангового или рейтингового распределения) и бинарных (между парами взаимно влияющих факторов) связей. Коэффициент функциональной связности равен сумме коэффициентов корреляции, разделенной на произведение числа строк на количество столбцов. Этот статистический показатель для всей сети родников применим при сопоставлении разных территорий. Первое место как влияющий параметр занимает общее микробное число, а как зависимый показатель – цветность. Анализ всех 484 моделей показал, что высокой предсказательной силой обладают слабые и средние факторные связи. Они же зачастую приводят к научно-техническим решениям мировой новизны на уровне изобретений. ...

11 05 2026 16:19:28

ЧЕРНОСЛИВ В ПРОИЗВОДСТВЕ СЫРОКОПЧЕНЫХ КОЛБАС

ЧЕРНОСЛИВ В ПРОИЗВОДСТВЕ СЫРОКОПЧЕНЫХ КОЛБАС Статья в формате PDF 244 KB...

10 05 2026 9:40:32

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РЫБОЛОВНЫХ СИСТЕМ

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РЫБОЛОВНЫХ СИСТЕМ Статья в формате PDF 106 KB...

09 05 2026 16:29:43

ОПЫТ ПРИМЕНЕНИЯ МЕТОДА АГРОСТЕПЕЙ ДЛЯ ВОССТАНОВЛЕНИЯ НАРУШЕННОЙ РАСТИТЕЛЬНОСТИ ДОЛИНЫ СРЕДНЕЙ ЛЕНЫ (ЦЕНТРАЛЬНАЯ ЯКУТИЯ)

ОПЫТ ПРИМЕНЕНИЯ МЕТОДА АГРОСТЕПЕЙ ДЛЯ ВОССТАНОВЛЕНИЯ НАРУШЕННОЙ РАСТИТЕЛЬНОСТИ ДОЛИНЫ СРЕДНЕЙ ЛЕНЫ (ЦЕНТРАЛЬНАЯ ЯКУТИЯ) Анализ опыта по восстановлению методом агростепей растительности на нарушенных кормовых угодьях долины средней Лены показал, что метод при соблюдении экологических условий и видового состава участков обеспечивает восстановление растительности, проявляющееся в повышении проективного покрытия и доминировании в травостое целинных видов. Соответствие экологических условий и видового состава травостоя при подборе участков обеспечивает восстановление растительности нарушенных участков до 70–75 % и доминирование в травостое целинных видов до 60–65 % в условиях нормального и сильного засоления. ...

06 05 2026 17:25:12

Еще:
Поддержать себя -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 ::