ПАРАЛЛЕЛЬНЫЙ АЛГОРИТМ (2,1)-МЕТОДА ПЕРЕМЕННОГО ШАГА
Предложен параллельный алгоритм переменного шага на основе (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.
Список литературы
- Новиков Е.А. Явные методы для жестких систем. - Новосибирск: Наука, 1997.
- Хайрер Э., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-алгебраические задачи. - М.: Мир, 1999.
- Исаев С.В., Малышев А.В., Шайдуров В.В. Развитие Красноярского центра параллельных вычислений // Вычислительные технологии. - 2006. - №11. - С. 28-33.
Статья в формате PDF 263 KB...
19 04 2024 2:24:34
Статья в формате PDF 110 KB...
18 04 2024 17:13:29
Статья в формате PDF 299 KB...
17 04 2024 1:11:47
В статье рассмотрены реакции 1,3-дегидроадамантана, относящегося к напряженным мостиковым [3.3.1]пропелланам, с диметилтрисульфидом. Установлено, что при взаимодействии образуются 1,3-бис(метилтио)адамантан, 1-(метилдитио)-3-(метилтио)адамантан и 1,3-бис(метилдитио)адамантан в соотношении 1:4,5:1. Структуры полученных соединений подтверждены методами хромато-масс-спектометрии и ЯМР1Н-спектроскопии. Выход целевого 1-(метилдитио)-3-(метилтио)адамантана составляет 50 %. Было предположено, что реакция протекает по радикальному механизму. Приведено описание эксперимента. ...
15 04 2024 5:23:26
Статья в формате PDF 258 KB...
14 04 2024 18:51:50
Статья в формате PDF 161 KB...
11 04 2024 22:45:36
Статья в формате PDF 142 KB...
10 04 2024 7:25:46
Статья в формате PDF 142 KB...
09 04 2024 9:19:42
Статья в формате PDF 138 KB...
08 04 2024 15:14:59
Статья в формате PDF 134 KB...
06 04 2024 21:19:45
Статья в формате PDF 300 KB...
05 04 2024 13:10:16
Статья в формате PDF 112 KB...
04 04 2024 19:34:48
02 04 2024 3:19:13
Статья в формате PDF 241 KB...
01 04 2024 17:52:56
Статья в формате PDF 141 KB...
31 03 2024 15:34:40
Статья в формате PDF 358 KB...
30 03 2024 6:37:21
Статья в формате PDF 120 KB...
29 03 2024 13:36:47
Статья в формате PDF 120 KB...
28 03 2024 0:55:19
Статья в формате PDF 144 KB...
27 03 2024 19:15:14
Статья в формате PDF 265 KB...
25 03 2024 14:26:31
Статья в формате PDF 123 KB...
24 03 2024 10:32:47
Статья в формате PDF 124 KB...
23 03 2024 19:21:10
Приведены результаты научных исследований сохранения и улучшения экологического состояния агроландшафтов Казахстана. Проведены экспериментальные работы с учетом дифференциации зональных систем земледелия. Исследования показали, что оценка в эрозионных агроландшафтах адаптивности основной обработки богарных светло-каштановых почв на уровне мезо – и микроландшафтных условий, вспашка более эффективна в северных и восточных экспозиций склонов, где плотность пахотного слоя была в среднем за вегетацию зерновых культур в основном на 0,02–0,04 г/см3 меньше по сравнению с плоскорезной обработкой. На склонах южной и западной экспозиций наоборот плоскорезная обработка способствовала снижению уплотненности почвы, на 0,03–0,05 г/см3 и повышению ее противоэрозионной устойчивости в 1,2–1,5 раза. На склонах северной и восточной экспозиции вспашка обеспечивает более эффективную борьбу с сорняками, а плоскорезная – на южных и западных склонах более высокое и равномерное накопление снега и рациональное использование влаги. Важнейшим звеном улучшения экологии почв является оптимизация севооборотов. В статье предлагается построить севооборот по количеству оставляемого в почве органического вещества, каждым предшественником. Для совершенствования севооборотов рекомендуется сидерация, уплотненные посевы, размещение многолетних и однолетних трав, применения органических удобрений и др. ...
22 03 2024 17:20:46
Статья в формате PDF 114 KB...
21 03 2024 4:18:21
Статья в формате PDF 146 KB...
20 03 2024 19:13:15
Статья в формате PDF 125 KB...
19 03 2024 0:37:17
Статья в формате PDF 128 KB...
18 03 2024 4:54:57
Настоящая работа посвящена экономико-математическому моделированию процесса кадрового обеспечения организации с учетом основных положений и методов индустриально-организационной психологии [1]. ...
17 03 2024 10:39:16
Статья в формате PDF 298 KB...
16 03 2024 10:36:15
Статья в формате PDF 120 KB...
14 03 2024 2:15:15
Статья в формате PDF 250 KB...
13 03 2024 21:40:57
Статья в формате PDF 170 KB...
11 03 2024 9:36: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 ::