ПАРАЛЛЕЛЬНЫЕ ЯВНЫЕ ОДНОШАГОВЫЕ МЕТОДЫ ДЛЯ ЧИСЛЕННОГО РЕШЕНИЯ ЖЕСТКИХ СИСТЕМ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
Предложены параллельные явные одношаговые методы первого, второго порядков, обеспечивающие возможность с минимальными вычислительными затратами интегрировать жесткие системы обыкновенных дифференциальных уравнений. В предлагаемых параллельных алгоритмах изменение величины шага построены на основе контроля точности и устойчивости численной схемы, а в неравенстве для контроля точности применяется оценка локальной ошибки метода.
В настоящее время одним из основных параметром, ха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 / hn |λmax | )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 % и ниже.
Список литературы
- Ващенко Г.В., Новиков Е.А. Параллельная реализация явных методов типа Рунге-Кутты // Вестник КрасГАУ. - 2010 - №2 - С. 14-18.
- Новиков Е.А. Явные методы для жестких систем. - Новосибирск: Наука, 1997.
- Хайрер Э., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-алгебраические задачи. - М.: Мир, 1999.
- Исаев С.В., Малышев А.В., Шайдуров В.В. Развитие Красноярского центра параллельных вычислений // Вычислительные технологии. - 2006. - №11. - С. 28-33.
Статья в формате PDF 100 KB...
28 03 2024 2:42:43
На биопсийном материале пяти первородящих женщин в возрасте от 20 до 38 лет с физиологической родовой деятельностью проводили количественное светооптическое изучение строения миометрия матки. Оценили тканевой состав, клеточный состав и число гладкомышечных клеток в поле зрения микроскопа. Показали, что основными компонентами миометрия были гладкомышечные волокна, элементы соединительной ткани и микрососудистого русла. Гладкомышечные клетки демонстрировали разное сродство к толуидиновому синему и были условно разделены на светлые, темные и промежуточные клетки. Выявлена внутригрупповая вариация всех оцененных количественных параметров. Полученные данные могут быть базовыми при оценке тех же параметров у рожениц с патологической родовой деятельностью. ...
27 03 2024 0:10:56
Статья в формате PDF 127 KB...
26 03 2024 12:46:26
Статья в формате PDF 274 KB...
25 03 2024 23:23:41
Статья в формате PDF 138 KB...
24 03 2024 22:46:14
Статья в формате PDF 112 KB...
23 03 2024 0:37:24
В статье рассматривается проблема изучения понятия «политическая дискурсия» в современной лингвистике. Развитие когнитивной парадигмы в лингвистике актуализировало изучение понятия «дикурс», исследование же политической дискурсии даёт возможность исследовать подробнее языковую личность политика. ...
22 03 2024 23:19:33
Рассматривается проблема организации продуктивной умственной деятельности учащихся общеобразовательных учреждений в системе дидактических принципов современной педагогики. Анализ принципов показывает, что отечественная дидактика в большой мере сохраняет черты традиционной модели обучения и недостаточно учитывает психологическую природу мышления и закономерности продуктивной умственной деятельности при разработке принципов обучения. Выделены основополагающие принципы организации продуктивной умственной деятельности на основе закономерностей развития знания и процесса познания, психологических закономерностей мышления. ...
21 03 2024 12:58:14
Статья в формате PDF 102 KB...
19 03 2024 21:42:24
Статья в формате PDF 114 KB...
18 03 2024 22:16:31
Статья в формате PDF 235 KB...
17 03 2024 21:30:51
Статья в формате PDF 100 KB...
16 03 2024 3:43:24
Статья в формате PDF 124 KB...
15 03 2024 3:52:28
Статья в формате PDF 113 KB...
14 03 2024 17:15:15
Статья в формате PDF 289 KB...
13 03 2024 17:20:32
Статья в формате PDF 250 KB...
12 03 2024 13:33:23
Статья в формате PDF 232 KB...
11 03 2024 5:57:56
Статья в формате PDF 108 KB...
10 03 2024 20:32:27
Статья в формате PDF 125 KB...
09 03 2024 2:24:46
Статья в формате PDF 113 KB...
08 03 2024 11:22:53
Статья в формате PDF 109 KB...
07 03 2024 14:52:33
Статья в формате PDF 103 KB...
06 03 2024 5:43:49
Статья в формате PDF 301 KB...
05 03 2024 1:47:44
Статья в формате PDF 100 KB...
04 03 2024 9:40:47
Статья в формате PDF 113 KB...
03 03 2024 11:52:43
Статья в формате PDF 244 KB...
02 03 2024 16:24:23
29 02 2024 21:58:40
Статья в формате PDF 116 KB...
28 02 2024 6:21:55
Статья в формате PDF 115 KB...
27 02 2024 19:24:38
Статья в формате PDF 112 KB...
26 02 2024 23:51:59
Статья в формате PDF 323 KB...
25 02 2024 4:20:33
Статья в формате PDF 110 KB...
24 02 2024 5:10:24
Статья в формате PDF 131 KB...
23 02 2024 0:15:23
Статья в формате PDF 113 KB...
21 02 2024 19:56:45
Статья в формате PDF 119 KB...
20 02 2024 7:46:36
Статья в формате PDF 140 KB...
19 02 2024 5:15:27
Статья в формате PDF 104 KB...
18 02 2024 18:15:31
Еще:
Поддержать себя -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 ::