|
Текст подпрограммы и версий ( Фортран ) aig0r.zip aig0d.zip aig0c.zip |
Тексты тестовых примеров ( Фортран ) taig0r.zip taig0d.zip taig0c.zip |
|
Текст подпрограммы и версий ( Си ) aig0r_c.zip aig0d_c.zip aig0c_c.zip |
Тексты тестовых примеров ( Си ) taig0r_c.zip taig0d_c.zip taig0c_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) aig0r_p.zip aig0e_p.zip aig0c_p.zip |
Тексты тестовых примеров ( Паскаль ) taig0r_p.zip taig0e_p.zip taig0c_p.zip |
Обращение вещественной матрицы методом Жордана без выбора ведущего элемента.
Заданная матрица А N - го порядка приводится к единичной матрице Е поледовательностью элементарных преобразований Жордана L1, L2, ..., LN так, что LNLN-1...L1*А = Е, при этом А-1 = LNLN-1 ... L1.
В.В.Воеводин, Р.В.Петрина, Комплекс алгоритмов, основанных на преобразованиях типа Гаусса, в пакете линейной алгебры. Сб. "Численный анализ на ФОPТPАНе", вып.3, Изд-во МГУ, 1973.
SUBROUTINE AIG0R (A, N)
Параметры
| A - | двумерный N на N массив, в котором задается исходная матрица; по окончании работы подпрограммы в массиве А запоминается вычисленная обратная матрица (тип: вещественный); |
| N - | заданный порядок матрицы (тип: целый). |
Версии
| AIG0D - | обращение матрицы, заданной с удвоенной точностью, методом Жордана без выбора ведущего элемента. |
| AIG0C - | обращение комплексной матрицы методом Жордана без выбора ведущего элемента. |
Вызываемые подпрограммы : нет
Замечания по использованию
| 1. |
В подпрограмме АIG0D массив А имеет тип DОUВLЕ РRЕСISIОN. | |
| 2. | В подпрограмме АIG0С массив А имеет тип СОМР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, 0.8, -3.2,
* -8.9, -7.2, 3.5, 9.3, 3.3/
CALL AIG0R (A, 4)
Результат:
| 0.0505 0.0542 0.0062 0.0350 |
| 0.0518 -0.0846 0.0721 -0.0003 |
A-1 = | -0.0084 0.0432 0.0202 -0.1211 |
| -0.0497 0.0279 0.0790 -0.0577 |