|
Текст подпрограммы и версий ( Фортран ) aet2c.zip , aet2p.zip |
Тексты тестовых примеров ( Фортран ) taet2c.zip , taet2p.zip |
|
Текст подпрограммы и версий ( Си ) aet2c_c.zip |
Тексты тестовых примеров ( Си ) taet2c_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) aet2c_p.zip , aet2z_p.zip |
Тексты тестовых примеров ( Паскаль ) taet2c_p.zip , taet2z_p.zip |
Вычисление всех собственных значений комплексной верхней матрицы Хессенберга QR - алгоpитмом со сдвигом.
Подпрограмма АЕТ2С вычисляет собственные значения верхней матрицы Хессенберга QR - алгоpитмом со сдвигом.
Предполагается, что исходная матрица была предварительно промасштабирована с помощью подпрограммы АМВ1С.
Уилкинсон Дж.Х. Алгебраическая проблема собственных значений. М.: "Наука", 1970.
SUBROUTINE AET2C (NM, N, LOW, IGH, HR, HI, WR, WI, IERR)
Параметры
| NM - | число строк двумерных массивов НR и НI, указанное при описании этих массивов в вызывающей подпрограмме (тип: целый); |
| N - | порядок исходной матрицы (тип: целый); |
|
LOW - IGH | выходные параметры подпрограммы АМВ1С (тип: целый); если матрица не масштабировалась, то можно взять LОW = 1, IGН = N; |
| HR, HI - | вещественные двумерные массивы размерности NМ на N, содержащие на входе в подпрограмму в своих первых N строках соответственно вещественную и мнимую части исходной матрицы Хессенберга; |
| WR, WI - | вещественные векторы длины N, содержащие на выходе из подпрограммы соответственно вещественные и мнимые части вычисленных собственных значений; |
| IERR - | целочисленная переменная, служащая для сообщения об ошибках, обнаруженных при работе подпрограммы; значение IЕRR полагается равным индексу собственного значения, для вычисления которого требуется более 30 итераций, при этом собственные значения с индексами IЕRR + 1, IЕRR + 2, ..., N правильны, если вычислены все собственные значения, то IЕRR = 0 . |
Версии
| AET2P - | вычисление всех собственных значений комплексной верхней матрицы Хессенберга, заданной с двойной точностью, QR - алгоpитмом со сдвигом. |
Вызываемые подпрограммы
| UTAE12 - | подпрограмма выдачи диагностических сообщений при работе подпрограмм АЕТ2С и АЕТ2Р. |
Замечания по использованию
| 1. |
В подпрограмме АЕТ2Р паpаметpы НR, НI, WR, WI имеют тип DОUВLЕ РRЕСISIОN. | |
| 2. | Подпрограмма АЕТ2R (АЕТ2Р) не сохраняет исходную матрицу. |
DIMENSION HR(4, 4), HI(4, 4), WR(4), WI(4)
DATA HR /4., 3., 0., 0., -2., -3., 2., 0., -2., -6., 6., 1., 2., 6., -2., 3./
DATA HI /1., 4*0., 1., 4*0., 1., 4*0., 1./
N = 4
CALL AET2C (N, N, 1, N, HR, HI, WR, WI, IERR)
Результаты:
IERR = 0
WR = (1., 2., 3., 4.)
WI = (1., 1., 1., 1.)