Тысяча полезных мелочей    

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

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

Ващенко Г.В. Статья в формате 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 240 KB...

13 07 2026 11:25:31

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

СОВРЕМЕННОЕ СОСТОЯНИЕ И ТЕНДЕНЦИИ РАЗВИТИЯ ГИДРОГРАФО – ГЕОДЕЗИЧЕСКОГО ОБЕСПЕЧЕНИЯ МОРСКИХ ГЕОЛОГОРАЗВЕДОЧНЫХ РАБОТ В статье обсуждаются последние достижения технологий гидрографических и геодезических съемок, таких как дифференциальная система GPS/ГЛОНАСС субметровой точности определения положения на поверхности моря, интегрированная DGPS с гидроакустической системой HPR для определения положения под водой, многолучевые эхолоты, гидролокаторы бокового обзора; морские датчики движения, специально разработанные для высокоточного измерения перемещений в море для пользователей, требующих высокой точности измерений дифферента, крена и перемещений по высоте. Аэролазерная батиметрия имеет значительный потенциал для замены эхолота при измерении глубин. Отмечено, что ROV (буксируемые подводные аппараты) и AUV (автономные подводные аппараты) становятся технически и экономически более выгодной платформой для съемки в специальных применениях и в будущем станут широко использующейся техникой. ...

04 07 2026 17:27:53

О ФИЗИОЛОГИИ РАЗВИТИЯ ЛИМФАТИЧЕСКОЙ СИСТЕМЫ

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

02 07 2026 17:43:46

СТАНДАРТИЗАЦИЯ СОЦИАЛЬНОЙ РЕКЛАМЫ

СТАНДАРТИЗАЦИЯ СОЦИАЛЬНОЙ РЕКЛАМЫ Статья в формате PDF 182 KB...

30 06 2026 9:12:28

ВОДНЫЕ РЕСУРСЫ КРУПНЫХ РЕК

ВОДНЫЕ РЕСУРСЫ КРУПНЫХ РЕК Приведены методы ранжирования и рангового моделирования гидрологических параметров у множества крупных рек Земли по примеру статистических данных из учебника. ...

28 06 2026 5:41:49

КОМПЬЮТЕРНЫЕ ТЕХНОЛОГИИ В ОБУЧЕНИИ ХИМИИ

КОМПЬЮТЕРНЫЕ ТЕХНОЛОГИИ В ОБУЧЕНИИ ХИМИИ Статья в формате PDF 251 KB...

27 06 2026 8:14:40

СОЦИАЛЬНАЯ РАБОТА КАК ТЕОРЕТИКО-МЕТОДОЛОГИЧЕСКИЙ ИНСТРУМЕНТ СОЦИАЛИЗАЦИИ В РОССИИ

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

24 06 2026 5:19:35

О ПОСЛЕДСТВИЯХ ОПУСТЫНИВАНИЯ В РАВНИННОМ ДАГЕСТАНЕ

О ПОСЛЕДСТВИЯХ ОПУСТЫНИВАНИЯ В РАВНИННОМ ДАГЕСТАНЕ Статья в формате PDF 117 KB...

22 06 2026 14:42:29

ОПРЕДЕЛЕНИЕ НАПРЯЖЕНИЙ В СЕЧЕНИЯХ ДЕТАЛЕЙ РВП

ОПРЕДЕЛЕНИЕ НАПРЯЖЕНИЙ В СЕЧЕНИЯХ ДЕТАЛЕЙ РВП Статья в формате PDF 589 KB...

21 06 2026 15:59:52

ПРОМЫШЛЕННОЕ ЗАГРЯЗНЕНИЕ ПОЧВ Г. О. ШУИ

ПРОМЫШЛЕННОЕ ЗАГРЯЗНЕНИЕ ПОЧВ Г. О. ШУИ Статья в формате PDF 250 KB...

20 06 2026 3:12:20

РАСЧЕТ И КОМПЕНСАЦИЯ МАГНИТНОЙ ДЕВИАЦИИ

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

15 06 2026 9:25:22

Алгоритм проведения дифференциальной диагностики

Алгоритм проведения дифференциальной диагностики Статья в формате PDF 104 KB...

12 06 2026 20:52:29

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