|
Текст подпрограммы и версий ( Фортран ) asc0c.zip |
Тексты тестовых примеров ( Фортран ) tasc0c.zip |
|
Текст подпрограммы и версий ( Си ) asc0c_c.zip |
Тексты тестовых примеров ( Си ) tasc0c_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) asc0c_p.zip |
Тексты тестовых примеров ( Паскаль ) tasc0c_p.zip |
Решение системы линейных алгебраических уравнений с циркулянтной матрицей.
Решается система А*x = b, где А - циркулянтная матрица порядка М. В основе алгорифма лежит подобное преобразование матрицы А к диагональной матрице на базе дискретного преобразования Фурье.
С.Н.Воеводина. Решение системы уравнений с клеточно - теплицевыми матрицами. В Сб. Вычислительные методы и программирование, вып. 24, Изд-во МГУ, 1975.
SUBROUTINE ASC0C ( A, X, R, M)
Параметры
| A - | комплексный одномерный массив длины М, в котором матрица системы задается своей первой строкой; |
| X - | комплексный одномерный массив длины М, в котором задается правая часть. На выходе из подпрограммы в Х получается решение; |
| R - | комплексный одномерный рабочий массив длины М; |
| M - | порядок матрицы (тип: целый). |
Версии: нет
Вызываемые подпрограммы: нет
Замечания по использованию: нет
COMPLEX A(4), X(4), R(4)
DATA A / (1., 1.), (2., 2.), (3., 3.), (4., 4.)/
DATA X /4* (10., 10.)/
M = 4
CALL ASC0C (A, X, R, M)
Результаты:
X = ( (1., 0.), (1., 0.), (1., 0.), (1., 0.) )