|
Текст подпрограммы и версий ( Фортран ) ask6c.zip |
Тексты тестовых примеров ( Фортран ) task6c.zip |
|
Текст подпрограммы и версий ( Си ) ask6c_c.zip |
Тексты тестовых примеров ( Си ) task6c_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) ask6c_p.zip |
Тексты тестовых примеров ( Паскаль ) task6c_p.zip |
Решение системы линейных алгебраических уравнений с клеточно - циркулянтной матрицей, все клетки которой суть клеточно - циркулянтные матрицы общего вида.
Решается система А*x = b с клеточно - циркулянтной матрицей А (см. [1]) клеточного размера K на К, с клетками размера (М*L) на (М*L); клетки матрицы А сами являются клеточно - циркулянтными матрицами клеточного размера L на L, с клетками общего вида размера М на М.
В основе алгорифма лежит подобное преобразование матрицы А к клеточно - диагональной матрице клеточного размера К на К, с клетками размера (М*L) на (М*L).
| 1. | В.В.Воеводин. Вычислительные основы линейной алгебры. М., Наука, 1977. |
| 2. | Е.Е.Тыртышников. О решении систем с матрицами типа теплицевых. В сб. "Численный анализ на ФОРТРАНе. Вычисл. Методы и инструментальные системы". Изд-во МГУ, 1978. |
SUBROUTINE ASK6C ( A, X, R, M, L, K)
Параметры
| A - | комплексный одномерный массив длины М2 на L на К, в котором задается матрица системы следующим образом. Массив А делится на К кусков длины М2 на L. В I - ом куске, I = 1, ..., К, задается I - ая клетка первой клеточной строки исходной матрицы. Задание I - ой клетки в своем куске осуществляется следующим образом: соответствующий кусок делится на L сегментов длины М2; в J - ом сегменте, J = 1, ..., L, записывается J - ая клетка первой клеточной строки I - ой клетки. Все внутренние клетки в своих сегментах записываются по столбцам; |
| X - | комплексный одномерный массив длины М на L на К, в котором задается правая часть; на выходе из подпрограммы в Х получается решение; |
| R - | комплексный одномерный рабочий массив длины max (Q, L, К), где Q - целая часть числа (М + 1)/2; |
| M - | порядок внутренней клетки (тип: целый); |
| L - | клеточный порядок клеток исходной матрицы (тип: целый); |
| K - | клеточный порядок исходной матрицы (тип: целый). |
Версии: нет
Вызываемые подпрограммы
| ASK3C - | решение системы линейных алгебраических уравнений с клеточно - циркулянтной матрицей все клетки которой суть матрицы общего вида. |
| ASG1C - | решение системы линейных алгебраических уравнений методом Гаусса с выбором по столбцу. |
Замечания по использованию
| В подпрограмме АSК6С используются служебные подпрограммы АR01С, АR02С. |
COMPLEX A(16), X(8), R(2)
DATA A / (1., 1.), (2., 2.), (2., 2.), (1., 1.), (2., 2.), (4., 4.),
* (4., 4.), (2., 2.), (2., 2.), (4., 4.), (4., 4.), (2., 2.),
* (4., 4.), (8., 8.), (8., 8.), (4., 4.) /
DATA X /8* (27., 27.)/
M = 2
L = 2
K = 2
CALL ASK6C (A, X, R, M, L, K)
Результат:
X = ( (1., 0.), (1., 0.), (1., 0.), (1., 0.), (1., 0.), (1., 0.),
(1., 0.), (1., 0.) )