|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://num-anal.srcc.msu.ru/lib_na/cat/af/afg0r.htm
Дата изменения: Fri Jun 20 10:59:10 2014 Дата индексирования: Sun Apr 10 01:09:12 2016 Кодировка: Windows-1251 |
|
Текст подпрограммы и версий ( Фортран ) afg0r.zip , afg0d.zip , afg0c.zip |
Тексты тестовых примеров ( Фортран ) tafg0r.zip , tafg0d.zip , tafg0c.zip |
|
Текст подпрограммы и версий ( Си ) afg0r_c.zip , afg0d_c.zip , afg0c_c.zip |
Тексты тестовых примеров ( Си ) tafg0r_c.zip , tafg0d_c.zip , tafg0c_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) afg0r_p.zip , afg0e_p.zip , afg0c_p.zip |
Тексты тестовых примеров ( Паскаль ) tafg0r_p.zip , tafg0e_p.zip , tafg0c_p.zip |
Треугольное разложение вещественной матрицы методом Гаусса без выбора ведущего элемента.
Выполняется факторизация матрицы А вида А = LU , где L - нижняя, U - верхняя треугольные матрицы. Матрица L имеет единичные диагональные элементы.
В.В.Воеводин, Р.В.Петрина, Комплекс алгоритмов, основанных на преобразованиях Гаусса, в пакете линейной алгебры, Сб. "Численный анализ на ФОРТРАНе", вып.3, Изд-во МГУ, 1973.
SUBROUTINE AFG0R (A, N)
Параметры
| A - | двумерный N*N массив, в котором задается исходная матрица (тип: RЕАL); по окончании работы подпрограммы в массиве А на соответствующих местах запоминаются элементы матрицы U и поддиагональные элементы матрицы L; |
| N - | заданный порядок исходной матрицы (тип: целый). |
Версии
| AFG0D - | треугольное разложение матрицы, заданной с удвоенной точностью, методом Гаусса без выбора ведущего элемента. |
| AFG0C - | треугольное разложение комплексной матрицы методом Гаусса без выбора ведущего элемента. |
Вызываемые подпрограммы: нет
Замечания по использованию
| 1. |
В подпрограмме АFG0D массив А имеет тип DОUВLЕ РRЕСISIОN. | |
| 2. | В подпрограмме АFG0С массив А имеет тип СОМРLЕХ. |
DIMENSION A(4, 4)
DATA A /7.9, 8.5, 4.3, 3.2, 5.6, -4.8, 4.2, -1.4, 5.7, .8, -3.2,
* -8.9, -7.2, 3.5, 9.3, 3.3/
CALL AFG0R (A, 4)
Результат:
| 7.900 5.600 5.700 -7.200 |
| 1.076 -10.825 -5.333 11.247 |
A = | 0.544 -0.106 -6.870 14.415 | ,
| 0.405 0.339 1.369 -17.323 |
это означает, что
| 1 0 0 0 |
| 1.076 1 0 0 |
L = | 0.544 -0.106 1 0 | ,
| 0.405 0.339 1.369 1 |
| 7.900 5.600 5.700 -7.200 |
| 0 -10.825 -5.333 11.247 |
U = | 0 0 -6.870 14.415 | .
| 0 0 0 -17.323 |